Plotting Pressure vs Volume graph for Otto Cycle and also computing its Thermal Efficiency

Objective: function [V] = engine_kinematics(bore,stroke,con_rod,cr,startc,endc) a=stroke/2; R=con_rod/a; Vs=pi/4*bore^2*stroke; theta=linspace(startc,endc,720) Vc=Vs/cr-1; t1=0.5*(cr-1); t2=R+1-cosd(theta); t3=(R^2-sind(theta).^2).^0.5; V=(1+t1*(t2-t3)).*Vc; end ***************************************** close all clear all…