P-V curve for Otto Cycle

Objective: function[v]=Engine_Kinematics(bore,stroke,conrod,r,start_crank,end_crank)a=stroke/2;R=conrod/a;vswept=(pi/4)*stroke*bore^2;vclear=vswept/(r-1);theta=linspace(start_crank,end_crank,720);term1=0.5*(r-1);term2=R+1-cosd(theta);term3=(R^2-sind(theta).^2).^0.5;v=(1+term1*(term2-term3)).*vclear;end clear allclose allclc%Otto…