[R-br] ajuda com integral

[DANIELA CUNHA DA SE]

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120413/938bb316/attachment-0001.html>
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: INTEGRAL.xlsx
Tipo: application/vnd.openxmlformats-officedocument.spreadsheetml.sheet
Tamanho: 10060 bytes
Descrição: não disponível
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120413/938bb316/attachment-0001.xlsx>