I have tried several things, but I just can't work out how to do this, so perhaps someone could help me here. This is an example taken from http://www.gnu.org/software/gsl/manual/html_node/ODE-Example-programs.html on how to carry out integrals.

``````#include <stdio.h>
#include <gsl/gsl_errno.h>
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_odeiv.h>

int
func (double t, const double y[], double f[],
void *params)
{
double mu = *(double *)params;
f = y;
f = -y - mu*y*(y*y - 1);
return GSL_SUCCESS;
}

int
jac (double t, const double y[], double *dfdy,
double dfdt[], void *params)
{
double mu = *(double *)params;
gsl_matrix_view dfdy_mat
= gsl_matrix_view_array (dfdy, 2, 2);
gsl_matrix * m = &dfdy_mat.matrix;
gsl_matrix_set (m, 0, 0, 0.0);
gsl_matrix_set (m, 0, 1, 1.0);
gsl_matrix_set (m, 1, 0, -2.0*mu*y*y - 1.0);
gsl_matrix_set (m, 1, 1, -mu*(y*y - 1.0));
dfdt = 0.0;
dfdt = 0.0;
return GSL_SUCCESS;
}

int
main (void)
{
const gsl_odeiv_step_type * T
= gsl_odeiv_step_rk8pd;

gsl_odeiv_step * s
= gsl_odeiv_step_alloc (T, 2);
gsl_odeiv_control * c
= gsl_odeiv_control_y_new (1e-6, 0.0);
gsl_odeiv_evolve * e
= gsl_odeiv_evolve_alloc (2);

double mu = 10;
gsl_odeiv_system sys = {func, jac, 2, &mu};

double t = 0.0, t1 = 100.0;
double h = 1e-6;
double y = { 1.0, 0.0 };

while (t < t1)
{
int status = gsl_odeiv_evolve_apply (e, c, s,
&sys,
&t, t1,
&h, y);

if (status != GSL_SUCCESS)
break;

printf ("%.5e %.5e %.5e\n", t, y, y);
}

gsl_odeiv_evolve_free (e);
gsl_odeiv_control_free (c);
gsl_odeiv_step_free (s);
return 0;
}``````

Now in the int main function there are these two lines:

``````double mu = 10;
gsl_odeiv_system sys = {func, jac, 2, &mu};``````

I made mu an array and then passed it into the gsl_odeiv_system, but then when I went to retrieve the values from the array in the int func function I tried a whole manner of things and just couldn't manage to get the values, does anyone know what I should do to be able to pass the array through and then retrieve the values in the func function?

Be a part of the DaniWeb community

We're a friendly, industry-focused community of 1.20 million developers, IT pros, digital marketers, and technology enthusiasts learning and sharing knowledge.