Question

在用fortran运行时出现forrt1:severe<161>:Program Exception - array bounds exceeded是怎么回事啊？

这是我写的程序：
program main
real::x
integer::n
print*,'请输入legendre多项式的阶n='
read*,n
print*,'请输入参数x='
read*,x
call my(n,x)
end

subroutine my(n,x)
dimension p(0:n)
if(x==1) then
do i=1,n
p(i)=x**i
enddo
print*,p(i)
else
p(0)=1.0
p(1)=x
do i=2,n
p(i)=(2.0-1.0/dble(i))*x*p(i-1)-(1.0-1.0/dble(i))*p(i-2)
enddo
print*,p(i)
endif
end

Answer 1

问题出在这：
print*,p(i)

假如n是5
上面的循环后
do i=2,n
p(i)=(2.0-1.0/dble(i))*x*p(i-1)-(1.0-1.0/dble(i))*p(i-2)
enddo

i为6
超过P的范围n =5.