Calcula o polinômio cúbico de ajuste para os pontos tabulados

// Calcula o polinômio cúbico de ajuste 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.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 cúbico de ajuste

[a,rmq]=num_ajuste_polinomial(3,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 p3

p3=num_avalia_polinomial(h,a)

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

// calculados (h,p3), em magenta

title('Polinômio cúbico de ajuste')

xlabel('Altitude (km)')

ylabel('Pressão (kPa)')

plot(alt,pressao,'rx',h,p3,'m--')

legend(['Pontos tabulados';'p3(h)']);