Mathematica 9 is now available
Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2012

[Date Index] [Thread Index] [Author Index]

Search the Archive

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



  • Prev by Date: Re: NDSolve output without interpolating function
  • Next by Date: Re: How to simplify this code with variable number of
  • Previous by thread: InverseFunction: how to manage?
  • Next by thread: Re: InverseFunction: how to manage?