Here's a simple program in Turbo C++ to perform various operations on matrices, including finding

the sum, product, transpose, trace, and inverse of matrices. We'll assume the matrices are 3x3 for
simplicity.

```cpp

#include <iostream>

#include <iomanip>

using namespace std;

const int N = 3;

// Function to input a matrix

void inputMatrix(double mat[N][N]) {

for (int i = 0; i < N; i++) {

for (int j = 0; j < N; j++) {

cin >> mat[i][j];

// Function to print a matrix

void printMatrix(double mat[N][N]) {

for (int i = 0; i < N; i++) {

for (int j = 0; j < N; j++) {

cout << setw(10) << mat[i][j] << " ";

cout << endl;

// Function to find the sum of two matrices

void sumMatrices(double mat1[N][N], double mat2[N][N], double result[N][N]) {

for (int j = 0; j < N; j++) {

result[i][j] = mat1[i][j] + mat2[i][j];

// Function to find the product of two matrices

void productMatrices(double mat1[N][N], double mat2[N][N], double result[N][N]) {

for (int i = 0; i < N; i++) {

for (int j = 0; j < N; j++) {

result[i][j] = 0;

for (int k = 0; k < N; k++) {

result[i][j] += mat1[i][k] * mat2[k][j];

// Function to find the transpose of a matrix

void transposeMatrix(double mat[N][N], double result[N][N]) {

for (int i = 0; i < N; i++) {

for (int j = 0; j < N; j++) {

result[j][i] = mat[i][j];

// Function to find the trace of a matrix

double traceMatrix(double mat[N][N]) {

double trace = 0;
 for (int i = 0; i < N; i++) {

trace += mat[i][i];

return trace;

// Function to find the inverse of a matrix

bool inverseMatrix(double mat[N][N], double result[N][N]) {

double det = mat[0][0] * (mat[1][1] * mat[2][2] - mat[1][2] * mat[2][1]) -

mat[0][1] * (mat[1][0] * mat[2][2] - mat[1][2] * mat[2][0]) +

mat[0][2] * (mat[1][0] * mat[2][1] - mat[1][1] * mat[2][0]);

if (det == 0) {

cout << "Matrix is singular, inverse does not exist." << endl;

return false;

double invDet = 1.0 / det;

result[0][0] = (mat[1][1] * mat[2][2] - mat[1][2] * mat[2][1]) * invDet;

result[0][1] = (mat[0][2] * mat[2][1] - mat[0][1] * mat[2][2]) * invDet;

result[0][2] = (mat[0][1] * mat[1][2] - mat[0][2] * mat[1][1]) * invDet;

result[1][0] = (mat[1][2] * mat[2][0] - mat[1][0] * mat[2][2]) * invDet;

result[1][1] = (mat[0][0] * mat[2][2] - mat[0][2] * mat[2][0]) * invDet;

result[1][2] = (mat[0][2] * mat[1][0] - mat[0][0] * mat[1][2]) * invDet;

result[2][0] = (mat[1][0] * mat[2][1] - mat[1][1] * mat[2][0]) * invDet;

result[2][1] = (mat[0][1] * mat[2][0] - mat[0][0] * mat[2][1]) * invDet;

result[2][2] = (mat[0][0] * mat[1][1] - mat[0][1] * mat[1][0]) * invDet;

return true;

}
int main() {

double mat1[N][N], mat2[N][N];

double result[N][N];

double trace;

cout << "Enter the elements of the first 3x3 matrix (row-wise):" << endl;

inputMatrix(mat1);

cout << "Enter the elements of the second 3x3 matrix (row-wise):" << endl;

inputMatrix(mat2);

// Sum of matrices

sumMatrices(mat1, mat2, result);

cout << "Sum of matrices:" << endl;

printMatrix(result);

// Product of matrices

productMatrices(mat1, mat2, result);

cout << "Product of matrices:" << endl;

printMatrix(result);

// Transpose of the first matrix

transposeMatrix(mat1, result);

cout << "Transpose of the first matrix:" << endl;

printMatrix(result);

// Trace of the first matrix

trace = traceMatrix(mat1);

cout << "Trace of the first matrix: " << trace << endl;
 // Inverse of the first matrix

if (inverseMatrix(mat1, result)) {

cout << "Inverse of the first matrix:" << endl;

printMatrix(result);

return 0;

```

### Explanation:

1. **Function Definitions**:

- `inputMatrix`: Takes input for a 3x3 matrix.

- `printMatrix`: Prints a 3x3 matrix.

- `sumMatrices`: Calculates the sum of two 3x3 matrices.

- `productMatrices`: Calculates the product of two 3x3 matrices.

- `transposeMatrix`: Calculates the transpose of a 3x3 matrix.

