clear, close all;

[p,e,t] = initmesh('geom');
pdemesh(p,e,t), axis equal


b=[
 1   1   1   1   1   1   1   1   1   1   1   1  %N
 1   0   0   0   0   1   1   0   1   1   1   1  %M
 1   1   1   1   1   1   1   1   1   1   1   1  %Q
 1   1   1   1   1   1   1   1   1   1   1   1  %G
 1   1   1   1   1   1   1   1   1   1   1   1  %H
 1   1   1   1   1   1   1   1   1   1   1   1  %R
'0' '0' '0' '0' '0' '0' '0' '0' '0' '0' '0' '0' %Q
'0' '0' '0' '0' '0' '0' '0' '0' '0' '0' '0' '0' %G
'1' '0' '0' '0' '0' '1' '1' '0' '1' '1' '1' '1' %H
'0' '0' '0' '0' '0' '0' '0' '0' '9' '9' '9' '9' %R 
];

[q,g,h,r]=assemb(b,p,e);