Diesel cycle-PV Diagram visualizer

Objective: %engine geometry% d=0.1; L=0.1; conrodlength=0.15; crankradius=L/2; R=conrodlength/crankradius; rc=17.5; cutoff=2; vs=(pi/4)*d*d*L; vc=vs/(rc-1); gamma=1.4; %state 1% p1=101325; T1=500; v1=vs+vc; %state 2% v2=vc; p2=p1*(rc)^gamma; T2=T1*(rc)^(gamma-1); theta_comp=linspace(180,360,180); term1_comp=0.5*(rc-1); term2_comp=R+1-cosd(theta_comp);…