hermit 曲线算法(待证)

void hermit(){

        double x[500],  y[500];

        int k=0;

        for(double u=0.0;u<1.0;u+=0.002){

              x[k]=-10*u*u*u+15*u*u+u+1;

             y[k]=4*u*u*u+4*u*u+2*u+6;

             k++;

       } 

       for(int j=0;j<k;j++){

             x[j]*=50;

            y[j]*=20;

        } 

        for(int i=0;i<k-1;i++){

             pDC->MoveTo(x[i],y[k]);

             pDC->LineTo(x[i+1],y[i+1]);

        }

}