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 ex01h=l/(n-1) = 0
注意 l 和 n 是整数.
另外你的编程风格非常不好,建议这样写.