Mercurial > repos > public > sbplib

diff +util/matrixborrow.m @ 0:48b6fb693025

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Initial commit.
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 17 Sep 2015 10:12:50 +0200
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/+util/matrixborrow.m	Thu Sep 17 10:12:50 2015 +0200
@@ -0,0 +1,15 @@
+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