求解R中的非线性方程

时间:2018-02-14 21:15:44

标签: r nonlinear-functions

我需要为P和K解决以下函数:

enter image description here

我想找到满足等式的所有(或少数)P和K.

我尝试过使用R的nleqslv软件包,但出了点问题。

MPK<-function(X){
  Y=numeric(2)
  Y[1] = 4.34783*((.3*(X[1]^.23)+.7*(X[2]^.23))^3.34783)*0.069*(X[1]^-.77)
  Y[2] = 0.3*((1-X[2])/(1-X[1]))^.7
  Y
} 

#solve for K, P
xstart = c(.5,.5)
nleqslv(x = xstart,fn = MPK)

我得到的是以下内容:

$x
[1] 1.214578e+10 1.006411e+00

$fvec
[1] 5.531138e-03 7.636165e-10

$termcd
[1] 5

$message
[1] "Jacobian is too ill-conditioned (1/condition=8.9e-013) (see 
allowSingular option)"

$scalex
[1] 1 1

$nfcnt
[1] 142

$njcnt
[1] 7

$iter
[1] 70

如何指定我需要Y 1等于Y [2]的解决方案?

0 个答案:

没有答案