Kalkulus Dengan EMT
Kalkulus Dengan EMT
Kalkulus Dengan EMT
EMT (bersama Maxima) dapat digunakan untuk melakukan semua perhitungan di dalam kalkulus, baik secara numerik maupun
analitik (eksak).
Mendefinisikan Fungsi
Setiap format harus diawali dengan perintah function (bukan sebagai ekspresi).
1
4.31977682472
20.7392088022
>g(0)
>g(1)
2.20920171961
>g(f(5))
0.950898070639
file:///G:/Euler/EMT4Kalkulus.html 1/13
10/13/2020 Kalkulus dengan EMT
3
x x > 0
f (x) = {
2
x x ≤ 0.
Fungsi tersebut tidak dapat didefinisikan sebagai fungsi numerik secara "inline" menggunakan format :=, melainkan
didefinisikan sebagai program. Perhatikan, kata "map" digunakan agar fungsi dapat menerima vektor sebagai input, dan
hasilnya berupa vektor. Jika tanpa kata "map" fungsinya hanya dapat menerima input satu nilai.
>f(1)
>f(-2)
>f(-5:5)
>aspect(1.5); plot2d("f(x)",-5,5):
x
2 E
3 x + 1
file:///G:/Euler/EMT4Kalkulus.html 2/13
10/13/2020 Kalkulus dengan EMT
3 x + 1
2 E
Latihan
Bukalah buku Kalkulus. Cari dan pilih beberapa (paling sedikit 5 fungsi berbeda tipe/bentuk/jenis) fungsi dari buku
tersebut, kemudian definisikan di EMT pada baris-baris perintah berikut (jika perlu tambahkan lagi). Untuk setiap
fungsi, hitung beberapa nilainya, baik untuk satu nilai maupun vektor. Gambar grafik tersebut.
Juga, carilah fungsi beberapa (dua) variabel. Lakukan hal sama seperti di atas.
Menghitung Limit
Perhitungan limit pada EMT dapat dilakukan dengan menggunakan fungsi Maxima, yakni "limit". Fungsi "limit" dapat
digunakan untuk menghitung limit fungsi dalam bentuk ekspresi maupun fungsi yang sudah didefinisikan sebelumnya. Nilai
limit dapat dihitung pada sebarang nilai atau pada tak hingga (-inf, minf, dan inf). Limit kiri dan limit kanan juga
dapat dihitung, dengan cara memberi opsi "plus" atau "minus". Hasil limit dapat berupa nilai, "und' (tak definisi),
"ind" (tak tentu namun terbatas), "infinity" (kompleks tak hingga).
Perhatikan beberapa contoh berikut. Perhatikan cara menampilkan perhitungan secara lengkap, tidak hanya menampilkan
hasilnya saja.
>$showev('limit(1/(2*x-1),x,0))
1
lim = −1
x→0 2x − 1
>$showev('limit((x^2-3*x-10)/(x-5),x,5))
2
x − 3 x − 10
lim = 7
x→5 x − 5
>$showev('limit(sin(x)/x,x,0))
sin x
lim = 1
x→0 x
>plot2d("sin(x)/x",-pi,pi):
>$showev('limit(sin(x^3)/x,x,0))
3
sin x
lim = 0
x→0 x
>$showev('limit(log(x), x, minf))
file:///G:/Euler/EMT4Kalkulus.html 3/13
10/13/2020 Kalkulus dengan EMT
>$showev('limit((-2)^x,x, inf))
x
lim (−2) = infinity
x→∞
>$showev('limit(t-sqrt(2-t),t,2,minus))
−−−−
lim t − √2 − t = 2
t↑2
−−−− –
lim t − √2 − t = 5 − √3 i
t↓5
>plot2d("x-sqrt(2-x)",-2,5):
>$showev('limit((x^2-9)/(2*x^2-5*x-3),x,3))
2
x − 9 6
lim =
2
x→3 2x − 5x − 3 7
>$showev('limit((1-cos(x))/x,x,0))
1 − cos x
lim = 0
x→0 x
>$showev('limit((x^2+abs(x))/(x^2-abs(x)),x,0))
2
|x| + x
lim = −1
2
x→0 x − |x|
>$showev('limit((1+1/x)^x,x,inf))
x
1
lim ( + 1) = e
x→∞ x
>$showev('limit((1+k/x)^x,x,inf))
x
k
k
lim ( + 1) = e
x→∞ x
>$showev('limit((1+x)^(1/x),x,0))
lim (x + 1) x
= e
x→0
file:///G:/Euler/EMT4Kalkulus.html 4/13
10/13/2020 Kalkulus dengan EMT
>$showev('limit((x/(x+k))^x,x,inf))
x
x
−k
lim ( ) = e
x→∞ x + k
>$showev('limit(sin(1/x),x,0))
1
lim sin( ) = ind
x→0 x
>$showev('limit(sin(1/x),x,inf))
1
lim sin( ) = 0
x→∞ x
>plot2d("sin(1/x)",-5,5):
Latihan
Bukalah buku Kalkulus. Cari dan pilih beberapa (paling sedikit 5 fungsi berbeda tipe/bentuk/jenis) fungsi dari buku
tersebut, kemudian definisikan di EMT pada baris-baris perintah berikut (jika perlu tambahkan lagi). Untuk setiap
fungsi, hitung nilai limit fungsi tersebut di beberapa nilai dan di tak hingga. Gambar grafik fungsi tersebut untuk
mengkonfirmasi nilai-nilai limit tersebut.
Turunan Fungsi
Definisi turunan:
f (x + h) − f (x)
′
f (x) = lim
h→0 h
Berikut adalah contoh-contoh menentukan turunan fungsi dengan menggunakan definisi turunan (limit).
n n
(x + h) − x
n−1
lim = nx
h→0 h
Mengapa hasilnya seperti itu? Tuliskan atau tunjukkan bahwa hasil limit tersebut benar, sehingga benar turunan
fungsinya benar. Tulis penjelasan Anda di komentar ini.
sin(x + h) − sin x
lim = cos x
h→0 h
file:///G:/Euler/EMT4Kalkulus.html 5/13
10/13/2020 Kalkulus dengan EMT
Mengapa hasilnya seperti itu? Tuliskan atau tunjukkan bahwa hasil limit tersebut
benar, sehingga benar turunan fungsinya benar. Tulis penjelasan Anda di komentar ini.
Sebagai petunjuk, ekspansikan sin(x+h) dengan menggunakan rumus jumlah dua sudut.
log(x + h) − log x 1
lim =
h→0 h x
Mengapa hasilnya seperti itu? Tuliskan atau tunjukkan bahwa hasil limit tersebut
benar, sehingga benar turunan fungsinya benar. Tulis penjelasan Anda di komentar ini.
Sebagai petunjuk, gunakan sifat-sifat logaritma dan hasil limit pada bagian sebelumnya di atas.
1 1
−
x+h x 1
lim = −
2
h→0 h x
Use assume!
Error in:
$showev('limit((E^(x+h)-E^x)/h,h,0)) // turunan f(x)=e^x ...
^
x+h x
e − e
lim .
h→0 h
>$showev('limit((E^h-1)/h,h,0))
h
e − 1
lim = 1
h→0 h
>$factor(E^(x+h)-E^x)
h x
(e − 1) e
h
e − 1
x x
( lim ) e = e
h→0 h
x
x
x+h x
(x + h) − x
lim = infinity
h→0 h
file:///G:/Euler/EMT4Kalkulus.html 6/13
10/13/2020 Kalkulus dengan EMT
x+h x
(x + h) − x
lim .
h→0 h
x+h x
(x + h) − x
x
lim = x (log x + 1)
h→0 h
[x > 0]
>&forget(x<0)
[x < 0]
>&facts()
[]
arcsin(x + h) − arcsin x 1
lim =
−−−−−
h 2
h→0 √1 − x
tan(x + h) − tan x 1
lim =
2
h→0 h cos x
sinh(x)
−x 2 x
e (e + 1)
x −x
e + e
= cosh(x).
2
>plot2d(["f(x)","df(x)"],-pi,pi,color=[blue,red]):
file:///G:/Euler/EMT4Kalkulus.html 7/13
10/13/2020 Kalkulus dengan EMT
Latihan
Bukalah buku Kalkulus. Cari dan pilih beberapa (paling sedikit 5 fungsi berbeda tipe/bentuk/jenis) fungsi dari buku
tersebut, kemudian definisikan di EMT pada baris-baris perintah berikut (jika perlu tambahkan lagi). Untuk setiap
fungsi, tentukan turunannya dengan menggunakan definisi turunan (limit), seperti contoh-contoh tersebut. Gambar grafik
fungsi asli dan fungsi turunannya pada sumbu koordinat yang sama.
Integral
EMT dapat digunakan untuk menghitung integral, baik integral tak tentu maupun integral tentu. Untuk integral tak tentu
(simbolik) sudah tentu EMT menggunakan Maxima, sedangkan untuk perhitungan integral tentu EMT sudah menyediakan
beberapa fungsi yang mengimplementasikan algoritma kuadratur (perhitungan integral tentu menggunakan metode numerik).
Pada notebook ini akan ditunjukkan perhitungan integral tentu dengan menggunakan Teorema Dasar Kalkulus:
b
′
∫ f (x) dx = F (b) − F (a), dengan F (x) = f (x).
a
Fungsi untuk menentukan integral adalah integrate. Fungsi ini dapat digunakan untuk menentukan, baik integral tentu
maupun tak tentu (jika fungsinya memiliki antiderivatif). Untuk perhitungan integral tentu fungsi integrate
menggunakan metode numerik (kecuali fungsinya tidak integrabel, kita tidak akan menggunakan metode ini).
>$showev('integrate(x^n,x))
n+1
x
n
∫ x dx =
n + 1
>$showev('integrate(1/(1+x),x))
1
∫ dx = log(x + 1)
x + 1
>$showev('integrate(1/(1+x^2),x))
file:///G:/Euler/EMT4Kalkulus.html 8/13
10/13/2020 Kalkulus dengan EMT
1
∫ dx = arctan x
2
x + 1
>$showev('integrate(1/sqrt(1-x^2),x))
1
∫ dx = arcsin x
−−−−−
√1 − x2
>$showev('integrate(sin(x),x,0,pi))
∫ sin x dx = 2
0
>$showev('integrate(sin(x),x,a,b))
>$showev('integrate(x^n,x,a,b))
b n+1 n+1
b a
n
∫ x dx = −
n + 1 n + 1
a
>$showev('integrate(x^2*sqrt(2*x+1),x))
7 5 3
(2 x + 1) 2
(2 x + 1) 2
(2 x + 1) 2
2 −−−−−−
∫ x √2 x + 1 dx = − +
28 10 12
>$showev('integrate(x^2*sqrt(2*x+1),x,0,2))
5
2
−−−−−− 252 2
2
∫ x √2 x + 1 dx = −
21 105
0
>$ratsimp(%)
7
2
−−−−−− 252 − 2
2
∫ x √2 x + 1 dx =
105
0
>$showev('integrate((sin(sqrt(x)+a)*E^sqrt(x))/sqrt(x),x,0,pi^2))
π
2
−
− √x
sin(√x + a) e
π π
∫ dx = (−e − 1) sin a + (e + 1) cos a
−
−
0 √x
>$factor(%)
π
2
−
− √x
sin(√x + a) e
π
∫ dx = (−e − 1) (sin a − cos a)
−
−
0 √x
2
- x
E
>$showev('integrate(f(x),x))
−
−
2 √π erf (x)
−x
∫ e dx =
2
file:///G:/Euler/EMT4Kalkulus.html 9/13
10/13/2020 Kalkulus dengan EMT
2
−x
e
erf (x) = ∫ dx.
−
−
√π
Kita tidak dapat menggunakan teorema Dasar kalkulus untuk menghitung integral tentu fungsi tersebut jika semua
batasnya berhingga. Dalam hal ini dapat digunakan metode numerik (rumus kuadratur).
Integral tentu
π
2
−x
∫ e dx
0
dapat dihampiri dengan jumlah luas persegi-persegi panjang di bawah kurva y=f(x) tersebut. Langkah-langkahnya adalah
sebagai berikut.
Hasilnya adalah:
π
x
∫ x dx = 0.8362196102528469
0
Jumlah tersebut diperoleh dari hasil kali lebar sub-subinterval (=0.1) dan jumlah nilai-nilai f(x) untuk x = 0.1, 0.2,
0.3, ..., 3.2.
0.836219610253
Untuk mendapatkan nilai integral tentu yang mendekati nilai sebenarnya, lebar sub-intervalnya dapat diperkecil lagi,
sehingga daerah di bawah kurva tertutup semuanya, misalnya dapat digunakan lebar subinterval 0.001. (Silakan dicoba!)
Meskipun Maxima tidak dapat menghitung integral tentu fungsi tersebut untuk batas-batas yang berhingga, namun integral
tersebut dapat dihitung secara eksak jika batas-batasnya tak hingga. Ini adalah salah satu keajaiban di dalam
file:///G:/Euler/EMT4Kalkulus.html 10/13
10/13/2020 Kalkulus dengan EMT
matematika, yang terbatas tidak dapat dihitung secara eksak, namun yang tak hingga malah dapat dihitung secara eksak.
>$showev('integrate(f(x),x,0,inf))
∞ −
−
2 √π
−x
∫ e dx =
2
0
Berikut adalah contoh lain fungsi yang tidak memiliki antiderivatif, sehingga integral tentunya hanya dapat dihitung
dengan metode numerik.
x
x
>$showev('integrate(f(x),x,0,1))
1 1
x x
∫ x dx = ∫ x dx
0 0
Maxima gagal menghitung integral tentu tersebut secara langsung menggunakan perintah integrate. Berikut kita lakukan
seperti contoh sebelumnya untuk mendapat hasil atau pendekatan nilai integral tentu tersebut.
1
x
∫ x dx = 0.7834935879025506
0
Latihan
Bukalah buku Kalkulus.
Cari dan pilih beberapa (paling sedikit 5 fungsi berbeda tipe/bentuk/jenis) fungsi dari buku tersebut,
kemudian definisikan di EMT pada baris-baris perintah berikut (jika perlu tambahkan lagi).
file:///G:/Euler/EMT4Kalkulus.html 11/13
10/13/2020 Kalkulus dengan EMT
Untuk setiap fungsi, tentukan anti turunannya (jika ada), hitunglah integral tentu dengan batas-batas yang
menarik (Anda tentukan sendiri), seperti contoh-contoh tersebut.
Lakukan hal yang sama untuk fungsi-fungsi yang tidak dapat diintegralkan (cari sedikitnya 3 fungsi).
Gambar grafik fungsi dan daerah integrasinya pada sumbu koordinat yang sama.
Gunakan integral tentu untuk mencari luas daerah yang dibatasi oleh dua kurva yang berpotongan di dua
titik. (Cari dan gambar kedua kurva dan arsir (warnai) daerah yang dibatasi oleh keduanya.)
Gunakan integral tentu untuk menghitung volume benda putar kurva y= f(x) yang diputar mengelilingi sumbu x
dari x=a sampai x=b, yakni
b
2
V = ∫ π(f (x) dx.
a
(Pilih fungsinya dan gambar kurva dan benda putar yang dihasilkan. Anda dapat mencari contoh-contoh bagaimana cara
menggambar benda hasil perputaran suatu kurva.)
- Gunakan integral tentu untuk menghitung panjang kurva y=f(x) dari x=a sampai x=b dengan menggunakan rumus:
b −−−−−−−−−−
′ 2
S = ∫ √ 1 + (f (x)) dx.
a
dengan cara yang sama seperti mendefinisikan vektor dengan elemen-elemen beraturan (menggunakan titik dua
":");
menggunakan perintah "sequence" dan rumus barisan (suku ke -n);
menggunakan perintah "iterate" atau "niterate";
menggunakan fungsi Maxima "create_list" atau "makelist" untuk menghasilkan barisan simbolik;
menggunakan fungsi biasa yang inputnya vektor atau barisan;
menggunakan fungsi rekursif.
EMT juga dapat digunakan untuk menghitung jumlah deret berhingga maupun deret tak hingga, dengan menggunakan perintah
(fungsi) "sum". Perhitungan dapat dilakukan secara numerik maupun simbolik dan eksak.
Berikut adalah beberapa contoh perhitungan barisan dan deret menggunakan EMT.
[1, 2, 3, 4, 5, 6, 7, 8, 9, 10]
>1:2:30
[1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29]
>sum(1:2:30), sum(1/(1:2:30))
225
2.33587263431
n
n (n + 1)
∑k =
2
k=1
∞ ∞
1 1
∑ = ∑
k k
3 + k 3 + k
k=0 k=0
∞ 2
1 π
∑ =
2
x 6
x=1
∞ k−1 ∞ x
(−1) (−1)
∑ = −∑
k x
k=1 x=1
∞ k ∞ k
(−1) (−1)
∑ = ∑
2k − 1 2k − 1
k=1 k=1
>$ev(sum(1/n!, n, 0, inf),simpsum=true)
∞
1
∑
n!
n=0
∞
a
k
a ∑x =
1 − x
k=0
Deret Taylor
Deret Taylor suatu fungsi f yang diferensiabel sampai tak hingga di sekitar x=a adalah:
∞ k (k)
(x − a) f (a)
f (x) = ∑ .
k!
k=0
>$'e^x =taylor(exp(x),x,0,10) // deret Taylor e^x di sekitar x=0, sampai suku ke-11
10 9 8 7 6 5 4 3 2
x x x x x x x x x
x
e = + + + + + + + + + x + 1
3628800 362880 40320 5040 720 120 24 6 2
10 9 8 7 6 5 4 3 2
(x − 1) (x − 1) (x − 1) (x − 1) (x − 1) (x − 1) (x − 1) (x − 1) (x − 1)
log x = x − + − + − + − + − − 1
10 9 8 7 6 5 4 3 2
file:///G:/Euler/EMT4Kalkulus.html 13/13