function Boundarys

global Node BNDY1 BNDY2 BNDY3 BNDY4 BC1 BC2 BC3 BC4 


BNDY1=[3 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 2];
BNDY2=[4 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 1];
BNDY3=[1  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 2];
BNDY4=[4 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 3];

Sz1=size(BNDY1); Sz2=size(BNDY2); Sz3=size(BNDY3); Sz4=size(BNDY4);

for i=1:Sz1(2)
    BC1(i,:)=Node(BNDY1(i),:);
end
for i=1:Sz2(2)
    BC2(i,:)=Node(BNDY2(i),:);
end
for i=1:Sz3(2)
    BC3(i,:)=Node(BNDY3(i),:);
end
for i=1:Sz4(2)
    BC4(i,:)=Node(BNDY4(i),:);
end

disp('Boundary OK!')
