function [ x,y ] = Geom(bs,s)
nbs = 10;
if nargin==0
x=nbs;
return
end
dl=[ 0 1 2 3 4 5 6 7 8 9;
1 2 3 4 5 6 7 8 9 10;
0 0 0 0 0 0 0 0 1 1;
1 1 1 1 2 2 2 1 2 2];
if nargin==1
x=dl(:,bs);
return
end
x=zeros(size(s));
y=zeros(size(s));
[m,n]=size(bs);
if m==1 && n==1,
bs=bs*ones(size(s));
elseif m~=size(s,1) || n~=size(s,2),
error(message('pde:cardg:SizeBs'));
end
nn=size(bs,1)*size(bs,2);
if nn==0
x=[];
y=[];
return
end
for ibs=1:nn
switch bs(ibs)
case 1
xx=50 - (50-10)*s(ibs);
yy=-50 +0*s(ibs);
case 2
xx=10 - (10-5)*(s(ibs)-1);
yy=-50*(1-(s(ibs)-1));
case 3
xx=5*cos(pi/2*(s(ibs)-2));
yy=5*sin(pi/2*(s(ibs)-2));
case 4
xx=0*(s(ibs)-3);
yy=5 + ((20-12)-5)*(s(ibs)-3);
case 5
xx=0*(s(ibs)-4);
yy=(20-12) + (12)*(s(ibs)-4);
case 6
xx=50*(s(ibs)-5);
yy=20 + 0*(s(ibs)-5);
case 7
xx=50 + 0*(s(ibs)-6);
yy=20 - (12)*(s(ibs)-6);
case 8
xx=50 + 0*(s(ibs)-7);
yy=(20-12) - ((20-12)-(-50))*(s(ibs)-7);
case 9
xx=30*(s(ibs)-8);
yy=(20-12)+3*sin(2*pi/1*xx);
case 10
xx=30-(30-50)*(s(ibs)-9);
yy=20-12+0*(s(ibs)-9);
end
x(ibs)=xx;
y(ibs)=yy;
end
end
