下面这段程序总提示有错误,我找了很久都没找出来,请大家帮看看。
program Console1
implicit none
real R1,E1,dr,arf,pi,a0,a1,a2,ab,a
real fa,fb
integer n
n=0
pi=3.14
R1=0.0001
E1=1e9
dr=20.0
arf=0.00002
open(15,file='input.dat')
open(16,file='output.dat')
read(15,*)a0
fa(a)=a**2/R1-(2*pi*a*dr/E1)**(1.0/2.0)-arf
fb(a)=2*a/R1-1.0/2.0*(2*pi*dr/E1)**(1.0/2.0)*(a**(-1.0/2.0))
ab=a0+1.0
a1=a0
do while(abs(a1-ab)>1.0e-5)
a2=a1-fa(a1)/fb(a1)
ab=a1
a1=a2
n=n+1
end do
write(16,'(i4,f5.2)')n,a1
end program Console1
提示我的错误是:错误1 error #6410: This name has not been declared as an array or a function. [FB]
错误2 error #6410: This name has not been declared as an array or a function. [FA]
program Console1
implicit none
real R1,E1,dr,arf,pi,a0,a1,a2,ab,a
real fa,fb
integer n
n=0
pi=3.14
R1=0.0001
E1=1e9
dr=20.0
arf=0.00002
open(15,file='input.dat')
open(16,file='output.dat')
read(15,*)a0
fa(a)=a**2/R1-(2*pi*a*dr/E1)**(1.0/2.0)-arf
fb(a)=2*a/R1-1.0/2.0*(2*pi*dr/E1)**(1.0/2.0)*(a**(-1.0/2.0))
ab=a0+1.0
a1=a0
do while(abs(a1-ab)>1.0e-5)
a2=a1-fa(a1)/fb(a1)
ab=a1
a1=a2
n=n+1
end do
write(16,'(i4,f5.2)')n,a1
end program Console1
提示我的错误是:错误1 error #6410: This name has not been declared as an array or a function. [FB]
错误2 error #6410: This name has not been declared as an array or a function. [FA]