Question

fortran编程中遇到数组越界的问题： 运行到 y(i)=1/(1+x(i)*x(i))处程序停止，请问如何修改？

PROGRAM text1
implicit none
real x(7) !节点
real y(7) !函数值
real y1(7) !函数一阶导数
real q(7) !车比雪夫插值节点
real q1(7) !车比雪夫点处的函数值
real rt(100) !输出原函数函数值
real rt1(100) !输出拉格朗日插值函数值
real rt2(100) !输出三次样条插值函数值
real rt3(100) !输出车比雪夫点构成的拉格朗日插值函数值

real i,j,length
length=(5.0-(-5.0))/6

do i=1,7,1
x(i)=length*(i-1)-5.0
end do
do i=1,7,1
y(i)=1/(1+x(i)*x(i))
y1(i)=-2*(x(i))/((1+(x(i))**2)**2)
end do

Answer 1

real i,j,length
把这句换成：
real length
integer i, j
i 和毁键仿 j 两个变量应该用整型的亮猛。

试试看纤纤吧。

追问

试了一下，没有用。。