Newsgroup: it.comp.lang.c++
Titolo: aiuto algoritmo!!

cerco un algoritmo in c++ che mi risolva il metodo di integrazione
numerica di eulero in avanti....aiutatemi vi prego

Se mi spieghi in cosa consiste, magari una mano te la do....ho la teoria
matematica un'pò arrugginita!!
--------------------------------
Inviato via 
http://usenet.libero.it

 numerica di eulero in avanti....aiutatemi vi prego
"in avanti"??? vuoi forse dire con una precisione superiore al primo
ordine.
Se e' cosi' ci sono:
 - il metodo del punto mediano (secondo ordine)
 - il metodo di Runge-Kutta (e' molto diffuso quello del quarto ordine,
il punto mediano e' in realta' un RK2...)
-- 
Posted via Mailgate.ORG Server - 
http://www.Mailgate.ORG

 numerica di eulero in avanti....aiutatemi vi prego
non ho ben capito cosa voglia dire "in avanti"
cmq prova col metodo dei rettangoli, che usa eulero:
double F(double x), //la funzione che devi integrare
double eulero(double a, double b, double dx)
{
double area=0;
double x=a;
while (x<b)
{area=+F(x)*dx; x=+dx;}
return area;
}
double integrale(double a, double b)
{
double EPS=...; //lo scostamento max che vuoi tra le due aree
                         //o errore "epsilon"
short MAXP=...;//il num max di passaggi che vuoi fare
                            //vedi tu al precisione
double area_min, area_max;
double dx=(b-a)/2;
short pass=0;
area_min=(b-a)*F(a);
do
{area_max=area_min;
area_min= eulero(a,b,dx);
dx=dx/2;}
while ((area_max-area_min)<EPS && pass<MAXP);
return area_min;
}
non l'ho compilato, è quello che mi ricordo che ho studiato
per l'esame di fondamenti.. e te l'ho scritto di getto
quindi non volermene per gli eventuali errori.
il procedimento cmq dovrebbe essere giusto.
se vuoi algoritmi più efficienti prova ad usare quello dei trapezi,
che è più accurato (nel risultato intendo)

 il punto mediano e' in realta' un RK2...)
Devo risolvere l'equazione differenziale del primo ordine della caduta
di un corpo con il metodo di eulero in avanti (forward eulero) con una
precisione del secondo ordine.In pratica dovrebbe essere molto simile
a quello del punto intermedio solo che invece di prendere il punto
intermedio prende il punto successivo, comunque se hai solo
l'algoritmo del punto intemedio va bene lo stesso al limite lo
ritocco....ti ringrazio comunque per la cortesia..spero di ricevere
risposta.

 intermedio prende il punto successivo, comunque se hai solo
ok, scusa il ritardo...
http://guide.supereva.it/fortran/interventi/2001/08/61578.shtml
http://math.gmu.edu/~kblake/proj4/Euler.html
http://www.myphysicslab.com/runge_kutta.html
http://www.ap.univie.ac.at/users/ves/cp0102/dx/node69.html
http://www.csit.fsu.edu/~navon/RK4_C.pdf
http://www.ulib.org/webRoot/Books/Numerical_Recipes/bookcpdf/c16-1.pdf
http://www.cnbc.cmu.edu/~bard/numerics/node4.html
http://www.sali.freeservers.com/engineering/fortran_codes/nystrom.html
http://www.vislab.usyd.edu.au/users/manuals/maths/matlab5/Runge.html
ti bastano????


 wrote in message
c889939c6028cfdb672...@mygate.mailgate.org...

LOL
Ma perchč in italia lo chiamiamo algoritmo di eulero e all'estero lo
chiamano runge-kutta algorithm?
Se č lo stesso non dovrebbe essere conosciuto con lo stesso nome?
Avevo provato anche io a cercare questi link su internet per rispondere
all'op ma cercando Eulero avevo trovato solo il primo della lista (l'unico
in italiano), e non sapendo che all'estero lo chiamano in un altro modo mi
sono fermato lě......

 Se č lo stesso non dovrebbe essere conosciuto con lo stesso nome?
Giusta osservazione, ma non e' lo stesso! o meglio, l'algoritmo di
Eulero
e' una "specializzazione" dell'algoritmo di Runge-Kutta (che, come avrai
visto, ha molte varianti) fermandosi a considerare soltanto il primo
ordine (e da qui l'errore e' sul secondo).


 ordine (e da qui l'errore e' sul secondo).
Ho capito.
Grazie.