| 查看: 638 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有12人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有18人回复
-
2027博士申请-临床检验诊断学、病原生物学、微生物方向
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有17人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有10人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有22人回复
-
越来越担心硕士毕业后没有竞争力了
已经有14人回复
-
论文的挂名怎么评判贡献？基金申请中，排名第一的挂名文章
已经有3人回复
-
申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请
已经有3人回复
-
闲聊
已经有6人回复
找到一些相关的精华帖子，希望有用哦~
一元三次方程以及一元四次方程的求根公式
已经有11人回复
- 求助——求解三元二次多项式方程组
已经有9人回复
- 求助数学大神【两个多项式如何看其是否“相近”】
已经有3人回复
- 求助一道矩阵的题！好人一生平安！
已经有3人回复
- 请教一个关于矩阵多项式的题目
已经有6人回复
- 求助：如何用MATLAB程序提取多项式系数。。。
已经有6人回复
- 求助个多项式二次型的问题
已经有8人回复
- stdafx.h是MFC中处理矩阵的函数，能不能在DSP中调用。
已经有3人回复
- 【求助】用matlab对多项式函数求导
已经有3人回复
- 【求助】怎样用chebyshev多项式拟合XRPD 背景?
已经有6人回复
- 【求助】MATLAB中多项式除法问题
已经有4人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人