def dfdr(r,theta):
out = np.sin(theta)*(np.exp(-r/2.) - (1./2.)*r*np.exp(-r/2.))
out[-1] = 0
return out
微信
公众号
扫码关注公众号
微信
公众号