* sh +7), 70)); (1) = c (1) +0.001 * (max ((387.339 * sh ^ 2-666.222 * sh +293.476), 70));; end; pr == 2sh <0.37 c (1) = c (1) +0.01 * (max ((-679.328 * sh ^ 2 +502.703 * sh +7), 70)); (1) = c (1) +0.01 * (max ((387.339 * sh ^ 2-666.222 * sh +293.476), 70));; end; pr == 3sh <0.37 c (1) = c (1) +0.1 * (max ((-679.328 * sh ^ 2 +502.703 * sh +7), 70)); (1) = c (1) +0.1 * (max ((387.339 * sh ^ 2-666.222 * sh +293.476), 70));; end; pr == 4sh <0.37 c (1) = c (1) + (max ((-679.328 * sh ^ 2 +502.703 * sh +7), 70)); (1) = c (1) + (max ((387.339 * sh ^ 2-666.222 * sh +293.476), 70));; end; = [a (1), b
, 0, a (2), b (2), 0,0,-a (4),-b
...
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, a (9), b [C (1); c (2) - c (4); c (2) - c (5); 0; c (18) - c (26); c (18) - c (27); 0; c (28) - c (34); c (28) - c (35); 0; c (36) +9; c (37) - c (38); c (37) - c (39); 0 ; c (40) +9; c (41) - c (42); c (41) - c (43); 0; c (44) - c (46); c (44) - c (47); 0; c (48) +9; c (49) +9; c (45) +9; c (50) - c (52); c (50) - c (53); 0; c (54) + 9; c (55) +9; c (51) +9; c (29) - c (30); c (29) - c (31); 0; c (32) +9; c (33) + 9; c (19) - c (20); c (19) - c (21); 0; c (22) +9; c (23) - c (24); c (23) - c (25) ; 0; c (16) +9; c (17) +9; c (3) - c (6); c (3) - c (7); 0; c (10) +9; c (11) - c (12); c (11) - c (13); 0; c (14) +9; c (15) +9; c (8) +9; c (9) +9];
[n_l, n] = size (L);
% ************* right boundary condition ************************* * j = 1: Nvet, (j) = V (1, N +1, j); (j) = V (2, N +1, j); (j) = gamma (j)/(A0 (j ) * (sqrt (An (j)) + sqrt (A0 (j)))); (j) = (rho * Qn (j))/(2 * An (j) ^ 2); (j) = gamma (j)/(sqrt (An (j)) + sqrt (A0 (j)));; = [a (1), b (1),-a (2),-b (1), 0,0,-a (3),-b
...
% ************* assignment *************************** ********* _od = n-n_l;% number of linealy independent solutions homogeneous left boundary condition for one vessel_0 = zeros ([n, n_od, N +1]); _f = zeros ([n, N +1]); _0 = zeros ([n, n_od, N +1]); _f = zeros ([n, N +1]); = zeros ([n_od +1, n_od +1, N]);% orthogonalising matrices = zeros ([Nfun, N +1, Nvet]);% solution = zeros ([n_od +1, N +1]);% coefficients of return marching
