Mercurial > repos > public > sbplib

changeset 218:da058ce66876 feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merged with default.
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 27 Jun 2016 13:24:59 +0200
parents 39b7dcb2c724 (current diff) 8b4993d53663 (diff)
children f66513508c75 a2cb0d4f4a02
files
diffstat 1 files changed, 4 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/+util/calc_borrowing.m	Thu Jun 16 09:21:17 2016 +0200
+++ b/+util/calc_borrowing.m	Mon Jun 27 13:24:59 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')