[Date Index]
[Thread Index]
[Author Index]
Re: Math5.0 frustrated Compatibility Problem
*To*: mathgroup at smc.vnet.net
*Subject*: [mg56958] Re: [mg56932] Math5.0 frustrated Compatibility Problem
*From*: DrBob <drbob at bigfoot.com>
*Date*: Wed, 11 May 2005 05:24:17 -0400 (EDT)
*References*: <200505100742.DAA08390@smc.vnet.net>
*Reply-to*: drbob at bigfoot.com
*Sender*: owner-wri-mathgroup at wolfram.com
I suggest you convert that to InputForm and post again. It's useless in that form.
Bobby
On Tue, 10 May 2005 03:42:48 -0400 (EDT), 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) + (ϵ/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;
> ϵ = 13.1*8.854*10^-14;
> db = 2.5*10^-6;
> Cb = ϵ/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.
> Graduate Assistance
>
> _________________________________________________________________
> FREE pop-up blocking with the new MSN Toolbar â?? get it now!
> http://toolbar.msn.click-url.com/go/onm00200415ave/direct/01/
>
>
>
>
--
DrBob at bigfoot.com
Prev by Date:
**Re: Re: FilledPlot: Curves->Back option and Epilog not working?**
Next by Date:
**Re: Re: New Mathematica menu: New Package**
Previous by thread:
**Math5.0 frustrated Compatibility Problem**
Next by thread:
**Re: Math5.0 frustrated Compatibility Problem**
| |