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