如何用maple求解贝塞尔函数零点
答案:2 悬赏:10 手机版
解决时间 2021-03-23 03:57
- 提问者网友：最美的风景
- 2021-03-22 15:21
如何用maple求解贝塞尔函数零点
最佳答案
- 五星知识达人网友：风格不统一
- 2021-03-22 16:27
FUNCTION bessj0(x)
REAL bessj0,x
REAL ax,xx,z
DOUBLE PRECISION p1,p2,p3,p4,p5,q1,q2,q3,q4,q5,
r1,r2,r3,r4,r5,r6,s1,s2,s3,s4,s5,s6,y
SAVE p1,p2,p3,p4,p5,q1,q2,q3,q4,q5,r1,r2,r3,r4,r5,r6,
s1,s2,s3,s4,s5,s6
DATA p1,p2,p3,p4,p5/1.d0,-.d-2,.d-4,
-.d-5,.d-6/, q1,q2,q3,q4,q5
/-.d-1,.d-3,-.d-5,
.d-6,-.d-7/
DATA r1,r2,r3,r4,r5,r6/.d0,-.d0,
.7d0,-.18d0,77392.33017d0,
-184.d0/,s1,s2,s3,s4,s5,s6/.d0,
.d0,.718d0,59272.64853d0,
267.d0,1.d0/
if(abs(x)<8.) then
y=x**2
bessj0=(r1+y*(r2+y*(r3+y*(r4+y*(r5+y*r6)))))/
(s1+y*(s2+y*(s3+y*(s4+y*(s5+y*s6)))))
else
ax=abs(x)
z=8./ax
y=z**2
xx=ax-.
bessj0=sqrt(./ax)*(cos(xx)*(p1+y*(p2+y*
(p3+y*(p4+y*p5))))-z*sin(xx)*(q1+y*(q2+y*
(q3+y*(q4+y*q5)))))
endif
END FUNCTION bessj0
REAL bessj0,x
REAL ax,xx,z
DOUBLE PRECISION p1,p2,p3,p4,p5,q1,q2,q3,q4,q5,
r1,r2,r3,r4,r5,r6,s1,s2,s3,s4,s5,s6,y
SAVE p1,p2,p3,p4,p5,q1,q2,q3,q4,q5,r1,r2,r3,r4,r5,r6,
s1,s2,s3,s4,s5,s6
DATA p1,p2,p3,p4,p5/1.d0,-.d-2,.d-4,
-.d-5,.d-6/, q1,q2,q3,q4,q5
/-.d-1,.d-3,-.d-5,
.d-6,-.d-7/
DATA r1,r2,r3,r4,r5,r6/.d0,-.d0,
.7d0,-.18d0,77392.33017d0,
-184.d0/,s1,s2,s3,s4,s5,s6/.d0,
.d0,.718d0,59272.64853d0,
267.d0,1.d0/
if(abs(x)<8.) then
y=x**2
bessj0=(r1+y*(r2+y*(r3+y*(r4+y*(r5+y*r6)))))/
(s1+y*(s2+y*(s3+y*(s4+y*(s5+y*s6)))))
else
ax=abs(x)
z=8./ax
y=z**2
xx=ax-.
bessj0=sqrt(./ax)*(cos(xx)*(p1+y*(p2+y*
(p3+y*(p4+y*p5))))-z*sin(xx)*(q1+y*(q2+y*
(q3+y*(q4+y*q5)))))
endif
END FUNCTION bessj0
全部回答
- 1楼网友：一把行者刀
- 2021-03-22 18:06
支持一下感觉挺不错的
我要举报
如以上问答信息为低俗、色情、不良、暴力、侵权、涉及违法等信息，可以点下面链接进行举报！
大家都在看
推荐资讯