On 11/7/05, Silke Diedenhofen <[EMAIL PROTECTED]> wrote:

> This is the function I like to integrate:
>
> double f (double omega, void * params) {
>       double Wg = *(double *) params;
>       double omega_0 = *(double *) params;
[snip]

Are you sure you're passing those parameters as intended?  I would
have expected to see something like

double f(double omega, void *vparams) {
    double *params = (double *) vparams;
    double Wg = params[0];
    double omega_0 = params[1];
    // etc.

Make sure that f() evaluates to the expected result for particular
values of omega and parameters.

-- mj


_______________________________________________
Help-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/help-gsl

Reply via email to