#include <stdio.h>
#include <math.h>
#define h 0.01
#define N 10000
double f(double x[],double dx,double y[],double dy,int i)
{
double xi,yi;
xi=x[i]+dx;
yi=y[i]+dy;
return (10*(yi-xi));
}
double p(double x[],double dx,double y[],double dy,double z[],double dz,int i)
{
double xi,yi,zi;
xi=x[i]+dx;
yi=y[i]+dy;
zi=z[i]+dz;
return (28*xi-xi*zi-yi);
}
double q(double x[],double dx,double y[],double dy,double z[],double dz,int i)
 {
  double xi,yi,zi;
  xi=x[i]+dx;
  yi=y[i]+dy;
  zi=z[i]+dz;
  return (xi*yi-8/3*zi);
  }
  
 main()
{
  double Kx[3][4],Ky[3][4],Kz[3][4];
  double x[3]={2,4,2},y[3]={5,7,4},z[3]={3,9,5};
  int i,j;      
  FILE *fp;
  fp=fopen("l_orenzrr.txt","w");
  for(j=1;j<=N;j++)
  {
    for(i=0;i<3;++i)
    {
     Kx[i][0]=f(x,0,y,0,i);
     Ky[i][0]=p(x,0,y,0,z,0,i);
     Kz[i][0]=q(x,0,y,0,z,0,i);
     
     Kx[i][1]=f(x,h/2*Kx[i][0],y,h/2*Ky[i][0],i);
     Ky[i][1]=p(x,h/2*Kx[i][0],y,h/2*Ky[i][0],z,h/2*Kz[i][0],i);
     Kz[i][1]=q(x,h/2*Kx[i][0],y,h/2*Ky[i][0],z,h/2*Kz[i][0],i);
   
     Kx[i][2]=f(x,h/2*Kx[i][1],y,h/2*Ky[i][1],i);
     Ky[i][2]=p(x,h/2*Kx[i][1],y,h/2*Ky[i][1],z,h/2*Kz[i][1],i);
     Kz[i][2]=q(x,h/2*Kx[i][1],y,h/2*Ky[i][1],z,h/2*Kz[i][1],i);
     
     Kx[i][3]=f(x,h*Kx[i][2],y,h*Ky[i][2],i);   
     Ky[i][3]=p(x,h*Kx[i][2],y,h*Ky[i][2],z,h*Kz[i][2],i);
     Kz[i][3]=q(x,h*Kx[i][2],y,h*Ky[i][2],z,h*Kz[i][2],i);
     
     x[i]=x[i]+(Kx[i][0]+2*Kx[i][1]+2*Kx[i][2]+Kx[i][3])/6*h;
     y[i]=y[i]+(Ky[i][0]+2*Ky[i][1]+2*Ky[i][2]+Ky[i][3])/6*h;
     z[i]=z[i]+(Kz[i][0]+2*Kz[i][1]+2*Kz[i][2]+Kz[i][3])/6*h;
     
     printf("  %8f ,  %8f ,  %8f",x[i],y[i],z[i]);
     fprintf(fp,"%8f , %8f, %8f \n",x[i],y[i],z[i]);
     printf("\n");   
   }
  }
}

就是平滑连接

这个程序出的图是这种的   点对咯但是连接有问题

求指点啊

???????