To write a program to find the solution of a matrix using Gaussian Elimination.
- Hardware โ PCs
- Anaconda โ Python 3.7 Installation / Moodle-Code Runner
- First,we want to import numpy,then import sys,assume a variable.
- For gaussian elimination method, we want to make 2nd and 3rd column zero.
- For that we want to make a range accorting to our program output.
- Then print the program with correct form then the output will display.
'''Program to solve a matrix using Gaussian elimination without partial pivoting.
Developed by: MANOJ KUMAR S
RegisterNumber: 23002959
'''
import sys
import numpy as np
n=int(input())
matrix=np.zeros((n,n+1))
x=np.zeros(n)
for i in range(n):
for j in range(n+1):
matrix[i][j]=int(input())
for i in range(n):
if matrix[i][i]==0.0:
sys.exit("Divide by zero error")
for j in range(i+1,n):
ratio=matrix[j][i]/matrix[i][i]
for k in range(n+1):
matrix[j][k]=matrix[j][k]-ratio*matrix[i][k]
x[n-1]=matrix[n-1][n]/matrix[n-1][n-1]
for i in range(n-2,-1,-1):
x[i]=matrix[i][n]
for j in range(i+1,n):
x[i]=x[i]-matrix[i][j]*x[j]
x[i]=x[i]/matrix[i][i]
for i in range(n):
print("X%d = %0.2f" %(i,x[i]),end=" ")
Thus the program to find the solution of a matrix using Gaussian Elimination is written and verified using python programming.