1

在我的代码中，我在 while 循环之前定义了一个步长“h”。当我尝试在循环中使用它时，它似乎会自行改变。如果我在循环内定义它似乎没问题，但我得到的数据似乎不正确，所以我猜这个问题可能是相关的。

即使我在这个位置打印它（参见代码中的 printf），输出也是 3 个值，我不知道为什么。如果您看到任何其他不相关但似乎错误的内容，请告诉我，正如我所说，我得到了意想不到的值（可能只是我的公式）。

int main()
{
    FILE *f1;
    f1 = fopen("Question2 part 3 solution.txt", "w");
   double r0=0.05;
   double dr0=-a*a*r0/sqrt(1+pow(a*a*r0,2)),h=0.01;
   double k[4][3],x[]={dr0,z0,T0,r0},x1[]={0,0,0,0}, s=0;
   int i,j;

   while(s<=1)
    {
        //Runge-Kutta
        for (j=0;j<4;j++)
        {
            for (i=0;i<4;i++)
            {
                if (j==0){k[i][0]=h*System(i,x[0],x[1],x[2],x[3]);}
                if (j==1){k[i][1]=h*System(i,x[0]+k[0][0]/2.0,x[1]+k[1][0]/2.0,x[2]+k[2][0]/2.0,x[3]+k[3][0]/2.0);}
                if (j==2){k[i][2]=h*System(i,x[0]+k[0][1]/2.0,x[1]+k[1][1]/2.0,x[2]+k[2][1]/2.0,x[3]+k[3][1]/2.0);}
                if (j==3){k[i][3]=h*System(i,x[0]+k[0][2],x[1]+k[1][2],x[2]+k[2][2],x[3]+k[3][2]);}
            }
        }

        for (i=0;i<4;i++)
        {
            x[i]=x[i]+(k[i][0]+2.0*k[i][1]+2.0*k[i][2]+k[i][3])/6.0;
        }
printf("%8.3lf",h);
        s+=h;
    }

   fclose(f1);
    return(0);
}


double System(int i,double dr, double z, double T, double r)
    {
    //printf("%e\t%e\t%e\t%e\n",dr,z,T,r);
    if (T==T0 && z==z0 && i==0) {return  (-a*a*dr)*pow(1-dr*dr,3/2)/2.0;}
    if (i==0 && T!=0){return  (-a*a*r*(1-dr*dr)-dr*sqrt(1-dr*dr))/T;}
    if (i==1){return  (-sqrt(1-dr*dr));}
    if (i==2){return  (-a*a*r*dr+sqrt(1-dr*dr));}
    if (i==3){return  (dr);}
    //if (i==3){return  (-m2*l1*l2*B*theta1dt*theta2dt*sin(theta2-theta1)-l2*m2*g*B*sin(theta2));}
    }

提前致谢！

4

2 回答 2

3

见声明k：

double k[4][3]

然后看到这个说法

k[i][3]=...

在这里，您超出了数组的边界写入，导致未定义的行为。

于 2013-11-10T10:27:07.707 回答
3

您超出了内存，变量k定义为，但是double k[4][3]您正在更新k[i][3]j==3

于 2013-11-10T10:27:23.360 回答