Solving this equation is a good opportunity for me to getting to know ODE solvers: Euler, Runge-Kutta etc and their implementations in Matlab and Mathematica. I'm off to get a few books to read on them. The general plan is:
1. Double-check that the equation is right
2. Study numerical techniques theoretically
3. Write code in Matlab and/or Mathematica
My only concern is how am I going to implement my external function/cause which is a weird bunch of gaussians modulated by another gaussian and then a decreasing exponential.
