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[], char *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 ( );