Question

高手帮忙解决一下啊!!!关于fortran的一个错误

我的程序是PROGRAM aero
PARAMETER(max=100)
PARAMETER(pi=3.1415926)
PARAMETER(E=2.71828)
REAL :: a(max),b(max),I,c(max)
OPEN(1,FILE='aero_in.txt')
OPEN(2,FILE='aero_ou.txt')
DO I=1,max
READ(1,*) a(I),b(I)
if(a(I)/=0) then
c(I)=b(I)+F(a(I),b(I))
write(2,*) a(I),c(I),0
else
write(2,*) a(I),b(I),0
endif
ENDDO
CONTAINS
FUNCTION F(m,n)
real :: m,n,e(7),g(7),x(14),F,sum
OPEN(3,FILE='ek.txt')
OPEN(4,FILE='aero.txt')
read(3,*) (e(I),I=2,7)
read(4,*) (x(14),I=1,14)
do I=1,7
IF(I/=1) THEN
g(I)=(sin(pi*(m**e(I))))**3
ELSE
w=(m**0.25)*(1-m)*(E**(-20*m))
ENDIF
enddo
if(n>0) then
sum=0
do I=1,7
sum=sum+g(I)*x(I)
enddo
else
sum=0
do I=1,7
sum=sum+g(I)*x(I+7)
enddo
F=sum
endif
end FUNCTION F
END PROGRAM aero
错误是
--------------------Configuration: Text1 - Win32 Debug--------------------
Compiling Fortran...
G:\wx\Text1.f90
G:\wx\Text1.f90(28) : Error: The shapes of the array expressions do not conform. [W]
w=(m**0.25)*(1-m)*(E**(-20*m))
^
Error executing df.exe.

Text1.obj - 1 error(s), 0 warning(s)

Answer 1

首先，你的 w 没有定义，于是系统认为是一个单变量。

然后，你的 E 是一个 7 长度的数组。

w=(m**0.25)*(1-m)*(E**(-20*m))

这一句从逻辑上就通不过，因为对一个数组（E）进行操作，结果必然是一个数组，不能是一个单变量（W）

Answer 2

说句实话，我看不太懂，你看你把你的原题拿出来分享一下，帮你看看。

好像你定义的i就有问题。写成：integer::i。

然后后面f=sum应该放外面吧。

你那个w是做什么的？