4 const gsl_odeiv_step_type * T
8 = gsl_odeiv_step_alloc (T, 2);
11 gsl_odeiv_system sys = {func, jac, 2, &mu};
13 double t = 0.0, t1 = 100.0;
15 double y[2] = { 1.0, 0.0 }, y_err[2];
16 double dydt_in[2], dydt_out[2];
18 /* initialise dydt_in from system parameters */
19 GSL_ODEIV_FN_EVAL(&sys, t, y, dydt_in);
23 int status = gsl_odeiv_step_apply (s, t, h,
29 if (status != GSL_SUCCESS)
32 dydt_in[0] = dydt_out[0];
33 dydt_in[1] = dydt_out[1];
37 printf ("%.5e %.5e %.5e\n", t, y[0], y[1]);
40 gsl_odeiv_step_free (s);