secs1d pin diode help

[tkroezen]

Hi,

I am having some trouble trying to simulate a PIN diode in secs1d. My
approach was to alter the PN diode code included as the example in the
documentation file, but unfortunately this hasn't produced any results yet.

The only changes I made to the code are listed below.

change lines 63 until 68 to:
 /Na = 1e23 * (x <= x_pi_jun);
 N_0 = 1e16 .* (x > x_pi_jun) .* (x <= x_in_jun);
 Nd = 1e23 * (x > x_in_jun);
 D  = Nd - Na + N_0;  
/
change lines 70 until 83 to:
/V_p = -1;
V_n =  0;
 
 Fp = V_p * (x <= x_pi_jun);
 Fn = Fp;
 
 p = Na .* (x <= x_pi_jun) + ni^2/Nd .* (x >= x_in_jun);
 n = ni^2/Na .* (x <= x_pi_jun) + Nd.* (x >= x_in_jun);
 V = Fn + Vth * log (n / ni);/

Unfortunately this produces the following errors:
 
/error: bimu_bernoulli: invalid conversion from NaN to logical
error: called from
    bimu_bernoulli at line 45 column 11
    secs1d_dd_gummel_map at line 94 column 12
    secs1d_dd_gummel_map_example_PIN at line 147 column 49
error: evaluating argument list element number 1
error: called from
    bimu_bernoulli at line 45 column 11
    secs1d_dd_gummel_map at line 94 column 12
    secs1d_dd_gummel_map_example_PIN at line 147 column 49/

I do not know what produces this error, does anyone perhaps know what I did
wrong?

Many thanks in advance,

Tijmen



--
Sent from: http://octave.1599824.n4.nabble.com/Octave-General-f1599825.html