Re: InverseFunction: how to manage?
- To: mathgroup at smc.vnet.net
- Subject: [mg124298] Re: InverseFunction: how to manage?
- From: Murray Eisenberg <murray at math.umass.edu>
- Date: Mon, 16 Jan 2012 17:09:07 -0500 (EST)
- Delivered-to: l-mathgroup@mail-archive0.wolfram.com
- References: <201201142217.RAA13006@smc.vnet.net>
- Reply-to: murray at math.umass.edu
First, Mathematica indicates that you need to impose a condition in order that your paper-and-pencil integral gives the result you claim: Integrate[1/Sqrt[Sin[u]], {u, 0, f}] ConditionalExpression[ -2*EllipticF[(-2*f + Pi)/4, 2] + Sqrt[2]*EllipticK[1/2], Sin[f] >= 0] Thus take: t[f_] = Integrate[1/Sqrt[Sin[u]], {u, 0, f}, Assumptions -> 0 <= f <= Pi] -2*EllipticF[(-2*f + Pi)/4, 2] + Sqrt[2]*EllipticK[1/2] Is the following of any use? f /. First@Solve[t[f] == tau, f] (Pi - 4*JacobiAmplitude[(-tau + Sqrt[2]*EllipticK[1/2])/2, 2])/2 (You'll get a warning that inverse functions are being used by Solve and so that some solutions may not be found.) On 1/14/12 5:17 PM, Dr. Wolfgang Hintze wrote: > I apologize for asking this very elementary question but how do I > manage InverseFunction? > > Here is an example > > When I solve the equation of motion for a pendulum > > f''[t] == Cos[f[t]], f[0]== 0, f'[0] == 0 > > I get (with paper and pencil) the time t as a function of the angle f > thus > > t[f_] = Integrate[1/Sqrt[Sin[u]], {u, 0, f}] > Out[22]= > 2*(EllipticF[f/2 - Pi/4, 2] + EllipticF[Pi/4, 2]) > > Now I want the the angle as a function of time (f[t]) like this > > "f[t_] = InverseFunction[t[f]]" > > But this does not work. I also tried to define t as a pure function > > t = 2*(EllipticF[#1/2 - Pi/4, 2] + EllipticF[Pi/4, 2])& > > but again, I have not seen a way to invert this, and for instance carry > out Plot[f,{t,0,2 Pi}]. > > Thanks in advance for any hints. > > Best regards, > Wolfgang > > > > -- Murray Eisenberg murray at math.umass.edu Mathematics & Statistics Dept. Lederle Graduate Research Tower phone 413 549-1020 (H) University of Massachusetts 413 545-2859 (W) 710 North Pleasant Street fax 413 545-1801 Amherst, MA 01003-9305
- References:
- InverseFunction: how to manage?
- From: "Dr. Wolfgang Hintze" <weh@snafu.de>
- InverseFunction: how to manage?