Mercurial > repos > public > sbplib
diff +util/calc_borrowing.m @ 212:8b4993d53663
Fixed a mistake in util.calc_borrowing.
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Mon, 20 Jun 2016 13:03:35 +0200 |
| parents | 03151dd46534 |
| children | d24869abc7cd |
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')