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