%let n = 80; 
goptions reset=all; 
goptions cback="blue" ctext="white"; 
Data fabric; 
if _n_=1 then do;  X0=0; R1=0; pi = 4*atan(1); end;
input R;  X0=X0+(R+R1)+ R1*.05 ; 
do wheel = 1 to 2; Y0=R*((-1)**wheel); id+1;
do i = 0 to &n; 
S = R*sin(2*pi*i/&n); C=R*cos(2*pi*i/&n); 
X = X0+C; Y=Y0+S; output; keep X Y X0 Y0 id; 
end; end; R1=R; retain;
cards; 
11 
 8 
 5 
16 
6  
20 
 ; 
 proc means; 
 axis1 length = 3.4 inches offset = (.2 inches,.2 inches);
 axis2 length = 5.7 inches offset = (.2 inches,.2 inches); 
 Title "Rollers"; 
 proc gplot; 
 plot Y*X= id /haxis=axis2 vaxis=axis1 nolegend vref=0 cvref=white; 
 symbol1 v=none i=spline r=12 c=white w=2; 
 run;