Mercurial > repos > public > sbplib

view +util/matrixborrow.m @ 1017:2d7c1333bd6c feature/advectionRV

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add support for using the ODE to approximate the time derivative in the residual
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Tue, 11 Dec 2018 16:29:21 +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