Mira, estuve mirando muy por arriba la pagina y te digo que si queres resolver ecuaciones diferenciales en C, C++, Pascal, lo que sea, se me va de las manos. Pero resolver ecuaciones diferenciales por metodos numericos, o sea, por ecuaciones en diferencia, es cuestion de que entiendas un poco de matematica basica, y luego tendras que escribir un metodo para calcularlo.
En este momento me acuerdo del metodo de Euler para discretizar esta ecuacion:
dy/dt=f(u,t)
entonces Euler propone: Un+1 = Un + k.f(Un,Tn)
donde k es el paso de calculo y Un sera la solucion numerica de la ecuacion.
Ahora, si miras este metodo, no tiene complejidad para programarlo. Primero el usuario tiene que ingresar Uo y luego podes empezar a iterar.
double Solucion = Uo+k.( funcion(Uo, T) )
y listo y si queres los podes ir cargando a un vector de double:
VecSolucion=Solucion;
Espero que sirva
P.D.: el metodo de Euler no es el mejor que existe, pero si te asegura convergencia a la solucion real, pero el orden de error que provee es muy grande.
Salu2!!!!!
