Question

为啥我写的fortran结果算出来都是零，程序如下。

program ex01
use imsl
parameter(n=50,l=1,p=1,h=l/(n-1))
parameter(pi=3.1415926535)
dimension a(n-2,n),b(n-2,1),y(n,1),q(n,1)

!open(1,file='naodu.dat')

y(:,:)=0

do i=1,n
q(i,1)=p*sin((i-1)*h*pi/l)
end do

do i=1,n
write(*,*) q(i,1)
end do

end

Answer 1

program ex01
  implicit none
  integer , parameter :: n=50 , l=1 , p=1 
  real , parameter :: h= 1.0* l / (n-1)
  Real , parameter :: pi=3.1415926535
  Real :: a(n-2,n),b(n-2,1),y(n,1),q(n,1)
  integer :: i
  y(:,:)=0
  do i=1,n
    q(i,1)=p*sin((i-1)*h*pi/l)
  end do
  do i=1,n
    write(*,*) q(i,1)
  end do
end program ex01

h=l/(n-1) = 0

注意 l 和 n 是整数。

另外你的编程风格非常不好，建议这样写。