area<-function(m,t,na.rm=TRUE) {
#na.rm não funciona
u <- matrix(nrow=length(m),ncol=3)
colnames(u)<-c("Method1","Method2","Method3")
for(i in 1:length(m)){
for (j in 1:3){
f1<-function(x, m, t) {m*x+t} # m e t variáveis locais
mt<-function(x) {b[j]*x+s[j]} # b e s variáveis globais
dif<-function(x, m,t) {abs(f1(x,m,t)-mt(x))}
int[i]<-integrate(dif, -1,1,m = m[i], t = t[i])$value
u[i,j]<-int[i]
}}
return(u)}