Shirp[0,{i1_,j1_}]:={Rectangle[{i1,j1}]}; Shirp[n_,{i2_,j2_}]:=Module[{s={}},Do[If[Not[u==v==2], s=Union[s,Shirp[n-1,{u*3^(n-1)+i2,v*3^(n-1)+j2}]]],{u,3},{v,3}];s]; Show[Graphics[Shirp[5,{0,0}]],AspectRatio->1];