Drexel dragonThe Math ForumDonate to the Math Forum



Search All of the Math Forum:

Views expressed in these public forums are not endorsed by Drexel University or The Math Forum.


Math Forum » Discussions » Software » comp.soft-sys.math.mathematica

Topic: another query
Replies: 0  

Advanced Search

Back to Topic List Back to Topic List  
Trichy V. Krishnan

Posts: 8
Registered: 11/22/12
another query
Posted: Jun 23, 2013 10:50 PM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

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



Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© Drexel University 1994-2014. All Rights Reserved.
The Math Forum is a research and educational enterprise of the Drexel University School of Education.