Mercurial > forge
changeset 11744:d995950329f1 octave-forge
compute displacement current
| author | cdf |
|---|---|
| date | Wed, 05 Jun 2013 05:23:23 +0000 |
| parents | b99bf9d12345 |
| children | a6f6a1deefbb |
| files | extra/secs1d/inst/secs1d_newton.m |
| diffstat | 1 files changed, 8 insertions(+), 8 deletions(-) [+] |
line wrap: on
line diff
--- a/extra/secs1d/inst/secs1d_newton.m Wed Jun 05 00:12:39 2013 +0000 +++ b/extra/secs1d/inst/secs1d_newton.m Wed Jun 05 05:23:23 2013 +0000 @@ -150,15 +150,15 @@ Itot(:, tstep) = Jn([1 end], tstep) + Jp([1 end], tstep); - # Itot(1, tstep) += constants.e0 * material.esir * ... - # ((V(2, tstep) - V(1, tstep)) - - # (V(2, tstep-1) - V(1, tstep-1))) / ... - # (device.x(2) - device.x(1)); + Itot(1, tstep) += constants.e0 * material.esir * ... + ((V(2, tstep) - V(1, tstep)) - + (V(2, tstep-1) - V(1, tstep-1))) / ... + (device.x(2) - device.x(1)); - # Itot(2, tstep) += constants.e0 * material.esir * ... - # ((V(end, tstep) - V(end-1, tstep)) - - # (V(end, tstep-1) - V(end-1, tstep-1))) / ... - # (device.x(end) - device.x(end-1)); + Itot(2, tstep) += constants.e0 * material.esir * ... + ((V(end, tstep) - V(end-1, tstep)) - + (V(end, tstep-1) - V(end-1, tstep-1))) / ... + (device.x(end) - device.x(end-1)); Itot(:, tstep) *= device.W;