function c=compilemaxpower(xl,xr,yb,yt,M,N)
c=[];
format long
for K=1:0.1:5;
    a=K+3;
    b=a+2;
    p=findmaxpower(xl,xr,yb,yt,M,N,K,a,b);
    P=p(end,1);
    max_w=p(end,2);
    k=[K P max_w];
    c=cat(1,c,k);
end