       Re: NSolve Vs. Elliptic Integral

• To: mathgroup at smc.vnet.net
• Subject: [mg62605] Re: NSolve Vs. Elliptic Integral
• From: "nilaakash at gmail.com" <nilaakash at gmail.com>
• Date: Tue, 29 Nov 2005 04:45:21 -0500 (EST)
• References: <dme6pg\$i69\$1@smc.vnet.net>
• Sender: owner-wri-mathgroup at wolfram.com

```Dear All,
Thanks everybody to answer my question. Though I have
got my answer still there are some prblems.
Basically I wanted to get a particulam "m" value for
fixed g[m] value. Mr. Valeri and Mr. Bob both have given nice answer.
Now if I change my g[m] value, means if I increase
g[m], the "m" value reaches close Pi/2 and the integral diverges there.
As my integral function is symmetric about Pi/2 axis, when I ask to
evaluate "m" value near Pi/2 region suddenly it jumps from left side
curve to right side curve, and I get jump.
I am giving the code and please suggest me how to get
those "m" values which are belong to curve which is only extended
from 0.003 to Pi/2.

Clear[f,g];

f[x_,m_]:=Sqrt[(1+0.176*Sin[m]^2*Sin[x]^2)*
(1+1.01769*Sin[m]^2*Sin[x]^2)/(1-Sin[m]^2*Sin[x]^2)];

g[m_]:=W*46381*10^(-6)*Sqrt[1/(-0.0012245+Sin[m]^2)];

h[m_?NumericQ] := NIntegrate[f[x, m], {x, ArcSin[0.002/Sin[m]],Pi/2}]

t={};For[W = 1, W <= 400, W = W + 1,
a = FindRoot[h[m]==g[m], {m, 0.1}];
num = m /. a;
t=Join[t,Table[{num,W},{i,1,1}]];]
MatrixForm[t]

ListPlot[t, PlotRange -> All, PlotRange -> All, PlotJoined -> False,
Frame -> True, PlotStyle -> {PointSize[0.008], Hue[0.6]}];

Thanks.

nilaakash

```

• Prev by Date: Re: references to numbered equations in arXiv style sheet
• Next by Date: Re: MatrixForm
• Previous by thread: Re: Re: Re: NSolve Vs. Elliptic Integral
• Next by thread: Notation package in batch mode