       Re: Help with solving ODE

• To: mathgroup at smc.vnet.net
• Subject: [mg81789] Re: Help with solving ODE
• From: Norbert Marxer <marxer at mec.li>
• Date: Wed, 3 Oct 2007 06:10:59 -0400 (EDT)
• References: <fdvdq7\$sul\$1@smc.vnet.net>

```On 3 Okt., 08:45, Pioneer1 <1pione... at gmail.com> wrote:
> Hi,
>
> Can anyone help solve this linearized differential equation:
>
> Iy'' + ky' = 2GMmd/a^2
>
> Primes are time derivates of y (=theta=excursion angle). Is it
> possible to solve this for the initial conditions y(0)=0 and y'(0)=0?
>
> I got the solution at sci.math for the non-linear version and I want
>
>
> Further information is also available at sci.physics.research
>
>
> Parameters are:
>
> > y = theta = excursion angle in radians
> > A = I = moment of inertia = 13,138,117.34 g cm^2
> > B = R = damping = for now I assume this to be zero
> > C = k = torsion constant  = 724.68 g cm^2 sec^-2
> > d = moment arm = 93.09 cm
> > D = 2GMmd = 2 * 6.67*10^-8 * 158100 * 729.8 * 93.09 = 1432.82
> > a = distance between weights = 22.10 cm
>
> I would truly appreciate help with this. Thanks

Hello

With your data and units you can use

sol=DSolve[{13.138 10^6 y''[t]+724.68
y[t]==1432.82/22.10^2,y'==0,y==0},y[t],t];
Plot[sol[[1,1,2]],{t,0,3000}]

I hope this helps.

Best Regards
Norbert Marxer

```

• Prev by Date: Re: Equivalent functionality to colorbar in Mathematica?
• Next by Date: Re: Re: Tooltips in ContourPlot
• Previous by thread: Re: Help with solving ODE
• Next by thread: Re: Help with solving ODE