#include <complex.h>

extern void matprint(int m,int n,double A[m][n]);
extern void vecprint(int m,double X[m]);
extern void multAx(int m,int n,double A[m][n],
    double X[n],double B[m]);
extern void cmatprint(int m,int n,complex A[m][n]);
extern void cvecprint(int m,complex X[m]);
extern void cmultAx(int m,int n,complex A[m][n],
    complex X[n],complex B[m]);
void LUfact(int n,double A[n][n]);
void LUsolve(int n,double LU[n][n],
    double x[n],double b[n]);
void PLUsolve(int n,double LU[n][n],double *P[n],
    double x[n],double b[n]);
void PLUfact(int n,double A[n][n],double *P[n]);