Câmbio para Gretl

dd=[1 3 4 5 12 16 17 23 25 29 32 38 40 43 46 47 53 54 57 60 63]

y=[2.2456 2.2362 2.2434 2.2396 2.3351 2.3286 2.3572 2.4070 2.4299 2.4557 2.4918 2.3860 2.4236 2.3927 2.4645 2.4325 2.5137 2.4570 2.5229 2.4079 2.5005]

plot(dd,y,'rx')

x=1:1:93

[a0,a1,rmq]=num_ajuste_exponencial(dd,y)

gexp=num_avalia_exponencial(x,a0,a1)

mprintf("Ajuste exponencial: a0=%g a1=%g rmq=%g\n  valor estimado=%g\n",a0,a1,rmq,gexp(93))

plot(x,gexp,'b-')

[a0,a1,rmq]=num_ajuste_potencial(dd,y)

gpot=num_avalia_potencial(x,a0,a1)

mprintf("Ajuste potencial: a0=%g a1=%g rmq=%g\n  valor estimado=%g\n",a0,a1,rmq,gpot(93))

plot(x,gpot,'g-')

[a,rmq]=num_ajuste_polinomial(3,dd,y)

gcub=num_avalia_polinomial(x,a)

mprintf("Ajuste polinomial: a0=%g a1=%g a2=%g a3=%g rmq=%g\n  valor estimado=%g\n",a(1),a(2),a(3),a(4),rmq,gcub(93))

plot(x,gcub,'c-')

[a,rmq]=num_ajuste_polinomial(2,dd,y)

gquad=num_avalia_polinomial(x,a)

mprintf("Ajuste polinomial: a0=%g a1=%g a2=%g rmq=%g\n  valor estimado=%g\n",a(1),a(2),a(3),rmq,gquad(93))

plot(x,gquad,'k-')