Pagina 1 di 2

Risoluzione eq differenziale in mathematica

MessaggioInviato: 23 ott 2013, 16:33
da Rabeluk
Come si risolve questa equazione numericamente in mathematica?

s''(t) \left(\frac{I_g}{R^2}+M\right)-g M \sin (\theta )+k s(t)=0

Re: risoluzione eq differenziale in mathematica

MessaggioInviato: 23 ott 2013, 18:15
da PietroBaima
Prova così:

Codice: Seleziona tutto
IG = 1;
R = 1;
k = 1;
M = 1;
g = 9.81;
Theta = Pi/3;
EXT = 20;
INIT0 = 0;
INIT1 = 0;
s = NDSolve[{y''[x]*(IG/R^2 + M) + k*y[x] - g*M*Sin[Theta] == 0,
   y[0] == INIT0, y'[0] == INIT1}, y, {x, 0, EXT}]
Plot[Evaluate[y[x] /. s], {x, 0, EXT}, PlotRange -> All]


La equazione differenziale verrà integrata per x\in (0,EXT).
Vedi tu come definire tutte le variabili, condizioni iniziali INIT0 e INIT1 comprese.

Ciao,
Pietro.

Re: risoluzione eq differenziale in mathematica

MessaggioInviato: 23 ott 2013, 18:23
da Rabeluk
grazie mille :D

non è che potresti indicarmi cosa scrivere per risolvere questo sistema???

\left\{-g M \sin (\theta )+k s(t)+p'(t)=0,s'(t)-\frac{p(t)}{\frac{I_g}{R^2}+M}=0\right\}

Re: risoluzione eq differenziale in mathematica

MessaggioInviato: 23 ott 2013, 18:24
da PietroBaima
sì, certo che potrei :mrgreen:

Re: risoluzione eq differenziale in mathematica

MessaggioInviato: 23 ott 2013, 18:28
da PietroBaima
...se però non definisci p(t), numericamente è un pochino... impossibile!

Re: risoluzione eq differenziale in mathematica

MessaggioInviato: 23 ott 2013, 18:30
da Rabeluk
p'[t]=s'(t) \left(\frac{I_g}{R^2}+M\right)

cosi?

Re: risoluzione eq differenziale in mathematica

MessaggioInviato: 23 ott 2013, 18:32
da PietroBaima
ah, hai cambiato anche la prima, ok.

Re: risoluzione eq differenziale in mathematica

MessaggioInviato: 23 ott 2013, 18:33
da Rabeluk
comunque ho sbagliato nn è p'[t] ma p[t] sorry #-o

Re: risoluzione eq differenziale in mathematica

MessaggioInviato: 23 ott 2013, 18:36
da PietroBaima
:-M
Mi riscrivi il sistema, per favore?

Re: risoluzione eq differenziale in mathematica

MessaggioInviato: 23 ott 2013, 18:37
da PietroBaima
aspetta, ma hai modificato di nuovo il post...
adesso è corretto?