Mercurial > repos > public > sbplib

changeset 212:8b4993d53663

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed a mistake in util.calc_borrowing.
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 20 Jun 2016 13:03:35 +0200
parents 4f7930d2d2c4
children da058ce66876 7f54cef985de
files +util/calc_borrowing.m
diffstat 1 files changed, 4 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/+util/calc_borrowing.m	Thu May 26 16:45:24 2016 +0200
+++ b/+util/calc_borrowing.m	Mon Jun 20 13:03:35 2016 +0200
@@ -1,5 +1,5 @@
 
-m = 30;
+m = 100;
 h = 1;
 
 
@@ -9,9 +9,11 @@
 S2 = S2_1*S2_1' + S2_m*S2_m';
 S3 = S3_1*S3_1' + S3_m*S3_m';
 
+alpha_I  = util.matrixborrow(M4, h^-1*S1  );
 alpha_II  = util.matrixborrow(M4, h*S2  );
 alpha_III = util.matrixborrow(M4, h^3*S3);
 fprintf('4th order non-compatible\n')
+fprintf('alpha_I1:  %.10f\n',alpha_I)
 fprintf('alpha_II:  %.10f\n',alpha_II)
 fprintf('alpha_III: %.10f\n',alpha_III)
 fprintf('\n')
@@ -69,7 +71,7 @@
 fprintf('6th order compatible\n')
 fprintf('alpha_II:  %.10f\n',alpha_II)
 fprintf('alpha_III: %.10f\n',alpha_III)
-fprintf('\n')3
+fprintf('\n')