Difference between revisions of "Problem set 2.1 No.6"
Dikaauliady (talk | contribs) (Created page with "File:Kuis.png Baris pertama ditukar dengan baris ketiga") |
Dikaauliady (talk | contribs) |
||
Line 1: | Line 1: | ||
− | [[File: | + | [[File:Hasil kuis.PNG]] |
− | + | import numpy as np | |
+ | |||
+ | #1. definisi matriks | ||
+ | A = np.array ([[1,2,0,-2,0], [0,1,0,2,-1], [0,0,2,1,2], [0,0,0,-1,1], [0,1,-1,1,-1]], float) | ||
+ | B = np.array ([1,1,-4,-2,-1], float) | ||
+ | |||
+ | n = len(A) | ||
+ | |||
+ | #2. eliminasi gauss | ||
+ | |||
+ | for k in range (0,n-1) : #k adalah baris pivot dimulai dari 0 sampai n-1 | ||
+ | for i in range (k+1, n) : #i adalah baris dibawah pivot sampai terakhir | ||
+ | if A[i,k]!=0 : | ||
+ | lam = A[i,k]/A[k,k] | ||
+ | A[i,k:n] = A[i,k:n]-(A[k,k:n]*lam) | ||
+ | B[i] = B[i]-(B[k]*lam) | ||
+ | |||
+ | print ('matriks A :', '\n', A) | ||
+ | |||
+ | #3. back substitution | ||
+ | x=np.zeros (n,float) #menyatakan kalo dia matrix | ||
+ | for m in range (n-1, -1, -1) : #m dari bawah ke atas | ||
+ | x[m] = (B[m]-np.dot(A[m, m+1:n], x[m+1:n]))/A[m,m] | ||
+ | print ('nilai X', m+1, '=', x[m]) |
Latest revision as of 21:48, 19 October 2019
import numpy as np
- 1. definisi matriks
A = np.array ([[1,2,0,-2,0], [0,1,0,2,-1], [0,0,2,1,2], [0,0,0,-1,1], [0,1,-1,1,-1]], float) B = np.array ([1,1,-4,-2,-1], float)
n = len(A)
- 2. eliminasi gauss
for k in range (0,n-1) : #k adalah baris pivot dimulai dari 0 sampai n-1
for i in range (k+1, n) : #i adalah baris dibawah pivot sampai terakhir if A[i,k]!=0 : lam = A[i,k]/A[k,k] A[i,k:n] = A[i,k:n]-(A[k,k:n]*lam) B[i] = B[i]-(B[k]*lam)
print ('matriks A :', '\n', A)
- 3. back substitution
x=np.zeros (n,float) #menyatakan kalo dia matrix for m in range (n-1, -1, -1) : #m dari bawah ke atas
x[m] = (B[m]-np.dot(A[m, m+1:n], x[m+1:n]))/A[m,m] print ('nilai X', m+1, '=', x[m])