Calcula os coeficientes do polinômio interpolador de Newton

// Calcula o polinômio interpolador para os pontos tabulados

// abaixo, que representam a pressão atmosférica medida

// em função da altitude, e exibe o gráfico correspondente

alt=[0 10 20 30 40 50 60 70 80 90 100]

//pressao=[101.325 101.20494 101.08503 100.96525 100.84562 100.72613 100.60678 100.48758 100.36851 100.24958 100.1308]

pressao=[101.39124 101.27758 101.10488 101.01968 100.86883 100.74925 100.62843 100.57591 100.43376 100.28035 100.2241]

// calcula os coeficientes do polinômio interpolador de Newton

c_newton=polinomio_interp_newton(alt,pressao)

// cria um arranjo contendo valores dentro do intervalo de valores em alt

h=(0:1:100)'

// avalia o polinômio nos valores em h, guardando em p

for i=1:max(size(h))

    p(i)=avalia_polinomio(c_newton,alt,h(i))

end

// exibe os pontos tabulados (alt,pressao), em vermelho, e os pontos

// calculados (h,p), em azul

title('Polinômio interpolador de Newton')

xlabel('Altitude (km)')

ylabel('Pressão (kPa)')

plot(alt,pressao,'rx',h,p,'b--')

legend(['Pontos tabulados';'Polinômio interpolador']);