main()
{
    int i,j;
    void rktlf(double,double[],int,double[]);
    double t,h,y[2],z[2][101];
    FILE * fp,*fp2;   
  fp=fopen("e1.txt","w");
    y[0]=0.0;y[1]=0.0;
    t=0.0;h=0.1;
    rkt1(t,y,2,h,11,z,rkt1f);
    printf ("\n");
    for (i=0;i<=10;i++)
    {
        t=i*h;
        printf("t=%5.2f\n",t);
        for (j=0;j<=0;j++)
        printf ("y(%d)=%0.7f\t",i,z[j][i]);
        printf("\n");
        fprintf(fp,"%10f   %10f\n",t,z[j][i]);

    }

    printf("finish");
    fclose(fp);
    getch();
}为什么输出的Ｚ【Ｊ】【Ｉ】和打印出来的Ｚ【Ｊ】【Ｉ】会不一样呢？

头文件？？
main()
{
    int i,j;
    void rktlf(double,double[],int,double[]);函数体呢？？
    double t,h,y[2],z[2][101];
    FILE * fp,*fp2;   
  fp=fopen("e1.txt","w");
    y[0]=0.0;y[1]=0.0;
    t=0.0;h=0.1;
    rkt1(t,y,2,h,11,z,rkt1f);
    printf ("\n");
    for (i=0;i<=10;i++)
    {
        t=i*h;
        printf("t=%5.2f\n",t);
        for (j=0;j<=0;j++)这行代码  有意义？？
        printf ("y(%d)=%0.7f\t",i,z[j][i]);
        printf("\n");
        fprintf(fp,"%10f   %10f\n",t,z[j][i]);你这个j是经过j++的值

    }

    printf("finish");
    fclose(fp);？？？
    getch();
}

[ 本帖最后由 A13433758072 于 2011-10-5 08:46 编辑 ]

#include "stdlib.h"
#include "math.h"
#include "stdio.h"
void rkt1(t,y,n,h,k,z,f)
void (*f)();
int n,k;
double t,h,y[],z[];
{
    int i,j,l;
    double a[4],tt,*b,*d;
    b=malloc(n*sizeof(double));
    d=malloc(n*sizeof(double));
    a[0]=h/2.0;a[1]=a[0];
    a[2]=h;a[3]=h;
    for (i=0;i<=n-1;i++) z[i*k]=y[i];
    for (l=1;l<=k-1;l++)
      {
       (*f)(t,y,n,d);
       for(i=0;i<=n-1;i++) b[i]=y[i];
       for(j=0;j<=2;j++)
        {
            for(i=0;i<=n-1;i++)
           {
            y[i]=z[i*k+l-1]+a[j]*d[i];
            b[i]=b[i]+a[j+1]*d[i]/3.0;
           }
          tt=t+a[j];
         (*f)(tt,y,n,d) ;
        }
      for (i=0;i<=n-1;i++)
        z[i*k+l]=y[i];
      t=t+h;
     }
   free(b);free(d);
}
void rkt1f(t,y,n,d)
int n;
double t,y[],d[];
{
    t=t;n=n;
    d[0]=y[1];
    d[1]=-y[0]+2*y[1]+t*exp(t)-t;
}
main()
{
    int i,j;
    void rktlf(double,double[],int,double[]);
    double t,h,y[2],z[2][101];
    FILE * fp,*fp2;   
  fp=fopen("e1.txt","w");
    y[0]=0.0;y[1]=0.0;
    t=0.0;h=0.1;
    rkt1(t,y,2,h,11,z,rkt1f);
    printf ("\n");
    for (i=0;i<=10;i++)
    {
        t=i*h;
        printf("t=%5.2f\n",t);
        j=0;
        printf ("y(%d)=%0.7f\t",i,z[j][i]);
        printf("\n");
        fprintf(fp,"%10f   %10f\n",t,z[j][i]);

    }

    printf("finish");
    fclose(fp);
    getch();
}麻烦再帮看下

我看晕了。。。

呵呵 ＳＯＲＲＹ

你写的什么啊，别把代码随便贴出来让大家改好不好？？？拜托，下面大括号什么意思？？？？

double t,h,y[],z[];
{

  好多的代码哦  我有点晕