Olá a todos,

Estou com o seguinte problema, quero calcular a integral de 0 até HT( valor da altura total dde cada arvore presente na base de dados em anexo).

Rotininha:

rm(list=ls(all=T))

setwd("C:/Documents and Settings/Daniela/Desktop/Ajustes no R")

mdados<-read.table("INTEGRAL.txt", h=T, dec=",");mdados

# funçao de afilamento
fx<-function(hi) 0.91-(1.28*(hi/HT)^1)+(2.98*(hi/HT)^2)-(4.50*(hi/HT)^3)+(2.58*(hi/HT)^4)-(0.69*(hi/HT)^5)

for (i in 1:nrow(mdados))##numero total de arvores do arquivo
{
integral<-integrate(fx,lower=0,upper=mdados$HT)
volume<-with(mdados, ((pi*(mdados$DAP^2))/40000)*integral)
mdados<-data.frame(mdados, volume)
}

Mas não está dando certo.

A equação (polinomio do 5 grau) gerada é: 0.91-(1.28*(hi/HT)^1)+(2.98*(hi/HT)^2)-(4.50*(hi/HT)^3)+(2.58*(hi/HT)^4)-(0.69*(hi/HT)^5)

Objetivo, calcular somente a integral para cada árvore (ex: 33 observações), para nao precisar faze-lo manualmente ao introduzir árvore a árvore (HT).

O arquivo está no excel, a integral é de 0 até HT. O DAP servirá somente para o calculo do Volume.

Alguém pode me ajudar?

Meu computador é Windows.

Agradeço a todos, desde já!

Daniela