Си. 2 интерполяционная формула Ньютона

[Моркковь!)]

Здравствуйте. Долго мучалась с первой формулой.
В общем, работает. Теперь думаю, как программу в обратную сторону повернуть для второй формулы Ньютона. Если есть какие-то идеи, прошу подсказать

[KobolD]

Код первой части в студию, и желательно теорию того как должна вторая часть работать.

[Моркковь!)]

#include <stdlib.h>
#include <conio.h>
#include <stdio.h>
void inst_r(float **r ,float *y, int m);
void tapping(float *x, float *y, float xtest);
int *In(int m);
main ()
{
int m; float *x,*y,xtest, newt;
newt=In(x, y, xtest);
tappint(newt, r);

instr_r(r, m);
}



int *In( int m)
{
int l; float **r;

printf ("\n input n \n");
scanf(" %d", &m);
r=(float**)malloc(m*sizeof(float*)) ;
for(l=0; l<m; l++)
r[l]=(float*)malloc(m*sizeof(float));
x=(float*)malloc(m*sizeof(float));
y=(float*)malloc(m*sizeof(float));
printf("\n input x`s\n ");
for(l=0; l<m; l++)
scanf( "%f", &x[l]);
printf("\n input y`s\n ");
for(l=0; l<m; l++)
scanf( "%f", &y[l]);
printf("\n input X\n ");
scanf( "%f", &xtest);
// tapping(x, y, xtest, m);
printf("\n y=%f",newt);
return newt;
getch();

}
void tapping(float *x, float *y, float xtest, int m)
{
float **r, h, newt, com; int l;
h=x[1]-x[0];
//inst_r (r, y,m);
newt=r[m-1][0]; com=1;
for (l=m-1;l>0;l--)
{
com=com*(xtest-x[l])/h/(m-l);
newt=newt+r[l-1][m-l]*com;
}

}

void inst_r(float **r, float *y, int m)
{
int i,j;
for(i=0;i<m;i++)
r[i][0]=y[i];
for (j=1;j<m;j++)
for (i=0;i<m-j;i++)
r[i][j]=r[i+1][j-1]-r[i][j-1];
}



это код второй части. Осталось ошибки исправить. проблемы только с передаваемыми и возвращаемыми параметрами))

[D*2]

Сообщение от Моркковь!)

Здравствуйте. Долго мучалась с первой формулой.
В общем, работает. Теперь думаю, как программу в обратную сторону повернуть для второй формулы Ньютона. Если есть какие-то идеи, прошу подсказать

Здравствуйте!!!
Скоро сдавать прогу, а я еще ничего не сделал. Помогите пожалуйста
Скинте пожалуйста код первой и второй интерполяционной формулы Ньютона
заранее спасибо