module bendy_bus

using EoM

include(joinpath("models", "input_ex_bendy_bus.jl"))

vpts = 0:800:100000
f(x) = input_ex_bendy_bus(; u=100 / 3.6, k=x, Î“=5000)
system = f.(vpts)
output = run_eom!.(system)
result = analyze.(output)

ss = :skip
impulse = :skip
vpt_name = ["k" "Angular stiffness" "Nm/rad"]
summarize(system, vpts, result; ss, impulse, vpt_name)

end

println("Done.")
