| 查看: 553 | 回复: 3 | |||
[求助]
请教高手看一下这是什么程序？Basic？ 已有1人参与
|
|
请问能告诉我它的运行顺序吗？谢谢！！ FUNCTION PLM(L,M,X) C LEGENDRE FUNCTIONS FOR (0.ge.L.le.8) IF(L.GT.0)GO TO 100 C L=0 PLM=0.7071068 RETURN 100 Y=X*X GO TO (1,2,3,4,5,6,7,8),L C L=1 1 GO TO (10,11) M+1 10 PLM=1.2247449*X RETURN 11 PLM=0.8660254*(1.-Y)**0.5 RETURN C L=2 2 GO TO (20,21,22) M+1 20 PLM=0.7905694*(-1.+3.*Y) RETURN 21 PLM=1.9364917*X*(1.-Y)**0.5 RETURN 22 PLM=0.9682458*(1.-Y) RETURN C L=3 3 GO TO (30,31,32,33) M+1 30 PLM=0.9354143*X*(-3.+5.*Y) RETURN 31 PLM=0.8100926*(1.-Y)**0.5*(-1.+5.*Y) RETURN 32 PLM=2.5617377*X*(1.-Y) RETURN 33 PLM=1.0458250*(1.-Y)**1.5 RETURN C L=4 4 GO TO (40,41,42,43,44) M+1 40 PLM=0.2651650*(3.+Y*(-30.+35.*Y)) RETURN 41 PLM=1.1858541*X*(1.-Y)**0.5*(-3.+7.*Y) RETURN 42 PLM=0.8385255*(1.-Y)*(-1.+7.*Y) RETURN 43 PLM=3.1374751*X*(1.-Y)**1.5 RETURN 44 PLM=1.1092650*(1.-Y)**2 RETURN C L=5 5 GO TO (50,51,52,53,54,55) M+1 50 PLM=0.2931510*X*(15.+Y*(-70.+63.*Y)) RETURN 51 PLM=0.8028270*(1.-Y)**0.5*(1.+Y*(-14.+21.*Y)) RETURN 52 PLM=4.2481614*X*(1.-Y)*(-1.+3.*Y) RETURN 53 PLM=0.8671523*(1.-Y)**1.5*(-1.+9.*Y) RETURN 54 PLM=3.6790157*X*(1.-Y)**2 RETURN 55 PLM=1.1634069*(1.-Y)**2.5 RETURN C L=6 6 GO TO (60,61,62,63,64,65,66) M+1 60 PLM=0.15934436*(-5.+Y*(105.+Y*(-315.+231.*Y))) RETURN 61 PLM=1.0326695*X*(1.-Y)**0.5*(5.+Y*(-30.+33.*Y)) RETURN 62 PLM=0.8163969*(1.-Y)*(1.+Y*(-18.+33.*Y)) RETURN 63 PLM=1.6327938*X*(1.-Y)**1.5*(-3.+11.*Y) RETURN 64 PLM=0.8943180*(1.-Y)**2*(-1.+11.*Y) RETURN 65 PLM=4.1947232*X*(1.-Y)**2.5 RETURN 66 PLM=1.2109123*(1.-Y)**3 RETURN C L=7 7 GO TO (70,71,72,73,74,75,76,77) M+1 70 PLM=0.1711633*X*(-35.+Y*(315.+Y*(-693.+429*Y))) RETURN 71 PLM=0.1601086*(1.-Y)**0.5*(-5.+Y*(135.+Y*(-495.+429*Y))) RETURN 72 PLM=0.3921844*X*(1.-Y)*(15.+Y*(-110.+143.*Y)) RETURN 73 PLM=0.2773162*(1-Y)**1.5*(3.+Y*(-66.+143.*Y)) RETURN 74 PLM=1.8395078*X*(1.-Y)**2*(-3.+13.*Y) RETURN 75 PLM=0.9197539*(1.-Y)**2.5*(-1.+13.*Y) RETURN 76 PLM=4.6898432*X*(1.-Y)**3 RETURN 77 PLM=1.2534133*(1.-Y)**3.5 RETURN C L=8 8 GO TO (80,81,82,83,84,85,86,87,88) M+1 80 PLM=0.022777156*(35.+Y*(-1260.+Y*(6930.+Y*(-12012.+6435.*Y)))) RETURN 81 PLM=0.1932706*X*(1.-Y)**0.5*(-35.+Y*(385.+Y*(-1001.+715.*Y))) RETURN 82 PLM=0.8085088*(1.-Y)*(-1.+Y*(33.+Y*(-143.+143.*Y))) RETURN 83 PLM=2.1894523*X*(1.-Y)**1.5*(3.+Y*(-26.+39.*Y)) RETURN 84 plm=0.8479712*(1.-y)**2*(1.+y*(-26.+65.*y)) RETURN 85 PLM=6.1148074*X*(1.-Y)**2.5*(-1.+5.*Y) RETURN 86 PLM=0.9435353*(1.-Y)**3*(-1.+15.*Y) RETURN 87 PLM=5.1679555*X*(1.-Y)**3.5 RETURN 88 PLM=1.2919889*(1.-Y)**4 RETURN END |
回复此楼» 猜你喜欢
[复试调剂]西南科技大学国防/材料导师推荐
已经有6人回复
-
化学工程321分求调剂
已经有12人回复
-
211本，11408一志愿中科院277分，曾在中科院自动化所实习
已经有4人回复
-
材料专硕326求调剂
已经有5人回复
-
东南大学364求调剂
已经有5人回复
-
国自科面上基金字体
已经有7人回复
-
药学383 求调剂
已经有4人回复
-
286求调剂
已经有5人回复
-
085601求调剂
已经有3人回复
-
302求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 杂化计算小疑问，请教高手
已经有3人回复
- Fortran有扩展名为n的文件吗？
已经有5人回复
- 请教高手 晶体是否同构
已经有6人回复
- 请教高手一个关于晶系的问题：英文rhombohedral 与 trigonal 分别是指什么晶系？
已经有4人回复
- visual studio 里面怎么运行fortran
已经有5人回复
- 请教高手这个核磁图谱与哪一个结构式（A，B,C,D）相对应
已经有11人回复
- 求高手帮我解决我的check cif
已经有3人回复
- 请教高手：comsol热应力模块中，某个物体或面满足什么条件才适合被设置为固定约束呢？
已经有8人回复
- 请教一下高手，邮箱取回密码这样显示，是什么意思？
已经有23人回复
- 优化结构finished，但是OUTCAR中没有信息。求助！！！！！
已经有15人回复
- 请教高手：单位同意报考的证明万一开不到怎么办？可不可以延后提交
已经有10人回复
- latex 页码和参考文献的问题，急！
已经有6人回复
- 看到文献中的柱形图画得很漂亮，请教高手，图中的柱形图是用什么软件画的？
已经有12人回复
- 请教高手，有知道醋酸水溶液中冰醋酸的回收有什么好办法的么 谢谢
已经有18人回复
- 请教PCR高手，使用SYBR Green I法进行荧光定量PCR对于产物大小有什么要求？
已经有21人回复
- 请教高手用chemdrew如何画COSY\nOESY弯曲的双箭头呢？请详细点，最好有图，谢谢
已经有5人回复
- 请教高手: LJ势中的指数只能是6-12吗？
已经有4人回复
- 请教高手：使用hypermesh如何对静力问题进行前处理
已经有6人回复
- 【求助】请教高手：有图--化学位移约5.2处的大包峰是什么？
已经有17人回复
- 【讨论】固溶处理后枝晶消除，形成细小等轴晶粒，如何解释？请教高手！
已经有45人回复
- 【求助】请教高手：关于sonogashira反应，到底是什么原因呢？
已经有12人回复
- 请教高手这种折线图怎么做？
已经有27人回复
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 应助: 454 (硕士)
- 金币: 20699.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.3小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
2楼2015-05-14 00:19:50
3楼2015-05-18 17:37:23
letdonkey
木虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 5456.2
- 红花: 1
- 帖子: 196
- 在线: 237.9小时
- 虫号: 1369479
- 注册: 2011-08-16
- 性别: GG
- 专业: 化学工程及工业化学
4楼2015-05-18 18:47:53
