#include <stdio.h>
#include <math.h>
#include "genlib.h"


/*stable covariance function*/
double stable(double h, double alpha)
{
  return(exp(-3.*(double)pow(h,alpha)));
}