Program Metode Tabulasi Menggunakan Python

Assalamualaikum Wr. Wb. Pada kesempatan kali ini saya akan mencoba berbagi hasil pekerjaan saya ketika mencoba memecahkan soal pada postingan penyelesain persamaan tak linear menggunakan metode tabulasi. Jadi postingan kali ini membahas penggunaan bahasa pemrograman python untuk kasus metode numerik di bidang penyelesaian persamaan tak linear menggunakan metode tabulasi. Poin-poin yang dibahas adalah soal, code, dan hasil running program.

Soal

Cari satu nilai x pada f(x)=x3 – 7x + 1, sehingga f(x) mendekati nilai 0, 

• batas presisi 10-7, dimana x berada diantara 0.1 dan 0.9.
• buatlah programnya dan tampilkan hasilnya.

Code

 

# Sufyan Saori, Numeric Method's, 
# "Tabulation Method" for Problem solving Nonlinear equations 
# Algotihm : input (x1 and x2). check if f(x1)*f(x2)<0. 
# email : xsufyan@gmail.com, blogspot:sufyan97
import os

def f(x):
 # function for f(x). it regulary change when case change too.
 y = x**3 - 7*x + 1
 return y

def inputFirst(x1,x2): 
 # function for input.
 x1 = float(input('Enter x1 : '))
 x2 = float(input('Enter x2 : '))
 return (x1,x2)

def check1(x1,x2):
 # function to check if x1 and x2 are qualified to continue.
 times = f(x1)*f(x2)
 if(times<0):
  return (True)
 else:
  os.system('cls')
  print("Nothing found when x1=",x1," and x2= ",x2);
  return(False)
 

# def tesRun1(x1,x2):
#  check1(x1,x2)
# x1 = 0.0
# x2 = 0.0
# x1,x2 = inputFirst(x1,x2)
# print (f(x1),f(x2))
# tesRun1(x1,x2)


def update1(x1,x2):
 step = (x2-x1)/10.0
 x2 = x1+step
 if(check1):
  while((f(x2)*f(x1))>0):
   temp = x2
   x2 = x2+step
   x1 = temp
  return x1,x2
 

def checkStop(x1,x2,tolerance):
 while (abs(f(x1))>tolerance and abs(f(x2))>tolerance):
  x1,x2 = update1(x1,x2)
  ##print ("masuk sini")
 if(f(x1)>tolerance):
  return x2
 else:
  return x1

# print (update1(0,1))
# print (f(0),f(0.1),f(0.2))

def main():
 tolerance = 0.00000001
 x1 = 0.0
 x2 = 0.0
 x1,x2 = inputFirst(x1,x2)
 if(check1(x1,x2)):
  x1,x2=update1(x1,x2)
  result = checkStop(x1,x2,tolerance)
  print (result," -- ", f(result))
 else:
  print ("Doesn't work")
 
main()

Hasil Run Program

Apa yang harus dilakukan :

• jalankan program python. 
• buka folder dimana file disimpan.
• pada addressnya ketik cmd lalu enter. maka akan muncul command prompt.
  
• pada command prompt ketikan python namafile.py
• maka program akan jalan.
• Masukan inputan.
• x1 = 0.1
• x2 = 0.9
•  Maka akan muncul hasil, di sebelah kanan merupakan nilai x yang dicari. sebelah kiri merupakan nilai f(x) yang dicari.

Sekian kurangnya mohon maaf. Mohon hargai code orang lain sertakan sumber ketika hendak digunakan lagi. Terimakasih, Wassalamuaalaikum Wr. Wb.
"Keep Learning Alon-alon asal kelakon, Best regards : sufyan97".

Lokasi:

Berbagi :