Model I-V.

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`

to `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 `-inf`

to `-gap-e*v`

and `gap`

to `inf`

.

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?