Pseudocódigo - Interpolação Polinomial

Variáveis: n, x, y, valor, resultado;
Onde -> n = Número de pontos
        x = Abscissa
        y = Ordenada
        valor = Valor à ser interpolado
        resultado = Armazena o resultado

Algoritmo:
r <- 0
para i <- 1 até n faça
c <- 1; d<-1
   para j <- 1 até n faça
   se i != j então
   c <- c*(valor - x[i])
   d <- d*(x[j] - x[i])
   fimse
   fimpara

resultado <- resultado + y[i]*c/d
fimpara
fimalgoritmo