elr10.c

在c环境下

  #include "stdio.h"
  #include "elr1.h"
  main()
  { 
	int i,j;
    double y[3],z[3][11],t,h,x;
    void elr1(t,y,n,h,k,z,f);
    void elr1f(t,y,n,d);
       y[0]=-1.0; y[1]=0.0; y[2]=1.0;
    t=0.0; h=0.01;
    elr1(t,y,3,h,11,z,elr1f);
    printf("\n");
    printf("*------------------------------------*\n");
    printf("     实验1：elrl法的实现    \n");
    printf("  完成人：王煜鹏    学号：05370010\n");
    printf("*------------------------------------*\n");
    printf("原偏微方程组为：\n");
    printf("y0`=y1,y0(0)=-1.0\n");
    printf("y1`=-y0,y1(0)=0.0\n");
    printf("y2`=-y2,y2(0)=1.0\n");
    printf("用elrl法的结果为(h=0.01)：\n");
    for (i=0; i<=10; i++)
      { x=i*h;
        printf("t=%5.2f\n",x);
        for (j=0; j<=2; j++)
          printf("y(%d)=%e  ",j,z[j][i]);
        printf("\n");
      }
    printf("\n");
  }