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
Selesaikan y″ + 6y′ + 8y = 0 dengan y(0) = 3 dan y′(0) = 4. Cari y(2).
y″ + 6y′ + 8y = 0y(0) = 3y′(0) = 4
y(2) =
2.Persamaan diferensialSulit
Selesaikan masalah nilai awal y′ + 3y/x = 4x2 untuk x > 0, dengan y(1) = −3. Cari y(2).
y′ + 3y/x = 4x2y(1) = −3
y(2) =
3.IntegralSedang
Hitung I = ∫₁1,5 x·ln x dx.
∫₁1,5 x ln x dx
I =
4.IntegralSulit
Hitung I = ∫₀2 x²·eˣ dx.
∫₀2 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
Misalkan f(x) = x3 − 13,5x2 + 54x − 3. Tentukan koordinat x x dari minimum lokal.
f(x) = x3 − 13,5x2 + 54x − 3
x =
7.Metode numerikSedang
Hampiri ∫01 x³ dx dengan kaidah trapesium, n = 4, dan bandingkan dengan nilai eksak. Carilah galat mutlak |E| = |T − I|.
f(x) = x³[a, b] = [0, 1]n = 4
|E| =
8.Metode numerikSulit
Selesaikan y′ = −2ty, y(0) = −1 dengan metode Euler dan langkah h = 0,25 sampai t = 0,75. Solusi eksaknya y = y₀e−t². Carilah galat global |E| = |y(0,75) − y3|.
y′ = −2tyy(0) = −1h = 0,25
|E| =
3Larutan
1.Persamaan diferensial
Jawaban: y(2) = 0,1448
Persamaan linear orde dua: untuk y″ + b·y′ + c·y = 0 coba y = ert: ini memberi persamaan karakteristik r² + b·r + c = 0. Dua akar real r₁ ≠ r₂ memberi y = A·er₁t + B·er₂t (redaman lebih); akar ganda r memberi y = (A + Bt)·ert (redaman kritis); akar kompleks −α ± jω memberi y = e−αt(A cos ωt + B sin ωt) (redaman kurang, berosilasi). Nilai awal y(0) dan y′(0) menentukan A dan B. Ditulis sebagai y″ + 2ζωny′ + ωn²y = 0, ωn adalah frekuensi alami dan ζ rasio redaman: ζ < 1 berosilasi. Dengan masukan konstan, y″ + by′ + cy = c·K, penyelesaiannya adalah nilai tunak K ditambah respons bebas.
Akar karakteristik
(r + 2)(r + 4) = 0 → r = −2, −4
Konstanta dari nilai awal
A + B = 3, −2A − 4B = 4
A = 8, B = −5
Hitung nilainya
y(2) = 8·e−4 + (−5)·e−8 = 0,1448
Jawaban
y(2) = 0,1448
2.Persamaan diferensial
Jawaban: y(2) = 4,875
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∫3/x dx = x3
Integralkan
(x3y)′ = 4x5 → x3y = 4x6/6 + C
Kondisi awal
−3 = 4/6 + C → C = −3,667
Hitung nilainya
y(2) = (4·26/6 + (−3,667))/23 = 4,875
Jawaban
y(2) = 4,875
3.Integral
Jawaban: I = 0,1436
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 = ln x, du = dx/x; dv = x dx, v = x²/2
Antiturunan
∫x ln x dx = (x²/2)ln x − ∫(x/2) dx = (x²/2)ln x − x²/4
Limit
(1,125·ln 1,5 − 0,5625) − (0 − 0,25) = 0,1436
Jawaban
I = 0,1436
4.Integral
Jawaban: I = 12,78
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
7,389·(2) − 1·2 = 12,78
Jawaban
I = 12,78
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: x = 6
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.
Turunan
f′(x) = 3x2 − 27x + 54 = 3·(x2 − 9x + 18)
Selesaikan f′(x) = 0
x = 3 atau x = 6
Mana yang mana
f″(x) = 6x − 27: negatif di maksimum, positif di minimum.
f″(3) = −9, f″(6) = 9
x = 6
Jawaban
x = 6
7.Metode numerik
Jawaban: |E| = 0,01563
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.
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
y3 = −0,65625
Nilai eksak
y(0,75) = −0,569783
Galat global
|E| = |−0,569783 + 0,65625| = 0,08647
Grafik
Jawaban
|E| = 0,08647
Membagi dua h kira-kira membagi dua galat metode Euler.