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