int i4_log_10 ( int i );
int i4_max ( int i1, int i2 );
int i4_min ( int i1, int i2 );
int i4_power ( int i, int j );
int i4_uniform_ab ( int a, int b, int *seed );

double r8_uniform_01 ( int *seed );

double *r83_indicator ( int m, int n );
void r83_print ( int m, int n, double a[], char *title );
void r83_print_some ( int m, int n, double a[], int ilo, int jlo, int ihi, 
  int jhi, char *title );

void r83v_cg ( int n, double a1[], double a2[], double a3[], double b[], 
  double x[] );
void r83v_copy ( int m, int n, double a1[], double a2[], double a3[], 
  double b1[], double b2[], double b3[] );
double *r83v_cr_fa ( int n, double a1[], double a2[], double a3[] );
double *r83v_cr_sl ( int n, double a_cr[], double b[] );
double *r83v_cr_sls ( int n, double a_cr[], int nb, double b[] );
void r83v_dif2 ( int m, int n, double a[], double b[], double c[] );
double *r83v_fs ( int n, double a1[], double a2[], double a3[], double b[] );
void r83v_gs_sl ( int n, double a1[], double a2[], double a3[], double b[], 
  double x[], int it_max );
void r83v_indicator ( int m, int n, double a[], double b[], double c[] );
void r83v_jac_sl ( int n, double a1[], double a2[], double a3[], double b[], 
  double x[], int it_max );
double *r83v_mtv ( int m, int n, double a1[], double a2[], double a3[],
  double x[] );
double *r83v_mv ( int m, int n, double a[], double b[], double c[], double x[] );
void r83v_print ( int m, int n, double a[], double b[], double c[], 
  char *title );
void r83v_print_some ( int m, int n, double a[], double b[], double c[],
  int ilo, int jlo, int ihi, int jhi, char *title );
void r83v_random ( int m, int n, int *seed, double a[], double b[], double c[] );
double *r83v_res ( int m, int n, double a[], double b[], double c[],
  double x[], double ax[] );
double *r83v_to_r8ge ( int m, int n, double a1[], double a2[], double a3[] );
double *r83v_to_r8vec ( int m, int n, double a1[], double a2[], double a3[] );
void r83v_transpose ( int m, int n, double a1[], double a2[], double a3[], 
  double b1[], double b2[], double b3[] );
void r83v_zeros ( int m, int n, double a[], double b[], double c[] );

double *r8ge_indicator ( int m, int n );
double *r8ge_mtv ( int m, int n, double a[], double x[] );
double *r8ge_mv ( int m, int n, double a[], double x[] );
void r8ge_print ( int m, int n, double a[], char *title );
void r8ge_print_some ( int m, int n, double a[], int ilo, int jlo, int ihi, 
  int jhi, char *title );
void r8ge_to_r83v ( int m, int n, double a[], double a1[], double a2[], 
  double a3[] );

double r8vec_dot_product ( int n, double a1[], double a2[] );
double *r8vec_indicator1_new ( int n );
double r8vec_norm ( int n, double a[] );
double r8vec_norm_affine ( int n, double v0[], double v1[] );
void r8vec_print ( int n, double a[], char *title );
void r8vec_to_r83v ( int m, int n, double a[], double a1[], double a2[], 
  double a3[] );
void r8vec_uniform_01 ( int n, int *seed, double x[] );
double *r8vec_uniform_01_new ( int n, int *seed );

void r8vec2_print ( int n, double a1[], double a2[], char *title );