求助，帮忙看看这个代码运行没有数据出来
Module mod_xyzImplicit none
real,protected::pi=acos(-1.0)
real,protected::L=1
real,protected::N=11
Contains
Subroutine Wik_(G,W,M,x,y,z,H)
Implicit none
real*8::G(1),M(1),W(1),x(11),y(11),z(11),H
Integer::i=1,k=0,a,b,c
CALL xyz_(x,y,z)
do
k=1+k
if (i==k)then
call G_Mii(G,W,M,x,y,z,a,b)
write(*,*)"G=",G
write(*,*)"W=",W
else
call G_Mik(G,W,M,x,y,z,a,b)
write(*,*)"z=",z(b)
write(*,*)"x=",x(a)
write(*,*)"G=",G
write(*,*)"M=",M
write(*,*)"H=",H
write(*,*)"W=",W
end if
end do
end Subroutine Wik_
Subroutinexyz_(x,y,z)
real*8::x(11),y(11),z(11)
Integer::a,b,c
Do a=1,N,1
x(a)=(L/2)*(1-cos(((a-1.0)/(N-1.0))*pi))
END Do
Do b=1,N,1
z(b)=(L/2)*(1-cos(((b-1.0)/(N-1.0))*pi))
END Do
Do c=1,N,1
y(c)=(L/2)*(1-cos(((c-1.0)/(N-1.0))*pi))
END Do
end Subroutine xyz_
Subroutine G_Mii(G,W,M,x,y,z,a,b)
Implicit none
real*8::G(1),M(1),W(1),x(11),y(11),z(11)
Integer::i,k,a,b
a=i
b=k
G=product(x(a)-y,y/=x(a))
M=product(z(b)-y,y/=z(b))
W=1/G
end Subroutine G_Mii
Subroutine G_Mik(G,W,M,x,y,z,a,b)
Implicit none
real*8::G(1),M(1),W(1),x(11),y(11),z(11)
Integer::i,j,k,H,N,a,b
a=i
b=k
G=product(x(a)-y,y/=x(a))
M=product(z(b)-y,y/=z(b))
H=1/(x(i)-x(k))
W=H*G/M
end Subroutine G_Mik
end Module mod_xyz
program main_mik
use mod_xyz
end program main_mik
在 program main_mil 中没有执行语句。
页:
[1]