Mercurial > repos > public > sbplib

changeset 219:f66513508c75 feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merged with feature/grids.
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 27 Jun 2016 13:26:02 +0200
parents fc07ebc49412 (current diff) da058ce66876 (diff)
children 5df8d20281fe
files
diffstat 1 files changed, 4 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
diff -r fc07ebc49412 -r f66513508c75 +util/calc_borrowing.m
--- a/+util/calc_borrowing.m	Mon Jun 27 13:23:46 2016 +0200
+++ b/+util/calc_borrowing.m	Mon Jun 27 13:26:02 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')