Re: NDSolve/InterpolatingFunction and vectors
- To: mathgroup at smc.vnet.net
- Subject: [mg53810] Re: NDSolve/InterpolatingFunction and vectors
- From: D Herring <dherring at at.uiuc.dot.edu>
- Date: Fri, 28 Jan 2005 02:43:56 -0500 (EST)
- References: <csqrsc$1sl$1@smc.vnet.net> <ct55s7$e6s$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Hi Jens,
Thanks for looking at my problem. Your solution helped a lot and gives
nearly the result I was wanting. Its probably the best option
available. I'm still disappointed that Mathematica doesn't properly
handle the "dimension" of InterpolatingFunctions.
Later,
Daniel
Jens-Peer Kuska wrote:
> Hi,
>
> soln = NDSolve[{
> xx'[t] == {{0, 1}, {-1, 0}}.xx[t], xx[0] == {{0}, {1}}
> }, xx, {t, 0, 10}]
>
>
>
> Needs["DifferentialEquations`NDSolveUtilities`"]
>
> and
>
> time = DifferentialEquations`NDSolveUtilities`Private`GetTimeData[soln];
> grid = First[
> DifferentialEquations`NDSolveUtilities`Private`GetGridData[soln]];
> ip1 = Interpolation[Transpose[{time, #[[1, 1, 1]] & /@ Transpose[grid]}]];
> ip2 = Interpolation[Transpose[{time, #[[1, 2, 1]] & /@ Transpose[grid]}]]
>
> helps not ?
>
> at least
>
> Plot[{ip1[t], ip2[t]}, {t, 0, 10}]
>
> work fine.
>
> Regards
>
> Jens