otto_cycle

Objective: function [v] = otto_Cycle_function1(bore,stroke,con_rod,cr,start_crank,end_crank) a = stroke/2; R = con_rod/a; v_s = (pi/4)*bore^2*stroke; v_c = v_s/(cr-1); theta = linspace(start_crank,end_crank,180); term1 = 0.5*(cr-1); term2 = R+1-cosd(theta); term3 = (R^2-sind(theta).^2).^0.5; v = v_c.*(1+term1*(term2-term3)); end…