PLEASE EXPLAIN/DISCUSS THE MATLAB CODES BELOW% Let's start the program% start the program % request the user to input the desired option with regard the unit% if the user input is ‘MKS’, then the units will be in meters and seconds% If the user input is ‘FPS’, then the units will be in feet and secondsunit = input("Please choose a unit: MKS (meters and seconds) or FPS (feet and seconds): ", "s"); % conditional statements enable you to select at run time which block of code to execute% in the situation where the user has input “MKS”if (unit == "MKS" || unit == "mks")% gravity is equivalent to -9.81gravity = -9.81;% in the situation where the user has input “FPS”elseif (unit == "FPS" || unit == "fps")% gravity is equivalent to -32.174gravity = -32.174;else% in the situation where the user has input none of the provided optionserror = " Please enter an input only within the provided options ";% disp( ) displays the value of variable without printing the variable name% display errordisp(error)% return control to the invoking program before it reaches the end of the script or functionreturn% to terminate statementsend% for the variables needed in accordance with the formula% request the user to input the desired value for the magnitude of the initial velocityVo = input("Enter magnitude of the initial velocity: ");% request the user to input the desired value for the angletheta = input("Enter angle: ");% request the user to input the desired value for the heighth = input("Enter height: ");% request the user to input the desired value for the accelerationax = input("Enter acceleration: "); % mathematical formula equations% formula equation for timetime = (Vo.*sind(theta) + sqrt((Vo.*sind(theta)).^2-2*gravity*h)) / (-gravity);t = linspace(0,time); % formula equation for velocity (x)velocityx = Vo*(cosd(theta)); % formula equation for velocity (y)velocityy = Vo*(sind(theta)); % formula equation for displacement (x)displacementx = (velocityx.*t) + ((1/2)*ax*(t.^2)); % formula equation for displacement (y)displacementy = (velocityy.*t) + ((1/2)*gravity*(t.^2)); % formula equation for final velocity (x)finalvelocityx = sqrt(velocityx.^2 + 2*ax*displacementx); % formula equation for final velocity (y)finalvelocityy = sqrt(velocityy.^2 + 2*gravity*displacementy); % formula equation for magnitudemagnitude = sqrt((finalvelocityx.^2)+(finalvelocityy.^2)); % For Plotting% for figure 1figure(1)% plot function is used to create a graphical representation% linewidth specifies the width of the line% -k = solid blackplot (t, magnitude,'-k','linewidth',3)% title adds the specified title to the current axes or standalone visualizationtitle('Magnitude of Velocity Versus Time')% labels the x-axis of the current axes or standalone visualizationxlabel('Time')% labels the y-axis of the current axes or standalone visualizationylabel('Magnitude of Velocity')% displays the major grid lines for the current axesgrid on% for figure 2figure(2)% plot function is used to create a graphical representation% linewidth specifies the width of the line% -k = solid black plot (displacementx,displacementy,'-k','LineWidth',3)% hold on is utilized to add a second line plot without deleting the existing line plothold on% title adds the specified title to the current axes or standalone visualizationtitle('Trajectory of the body') % labels the x-axis of the current axes or standalone visualizationxlabel("Horizontal Displacement(m)")% labels the y-axis of the current axes or standalone visualizationylabel("Vertical Displacement(m)")% displays the major grid lines for the current axesgrid on% for figure 3figure(3);% hold on is utilized to add a second line plot without deleting the existing line plothold on% length(t) returns the length of the largest array dimension in tfor n = 1:length(t)% plot function is used to create a graphical representationplot (displacementx(1:n),displacementy(1:n), '--r','linewidth',3)% labels the x-axis of the current axes or standalone visualizationxlabel("Horizontal Displacement(m)")% labels the y-axis of the current axes or standalone visualizationylabel("Vertical Displacement(m)")% title adds the specified title to the current axes or standalone visualizationtitle("Trajectory")% temporarily stops MATLAB execution and waits for the user to press any keypause(0.1)% hold on is utilized to add a second line plot without deleting the existing line plothold on% for terminationend

Question

