Air standard cycle
Objective: function[V_compression]=engine_kinematics(bore,stroke,con_rod,cr,start_crank,end_crank) a=stroke/2; R=con_rod/a; v_s=(pi/4)*bore^2*stroke; v_c= v_s/(cr-1); theta=linspace(start_crank,end_crank,100); term1=0.5*(cr-1); term2=R+1-cosd(theta); term3=(R^2-sind(theta).^2).^0.5; V_compression =(1+term1*(term2-term3)).*v_c; end…