#include #include /* hermite.c : 1993 3 13 by Oguni */ #define M 21 #define N 11 main(){ double x[N], y[N], dy[N], xnew[M], ynew[M], yy[M]; int i, m=21, n=11; for (i=0; ix[i] && xnew[j]<=x[i+1]){ a=xnew[j]-x[i]; b=x[i+1]-x[i]; c=xnew[j]-x[i+1]; f0=a*a*(2.0*(-c)+b)/(b*b*b); f1=a*a*c/(b*b); f2=c*c*(2.0*a+b)/(b*b*b); f3=c*c*a/(b*b); ynew[j]=y[i+1]*f0+dy[i+1]*f1+y[i]*f2+dy[i]*f3; } } } return; }