- `traceMatrix`: Calculates the trace (sum of diagonal elements) of a 3x3 matrix.

- `inverseMatrix`: Calculates the inverse of a 3x3 matrix using the formula for the inverse of a 3x3
matrix.

2. **Main Function**:

- Prompts the user to input the elements of two 3x3 matrices.

- Calculates and prints the sum, product, transpose, and trace of the matrices.

- Attempts to calculate and print the inverse of the first matrix, checking if it is singular.

This program covers basic matrix operations for a 3x3 matrix. For larger matrices or more complex
operations, consider using a library designed for linear algebra.
#include <iostream.h>

#include <iomanip.h>

#include <conio.h>

const int N = 3;

// Function to input a matrix

void inputMatrix(double mat[N][N]) {

for (int i = 0; i < N; i++) {

for (int j = 0; j < N; j++) {

cin >> mat[i][j];

// Function to print a matrix

void printMatrix(double mat[N][N]) {

for (int i = 0; i < N; i++) {

for (int j = 0; j < N; j++) {

cout << setw(10) << mat[i][j] << " ";

cout << endl;

// Function to find the sum of two matrices

void sumMatrices(double mat1[N][N], double mat2[N][N], double result[N][N]) {

for (int i = 0; i < N; i++) {

for (int j = 0; j < N; j++) {

result[i][j] = mat1[i][j] + mat2[i][j];

}
 }

// Function to find the product of two matrices

void productMatrices(double mat1[N][N], double mat2[N][N], double result[N][N]) {

for (int i = 0; i < N; i++) {

for (int j = 0; j < N; j++) {

result[i][j] = 0;

for (int k = 0; k < N; k++) {

result[i][j] += mat1[i][k] * mat2[k][j];

// Function to find the transpose of a matrix

void transposeMatrix(double mat[N][N], double result[N][N]) {

for (int i = 0; i < N; i++) {

for (int j = 0; j < N; j++) {

result[j][i] = mat[i][j];

// Function to find the trace of a matrix

double traceMatrix(double mat[N][N]) {

double trace = 0;

for (int i = 0; i < N; i++) {

trace += mat[i][i];

return trace;
}

// Function to find the inverse of a matrix

bool inverseMatrix(double mat[N][N], double result[N][N]) {

double det = mat[0][0] * (mat[1][1] * mat[2][2] - mat[1][2] * mat[2][1]) -

mat[0][1] * (mat[1][0] * mat[2][2] - mat[1][2] * mat[2][0]) +

mat[0][2] * (mat[1][0] * mat[2][1] - mat[1][1] * mat[2][0]);

if (det == 0) {

cout << "Matrix is singular, inverse does not exist." << endl;

return false;

double invDet = 1.0 / det;

result[0][0] = (mat[1][1] * mat[2][2] - mat[1][2] * mat[2][1]) * invDet;

result[0][1] = (mat[0][2] * mat[2][1] - mat[0][1] * mat[2][2]) * invDet;

result[0][2] = (mat[0][1] * mat[1][2] - mat[0][2] * mat[1][1]) * invDet;

result[1][0] = (mat[1][2] * mat[2][0] - mat[1][0] * mat[2][2]) * invDet;

result[1][1] = (mat[0][0] * mat[2][2] - mat[0][2] * mat[2][0]) * invDet;

result[1][2] = (mat[0][2] * mat[1][0] - mat[0][0] * mat[1][2]) * invDet;

result[2][0] = (mat[1][0] * mat[2][1] - mat[1][1] * mat[2][0]) * invDet;

result[2][1] = (mat[0][1] * mat[2][0] - mat[0][0] * mat[2][1]) * invDet;

result[2][2] = (mat[0][0] * mat[1][1] - mat[0][1] * mat[1][0]) * invDet;

return true;

int main() {

double mat1[N][N], mat2[N][N];

double result[N][N];
double trace;

cout << "Enter the elements of the first 3x3 matrix (row-wise):" << endl;

inputMatrix(mat1);

cout << "Enter the elements of the second 3x3 matrix (row-wise):" << endl;

inputMatrix(mat2);

// Sum of matrices

sumMatrices(mat1, mat2, result);

cout << "Sum of matrices:" << endl;

printMatrix(result);

// Product of matrices

productMatrices(mat1, mat2, result);

cout << "Product of matrices:" << endl;

printMatrix(result);

// Transpose of the first matrix

transposeMatrix(mat1, result);

cout << "Transpose of the first matrix:" << endl;

printMatrix(result);

// Trace of the first matrix

trace = traceMatrix(mat1);

cout << "Trace of the first matrix: " << trace << endl;

// Inverse of the first matrix

if (inverseMatrix(mat1, result)) {

cout << "Inverse of the first matrix:" << endl;

printMatrix(result);
 }

getch(); // To keep the console window open

return 0;