10wity0.c

来自「常用算法的C程序。主要内容包括多项式的计算、复数运算、随机数的产生、矩阵运算、矩」· C语言 代码 · 共 28 行

  #include "stdio.h"
  #include "10wity.c"
  main()
  { int i,j;
    void  wityf(double,double [],int,double []);
    double t,h,x,y[3],z[3][11];
    y[0]=-1.0; y[1]=0.0; y[2]=1.0;
    t=0.0; h=0.1;
    wity(t,y,3,h,11,z,wityf);
    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");
      }
  }

  void wityf(t,y,n,d)
  int n;
  double t,y[],d[];
  { t=t; n=n;
    d[0]=y[1]; d[1]=-y[0]; d[2]=-y[2];
    return;
  }