Mercurial > repos > public > sbplib

view +util/matrixborrow.m @ 975:47b48a97c675 feature/getBoundaryOperator

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Stupid merge makes this branch dangerous. Merging with other branches will result in unexpected changes. Never merge this branch with anything. Closing branch.
author Martin Almquist <malmquist@stanford.edu>
date Thu, 03 Jan 2019 13:24:39 +0100
parents 48b6fb693025
children
line wrap: on
line source

function alpha = matrixborrow(A,B)

    function y = neg(x)
        y = min(eig(A-x*B));
    end


    x0 = 1;
    while neg(x0) > -1
        x0 = 2*x0;
    end

    opt = optimset('Display','none');
    alpha = fsolve(@neg, x0,opt);
end