% Loops and Basic Programming
% Dr Owen Williams
% University of Wales TSD

clear;clc

% Calculate the Real Cylinder Volume

Bore=98;                %Cylinder Bore (mm)
Stroke=66;              %Piston Stroke (mm)
Rod=120;                %Con Rod Length (mm)
Clearance=58.8;         %Clearance Volume (cc)

Vc=Clearance*1000                                       %Clearance Volume (mm^3)
Rad=Stroke/2                                            %Crank Radius (mm)
Comp=((Stroke*((pi)*(Bore.^2))/4)+Vc)/Vc                %Compression Ratio
Crank=-180:1:180;                                       %Crank Angle Range
s=Rad*cosd(Crank)+((Stroke.^2)-(Rad^2)*sind(Crank).*sind(Crank)).^(1/2);  %s  
Vol=(Vc+((((pi)*(Bore.^2))/4))*(Stroke+Rad-s))/1000;    %Cylinder Volume (cc)

plot(Crank,Vol)