htswanalysis/MACS/lib/gsl/gsl-1.11/doc/examples/odefixed.c

   1 int
   2 main (void)
   3 {
   4   const gsl_odeiv_step_type * T
   5     = gsl_odeiv_step_rk4;
   6
   7   gsl_odeiv_step * s
   8     = gsl_odeiv_step_alloc (T, 2);
   9
  10   double mu = 10;
  11   gsl_odeiv_system sys = {func, jac, 2, &mu};
  12
  13   double t = 0.0, t1 = 100.0;
  14   double h = 1e-2;
  15   double y[2] = { 1.0, 0.0 }, y_err[2];
  16   double dydt_in[2], dydt_out[2];
  17
  18   /* initialise dydt_in from system parameters */
  19   GSL_ODEIV_FN_EVAL(&sys, t, y, dydt_in);
  20
  21   while (t < t1)
  22     {
  23       int status = gsl_odeiv_step_apply (s, t, h,
  24                                          y, y_err,
  25                                          dydt_in,
  26                                          dydt_out,
  27                                          &sys);
  28
  29       if (status != GSL_SUCCESS)
  30           break;
  31
  32       dydt_in[0] = dydt_out[0];
  33       dydt_in[1] = dydt_out[1];
  34
  35       t += h;
  36
  37       printf ("%.5e %.5e %.5e\n", t, y[0], y[1]);
  38     }
  39
  40   gsl_odeiv_step_free (s);
  41   return 0;
  42 }