Difference between revisions of "Kelompok 4"
(→TUGAS 2: Gauss Jordan dalam Penyelesaian Statika Struktur) |
|||
Line 29: | Line 29: | ||
Pada tugas ini kami menyelesaikan suatu persamaan Statika Struktur dengan metode Gauss Jordan menggunakan python. | Pada tugas ini kami menyelesaikan suatu persamaan Statika Struktur dengan metode Gauss Jordan menggunakan python. | ||
+ | |||
+ | Kami mengambil contoh soal | ||
+ | |||
+ | |||
+ | |||
+ | import numpy as np | ||
+ | |||
+ | a=np.array ([[1,1], | ||
+ | [7,-5]]) | ||
+ | |||
+ | b=np.array ([[1500], | ||
+ | [0]]) | ||
+ | |||
+ | c=[1,1,1,1,1,1,1,1] | ||
+ | n=len(b) | ||
+ | |||
+ | for k in range(0,n-1): | ||
+ | for i in range(n-1,k,-1): | ||
+ | if a[i][k]!=0.0: | ||
+ | op=a[i][k]/a[i-1][k] | ||
+ | b[i][0]=b[i][0]-op*b[i-1][0] | ||
+ | for f in range (0,n): | ||
+ | a[i][f]=a[i][f]-op*a[i-1][f] | ||
+ | |||
+ | |||
+ | if a[n-1][n-1]!=0: | ||
+ | c[n-1]=b[n-1][0]/a[n-1][n-1] | ||
+ | else: | ||
+ | c[n-1]=0 | ||
+ | for i in range(n-2,-1,-1): | ||
+ | sigma=0 | ||
+ | for k in range(0,n,1): | ||
+ | y=a[i][k]*c[k] | ||
+ | sigma=sigma+y | ||
+ | sigma=sigma-a[i][i] | ||
+ | if a[i][i]!=0: | ||
+ | c[i]=(b[i][0]-sigma)/a[i][i] | ||
+ | else: | ||
+ | a[i][i]=0 | ||
+ | for v in range(0,n,1): | ||
+ | print(a[v]) | ||
+ | print(c[0:n]) |
Revision as of 21:41, 15 October 2019
TUGAS KELOMPOK 4
Bismillahirrahmanirrahim.
PERKENALAN
Perkenalkan kami dari kelompok 4 yang beranggotakan:
1. Satria Nur Muhammad
2. Stefanus Nandiwardhana
3. Rafiq Ali Abdillah
TUGAS 1: Gauss Jordan dalam Penyelesaian Persamaan Aljabar Linear
Pada tugas ini kami menyelesaikan suatu persamaan Aljabar Linear dengan metode Gauss Jordan menggunakan python.
TUGAS 2: Gauss Jordan dalam Penyelesaian Statika Struktur
Pada tugas ini kami menyelesaikan suatu persamaan Statika Struktur dengan metode Gauss Jordan menggunakan python.
Kami mengambil contoh soal
import numpy as np a=np.array ([[1,1], [7,-5]]) b=np.array ([[1500], [0]]) c=[1,1,1,1,1,1,1,1] n=len(b) for k in range(0,n-1): for i in range(n-1,k,-1): if a[i][k]!=0.0: op=a[i][k]/a[i-1][k] b[i][0]=b[i][0]-op*b[i-1][0] for f in range (0,n): a[i][f]=a[i][f]-op*a[i-1][f]
if a[n-1][n-1]!=0: c[n-1]=b[n-1][0]/a[n-1][n-1] else: c[n-1]=0 for i in range(n-2,-1,-1): sigma=0 for k in range(0,n,1): y=a[i][k]*c[k] sigma=sigma+y sigma=sigma-a[i][i] if a[i][i]!=0: c[i]=(b[i][0]-sigma)/a[i][i] else: a[i][i]=0 for v in range(0,n,1): print(a[v]) print(c[0:n])