def BSFwdNormalDelta( IsCall, Fwd, Strike, Vol, Texp, Rd, Rf = 0.0 ): d1 = (Fwd - Strike)/Vol/numpy.sqrt(Texp) return exp( -Rd * Texp ) * cnorm(d1)