       Re: Math5.0 frustrated Compatibility Problem

• To: mathgroup at smc.vnet.net
• Subject: [mg56961] Re: [mg56932] Math5.0 frustrated Compatibility Problem
• From: Chris Chiasson <chris.chiasson at gmail.com>
• Date: Wed, 11 May 2005 05:24:27 -0400 (EDT)
• References: <200505100742.DAA08390@smc.vnet.net>
• Reply-to: Chris Chiasson <chris.chiasson at gmail.com>
• Sender: owner-wri-mathgroup at wolfram.com

```You should upload your notebook somewhere. The code was trashed by my
mail program.

On 5/10/05, Jingming Yao <jingming_yao at hotmail.com> wrote:
> Hello Everyone
>   There is a very interesting problem. I have a 4 equations 4 knowns
> nonlinear algebra equation system. It CAN be solved by Math 4.1 or 4.2 but
> it can NOT be solved by Math 5.0. It is the same program. And the error
> message in Math 5.0 is "FindRoot:: jsing: Encountered a singular Jacobian at
> the point {p02n, Æ?Â¢Vb2n, n02n, JCn, Æ?Ã?n}={2.9*10^18, 0.61, 3.7*10^11, 0.5}.
> Try perturbing the initial point(s)." My program is as follows.
>
> u2 = d0 - ((Cb*(
>       zVb2 - Æ?Â¢Vb2n))/(q*Ns02) + q*Lz*(p02n)/(q*Ns02) + (&#1013;/d)*(zÆ?Â³s +
> Æ?Â¢Æ?Â³sn)/(q*
>             Ns02));
>
> u4 = JBp - Jpe0*((p02n/zp02)*
>           Exp[-Æ?Â¢Æ?Â³sn/VT] - 1) - Jpb0*((
>                   p02n/zp02)*Exp[(Æ?Â¢Vb2n)/VT] - 1) - Jrc0*(Â?Ã£(p02n/
>     zp02)*Exp[-Æ?Â¢Æ?Â³sn/(2*VT)] -
>         1) - Jrb0*(Â?Ã£(p02n/zp02)*Exp[(Æ?Â¢Vb2n)/(2*VT)] - 1) -
> q*A*Lz*p02n*n02n;
>
> u6 = Æ?Â¢Æ?Â³sn - (1 - Æ?Â¢Vb2n);
>
> u12 = n02n*q*Vn*Exp[-Æ?Â³emc2/VT] + n02n*q*Vn*Exp[-Æ?Â³eme2/VT] +
>         q*A*Lz*p02n*n02n + (n02n*p02n/(n02n + p02n))*(q*Lz)/Æ?Ã?NR - (Æ?Ã?nc*
>             Jnc0*(Exp[(Æ?Â¢Vb2n)/VT] - 1) + Æ?Ã?ne*Jne0*(Exp[-Æ?Â¢Æ?Â³sn/VT] - 1));
>
> JBp = 0.001;
>
> Answer[i] = FindRoot[{u2, u4, u6, u12}, {p02n, 2.9*10^18}, {Æ?Â¢Vb2n, 0.61}, {
>     n02n, 3.7*10^11}, {Æ?Â¢Æ?Â³sn, 0.5}]
>
> The starting values I use above are close to answers Math 4.1 gives me. All
> the others are constants with values like
>
> d0 = 8*10^-7;
> d = 5*10^-5;
> &#1013; = 13.1*8.854*10^-14;
> db = 2.5*10^-6;
> Cb = &#1013;/db;
> zVb2 = 0.843;
> q = 1.60218*10^-19;
> Ns02 = 3.2*10^18;
> Lz = 6*10^-7;
> Jpe0 = 6.29*10^-19;
> zp02 = 4*10^9;
> zÆ?Â³s = 0.806;
> VT = 0.0259;
> Jpb0 = 1.6*10^-29;
> Jrc0 = 1.065*10^-10;
> Jrb0 = 5.33*10^-12;
> A = 10^-10;
> Vn = 10^7;
> Æ?Â³emc2 = 0.43;
> Æ?Â³eme2 = 0.446;
> Æ?Ã?NR = 10^-9;
> Æ?Ã?nc = 0.001;
> Æ?Ã?ne = 0.001;
> Jnc0 = 2.57*10^-10;
> Jne0 = 2.57*10^-10;
> Æ?Â³sn = zÆ?Â³s + Æ?Â¢Æ?Â³sn;
>
> Anybody has some idea? Thnaks a lot.
>
> Jingming
> Uconn, ECE Dept.
>
> _________________________________________________________________
> FREE pop-up blocking with the new MSN Toolbar â?? get it now!
> http://toolbar.msn.click-url.com/go/onm00200415ave/direct/01/
>
>

--
Chris Chiasson
http://chrischiasson.com/
1 (810) 265-3161

```

• Prev by Date: Re: Re: New Mathematica menu: New Package
• Next by Date: Re: Representation and Simulation of Dynamic Systems
• Previous by thread: Re: Math5.0 frustrated Compatibility Problem