exercise_batterypower;

m_new_batteries = [1:3000];
E_unit_batteries = E_batteries_total / m_batteries;
E_new_batteries = m_new_batteries * E_unit_batteries;

m_segway = m_segway_nobat + m_new_batteries;
Torque = (d_wheel / 2) * F_total(fsa(h_rider, m_rider) + A_segway, v_kmh2msec(v_const), c_D, m_segway + m_rider, alpha);
P_mech_effective = Torque * v_msd2rad(v_kmh2msec(v_const), d_wheel);
P_from_batteries = P_mech_effective * (1 / eff_overall);
t_travelled = E_new_batteries ./ P_from_batteries;

dist_travelled = t_travelled * v_kmh2msec(v_const);
dist_travelled_km = dist_travelled / 1000;

plot(m_new_batteries, dist_travelled_km), xlabel('Weight of Batteries in kg'), ylabel('Distance travelled in km'), title('Weight of batteries vs. distance travelled at constant velocity of 20km/h'), grid;