#ifndef _LEGENDRE_RULE_HPP_
#define _LEGENDRE_RULE_HPP_

# include <cstring>

void cdgqf ( int nt, int kind, double alpha, double beta, double t[],
  double wts[] );
void cgqf ( int nt, int kind, double alpha, double beta, double a, double b,
  double t[], double wts[] );
double class_matrix ( int kind, int m, double alpha, double beta, double aj[],
  double bj[] );
void imtqlx ( int n, double d[], double e[], double z[] );
void parchk ( int kind, int m, double alpha, double beta );
double r8_abs ( double x );
double r8_epsilon ( );
double r8_sign ( double x );
void r8mat_write ( std::string output_filename, int m, int n, double table[] );
void rule_write ( int order, std::string filename, double x[], double w[],
  double r[] );
void scqf ( int nt, double t[], int mlt[], double wts[], int nwts, int ndx[],
  double swts[], double st[], int kind, double alpha, double beta, double a,
  double b );
void sgqf ( int nt, double aj[], double bj[], double zemu, double t[],
  double wts[] );
void timestamp ( );

#endif