clear;

degmin=5;degmax=5;step=1;
t=linspace(-1,1,50);
[x,y,z]=meshgrid(t);
cmesh=[x(:) y(:) z(:)];

for n=degmin:step:degmax    
% controlmesh=wamlissa(4*n);
[wam,theta]=wamlissa(n);
[fek,ind]=dexsets(n,wam,0,[-1 1 -1 1 -1 1]);
lebc(n)=wamleb(n,wam,cmesh,[-1 1 -1 1 -1 1])
end;


