Solve numerical system of nonlinear equations? I need to solve a nonlinear system of equations that looks like this
$$\frac{1}{n_4}(\frac{c_4}{b_4R}-\ln(1-x_1-x_2-x_3))=\frac{1}{n_3}(\frac{c_3}{b_3R}-\ln(x_3))=\frac{1}{n_2}(\frac{c_2}{b_2R}-\ln(x_2))=\frac{1}{n_1}(\frac{c_1}{b_1R}-\ln(x_1))
$$
The numerical equivalent of what i'm trying to solve right now is
$$\frac{4100}{4100/297-1.987\ln(x_1)}=\frac{6200}{6200/303-1.987\ln(x_2)}=\frac{7000}{7000/327.5-1.987\ln(x_3)}=\frac{4100}{4100/404-1.987\ln(1-x_1-x_2-x_3)}
$$
What algorithm should i use to solve this? If anyone can solve this, could you post how you did it? I tried using a multivariable newton's algorithm but my jacobian exploded
 A: Hint
I think that you could reduce your problem to a single equation of a single variable.
Suppose that you define $$A=\frac{1}{n_4}(\frac{c_4}{b_4R}-ln(1-x_1-x_2-x_3))$$ So, you have $$A=\frac{1}{n_3}(\frac{c_3}{b_3R}-ln(x_3))=\frac{1}{n_2}(\frac{c_2}{b_2R}-ln(x_2))=\frac{1}{n_1}(\frac{c_1}{b_1R}-ln(x_1))$$ So $x_1,x_2,x_3$ can be expressed as a function of $A$ $$x_i(A)=\exp({\frac{c_i}{b_i~ R}-A~ n_i})$$ and you end with a single nonlinear equation $A=f(A)$ where only $A$ appears. 
$$A=\frac{1}{n_4}(\frac{c_4}{b_4R}-ln(1-x_1(A)-x_2(A)-x_3(A)))$$ 
You could notice that your equation in $A$ also write
$$F(A)=\sum _{i=1}^{i=4} \exp({\frac{c_i}{b_i~ R}-A~ n_i})=1$$ which, more than likely, would much better conditionned if written
$$G(A)=\log \Big(\sum _{i=1}^{i=4} \exp({\frac{c_i}{b_i~ R}-A~ n_i}) \Big )=0$$ Assuming  all $n_i>0$, this shows that an upper bound of $A$ corresponds to the maximum,over $i$, of the $(\frac {c_i}{n_i~b_i~R})$ terms. This gives a starting value for unknown $A$.
Solve it for $A$ (Newton method seems to be very appropriate) and go backwards to get $x_1,x_2,x_3$. 
In such a manner, you could notice that you could have as many equations as you could wish and apply the same procedure.
Using the data from your problem, we can easily find that the solution is $A=271.069$ from which we get $x_1=0.514479,x_2=0.297290,x_3=0.106527$ which are your answers.
