disp('go_laplace')

% create n1=10; n2=5; h1=0.5;h2=2;
a=-1; b=5;c=3;

disp('creating local matrices')
jmf=model_laplace(n1,n2,h1,h2);

disp('assembling')
K=assemble(jmf.Kloc,jmf.ien,jmf.id);

ndof=length(K);
f=rand(ndof,1);

disp('applying point constraints')
cnodes=jmf.edge{1}; % nodes
cdofs=unique(nozeros(id(:,cnodes))); % constrained dofs