26 Jan
2012
26 Jan
'12
19:13
Com ?nls.lm você resolveria seu problema. No entanto, segue um exemplo: f <-function(luz,a,b){ (((0.01940*luz+a+0.60375)-(((0.01940*luz+a+0.60375)^2)-4*0.01940*luz*b*(a+0.60375))^0.5)/2*b)-0.60375 } preditos <- do.call('f',c(list(luz=luz),lin1$par)) plot(luz,Pn,ylim=c(-1,5),type='l') par(new=TRUE) plot(luz,preditos,ylim=c(-1,5),col='red',xlab='',ylab='',type='l') (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}