MathGroup Archive 2002

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

Search the Archive

Re: Inv.Interpol.Function

  • To: mathgroup at smc.vnet.net
  • Subject: [mg36720] Re: Inv.Interpol.Function
  • From: Erich Mueller <emueller at mps.ohio-state.edu>
  • Date: Fri, 20 Sep 2002 04:16:40 -0400 (EDT)
  • Organization: Ohio State University
  • References: <am95s4$k3d$1@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

A slick way to solve this particular problem is to use NDSolve.  The
inverse function t[y] sattisfies the differential equation

t'[y]==1/R[t[y]]^4

with the inital condition

t[0]==0

So the following command should give an interpolating function for t[y]
(you should change the 1 in {y,0,1} to whatever the appropriate value is)

NDSolve[{t'[y]==1/R[t[y]]^4,t[0]==0},t,{y,0,1}]

Erich

On Wed, 18 Sep 2002, martin skogstad wrote:

> Dear NG
>
>
> Mabye this is too simple, but I cant just figure it out
>
> I want to get an inverse function for y[t] where
>
> y[t_]:=NIntegrate[R[x]^4,{x,0,t}] /. ndsolution[[1]]
> and
> R[t] is an interpolatingfunction(R>0 from NDSolve) on the interval 0=<t=<T
>
> after that I hope to be able to calculate the integral :
>
> a = (1/y[T])*
>         NIntegrate[R[InvFunction[y[t]]*Cos[y[t]], {y[t], 0, y[T]}]
>
> it works with a constant instead of R[t].
>
> Hope you can help, and that the above is understandable.
>
>
> Martin Skogstad
>
>



  • Prev by Date: Re: Front end proble in Linux
  • Next by Date: RE: Drawing an ellipse
  • Previous by thread: Inv.Interpol.Function
  • Next by thread: Re: problems with the definition of a function