求助]各位高手看看我的程序到底那里错了
n=3real v(3,3),a(3,3)
u=1e-6
open(1,file='v.dat')
open(2,file='ev.dat')
write(2,*)'input stiffness information'
read(1,*)((a(i,j),j=1,n),i=1,n)
write(2,100)((a(i,j),j=1,n),i=1,n)
write(2,*)'input vector information'
do 10 i=1,n
v(i,i)=1
do 15 j=1,n
if(i.ne.j) v(i,j)=0
15 continue
10 continue
write(2,100)((v(i,j),j=1,n),i=1,n)
call eig(v,a,n,u)
write(2,*)'acqure vector information'
write(2,100)((v(i),j=1,n),i=1,n)
close(1)
close(2)
100 format(1x,3f6.2)
end
subroutine eig(v,a,n,u)
real v(n,n),a(n,n)
integer p,q
66 fm=0.0
do20 i=1,n
do20 j=1,n
if((i.ne.j).and.(abs(a(i,j)).gt.fm))then
fm=abs(a(i,j))
p=i
q=j
end if
20 continue
if(fm.le.u)then
stop
endif
x=-a(p,q)
y=(a(q,q)-a(p,p))/2.0
fai=x/sqrt(x*x+y*y)
if(y.lt.0)fai=-fai
sn=1.0+sqrt(1.0-fai*fai)
sn=fai/sqrt(2*sn)
cn=sqrt(1-sn*sn)
fm=a(p,p)
a(p,p)=fm*cn*cn+a(q,q)*sn*sn+a(p,q)*fai
a(q,q)=fm*sn*sn+a(q,q)*cn*cn-a(p,q)*fai
a(p,q)=0.0
a(q,p)=0.0
do30 i=1,n
fm=v(i,p)
v(i,p)=fm*cn+v(i,q)*sn
v(i,q)=-fm*sn+v(i,q)*cn
30 contiue
goto 66
end
主要提示我第二行，和57行，我很想不明白，为什么会错 ，没有语法错误的。这是个雅克比程序！先谢谢大家了 渴望得到大家的指点！！ 这是一个fortran程序吧，这个到编程办应该比较容易得到答案
第二行错误很明显，变量的定义要在程序的顶部，把第一句和第二句换一下就好了 19行业有错误，改为write(2,100)((v(i,j),j=1,n),i=1,n) 57行单词写错了，改为30 continue 这些错误都是很基本的，建议搂主仔细检查一下程序，不难发现
页:
[1]