Hi all,

I am trying to solve the following system of two equations in order two obtain a value for the two parameters I am interested in.

The two equations are the Real and the Imaginary part of the same function which I called pol, and both of the equations should be equal to zero.

pol[M_, T_] =0.801526 - (M - (I*T)/2)^2 -0.198833 I (M - (I*T)/2)^2 ((1/((M - (I*T)/2)^6))0.0529229 ((-2.10673 + (M - (I*T)/2)^2) (-0.215392 + (M - (I*T)/2)^2))^(3/2)HeavisideTheta[-2.10673+(M-(I*T)/2)^2]^3+(1/((M - (I*T)/2)^6))0.947077 ((-1.08478 + (M - (I*T)/2)^2) (-0.00293504 + (M - (I*T)/2)^2))^(3/2)HeavisideTheta[-1.08478 + (M - (I*T)/2)^2]^3+(1/((M - (I*T)/2)^6))((-0.406049+(M-(I*T)/2)^2) (-0.128221 + (M - (I*T)/2)^2))^(3/2)HeavisideTheta[-0.406049 + (M - (I*T)/2)^2]^3)0.198833 I (M - (I*T)/2)^2 ((1/((M - (I*T)/2)^6))0.0529229 ((-2.10673 + (M - (I*T)/2)^2) (-0.215392 + (M - (I*T)/2)^2))^(3/2)HeavisideTheta[-2.10673 + (M - (I*T)/2)^2]^3 + (1/((M - (I*T)/2)^6))0.947077 ((-1.08478 + (M - (I*T)/2)^2) (-0.00293504 + (M - (I*T)/2)^2))^(3/2)HeavisideTheta[-1.08478 + (M - (I*T)/2)^2]^3 + (1/((M - (I*T)/2)^6))((-0.406049 + (M - (I*T)/2)^2) (-0.128221 + (M - (I*T)/2)^2))^(3/2)HeavisideTheta[-0.406049 + (M - (I*T)/2)^2]^3)

NSolve[{ComplexExpand[Re[POL[M,T]]]==0,ComplexExpand[Im[POL[M,T]]]==0},{M,T}]

Both equations depends on the parameters M and T whose values is what I am interested to obtain by solving the system of equations through the NSolve mathematica command.

After waiting lot of time trying to get an answer from mathematica I abort the computation because I think I am doing something wrong because to solve a system of two equations can not take such a long time.

Notice that some HeavisideTheta of complex numbers appears and maybe I have to treat it in a specific manner I don't know by the way.

Any help to solve it will be very acknowledged.

All my best