Re: Help with plotting and iterations

*To*: mathgroup at smc.vnet.net*Subject*: [mg72821] Re: [mg72808] Help with plotting and iterations*From*: Bob Hanlon <hanlonr at cox.net>*Date*: Sun, 21 Jan 2007 05:33:29 -0500 (EST)*Reply-to*: hanlonr at cox.net

result[p_]:=Module[{a=1/2,b=1/4,h=1/4, m,m1,z,t,t1,t2,b1,c,x,aa,bb}, m=b/a; m1=Sqrt[1-m^2]; z=(p-I*h); t=ArcSin[z]; t1=Re[t]; t2=-Im[t]; b1=Cot[t1]^2+m^2*Sinh[t2]-m1; c=m1^2*Cot[t1]^2; x=Reduce[{x^2 - b1*x - c==0,x>=0},x]// Last//N; aa=ArcCot[Sqrt[x]]; bb=ArcTan[Sqrt[x/m1]]; {EllipticF[aa,m],EllipticF[bb,m1]}]; TableForm[Join[ Table[Prepend[result[p],p],{p,0.01,0.99,0.01}], Table[Prepend[result[p],p],{p,1,10}]], TableHeadings->{None,{"p","Faa","Fbb"}}] ParametricPlot[result[p], {p, 0.01, Pi}, AxesLabel->{"Faa","Fbb"},ImageSize->432]; Bob Hanlon ---- ashesh <ashesh.cb at gmail.com> wrote: > > Hello All, > > Need help in working with Mathematica to determine (result_aa, > result_bb) values over a range of variable "p". > > p = 0.01; > > a = 0.5; b = 0.25; h = 0.25; > m = b/a; > m1 = Sqrt[1-m m]; > z = (p-i h)/b; (z is complex and (i h) represents imaginary part) > t = ArcSin[z]; > t1 = Re[t]; > t2 = -Im[t]; > b1 = Cot[t1]^2 + m^2*Sinh[t2]-m1; > c = m1^2*Cot[t1]^2; > Solve[x^2 - b1*x - c,x]; > x = x/. % > x = Last[x]; {i need the positive root from the two} > aa = ArcCot[Sqrt[x]]; > bb = Arctan[Sqrt[x/m1]]; > Faa = EllipticF[aa,m]; > Fbb = EllipticF[bb,m1]; > > result_aa(count) = Faa; > result_bb(count) = Fbb; > count = count+1 > > I would like to determine the values of Faa and Fbb over a range of p = > 0.01 to 10 with increments of 0.01 till 1.0 and then an increment of 1, > that is [0.01:0.01:1 2:1:10] > > Finally, I would like to plot with result_aa along x-axis and result_bb > along y-axis. > > Hope some one can solve the above problem. > > Ashesh >