From 7a88dbdd85d814df4faf950d544dfe3c3e3793d5 Mon Sep 17 00:00:00 2001
From: premsuryas <147473858+premsuryas@users.noreply.github.com>
Date: Tue, 26 Dec 2023 14:11:13 +0530
Subject: [PATCH] Update README.md

---
 README.md | 32 ++++++++++++++++++++++++++++----
 1 file changed, 28 insertions(+), 4 deletions(-)

diff --git a/README.md b/README.md
index 8852835d..63b4347a 100644
--- a/README.md
+++ b/README.md
@@ -14,16 +14,40 @@ To write a program to find the solution of a matrix using Gaussian Elimination.
 4. 
 
 ## Program:
-```
+
 /*
 Program to find the solution of a matrix using Gaussian Elimination.
-Developed by: 
-RegisterNumber: 
+Developed by: PREM KUMAR
+RegisterNumber: 23013598
 */
 ```
+import numpy as np
+import sys
+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][j]==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=' ')
 
+```
 ## Output:
-![gaussian elimination]()
+![output](/home/sec/Gaussian/Screenshot from 2023-12-26 13-57-34.png)
 
 
 ## Result: