clc
clear

r=10
phi=linspace(0,4*pi,1000);

x=r*(phi-sin(phi));
y=r*(1-cos(phi));
plot(x,y,'r-')
title('Cycloid Function for r = 10 inches')
xlabel('x (inches)')
ylabel('y (inches)')
grid on