Mercurial > repos > public > sbplib

comparison +sbp/higher_variable4.m @ 247:c2ca9717db4d feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed function name.
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 31 Aug 2016 17:31:25 +0200
parents 369c643b60c3
children 02423f9323c6
comparison
equal deleted inserted replaced
246:fe26791489e0 247:c2ca9717db4d
1 function [H, HI, D2, D4, e_1, e_m, M4, S2_1, S2_m, S3_1, S3_m, S_1, S_m] = higher4_compatible_halfvariable(m,h) 1 function [H, HI, D2, D4, e_1, e_m, M4, S2_1, S2_m, S3_1, S3_m, S_1, S_m] = higher_variable4(m,h)
2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3 %%% 4:de ordn. SBP Finita differens %%% 3 %%% 4:de ordn. SBP Finita differens %%%
4 %%% %%% 4 %%% %%%
5 %%% H (Normen) %%% 5 %%% H (Normen) %%%
6 %%% D1=H^(-1)Q (approx f?rsta derivatan) %%% 6 %%% D1=H^(-1)Q (approx f?rsta derivatan) %%%
10 10
11 %m=20; %problemstorlek 11 %m=20; %problemstorlek
12 %h=1/(m-1); 12 %h=1/(m-1);
13 %h=1; 13 %h=1;
14 14
15 c=ones(m,1); 15 c = ones(m,1);
16 16
17 17
18 H=diag(ones(m,1),0); 18 H=diag(ones(m,1),0);
19 H(1:4,1:4)=diag([17/48 59/48 43/48 49/48]); 19 H(1:4,1:4)=diag([17/48 59/48 43/48 49/48]);
20 H(m-3:m,m-3:m)=fliplr(flipud(diag([17/48 59/48 43/48 49/48]))); 20 H(m-3:m,m-3:m)=fliplr(flipud(diag([17/48 59/48 43/48 49/48])));