为啥我写的fortran结果算出来都是零,程序如下.programex01useimslparameter(n=50,l=1,p=1,h=l/(n-1))parameter(pi=3.1415926535)dimensiona(n-2,n),b(n-2,1),y(n,1),q(n,1)open(1,file='naodu.dat')y(:,:)=0doi=1,nq(i,1)=p*sin((i-1)*h*pi/l)end

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

h=l/(n-1) = 0

注意 l 和 n 是整数.

注意 l 和 n 是整数.

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

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