| 查看: 614 | 回复: 0 | ||
[求助]
求助厄米多项式
|
|
下面是一个关于厄米多项式的小程序，可是在运行中发现XA(I)和XB(I)输出的数值不一样，XB(I)都是1，如果把数组XB（I）去掉那么XA（I）的输出也变成了1......哪位高手帮忙看下~~ IMPLICIT DOUBLE PRECISION(A-H,O-Z) PARAMETER (NX=256) PARAMETER (NA=30) C C.....Declare arrays. DIMENSION XA(NX),XB(NX) DIMENSION FAT(NX),HER(NX,NA) C open(20,file='GRID.dat',status='unknown') open(30,file='HERMI.dat',status='unknown') C.....Variable data input. DATA ZMU/1/ DATA W/1/ DATA R0/-5/ DATA Rf/5/ DATA PI/3.141592653589793D0/ C C....Test that NX is even ITEST=MOD(NX,2) IF(ITEST.NE.0) THEN WRITE(6,*)' **** NX MUST BE EVEN-FATAL ERROR ****' STOP END IF C C.....Set up grid WRITE(6,*)'Grid paremeters:' WRITE(6,*)' Number of grid poins = ',NX RMIN=R0 RMAX=Rf ZL=(RMAX-RMIN) WRITE(6,*)' Grid length = ',ZL DX=ZL/DFLOAT(NX-1) WRITE(6,*)' Grid spacings = ',DX C C.....COMPUTE HERMI: X=RMIN DO 003 I=1,NX XA(I)=X XB(I)=X DO 004 J=0,NA IF(J.EQ.0)THEN HER(I,J)=1 ELSEIF(J.EQ.1)THEN HER(I,J)=2*XA(I) ELSEIF(J.EQ.2)THEN HER(I,J)=4*XA(I)**2-2 ELSEIF(J.GE.3)THEN HER(I,J)=2*XA(I)*HER(I,J-1)-2*(J-1)*HER(I,J-2) ENDIF 004 CONTINUE WRITE(20,*) XA(I),XB(I) X=X+DX 003 CONTINUE DO I=1,NX WRITE(30,*)I,HER(I,1) ENDDO STOP END |
回复此楼» 猜你喜欢
2026年国自然面上资助率
已经有10人回复
-
微信指数没变化，科研之友没阅读
已经有18人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有8人回复
-
HXDI做水性聚氨酯乳液，是不是特别容易出渣
已经有3人回复
-
小木虫看见有人已经知道结果了
已经有19人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有15人回复
-
申请2027年材料/化学类博士
已经有5人回复
-
你们的时间戳变了吗
已经有4人回复
-
今年的WR进展到哪一步了？
已经有9人回复
-
产物和副产物价值比较
已经有5人回复
找到一些相关的精华帖子，希望有用哦~
一元三次方程以及一元四次方程的求根公式
已经有11人回复
- 求助——求解三元二次多项式方程组
已经有9人回复
- 求助数学大神【两个多项式如何看其是否“相近”】
已经有3人回复
- 求助一道矩阵的题！好人一生平安！
已经有3人回复
- 请教一个关于矩阵多项式的题目
已经有6人回复
- 求助：如何用MATLAB程序提取多项式系数。。。
已经有6人回复
- 求助个多项式二次型的问题
已经有8人回复
- stdafx.h是MFC中处理矩阵的函数，能不能在DSP中调用。
已经有3人回复
- 【求助】用matlab对多项式函数求导
已经有3人回复
- 【求助】怎样用chebyshev多项式拟合XRPD 背景?
已经有6人回复
- 【求助】MATLAB中多项式除法问题
已经有4人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人