lagrange

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define comprimento 4 //comprimento dos vetores

int main() {
    float x[comprimento]={-1,0,1,2}; //preencha este vetor para atribuir valor em x
    float y[comprimento]={1,3,1,1}; //preencha este valor para atribuir valor em y
    int i, j;
    float p=0, s=0;
    float u=1.5; //ponto que deverá ser interpolado
    for(i=0; i<comprimento; i++){
        p = y[0];
        for(j=0; j<comprimento; j++){
            if(i!=j){
                p=p*((u-x[j])/(x[i]-x[j]));
            }
        }
        s=s+p;
    }
    printf("s= %f", s);
    return 0;
}