Otto Cycle PV plot and Thermal Efficiency

Objective: % Script 1 - Volume trace function for compression and expansion % splitting the formula for change of volume with rotation of crank % v = v_c(1+0.5*(cr-1)*(R +1 -cosd(theta)-((R^2-(sind(theta))^2))^0.5)) function [v] = engine_kinematics(bore, stroke, con_rod, cr, start_crank, end_crank) r = stroke/2; R = con_rod/r; %…