Air standard cycle

close all
clear all
clc
%Inputs
gamma=1.4;
t3=2300;
%state variables
p1=101325;
t1=500;

%Engine geometric parameters
bore=0.1;
stroke=0.1;
con_rod=0.15;
cr=12; %compression ratio

%calculting swept and clearance volume
v_swept=(pi/4)*bore^2*stroke;
v_clearance=v_swept/(cr-1);

v1=v_swept+v_clearance;
v2=v_clearance;

%state variables at point 2
%p2*v2^gamma=p1*v1^gamma
%p2=p1*cr^gamma since cr=(v1/v2)^gamma
p2=p1*cr^gamma;
%(p1*v1)/t1=(p2*v2)/t2 | t2= (p2*v2*t1)/p1*v1
t2= p2*v2*t1/(p1*v1);
constant_c= p1*v1^gamma;
V_compression = engine_kinematics(bore,stroke,con_rod,cr,180,0);
P_compression= constant_c./V_compression.^gamma; 

%state variables at point 3
v3=v2;
%(p3*v3)/t3=(p2*v2)/t2 
p3= p2*t3/t2;

constant_c= p3*v3^gamma;
V_expansion = engine_kinematics(bore,stroke,con_rod,cr,0,180);
P_expansion= constant_c./V_expansion.^gamma; 

%state variables at point 4
v4=v1;
%p3*v3^gamma=p4*v4^gamma
%p4=p3*(v3/v4)^gamma
p4=p3*(v3/v4)^gamma;

thermal_efficiency = 1-((1/cr)^(gamma-1));

figure(1)
hold on
plot(v1,p1,\'*\',\'color\',\'g\')
plot(V_compression,P_compression)
plot(v2,p2,\'*\',\'color\',\'g\')
plot(v3,p3,\'*\',\'color\',\'g\')
plot(V_expansion,P_expansion)
plot(v4,p4,\'*\',\'color\',\'g\')
plot([v2 v3],[p2 p3],\'color\',\'b\')
plot([v1 v4],[p1 p4],\'color\',\'b\')

I have used a call in fuction to get the Pv diagram for otto cycle and have also joined the plots. Also, found the thermal efficiency for an otto cycle.

I kindly request you to provide with a feedback, since when I asked the support team, they told that I have to use for loop to get the PV diagram, but when I asked in the whatsap group, they told its not required to use a for loop.