Method: Perform an integral, as a function of E, which outputs Current for each Voltage value used. This is repeated for an array of v_values. The equation can be found below.
Although the limits in this equation range from
inf, the limits must be restricted so that (E+eV)^2-\Delta^2>0 and E^2-\Delta^2>0, to avoid poles. (\Delta_1 = \Delta_2). Therefore there are currently two integrals, with limits from
However, I keep returning a
math range error although I believe I have excluded the troublesome E values by using the limits stated above.
Apologies for the vagueness of this question. But, can anybody see obvious mistakes or code misuse?