[Fortran] 纯文本查看 复制代码 Subroutine Wii_()
Implicit none
real*8::x(11),y(11),pi=3.141592653589793,L
Integer::N,i,j
L=1
N=11
Do i=1,N,1
x(i)=(L/2)*(1-cos(((i-1)/(N-1))*pi))
Do j=1,N,1
y(j)=(L/2)*(1-cos(((j-1)/(N-1))*pi))
END Do
END Do
end Subroutine Wii_
program main_mik
Implicit real*8(a-z)
real*8::G(1),W(1),x(11),y(11)
Integer::k
k=1
call Wii_()
G(k)=product(x(i)-y(j),y(j)/=x(i))
W(k)=1/G(k)
write(*,*)"W(k)=W(k)"
end program main_mik
|