Saturday, 17 August 2013

NSolve cannot solve system of equations in Mathematica

NSolve cannot solve system of equations in Mathematica

pI have the following two equations:/p precodeIeq1=-1 - (0.763546 I1
(1/(-(-8 + 10/I2)^2 + (1 - 9/I1 + 10/I2)^2) - ( 18 (1 - 9/I1 + 10/I2))/(
I1 (-(-8 + 10/I2)^2 + (1 - 9/I1 + 10/I2)^2)^2)) X^1.3)/((1.36- I1)
(I1/(-(-8 + 10/I2)^2 + (1 - 9/I1 + 10/I2)^2))^2.3) + ( 0.587343
X^1.3)/((1.36- I1) (I1/(-(-8 + 10/I2)^2 + (1 - 9/I1 + 10/I2)^2))^1.3) + (
0.587343 I1 X^1.3)/((1.36- I1)^2 (I1/(-(-8 + 10/I2)^2 + (1 - 9/I1 +
10/I2)^2))^1.3); Ieq2=1 + (0.763546 I1^2 ((20 (-8 + 10/I2))/I2^2 - (20 (1
- 9/I1 + 10/I2))/ I2^2) (-8 + (1 - 9/I1 + 10/I4)^2) X^1.3)/((1.36- I1)
(I1/(-(-8 + 10/I2)^2 + (1 - 9/I1 + 10/ I2)^2))^2.3 (-(-8 + 10/I2)^2 + (1 -
9/I1 + 10/I2)^2)^3) - ( 0.587343 I1 ((20 (-8 + 10/I2))/I2^2 - (20 (1 -
9/I1 + 10/I2))/ I2^2) (-8 + (1 - 9/I1 + 10/I4)^2) X^1.3)/((1.36- I1)
(I1/(-(-8 + 10/I2)^2 + (1 - 9/I1 + 10/ I2)^2))^1.3 (-(-8 + 10/I2)^2 + (1 -
9/I1 + 10/I2)^2)^2); /code/pre pand then I use the NSolve function:/p
precode NSolve[{Ieq1 == 0, Ieq2 == 0}, {I1, I2}] /code/pre pbecause I want
to solve each one as a function of X. But unfortunately it just gets stuck
while running. Do you have any suggestions to overcome this problem?/p

No comments:

Post a Comment