1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
|
function chainap(yt)
// Copyright INRIA
[n1,n2]=size(yt);
x=ones(n1/2,n2);
y=ones(n1/2,n2);
x(1,:)=2*r(1)*cos(yt(1,:));
y(1,:)=2*r(1)*sin(yt(1,:));
for i=2:n1/2,
x(i,:)=x(i-1,:)+2*r(i)*cos(yt(i,:));
y(i,:)=y(i-1,:)+2*r(i)*sin(yt(i,:));
end
x=[0*ones(1,n2);x];
y=[0*ones(1,n2);y];
rr=r(1)/(n1+1);
rect=2*[-(n1/2+1),-(n1/2+1),n1/2+1,n1/2+1];
plot2d(0,0,[0],"012"," ",rect);
pix=xget('pixmap')
if driver()=='Pos' then st=3;xset('pixmap',1);else st=1;end
for j=1:st:n2,
xset("wwpc");
xpoly(rect([1 3 3 1]),rect([2,2,4,4]),'lines',1)
chaindp([x(:,j),y(:,j)],rr,rect);
xset("wshow");
end
xset('pixmap',pix)
function chaindp(p,r,rect)
//draw chain given sequence of points
//!
[n,ign]=size(p);
th=xget("thickness");
xset("thickness",1)
arcs=[ -r ,r,2*r,2*r,0,64*360].*.ones(n,1);
arcs=[ p, 0*ones(n,4)] + arcs;
xarcs(arcs');
xset("thickness",4);
colors=1:n;colors(8)=n+1;
xpolys([p(1:(n-1),1),p(2:n,1)]',[p(1:(n-1),2),p(2:n,2)]',colors);
xset("thickness",th);
|
