| 查看: 579 | 回复: 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 |
回复此楼» 猜你喜欢
本人考085602 化学工程 专硕
已经有15人回复
-
283求调剂
已经有3人回复
-
274求调剂
已经有5人回复
-
302求调剂
已经有6人回复
-
11408 一志愿西电，277分求调剂
已经有3人回复
-
化学工程321分求调剂
已经有13人回复
-
341求调剂
已经有5人回复
-
材料与化工304求B区调剂
已经有7人回复
-
0703化学336分求调剂
已经有4人回复
-
[复试调剂]西南科技大学国防/材料导师推荐
已经有6人回复
找到一些相关的精华帖子，希望有用哦~
一元三次方程以及一元四次方程的求根公式
已经有11人回复
- 求助——求解三元二次多项式方程组
已经有9人回复
- 求助数学大神【两个多项式如何看其是否“相近”】
已经有3人回复
- 求助一道矩阵的题！好人一生平安！
已经有3人回复
- 请教一个关于矩阵多项式的题目
已经有6人回复
- 求助：如何用MATLAB程序提取多项式系数。。。
已经有6人回复
- 求助个多项式二次型的问题
已经有8人回复
- stdafx.h是MFC中处理矩阵的函数，能不能在DSP中调用。
已经有3人回复
- 【求助】用matlab对多项式函数求导
已经有3人回复
- 【求助】怎样用chebyshev多项式拟合XRPD 背景?
已经有6人回复
- 【求助】MATLAB中多项式除法问题
已经有4人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
5