double *lobatto_polynomial_derivative ( int m, int n, double x[] );
void lobatto_polynomial_derivatives ( int &n_data, int &n, double &x, 
  double &fx );
void lobatto_polynomial_plot ( int ndx_num, int ndx[], string prefix );
double *lobatto_polynomial_value ( int m, int n, double x[] );
void lobatto_polynomial_values ( int &n_data, int &n, double &x, double &fx );
double *r8vec_linspace_new ( int n, double a, double b );
void timestamp ( );