Matrix inverse (inv) in SciPy - Step-by-Step Execution

Concept Flow - Matrix inverse (inv)

Start with square matrix A

↓

Check if A is invertible

↓

Error: No inverse

↓

Calculate inverse A_inv

↓

Verify: A * A_inv = Identity

↓

Use A_inv for solving or analysis

We start with a square matrix, check if it can be inverted, then compute its inverse, and optionally verify the result.

Execution Sample

SciPy

import numpy as np
from scipy.linalg import inv

A = np.array([[4, 7], [2, 6]])
A_inv = inv(A)
print(A_inv)

This code computes the inverse of a 2x2 matrix A using scipy.linalg.inv and prints the result.

Execution Table

Step	Action	Matrix A	Check invertible	Result	Output
1	Define matrix A	[[4, 7], [2, 6]]	N/A	N/A	Matrix A created
2	Calculate determinant	[[4, 7], [2, 6]]	det = 4*6 - 7*2 = 24 - 14 = 10	det != 0	Matrix is invertible
3	Compute inverse using inv(A)	[[4, 7], [2, 6]]	N/A	Inverse calculated	[[0.6, -0.7], [-0.2, 0.4]]
4	Verify A * A_inv	[[4, 7], [2, 6]]	N/A	Product is identity matrix	[[1, 0], [0, 1]]
5	End	N/A	N/A	N/A	Inverse matrix ready for use

💡 Matrix is invertible because determinant is non-zero; inverse computed successfully.

Variable Tracker

Variable	Start	After Step 1	After Step 3	After Step 4	Final
A	undefined	[[4, 7], [2, 6]]	[[4, 7], [2, 6]]	[[4, 7], [2, 6]]	[[4, 7], [2, 6]]
det	undefined	undefined	10	10	10
A_inv	undefined	undefined	[[0.6, -0.7], [-0.2, 0.4]]	[[0.6, -0.7], [-0.2, 0.4]]	[[0.6, -0.7], [-0.2, 0.4]]
product	undefined	undefined	undefined	[[1, 0], [0, 1]]	[[1, 0], [0, 1]]

Key Moments - 3 Insights

Why do we check the determinant before calculating the inverse?

Because if the determinant is zero (see Step 2 in execution_table), the matrix is not invertible and inv() will fail.

What does the product of A and A_inv represent?

It should be the identity matrix (Step 4), confirming that A_inv is the correct inverse.

Can we invert a non-square matrix using inv()?

No, inv() requires a square matrix; otherwise, it will raise an error.

Visual Quiz - 3 Questions

Look at the execution_table, what is the determinant of matrix A at Step 2?

A 24

B 10

C 0

D 14

💡 Hint

Check the 'Check invertible' column at Step 2 in execution_table.

At which step does the code verify that the inverse is correct?

A Step 4

B Step 3

C Step 2

D Step 5

💡 Hint

Look for the step where the product A * A_inv is checked in execution_table.

If the determinant was zero, what would happen in the execution flow?

A The inverse would be calculated anyway

B The identity matrix would be returned

C An error would occur, stopping the process

D The matrix would be transposed

💡 Hint

Refer to the concept_flow where a non-invertible matrix leads to an error.

Concept Snapshot

Matrix inverse (inv):
- Use scipy.linalg.inv(matrix) for square matrices
- Check determinant != 0 before inversion
- Result satisfies: matrix * inverse = identity
- Raises error if matrix is singular or not square
- Useful for solving linear systems and transformations

Full Transcript

This visual execution traces how to compute the inverse of a square matrix using scipy.linalg.inv. We start by defining matrix A, then check if it is invertible by calculating its determinant. Since the determinant is 10 (non-zero), we proceed to compute the inverse matrix. We verify correctness by multiplying A with its inverse, which yields the identity matrix. The inverse matrix is then ready for use. Key points include the necessity of a non-zero determinant and that only square matrices can be inverted. If the determinant were zero, the process would stop with an error.