Plotting PV diagram for otto cycle and calculating its thermal efficiency

Objective: Firstly converting the program for determining the volume change plot during compression and expansion, into a function. function [V]=engine_kinematics(bore,stroke,con_rod,cr,start_crank,end_crank) a=stroke/2 R=con_rod/a v_s=(pi/4)*bore^2*stroke; v_c=v_s/(cr-1); theta=linspace(start_crank,end_crank,100); term1=0.5*(cr-1);…