Quote:
Originariamente inviato da dario fgx
Grazie, adesso funziona.
Ora considerando che la mia espressione in realtà è:
-LS/(gΩ)y'' - 2y = (1/(2g))y'|y'|(Lλ/D)
come scrivo un codice veloce che mi permetta di fare una analisi al variare di λ, D, L, S, Ω ?
In matlab si poteva fare in modo che in input venissero chiesti questi parametri.
Ad esempio sarebbe utile un codice che mi permettesse di plottare la soluzione al variare di λ tra 0 e 100 (ad esempio).
Grazie mille!
|
Ho scritto al volo questo codice, ho fatto anche una prova con parametri inventati e da robe sensate
f[omega_, L_, d_, S_, lambda_, g_] :=
NDSolve[{-L S/(g omega) y''[x] - 2 y[x] == 1/(2 g) y'[x] Abs[y'[x]] L lambda/d , y[0] == 1, y'[0] == 0}, y[x], {x, 0, 10}];
(* fissi i valori di omega,L,d,S,g qui sotto*)
omega = 20. ;
L = 1. ;
d = 2.;
g = 0.2;
S = 4.;
a = 5; (* è il passo con cui vari lambda da 1 a 100*)
fun = Table[y[x] /. f[omega, L, d, S, lambda, g], {lambda, 1, 100, a}];
Plot[fun, {x, 0, 10}]