Code that can solve an Otto cycle and make plots for it

Objective: # otto cycle simulator import matplotlib.pyplot as plt import math def engine_kinemetics( bore,stroke, con_rod,cr,starting_crank,ending_crank): # geometric parameters a=stroke/2 R=con_rod/a num_values = 20 s_crank=math.radians(starting_crank) e_crank=math.radians(ending_crank) v_s = (math.pi/4) * pow(bore,2) * stroke v_c…