The underlying problem (with n = 2000 unknowns at each parameter value) is a discretization of a two-point boundary value problem for a system of ordinary differential equations (ODE). The  numerical solution of such a system is standard, but requires nontrivial CPU time.
 

 

http://lcvmwww.epfl.ch