Re: Attempt to generalize a constant
- To: mathgroup at smc.vnet.net
- Subject: [mg57705] Re: [mg57684] Attempt to generalize a constant
- From: Pratik Desai <pdesai1 at umbc.edu>
- Date: Sun, 5 Jun 2005 04:17:47 -0400 (EDT)
- References: <200506040704.DAA11779@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Narasimham wrote:
>F[t_,mu_]= mu*JacobiSN[t,mu^2] is a function between +/- 1 extreme
>limits.
>
>FindRoot[F[t,mu]==1, {mu,1,1.5},MaxIterations-> 100 ] works for given t
>= Pi/4, solution is (mu-> 2.1236).
>
>But when it is attempted to make t as a variable to plot solutions for
>a range of t values in:
>
>ImplicitPlot[F[t,mu]==1,{mu,1,1.5},{t,Pi/6,Pi}], it does not compute
>fully. TIA for any tips.
>
>
>
Hi
I may be wrong about this, but the elliptic integral (functions)
have no closed form solution, as a result the idea of using t as a
explicit variable may not work here.
I have tried using Table where I change the values of t using the
iterator, Anyway here is my attempt, actually you can control the number
of iteration in order to get a smoother plot by changing n (see below).
Also I have used Mathematica syntax for the function
as
JacobiSN[u,m] instead of JacobiSN[t,mu^2], I think that is the reason
why I get values that are squared of what you have
Here is the code
Clear[F, t, m, sol1]
n = 2
sol1 = Table[First[FindRoot[Sqrt[m]*JacobiSN[u, m] ==1, {m, 1, 1.5},
MaxIterations -> 100]], {u, Ï?/4, Ï?, Ï?/(n*12)}]
l = Length[sol1]
m1 = Range[Ï?/4, Ï?, Ï?/(n*12)] // N
dat = Table[{m1[[s]], m /. sol1[[s]]}, {s, 1, l}]
ListPlot[dat, PlotJoined -> True]
Hope this helps
Pratik
--
Pratik Desai
Graduate Student
UMBC
Department of Mechanical Engineering
Phone: 410 455 8134
- References:
- Attempt to generalize a constant
- From: "Narasimham" <mathma18@hotmail.com>
- Attempt to generalize a constant