animalgeo

 animal geometry, [x,y] generated via [x,y]=ginput and clicking on picture

0001 function [geo, bc]=animalgeo(l)
0002  % animal geometry, [x,y] generated via [x,y]=ginput and clicking on picture
0003  x=[70.0819   62.4484   41.3096   10.7758   10.1886   40.1352   55.9893   55.4021   11.3630 ...
0004    12.5374   74.7794   74.1922   83.5872   90.6335  150.5267  151.1139  105.3132  104.1388 ...
0005   117.6441  148.1779  148.1779  118.8185   95.9181   90.0463];
0006 y=[17.0716   57.7190   62.0964   62.0964   72.7273   72.7273   84.6088  140.8898  154.6474 ...
0007   164.6529  175.9091  218.4325  218.4325  176.5344  164.6529  154.0220  142.7658   84.6088 ...
0008    72.7273   72.1019   63.3471   62.0964   57.0937   18.3223];
0009 
0010 r1=x'*l;r2=y'*(-l); geo=polygong(r1,r2);
0011  
0012  bbb=[[  2     2     2     2]; % Neumann BC, exported from gui for square
0013      [0     0     0     0];
0014      [1     1     1     1];
0015      [1     1     1     1];
0016      [1     1     1     1];
0017      [1     1     1     1];
0018      [1     1     1     1];
0019      [1     1     1     1];
0020     [48    48    48    48];
0021     [48    48    48    48];
0022     [48    48    48    48];
0023     [48    48    48    48];
0024     [48    48    48    48];
0025     [48    48    48    48];
0026     [48    48    48    48];
0027     [48    48    48    48];
0028     [48    48    48    48];
0029     [48    48    48    48];
0030     [48    48    48    48];
0031     [48    48    48    48];
0032     [49    49    49    49];
0033     [48    48    48    48];
0034     [48    48    48    48];
0035     [49    49    49    49];
0036     [48    48    48    48];
0037     [48    48    48    48]];
0038 bb=bbb(:,1);sgeo=size(geo);n=sgeo(2);bc=zeros(26,n);
0039 for i=1:n; bc(:,i)=bb; end