PLEASE EXPLAIN/DISCUSS THE MATLAB CODES BELOW% Let's start the program% start the program % request the user to input the desired option with regard the unit% if the user input is ‘MKS’, then the units will be in meters and seconds% If the user input is ‘FPS’, then the units will be in feet and secondsunit = input(&#34;Please choose a unit: MKS (meters and seconds) or FPS (feet and seconds): &#34;, &#34;s&#34;); % conditional statements enable you to select at run time which block of code to execute% in the situation where the user has input “MKS”if (unit == &#34;MKS&#34; || unit == &#34;mks&#34;)% gravity is equivalent to -9.81gravity = -9.81;% in the situation where the user has input “FPS”elseif (unit == &#34;FPS&#34; || unit == &#34;fps&#34;)% gravity is equivalent to -32.174gravity = -32.174;else% in the situation where the user has input none of the provided optionserror = &#34; Please enter an input only within the provided options &#34;;% disp( ) displays the value of variable without printing the variable name% display errordisp(error)% return control to the invoking program before it reaches the end of the script or functionreturn% to terminate statementsend% for the variables needed in accordance with the formula% request the user to input the desired value for the magnitude of the initial velocityVo = input(&#34;Enter magnitude of the initial velocity: &#34;);% request the user to input the desired value for the angletheta = input(&#34;Enter angle: &#34;);% request the user to input the desired value for the heighth = input(&#34;Enter height: &#34;);% request the user to input the desired value for the accelerationax = input(&#34;Enter acceleration: &#34;); % mathematical formula equations% formula equation for timetime = (Vo.*sind(theta) + sqrt((Vo.*sind(theta)).^2-2*gravity*h)) / (-gravity);t = linspace(0,time); % formula equation for velocity (x)velocityx = Vo*(cosd(theta)); % formula equation for velocity (y)velocityy = Vo*(sind(theta)); % formula equation for displacement (x)displacementx = (velocityx.*t) + ((1/2)*ax*(t.^2)); % formula equation for displacement (y)displacementy = (velocityy.*t) + ((1/2)*gravity*(t.^2)); % formula equation for final velocity (x)finalvelocityx = sqrt(velocityx.^2 + 2*ax*displacementx); % formula equation for final velocity (y)finalvelocityy = sqrt(velocityy.^2 + 2*gravity*displacementy); % formula equation for magnitudemagnitude = sqrt((finalvelocityx.^2)+(finalvelocityy.^2)); % For Plotting% for figure 1figure(1)% plot function is used to create a graphical representation% linewidth specifies the width of the line% -k = solid blackplot (t, magnitude,'-k','linewidth',3)% title adds the specified title to the current axes or standalone visualizationtitle('Magnitude of Velocity Versus Time')% labels the x-axis of the current axes or standalone visualizationxlabel('Time')% labels the y-axis of the current axes or standalone visualizationylabel('Magnitude of Velocity')% displays the major grid lines for the current axesgrid on% for figure 2figure(2)% plot function is used to create a graphical representation% linewidth specifies the width of the line% -k = solid black plot (displacementx,displacementy,'-k','LineWidth',3)% hold on is utilized to add a second line plot without deleting the existing line plothold on% title adds the specified title to the current axes or standalone visualizationtitle('Trajectory of the body') % labels the x-axis of the current axes or standalone visualizationxlabel(&#34;Horizontal Displacement(m)&#34;)% labels the y-axis of the current axes or standalone visualizationylabel(&#34;Vertical Displacement(m)&#34;)% displays the major grid lines for the current axesgrid on% for figure 3figure(3);% hold on is utilized to add a second line plot without deleting the existing line plothold on% length(t) returns the length of the largest array dimension in tfor n = 1:length(t)% plot function is used to create a graphical representationplot (displacementx(1:n),displacementy(1:n), '--r','linewidth',3)% labels the x-axis of the current axes or standalone visualizationxlabel(&#34;Horizontal Displacement(m)&#34;)% labels the y-axis of the current axes or standalone visualizationylabel(&#34;Vertical Displacement(m)&#34;)% title adds the specified title to the current axes or standalone visualizationtitle(&#34;Trajectory&#34;)% temporarily stops MATLAB execution and waits for the user to press any keypause(0.1)% hold on is utilized to add a second line plot without deleting the existing line plothold on% for terminationend

Accepted Answer

The given code is used to collect the information such as  unit used in  the calculation(either mks…