2

I have a system of equations to solve for x1 and x2: 

2*erf(A*x2)*exp(B*x1^2+C*x2^2)-D*erf(x1)=0

2*erf(D*x2)*exp(E*x1^2+F*x2^2)-G*erf(x1)-erf(H*x2)*exp(I*x1^2)=0

where B and C are negative while rest of the constants are positive.

I tried to solve it in C using multiroot function from the GSL library. It seems that multiroot was not able to solve the system. I tried the same implementation of multiroot for non-linear simultaneous polynomial functions (no exponential part in it) and it worked well and gave correct results. Any suggestions on how to solve such a system? Or any other powerful library which can be used for such equations?

Lee Taylor
  • 7,761
  • 16
  • 33
  • 49
Deepak
  • 21
  • 3

1 Answers1

0

Look for a numerical method, for example the "Numerical recipes" book is higly regarded. Solving nonalgebraic equations is a vast, complex area. Systems of such equations doubly so.

vonbrand
  • 11,412
  • 8
  • 32
  • 52