I am new to mathematica, it is blocked at the line "nen[te_, n2_]= ne/.Nsolve[nq1ne(te, n2] == 0, ne]".
Can I debug this line? Do I need to use some evaluate functions of mathematica?

Thanks for the reply. Now I put the code in https://www.dropbox.com/sh/ulu3sprw6spycsu/AADL9ExHzs1_wCt5tCNQYZbLa?dl=0. "1Plasma chmaber.nb" runs firstly, then "5Rate Constants.nb", then "61Patical Balance Equations.nb"´and run "62Patical Balance Equations.nb" at last. The Problem ist that it is bloc...

The following code run more than one day without a result. It was blocked at step3. How to improve the code and get a result in a limited peroid? Is my way Okay? Thanks in advance. Step1: q2n1[te,ne,n2_]=n1/.Solve[q2[te,ne,n2,n1]==0,n1] q2n1[te,ne,n2][[3]] Output: {0. -(8.248081032+1.428611033 I) te...