Help with plotting and iterations
- To: mathgroup at smc.vnet.net
- Subject: [mg72808] Help with plotting and iterations
- From: "ashesh" <ashesh.cb at gmail.com>
- Date: Sat, 20 Jan 2007 03:48:50 -0500 (EST)
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