Difference between revisions of "Problem set 2.1 No.6"

From ccitonlinewiki
Jump to: navigation, search
(Created page with "File:Kuis.png Baris pertama ditukar dengan baris ketiga")
 
 
Line 1: Line 1:
[[File:Kuis.png]]
+
[[File:Hasil kuis.PNG]]
  
Baris pertama ditukar dengan baris ketiga
+
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

Hasil kuis.PNG

import numpy as np

  1. 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)

  1. 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)

  1. 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])