Tugas rungekutta MFS

Koding

import numpy as np m = eval(input('massa benda: ')) g = eval(input('percepatan gravitasi: ')) u = eval(input('koefisien gesek bidang: ')) cd = eval(input('koefisien gesek udara: ')) teta = eval(input('sudut kemiringan bukit: ')) * (np.pi / 180) sint = np.sin(teta) cost = np.cos(teta) v0 = 0 t0 = 0 dt = 1 error = 100 tries = [] p = g * (sint - u * cost) q = cd / m def dvdt(t0, v0): return p - q * (v0 ** 1.5) while error > 0.005: k1 = dvdt(t0, v0) k2 = dvdt(t0 + 0.5 * dt, v0 + 0.5 * dt * k1) k3 = dvdt(t0 + 0.5 * dt, v0 + 0.5 * dt * k2) k4 = dvdt(t0 + dt, v0 + dt * k3) v1 = v0 + (1.0 / 6.0) * (k1 + 2 * k2 + 2 * k3 + k4) t0 = t0 + dt error = ((v1 - v0) / v1) * 100 v0 = v1 tries.append(v1) t1 = len(tries) print(f'kecepatan maksimum yang dapat dicapai adalah {round(v1, 2)} m/s dengan waktu {t1 + 1} s')