求解指数等式 非线性方程

zhou1_20 发表于 2021-2-9 09:55
你参考下

把参数结果带回原方程，得不到对应的P3，P10耶，请问哪里出错了，谢谢！
theta = optim(c(1,1), model, control = list(reltol=1e-99))
theta
percent <- function(P) {
  L=theta$par[1];
  S=theta$par[2];
  P50=3455;
P50*(1+L*S*qt(P,1379))**(1/L)
}
percent(c(0.0228,0.03,0.1,0.25,0.5,0.75,0.9,0.97))
结果：
[1]        NaN        NaN        NaN   943.6284  3455.0000  5971.4147  8239.6427 10482.0722