>>653
#include <stdio.h>
#include <math.h>

#define RAD(x) ((x) * M_PI / 180.)

struct sin_table {
double x;
double sinx;
};

double lagrange_interpolation(double x, int n, struct sin_table st[])
{
int i, j;
double t, y = 0.;

for (i = 0; i < n; i++) {
t = st[i].sinx;
for (j = 0; j < n; j++)
if (i != j)
t *= (x - st[j].x) / (st[i].x - st[j].x);
y += t;
}
return y;
}

int main(void)
{ };