function tau = external_torque(t)
% external torque
    tau = 0;
end
