main.c

一类kdv方程数值解法,kdv方程的数值解法

  #include "stdio.h"
  #include "part1.c"
  #include "math.h"
  #include "part2.c"
  main()
  { int i,j;
    double y[50],z[50][314],t,h,x;
    double n1,n2,n3,n4;
    printf("input n1 n2 n3 n4:");
    printf("\n");
    printf("n1 n2 n3 n4  aspectly stand for  a  b  arphar  betar");
    scanf("%lf%lf%lf%lf",&n1,&n2,&n3,&n4);
    t=0.00;h=0.02;
    for(i=0;i<=313;i++)
    y[i]=cos(n1*i*h);
    part1(t,y,315,h,50,z);
    printf("\n"); 
    for (i=0; i<=50; i++)
      { x=i*h;
        printf("t=%5.2f\n",x);
        for (j=0; j<=313; j++)
          printf("y(%d)=%e  ",j,z[j][i]);
        printf("\n");
        j=j-4;
      }
    printf("\n");
   part2();
  }