init.h = c(600, 700, 800, 950, 1100, 1300, 1500)
h.d = c(253, 337, 395, 451, 495, 534, 573)
plot(h.d~init.h)
mod1 <- lm(h.d~init.h)
mod2 <- update(mod1,.~. +I(init.h^2))
mod3<- update(mod2, .~.+I(init.h^3))
anova(mod1,mod2)
anova(mod2, mod3)

abline(mod1, col="red")
abline(mod2, col="blue")

cf.m2 <- coef(mod2)
cf.m3<- coef(mod3)

curve(cf.m2[1]+cf.m2[2]*x+cf.m2[3]*x^2, add=T, lty=2)

curve(cf.m3[1]+cf.m3[2]*x+cf.m3[3]*x^2+cf.m3[4]*x^3,  add=T, lty=2, col="green")

summary(mod2)
summary(mod3)
#O polinomio de terceiro grau é melhor modelo para descrever o experimento de Galileu, com a linha passando mais próximas aos pontos (visualizacao do grafica) e atraves da anova dos modelos.