Untuk akar: periksa pergantian tanda pada bisection, atau hitung f dan f′ untuk Newton, dan pertahankan enam angka atau lebih pada tabel iterasi.
Untuk nilai tabel: interpolasikan secara linear di antara dua titik yang bertetangga.
Untuk integral: cari h = (b − a)/n, daftarkan f di setiap simpul, terapkan bobot (1, 2, …, 2, 1 untuk trapesium; 1, 4, 2, …, 4, 1 untuk Simpson dengan n genap) dan kalikan dengan h/2 atau h/3.
Untuk PD: melangkahlah dari t₀ dengan yₙ₊₁ = yₙ + h·f(tₙ, yₙ) (Euler) atau kemiringan rata-rata (Heun), satu baris per langkah, sampai mencapai t yang diminta.
2Soal latihan
1.Persamaan diferensialSedang
Sebuah komponen baja tempa bersuhu 340 °C dibiarkan mendingin di bengkel bersuhu 23 °C. Setelah 14 menit suhunya 204 °C. Benda tersebut mengikuti hukum pendinginan Newton, dT/dt = −k·(T − Tₐ). Berapa lama setelah awal, t, yang diperlukan untuk mencapai 123 °C?
T0 = 340 °CTa = 23 °CT(14) = 204 °CT = 123 °C
t = min
2.Persamaan diferensialSulit
Selesaikan masalah nilai awal y′ + 2y/x = 3 untuk x > 0, dengan y(1) = 1. Cari y(3).
y′ + 2y/x = 3y(1) = 1
y(3) =
3.IntegralSedang
Hitung I = ∫₀1 x·cos x dx (x dalam radian).
∫₀1 x·cos x dx
I =
4.IntegralSulit
Hitung I = ∫₀1 x²·eˣ dx.
∫₀1 x²eˣ dx
I =
5.TurunanSedang
Misalkan f(x) = ln(2x2 + 4). Tentukan f′(−2).
f(x) = ln(2x2 + 4)x = −2
f′(−2) =
6.TurunanSulit
Garis singgung f(x) = 3x2 + 3x − 6 di x = 3 memotong sumbu y di y = b. Tentukan b.
f(x) = 3x2 + 3x − 6x = 3
b =
7.Metode numerikSedang
Hampiri ∫02 1/(1 + x³) dx dengan kaidah Simpson memakai n = 6 pita. Berikan S.
f(x) = 1/(1 + x³)[a, b] = [0, 2]n = 6
S =
8.Metode numerikSulit
Selesaikan y′ = y(2 − y), y(0) = 0,2 dengan metode Euler dan langkah h = 0,5 sampai t = 2,5. Solusi eksaknya y = 2y₀/(y₀ + (2 − y₀)e−2t). Carilah galat global |E| = |y(2,5) − y5|.
y′ = y(2 − y)y(0) = 0,2h = 0,5
|E| =
3Larutan
1.Persamaan diferensial
Jawaban: t = 28,82 min
Persamaan diferensial orde satu: dy/dt = k·y memiliki penyelesaian y = y₀·ekt (pertumbuhan untuk k > 0, peluruhan untuk k < 0); waktu paruh T memberikan k = ln 2/T. Hukum pendinginan Newton, dT/dt = −k·(T − Tₐ), memberikan T(t) = Tₐ + (T₀ − Tₐ)·e−kt. Persamaan linear y′ + p(x)·y = q(x) diselesaikan dengan faktor integrasi μ = e∫p dx: (μy)′ = μq, sehingga μy = ∫μq dx + C.
Konstanta pendinginan
k = ln(317/181)/14 = 0,04003 min⁻¹
Waktu untuk mencapai target
t = ln(317/100)/0,04003 = 28,82 min
Jawaban
t = 28,82 min
2.Persamaan diferensial
Jawaban: y(3) = 3
Persamaan diferensial orde satu: dy/dt = k·y memiliki penyelesaian y = y₀·ekt (pertumbuhan untuk k > 0, peluruhan untuk k < 0); waktu paruh T memberikan k = ln 2/T. Hukum pendinginan Newton, dT/dt = −k·(T − Tₐ), memberikan T(t) = Tₐ + (T₀ − Tₐ)·e−kt. Persamaan linear y′ + p(x)·y = q(x) diselesaikan dengan faktor integrasi μ = e∫p dx: (μy)′ = μq, sehingga μy = ∫μq dx + C.
Faktor integrasi
μ = e∫2/x dx = x2
Integralkan
(x2y)′ = 3x2 → x2y = 3x3/3 + C
Kondisi awal
1 = 3/3 + C → C = 0
Hitung nilainya
y(3) = (3·33/3 + 0)/32 = 3
Jawaban
y(3) = 3
3.Integral
Jawaban: I = 0,3818
Integrasi parsial: ∫u·dv = u·v − ∫v·du, dari aturan perkalian. Pilih u sebagai faktor yang menjadi lebih sederhana saat diturunkan (pangkat x, atau ln x), dan dv sebagai faktor yang dapat Anda integralkan (eˣ, sin x, cos x, atau x bila u = ln x). Untuk integral tentu, hitung u·v pada kedua batas: ∫ₐᵇ u dv = [u·v]ₐᵇ − ∫ₐᵇ v du. Faktor x² memerlukan dua putaran.
Pilih u dan dv
u = x, du = dx; dv = cos x dx, v = sin x
Antiturunan
∫x cos x dx = x sin x − ∫sin x dx = x sin x + cos x
Limit
(1·sin 1 + cos 1) − (0 + 1) = 0,3818
Jawaban
I = 0,3818
4.Integral
Jawaban: I = 0,7183
Integrasi parsial: ∫u·dv = u·v − ∫v·du, dari aturan perkalian. Pilih u sebagai faktor yang menjadi lebih sederhana saat diturunkan (pangkat x, atau ln x), dan dv sebagai faktor yang dapat Anda integralkan (eˣ, sin x, cos x, atau x bila u = ln x). Untuk integral tentu, hitung u·v pada kedua batas: ∫ₐᵇ u dv = [u·v]ₐᵇ − ∫ₐᵇ v du. Faktor x² memerlukan dua putaran.
Putaran pertama
∫x²eˣ dx = x²eˣ − ∫2x eˣ dx
Putaran kedua
∫2x eˣ dx = 2x eˣ − 2eˣ
∫x²eˣ dx = eˣ(x² − 2x + 2)
Limit
2,718·(1) − 1·2 = 0,7183
Jawaban
I = 0,7183
5.Turunan
Jawaban: f′(−2) = −0,6667
Diferensial: turunan f′(x) adalah kemiringan grafik di x. Aturan pangkat: (xⁿ)′ = n·xⁿ⁻¹. Aturan perkalian: (uv)′ = u′v + uv′. Aturan rantai: (f(g(x)))′ = f′(g(x))·g′(x). Di mana f′(x) = 0, grafik memiliki garis singgung mendatar: maksimum, minimum, atau titik pelana.
Aturan rantai untuk logaritma
f′(x) = 4x/(2x2 + 4)
Substitusikan
f′(−2) = −8/12 = −0,6667
Jawaban
f′(−2) = −0,6667
6.Turunan
Jawaban: b = −33
Diferensial: turunan f′(x) adalah kemiringan grafik di x. Aturan pangkat: (xⁿ)′ = n·xⁿ⁻¹. Aturan perkalian: (uv)′ = u′v + uv′. Aturan rantai: (f(g(x)))′ = f′(g(x))·g′(x). Di mana f′(x) = 0, grafik memiliki garis singgung mendatar: maksimum, minimum, atau titik pelana.
Titik dan kemiringan
f(3) = 30
f′(x) = 6x + 3 → f′(3) = 21
Garis singgung di x = 0
b = f(x₀) − f′(x₀)·x₀ = 30 − 21·3 = −33
Jawaban
b = −33
7.Metode numerik
Jawaban: S = 1,091
Integrasi numerik: bagi [a, b] menjadi n pita selebar h = (b − a)/n dengan simpul x₀ … xₙ. Kaidah trapesium: T = (h/2)[f₀ + 2f₁ + 2f₂ + … + 2fn−1 + fₙ]. Kaidah Simpson (n genap): S = (h/3)[f₀ + 4f₁ + 2f₂ + 4f₃ + … + 4fn−1 + fₙ]. Galat trapesium mengecil seperti h² (batas (b − a)h²·M/12 dengan M = maks|f″|), Simpson seperti h⁴, jadi Simpson biasanya jauh lebih akurat untuk simpul yang sama.
Integral ini tidak memiliki antiturunan elementer, dan justru pada saat seperti itulah kaidah numerik diperlukan.
8.Metode numerik
Jawaban: |E| = 0,04568
Penyelesaian numerik y′ = f(t, y): metode Euler melangkah menurut kemiringan di awal setiap langkah: yn+1 = yₙ + h·f(tₙ, yₙ). Metode Heun (Euler yang diperbaiki) merata-ratakan kemiringan di awal dan pada prediksi Euler: k₁ = f(tₙ, yₙ), k₂ = f(tₙ + h, yₙ + h·k₁), yn+1 = yₙ + (h/2)(k₁ + k₂). Metode Runge–Kutta klasik (RK4) menggabungkan empat kemiringan: yn+1 = yₙ + (h/6)(k₁ + 2k₂ + 2k₃ + k₄). Galat global Euler sebanding dengan h, Heun dengan h², RK4 dengan h⁴.
Langkah Euler
y5 = 1,931326
Nilai eksak
y(2,5) = 1,885651
Galat global
|E| = |1,885651 − 1,931326| = 0,04568
Grafik
Jawaban
|E| = 0,04568
Membagi dua h kira-kira membagi dua galat metode Euler.