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
- References:
- Math5.0 frustrated Compatibility Problem
- From: "Jingming Yao" <jingming_yao@hotmail.com>
- Math5.0 frustrated Compatibility Problem