// Author:  Douglas Wilhelm Harder
// Copyright (c) 2009 by Douglas Wilhelm Harder.  All rights reserved.

#include "Sparse.h"
#include <iostream>
using namespace std;

int main() {
	Matrix<12, 10> A;

	for ( int i = 0; i < 12; ++i ) {
		for ( int j = 0; j < 10; ++j ) {
			if ( (i + j) % 11 == 0 ) {
				A.set( i, j, -1e10/double(i*i*i*i + j*j*j*j*j*j*j*j*j*j + 1) );
			} else {
				A.set( i, j,  1e10/double(i*i*i*i + j*j*j*j*j*j*j*j*j*j + 1) );
			}
		}
	}

	A.set(2,2,0);

	cout << A << endl;
	Matrix<10, 12> B = transpose( A );

	cout << B << endl;

	return 0;

	
}