function [bnds,pert] = skewmu(matin,blk,fixnorm)

    [nblk,dum] = size(blk);
    blkp = ptrs(blk);
    if dum~=2
        error('Invalid Block structure');
        return
    end

    fixnorm = fixnorm(:);
    [nr,nc] = size(fixnorm);
    if nr~=nblk & nc~=1
        error('Invalid Fix-Norm Specification');
        return
    end
    
    [mtype,mrows,mcols,mnum] = minfo(matin);
    leftone = [];
    lefttwo = ones(mrows,1);
    rightone = [];
    righttwo = ones(mcols,1);
    for i=1:nblk
        ldim = abs(blk(i,1));
        if blk(i,2)==0
            rdim = ldim;
        else
            rdim = blk(i,2);
        end
        if isinf(fixnorm(i))
            leftone = [leftone blkp(i,2):blkp(i+1,2)-1];
            rightone = [rightone blkp(i,1):blkp(i+1,1)-1];
            lefttwo(blkp(i,2):blkp(i+1,2)-1) = ones(1,rdim);
            righttwo(blkp(i,1):blkp(i+1,1)-1) = ones(1,ldim);
        else
            lefttwo(blkp(i,2):blkp(i+1,2)-1) = sqrt(fixnorm(i))*ones(1,rdim);
            righttwo(blkp(i,1):blkp(i+1,1)-1) = sqrt(fixnorm(i))*ones(1,ldim);
        end    
    end        
    am = diag(lefttwo)*m*diag(righttwo);    
    while 
        sm = am;
        sm(leftone,:) = sqrt(alpha)*sm(leftone,:);
        sm(:,rightone) = sqrt(alpha)*sm(:,rightone);
        bb = mu(sm,blk);
        xdata = [xdata;alpha];
        ylow = [ylow;bb(2)];
        yupp = [yupp;bb(1)];        
        if bb(1)<1
            if alpha>alphau
                alphau = alpha;
            end
        end
        if bb(2)>=1
            if alpha<alphal
                alphal = alpha;
            end
        end
    end
            
