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')