Merge pull request Python-World#527 from Sergio-Torres/Solver_linear_equations

Solver linear equations

Author: chavarera

projects/Solver_linear_equations/README.md

@@ -0,0 +1,25 @@
+<h1>Solver linear equations</h1>
+
+<p>Enter values through spaces for each of the three equations, return the solution and a 3D graph.</p>
+
+<h3>Prerequisites</h3>
+<ul>
+  <li>Numpy</li>
+  <li>Matplotlib</li>
+</ul>
+
+<h3>Install prerequisites</h3>
+<ul>
+  <li>pip install numpy</li>
+  <li>pip install matplotlib</li>
+</ul>
+
+<h3>How to run the script</h3>
+<p>python3 linearEquations.py</p>
+
+<h3>Screenshot of the script running</h3>
+
+![screenshot1](https://user-images.githubusercontent.com/76462037/141659586-13865df9-f853-4dd5-9bae-a2e32f9ec8ef.jpeg)
+![screenshot2](https://user-images.githubusercontent.com/76462037/141659827-c4bd9b90-5f45-4eba-9f3c-f5b679c07c75.png)
+
+<h3><b>Author Name: </b></h3> <a href="https://github.com/Sergio-Torres">Sergio-Torres</a>

projects/Solver_linear_equations/linearEquations.py

@@ -0,0 +1,44 @@
+import numpy as np
+import matplotlib.pyplot as plt
+from mpl_toolkits.mplot3d import Axes3D
+from matplotlib import cm
+
+print('The 3 equations are entered individually, each value of the equation is entered separated by a space, for example: \ninput = 6 5 -3 4 \nThis will be equal to 6x + 5y- 3z = 4')
+
+print('Enter values for equation 1: ')
+a, b, c, d = map(float, input().split()) 
+
+print('Enter values for equation 2: ')
+e, f, g, h = map(float, input().split())
+
+print('Enter values for equation 3: ')
+i, j, k, l = map(float, input().split())
+#solve the linear equation
+A = np.array([[a,b,c],[e,f,g],[i,j,k]])
+b_a = np.array([d,h,l])
+sol = np.linalg.solve(A,b_a)
+print(sol)
+
+
+x,y = np.linspace(0,10,10), np.linspace(0,10,10)
+X,Y = np.meshgrid(x,y)
+
+Z1 = (d-a*X-b*Y)/c
+Z2 =(h-e*X-f*Y)/g
+Z3 =(l+X-Y)/k
+
+#create 3d graphics
+
+fig = plt.figure()
+ax = fig.add_subplot(111,projection='3d')
+ax.plot_surface(X,Y,Z1,alpha=0.5,cmap=cm.Accent,rstride=100,cstride=100)
+ax.plot_surface(X,Y,Z2,alpha=0.5,cmap=cm.Paired,rstride=100,cstride=100)
+ax.plot_surface(X,Y,Z3,alpha=0.5,cmap=cm.Pastel1,rstride=100,cstride=100)
+ax.plot((sol[0],),(sol[1],),(sol[2],),lw=2,c='k', marker='o', markersize=7, markeredgecolor='g', markerfacecolor='white')
+ax.set_xlabel('X axis')
+ax.set_ylabel('Y axis')
+ax.set_zlabel('Z axis')
+
+plt.show()
+
+