MathGroup Archive 2013

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

Search the Archive

another query

  • To: mathgroup at smc.vnet.net
  • Subject: [mg131281] another query
  • From: "Trichy V. Krishnan" <biztvk at nus.edu.sg>
  • Date: Sun, 23 Jun 2013 22:57:38 -0400 (EDT)
  • Delivered-to: l-mathgroup@mail-archive0.wolfram.com
  • Delivered-to: l-mathgroup@wolfram.com
  • Delivered-to: mathgroup-outx@smc.vnet.net
  • Delivered-to: mathgroup-newsendx@smc.vnet.net

Hi:

I am trying to run the following program but it is not converging after more than 5 hrs. The thing is that it used to converge a week ago but it is stalling now. Any idea? I tried exiting the Mathematica and rebooting it, but the problem persists.

The following is the program.

With regards

Trichy


nc=.; qs=.;qL=.;rk=.; rm=.;  rr=.;pk=.;pm=.; al=.; be=.;xx=.; xx1=.;x1=.; x2=.; x3=.;rhs1=.;lhs2=.; lhs3=.;Piks1=.; PikL23=.; Piks3=.; Pik=.; Pim=.;Pidk=.; Pidm=.; Pidpk=.;Pidpm=.; dPikpk=.;
 s1k=.; s1m=.;L2k=.; L2m=.; soln=.; Ssk=.; SLk=.; Ssm=.; SLm=.; h=.;
s3k=.; L3k=.; L3m=.;csk=.; csm=.; cLk=.; cLm=.; soln=.; soln1=.; soln2=.;
 Pik1=.;Pim1=.;c=.;
 nc=1; qs=1; qL=3;

 rhs1=(2(1-xx)+qs rm)pm+(2+qs rk)(1-pm);
xx1=Solve[2 xx+qs rk?rhs1, xx][[1]] // Simplify;
 x1=xx/.xx1;
s1k=nc al (x1+(1-x1)(1-pm));

 xx=.; xx1=.;
lhs2=(2xx+qL rk)pk+(2+qL rm)(1-pk);
xx1=Solve[lhs2?2(1-xx)+qL rm, xx] [[1]]//Simplify;
x2=xx/.xx1;
L2k=nc (1-al) be x2 pk;
 L2m=nc (1-al) be((1-x2)+x2(1-pk));
 xx=.; xx1=.;
lhs3=(2xx+qL rk)pk+(2 xx+(qL/qs) qs rk)(1-pk);
xx1=Solve[lhs3?2(1-xx)+qL rm, xx][[1]]//Simplify;
 x3=xx/.xx1;
 L3k=nc (1-al) (1-be)x3 pk;
 s3k=nc (1-al) (1-be)x3 (1-pk);
 L3m=nc (1-al) (1-be)(1-x3); )

s1k+s3k+L2k+L3k+s1m+L2m+L3m // FullSimplify // Rationalize;(

 Ssk=s1k+s3k (qL/qs)//FullSimplify //Rationalize;
 SLk=L2k+L3k //FullSimplify;
 Ssm=s1m //FullSimplify;
SLm=L2m+L3m // FullSimplify;
 csk=1; csm=0; cLk=0; cLm=0;


 h=2;
 Pik=(Ssk qs rk-qs csk (Ssk^h))+(SLk qL rk-qL cLk (SLk ^h))// FullSimplify//Rationalize;
Pim=(Ssm qs rm-qs csm (Ssm^h))+(SLm qL rm-qL cLm (SLm ^h))// FullSimplify//Rationalize;
Pidk=D[Pik,rk]//FullSimplify//Rationalize;
Pidm=D[Pim,rm]//FullSimplify//Rationalize;
Soln=Solve[{Pidk==0,Pidm?0},{rk,rm}][[1]]//FullSimplify//Rationalize;
rk=rk/.Soln;rm=rm/.Soln;
 Pik1=Pik // FullSimplify//Rationalize;
Pim1=Pim // FullSimplify//Rationalize


  • Prev by Date: Position[] bug or feature ?
  • Next by Date: Re: bug or feature ?
  • Previous by thread: Re: bug or feature ?
  • Next by thread: Decoupling E and I