Air Standard cycle

% PROGRAM FOR OTTO CYCLE

clear all
close all
clc
%Input
gamma = 1.4;
t3=2300State Variables
p1=101325;
t1=500;
%Engine geometric parameters
bore =0.1;
stroke=0.1;
con_rod=0.15;
cr=12;
%calculating the swept Volume& clearance volume
v_swept = (pi*(bore/2)^2*stroke);
v_clearance= v_swept/(cr-1);
v1 = v_swept+v_clearance;
v2 = v_clearance;
%State Variables@ point 2
%p1v1^gamma=p2v2^gamma
%p2=p1*(v1/v2)^gamma=p1*cr^gamma
p2=p1*cr^gamma;
%p1v1/t1=p2v2/t2
t2=p2*v2*t1/(p1*v1);
constant=p1*v1^gamma;
v_com=engine_kinematics(bore,stroke,con_rod,cr,180,0);
p_com=constant./v_com.^gamma;
%State variable @ point 3
v3=v2;
%p3v3/t3=p2v2/t2
p3=p2*t3/t2;
constant=p3*v3^gamma;
v_exp=engine_kinematics(bore,stroke,con_rod,cr,0,180);
p_exp=constant./v_exp.^gamma;
%State variable @ point 4
v4=v1;
%p3v3^gamma=p4v4^gamma
p4 = p3*(v3/v4)^gamma;
%Thermal Efficiency
therm_eff=(1-(1/(cr)^(gamma-1)))*100;
x=[\'Thermal Efficiency of otto cycle is \',num2str(therm_eff),\' %\'];
disp(x)
%Plotting
figure(1)
hold on
plot (v1,p1,\'*\',\'color\',\'r\')
plot (v2,p2,\'*\',\'color\',\'r\')
plot (v4,p4,\'*\',\'color\',\'r\')
plot(v_com,p_com,\'color\',\'b\')
plot(v_exp,p_exp,\'color\',\'g\')
plot([v2 v3],[p2 p3],\'color\',\'k\')
plot([v4 v1],[p4 p1],\'color\',\'r\')
xlabel(\'VOLUME\')
ylabel(\'PRESSURE\')
title(\'P-V DIAGRAM\',\'color\',\'k\')

%Plot for compression ratio vs thermal efficiency

figure(2)
%inputs
comp_ratio=1:1:20;
thermal_eff=(1-(1./(comp_ratio).^(gamma-1)))*100;
plot(comp_ratio,thermal_eff)
xlabel(\'Compression ratio\')
ylabel(\'Thermal efficiency\')
title(\'Compression ratio Vs Thermal efficiency\',\'color\',\'k\')

PROJECT LINK https://projects.skill-lync.com/projects/Air-standard-cycle-26490

Screenshots of plot