Otto cycle using matlab

Objective: close all clear all clc gamma=1.4 t3=1000 p1=101325 t1=500 bore = 0.1; stroke = 0.1; con_rod=0.15; cr=12; v_s=(pi/4)*bore^2*stroke; v_c=v_s/(cr-1); v1=v_s+v_c v2=v_c p2=p1*cr^gamma; t2=p2*v2*t1/(p1*v1); a=stroke/2 r=con_rod/a v_s=pi/4*bore^2*stroke; v_c=v_s/(cr-1) theta=linspace(180,0); term1=0.5*(cr-1); term2=r+1-cosd(theta);…