Mercurial > repos > public > sbplib

view +util/matrixborrow.m @ 549:ae905a11e32c feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change away from matrix formulation in Rk4 second order form The new formulation can be up to 8 times faster
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 28 Aug 2017 11:06:00 +0200
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