void sierpinski_carpet_chaos ( int n );
void r8mat_mv ( int m, int n, double a[], double x[], double ax[] );