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}