NumberOfVessels * 2); Q1 =
new Matrix (NumberOfVessels * 2, NumberOfVessels * 2); R_left =
new Matrix (NumberOfVessels * 2, NumberOfVessels * 2); Rr =
new Matrix (NumberOfVessels, NumberOfVessels); Y_f =
new Matrix (NumberOfVessels * 2, NumberOfSegmentation + 1) ; Z_f =
new Matrix (NumberOfVessels * 2, NumberOfSegmentation + 1); z_f =
new Matrix (2, 1); z_0 =
new Matrix (2, NumberOfVessels); Y =
new Matrix (NumberOfVessels * 2, NumberOfVessels + 1); Y1 =
new Matrix (NumberOfVessels * 2, NumberOfVessels * 2); Qpr =
new span> Matrix (NumberOfVessels * 2, NumberOfVessels * 2); Rpr =
new Matrix ( NumberOfVessels + 1, NumberOfVessels + 1); f1 =
new Matrix (2, NumberOfVessels); f2 =
new Matrix (2, NumberOfVessels); ff =
new Matrix (2, 1); del_mat =
new Matrix (2,2); del_fun =
new Matrix (2, 1); A_ =
new Matrix (2,2); alpha =
new Matrix (NumberOfVessels, 1); Beta = < b align = "justify"> new Matrix (NumberOfVessels + 1, NumberOfSegmentation + 1); Z =
new Matrix (NumberOfVessels * 2, NumberOfSegmentation + 1); Z_ =
new Matrix (NumberOfVessels * 2, NumberOfVessels + 1); r =
new Matrix (NumberOfVessels, 1); l =
new Matrix (NumberOfVessels, 1); Zero =
new Matrix (NumberOfVessels * 2, NumberOfVessels + 1); ZeroY_f =
new Matrix (NumberOfVessels * 2, NumberOfSegmentation + 1);
for (i = 0; i
{
for (j = 0; j
{[j]. set (0, i, vesselAreas...