Python （ NumPy ）

/**********************/
/* クラスMatrixの定義 */
/**********************/
class Matrix {       /* ２次元行列 */
		int n;              // 行の数
		int m;              // 列の数
	public:
		double **mat;       // 行列本体
		Matrix(int, int, double **);   // コンストラクタ（引数 3）
		Matrix(int, int);   // コンストラクタ（引数 2）
		Matrix(const Matrix &);   // 初期化のためのコンストラクタ
		Matrix() {n = 0;}   // コンストラクタ（引数無し）
		~Matrix()   // デストラクタ
		{
			if (n > 0) {
				int i1;
				for (i1 = 0; i1 < n; i1++)
					delete [] mat[i1];
				delete [] mat;
			}
		}
		Matrix &operator= (const Matrix &);   // ＝のオーバーロード
		Matrix operator+ (const Matrix &);          // ＋のオーバーロード
		Matrix operator- (const Matrix &);          // －のオーバーロード
		Matrix operator* (const Matrix &);          // ＊のオーバーロード
	friend ostream &operator << (ostream &, Matrix);     // << のオーバーロード
	friend Matrix sqrt(const Matrix &);     // 各要素の平方根
};

/****************************/
/* コンストラクタ（引数 3） */
/****************************/
Matrix::Matrix(int n1, int m1, double **x)
{
	n    = n1;
	m    = m1;
	mat  = new double * [n];
	for (int i1 = 0; i1 < n; i1++) {
		mat[i1] = new double [m];
		for (int i2 = 0; i2 < m; i2++)
			mat[i1][i2] = x[i1][i2];
	}
}

/****************************/
/* コンストラクタ（引数 2） */
/****************************/
Matrix::Matrix(int n1, int m1)
{
	n    = n1;
	m    = m1;
	mat  = new double * [n];
	for (int i1 = 0; i1 < n; i1++)
		mat[i1] = new double [m];
}

/********************************/
/* 初期化のためのコンストラクタ */
/********************************/
Matrix::Matrix(const Matrix &A)
{
	n    = A.n;
	m    = A.m;
	mat  = new double * [n];
	for (int i1 = 0; i1 < n; i1++) {
		mat[i1] = new double [m];
		for (int i2 = 0; i2 < m; i2++)
			mat[i1][i2] = A.mat[i1][i2];
	}
}

/********************************/
/* 出力（ << のオーバーロード） */
/********************************/
ostream &operator << (ostream &stream, Matrix A)
{
	for (int i1 = 0; i1 < A.n; i1++) {
		for (int i2 = 0; i2 < A.m; i2++)
			stream << " " << A.mat[i1][i2];
		stream << "\n";
	}

	return stream;
}

/*******************************/
/* ＝のオーバーロード          */
/*      A = B (A.operator=(B)) */
/*******************************/
Matrix& Matrix::operator= (const Matrix &B)
{
	if (&B != this) {   // 自分自身への代入を防ぐ
		if (n > 0) {   // 代入する前のメモリを解放
			for (int i1 = 0; i1 < n; i1++)
				delete [] mat[i1];
			delete [] mat;
		}
		n        = B.n;
		m        = B.m;
		mat  = new double * [n];   // メモリの確保と代入
		for (int i1 = 0; i1 < n; i1++) {
			mat[i1] = new double [m];
			for (int i2 = 0; i2 < m; i2++)
				mat[i1][i2] = B.mat[i1][i2];
		}
	}

	return *this;
}

/**********************/
/* ＋のオーバーロード */
/**********************/
Matrix Matrix::operator + (const Matrix &B)
{
	if (n != B.n || m != B.m) {
		cout << "***error  invalid data\n";
		exit(1);
	}

	Matrix C(n, m);

	for (int i1 = 0; i1 < n; i1++) {
		for (int i2 = 0; i2 < m; i2++)
			C.mat[i1][i2] = mat[i1][i2] + B.mat[i1][i2];
	}

	return C;
}

/**********************/
/* －のオーバーロード */
/**********************/
Matrix Matrix::operator - (const Matrix &B)
{
	if (n != B.n || m != B.m) {
		cout << "***error  invalid data\n";
		exit(1);
	}

	Matrix C(n, m);

	for (int i1 = 0; i1 < n; i1++) {
		for (int i2 = 0; i2 < m; i2++)
			C.mat[i1][i2] = mat[i1][i2] - B.mat[i1][i2];
	}

	return C;
}

/**********************/
/* ＊のオーバーロード */
/**********************/
Matrix Matrix::operator* (const Matrix &B)
{
	if (m != B.n) {
		cout << "***error  invalid data\n";
		exit(1);
	}

	Matrix C(n, B.m);

	for (int i1 = 0; i1 < n; i1++) {
		for (int i2 = 0; i2 < B.m; i2++) {
			C.mat[i1][i2] = 0.0;
			for (int i3 = 0; i3 < m; i3++)
				C.mat[i1][i2] += mat[i1][i3] * B.mat[i3][i2];
		}
	}

	return C;
}

/******************/
/* 各要素の平方根 */
/******************/
Matrix sqrt(const Matrix &A)
{
	Matrix C(A.n, A.m);

	for (int i1 = 0; i1 < A.n; i1++) {
		for (int i2 = 0; i2 < A.m; i2++)
			C.mat[i1][i2] = sqrt(A.mat[i1][i2]);
	}

	return C;
}