% create model
t = linspace(0,2*pi,10);
[x,y] = pol2cart(t,1);
X = [y*0;y*0;y*0+2;y*0+2];
Y = [y*0;y;y;y*0];
Z = [x*0;x;x;x*0];
h1 = surf(X,Y,Z); % closed cylinder
[x0,y0] = meshgrid(-1:3,-1:10);
h2 = surface(x0,y0,x0*0-1); % plate
% animte
axis equal
for i = 1:10
rotate(h1,[1 0 0],pi/10,[0 0 0])% rotate wheel
y1 = get(h2,'ydata');
set(h2,'ydata', y1-pi/10) %move plate
pause(0.5)
end
