Re: [Help-gsl] gsl_linear_interp example

[Daniel J Farrell]

Sorry made silly mistake:

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <gsl/gsl_interp.h>

int main (void)
{
      //sample data
      double x[4] = { 1970, 1980, 1990, 2000 };
      double y[4] = {   12,   11,   14,   13 };

      //inialise and allocate the gsl objects
      gsl_interp *interpolation = gsl_interp_alloc(gsl_interp_linear, 
4);
      gsl_interp_init(interpolation, x, y, 4);
      gsl_interp_accel * accelerator =  gsl_interp_accel_alloc();

      //get interpolation for x = 1981
      double value = gsl_interp_eval(interpolation, x, y, 1981,  
accelerator);

      printf("\n%g",value);

      //output:
      //11.3

  return 0;
}

On 6 Jun 2007, at 23:39, Daniel J Farrell wrote:

> #include <stdio.h>
> #include <stdlib.h>
> #include <math.h>
> #include <gsl/gsl_interp.h>
>
> int main (void)
> {
>       //sample data
>       double x[4] = { 1970, 1980, 1990, 2000 };
>       double y[4] = {   12,   11,   14,   13 };
>
>       //inialise and allocate the gsl objects
>       gsl_interp *interpolation = gsl_interp_alloc 
> (gsl_interp_linear,n);
>       gsl_interp_init(interpolation, x, y, n);
>       gsl_interp_accel * accelerator =  gsl_interp_accel_alloc();
>
>       //get interpolation for x = 1981
>       double value = gsl_interp_eval(interpolation, x, y, 1981,  
> accelerator);
>
>       printf("\n%g",value);
>
>       //output:
>       //11.3
>
>   return 0;
> }