computing parametre of otto cycle through MATLAB

Objective: MATLAB CODE: % To view processes in otto cycle clc clear close all bore = .1; con_rod = .15; cr = 10; gamma = 1.4; l = .2; v_s = pi/4*(bore)^2*l v_c = v_s/(cr -1) a = l/2; R = con_rod/a; p1 = 10100 %pascal t1 = 300 %kelvin v1 = v_s +v_c; v2 = v_c; v3 = v2; v4 = v1; const = p1*v1^gamma; p2 = const/v_c^gamma; theta1 = linspace(0,180,300);…