To write a program to find the solution of a matrix using Gaussian Elimination.
- Hardware โ PCs
- Anaconda โ Python 3.7 Installation / Moodle-Code Runner
1.Import the numpy module to use the built-in functions for calculation.
2.Import the sys module to use the built in functions.
3.Get input from the user for number of rows and add it by 1 for number of columns. 4.Using np.zeros() set the matrix as null marix.
5.Using for loop get input from the user for each element in the matrix.
6.Using for loop we can perform the elementary row operations and find the final matrix.
7.Print the values by solving the matrix using for loop with 2 point precision.
8.End the program.
Program to find the solution of a matrix using Gaussian Elimination.
Developed by: K.SANTHANA LAKSHMI
RegisterNumber: 212222240091
import numpy as np
import sys
n=int(input())
a=np.zeros((n,n+1))
x=np.zeros(n)
for i in range(n):
for j in range (n+1):
a[i][j]=float(input())
#print(a)
for i in range(n):
if a[i][j]==0.0:
sys. exit('Divide by zero found!')
for j in range(i+1,n):
scalar=a[j][i]/a[i][i]
for k in range(n+1):
a[j][k]=a[j][k]-scalar*a[i][k]
x[n-1]=a[n-1][n]/a[n-1][n-1]
for i in range(n-2, -1, -1):
x[i]=a[i][n]
for j in range(i+1,n):
x[i]=x[i]-a[i][j]*x[j]
x[i]=x[i]/a[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.