/* simeq_plus.h    matrices are n columns of rows with n elements [A] */
/*                 vectors  are n elements  |X| */
#include <stdlib.h>
#include <string.h>
#include <stdio.h>
#include <math.h>

void simeq(int n, double A[], double Y[], double X[]);
           /* solve linear system of equations  [A] * |X| = |Y|  */
double determinant(int n, double A[]);
                   /* D = det([A]) */
void inverse(int n, double A[], double AI[]);
             /* [AI] = [A]^-1 */
void mat_mul(int n, double A[], double B[], double C[]);
             /* |C| = [A]*|B| */
void mat_vec_mul(int n, double A[], double B[], double C[]);
                 /* |C| = [A]*|B| */
void mat_add(int n, double A[], double B[], double C[]);
             /* [C] = [A]+[B] */
void mat_sub(int n, double A[], double B[], double C[]);
             /* [C] = [A]-[B] */
void vec_add(int n, double A[], double B[], double C[]);
             /* |C| = |A|+|B| */
void vec_sub(int n, double A[], double B[], double C[]);
             /* |C| = |A|-|B| */
void mat_put(int n, double A[]);
void vec_put(int n, double A[]);

/* end simeq_plus.h */