Mathematica 9 is now available
Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2005
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2005

[Date Index] [Thread Index] [Author Index]

Search the Archive

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.
> Graduate Assistance
> 
> _________________________________________________________________
> 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
  • Next by thread: Debugging