pv diagram for otto cycle

Objective: % ENGINE KINEMATICS FUNCTION FOR CALCULATING VOLUME TRACE function [v]= engine_k(bore,stroke,c_rod,rc,startcrank,endcrank) r=(2*c_rod)/stroke; vs=(pi*bore^2*stroke)/4; vc=vs/(rc-1); theta=linspace(startcrank,endcrank,100); k1=.5*(rc-1); k2= r+1-cosd(theta); k3= (r^2- sind(theta).^2).^(0.5); v=(1+k1*(k2-k3)).*vc; end %…