part2.c

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

  #include"stdio.h"
  #include"math.h"
  #include"part1.c"
  void part2()
  { 
   int i,j; 
   double cs[50];
   for(i=1;i<=50;i++)
   cs[i]=cos(n1*157+n2*i*0.02);
   extern void part1();
   printf("k","        ");
   printf("(x,t)","     ");
   printf("numerical value","     ");
   printf("accurate value","      ");
   printf("|NV-AV|");
   gelr1();
   for(i=1;i<=floor(k/5);i++)
    {
     printf("%-9d",i*5);
     printf("(3.14,%5.2f)",h*i*5);
     printf("%-20.6f",z[157][i*5]);
     printf("%-20.6f",cs[i*5]);
     printf("%-5.4e\n",fabs(z[157][i*5]-cs[i*5]));
     } 
     printf("%-9d",k);
     printf("(3.14,%5.2f)",h*k);
     printf("%-20.6f",z[157][k]);
     printf("%-20.6f",cs[k]);
     printf("%-5.4e\n",fabs(z[157][k]-cs[k]));
   }