Sunday, November 8, 2015

Integral menggunakan Metode Numerik di Python dan Plot Fungsi menggunakan Matplotlib.

 Kode berikut digunakan untuk mencari nilai integral dari a ke b sebuah fungsi f(x).

 Saya menggunakan modul Numpy untuk fungsi pangkat (power).

from pylab import *
import numpy as np

def f(x):
    y = np.power(x,3)+np.power(x,2)+x+1
    return y

dx = .1
a = 0.
b = 1.

x = a
c = 0.  # untuk menampung hasil integral
while x<=b:
    y  = f(x)
    c += f(x)*dx
    x += dx

print 'hasil integral dari ',a,' ke ',b,' adalah ' ,c

#untuk plotting
x = np.arange(-10., 10., .01)
y = f(x)

plot(x,y)
grid(True)
show()


..