Интерполирование

#include <stdio.h>
#include <conio.h>
#include <math.h>
void main()
{
int i,k,N,M;
float x[100], y[100], z[10];
float DY, D2Y, L,L1, L2, q,H;
float dxi, dxi1, dxi2, P1, P2, h;
scanf("%d", &N);
H=sqrt(0.08/N);
i=0;
x[i]=N;
do
{
    y[i]=pow(N,3)*sin(x[i]/N);
    printf("%f%f\n", x[i], y[i]);
    x[i+1]=x[i]+H;
    i=i+1;
}
while (x[i]<=N+1);
M=i-1;
for (k=0; k<=9; k++)
{
    z[k]=N+0.05+(k+0.1);
    printf("%f\n", z[k]);
}
 
for (k=0; k<=9; k++)
{
    for (i = 0; i <= M; i++) {
         if ((x[i]<=z[k])&&  (z[k]<=x[i+1])) break;
        
        dxi=z[k]-x[i];
        dxi1=z[k]-x[i+1];
        dxi2=z[k]-x[i+2];
        L=pow(N,3)*sin(x[i]/N);
        P1=y[i]*dxi1/(-N) +y[i+1]*dxi/N;
        P2=y[i]*dxi1*dxi2/(2*N*N)+y[i+1]*dxi*dxi2/(-N*N)+y[i+2]*dxi*dxi1/(2*N*N);
        q=(z[k]-x[i])/H;
        DY=y[i+1]-y[i];
        D2Y=y[i+2]-y[i+1]-DY;
        L=pow(N,3)*sin(z[k]/N);
        L1=y[i]+q*DY;
        L2=y[i]+q*DY+q*(q-1)*D2Y/2;
        printf("%f%f%f%f%f%f\n", z[k], L, L1, L2, P1, P2);
    }
    
    printf("x ne prinadlejit in bound ");
}
 
getch();
}