A new type of bird capable of self-propulsion is designed inorder to maximise its flying distance in the presence ofaerodynamic drag. This new bird, having the initial mass of mo,generates thrust by ejecting 1% of its initial mass (0.01mo) withthe velocity of 50m/s in the direction opposite to the bird’sinstantaneous velocity relative to the bird every 0.1s (similar tohow rockets generate thrust).

function [vx_f,vy_f,m_f] =mass_ejection(vx_i,vy_i,v_eb,m_i,m_e)

% Inputs:

% vx_i: Velocity component in x direction before ejection

% vy_i: Velocity component in y direction before ejection

% v_eb: Ejection velocity of the mass relative to the bird

% m_i : Mass of the bird before ejection

% m_e : Ejection mass

% Outputs: % vx_f: Velocity component in x direction afterejection

% vy_f: Velocity component in y direction after ejection

% m_f : Mass of the bird after ejection

(d) Using the function written incorporate the mass ejectionmechanism into the existing model with drag (remember that you needto run this function every 0.1s and use the output to update thevelocity and mass of the bird). Compute and plot the trajectory ofthe new bird and compare it with the trajectory of the originalbird with the same initial conditions. On your plot, indicate withsymbols where mass ejections happen. The bird is launched with theinitial velocity of 10m/s at the angle of 50? from the horizontal.The initial mass of the bird, mo is 1kg. The release point of thebird is 2m directly above the base of slingshot and the base ofslingshot can be considered as the origin (figure 5). The densityof air is 1.2 kg/m3 and use CD = 1. The bird can be treated as asphere with a radius of 0.15 m. on matlab