ottocycle with pv diagram

Objective: %Determinining a function function [V] = engine_kinematics(bore, stroke, con_rod, cr, start_crank, end_crank) %Data a = stroke/2; R = con_rod/a; theta = linspace(start_crank, end_crank, 100); v_swept = (pi/4)*bore^2*stroke; v_clearance = v_swept/(cr-1); %Determining V term1 = 0.5*(cr-1); term2 = R + 1 - cosd(theta); term3…