oop
{
double q, a_, q_k, a_k;
if (k2 == 0)
{= V [j]. get (1, i + 1); _ = V [j]. get (0, i + 1); _k = 0; _k = 0;
}
else
{= U0 [j]. get (1, i + 1); _ = U0 [j]. get (0, i + 1); _k = V [j]. get (1, i + 1); _k = V [j]. get (0, i + 1);
}
double c_ = (vesselBeta [j] * Math. sqrt (a_))/(2 * rho * vesselAreas [j]); [j]. set (0, 0, (- sign [j]/tDelta) * 2 * q/(q * q/a_ - c_ * a_)); [j]. set (0, 1, (- sign [j]/tDelta)/(- (q/a_) * (q/a_) + c_)); [j]. set (1, 0, - sign [j]/tDelta); [j]. set (1, 1, 0);
if (k2 == 0)
{. set (0, j, sign [j] * (q * (1/tDelta + (8 * Math. PI * nu)/a_ )/((q/a_) * (q/a_) - c_)));. set (1, j, sign [j] * a_/tDelta);
}
else
{. set (0, j, [j]
* ((1/tDelta) * (q_k + (8 * Math. PI * nu) * tDelta * q/a_)/((q/a_)
* (q/a_) - c_)));. set (1, j, sign [j] * a_k/tDelta);
}
// solution of Cauchy problem
double del = vesselLegths [j]/(NumberOfSegmentation * SegmentationForIntegrating); _mat. setMatrix (0, 1, 0, 1, (G2 [j]. minus (G1 [j])). times (1.0/SegmentationForIntegrating)); _fun. set (0, 0, (f2. get (0, j) - f1. get (0, j))/SegmentationForIntegrating); _fun. set (1, 0, (f2. get (1, j) - f1. get (1, j))/SegmentationForIntegrating); _0 = Z_0 [i]. getMatrix (j * 2, j * 2 + 1, 0, NumberOfVessels - 1); _f = Z_f. getMatrix (j * 2, j * 2 + 1, i, i);
for (jj = 0; jj
{_. setMatrix (0, 1, 0, 1, G1 [j]. plus (del_mat. times (jj + 0.5)));// intermediate
// matrix. setMatrix (0, 1, 0, 0, (f1. getMatrix (0, 1, j, j)). plus (del_fun. times (jj + 0.5)));// intermediate
// right
// part_0. plusEquals (A_. times (z_0. times (del))); _f. plusEquals (((A_. times (z_f)). plus (ff)). times (del));
}
// new value_0 [i + 1]. setMatrix (j * 2, j * 2 + 1, 0, NumberOfVessels - 1, z_0); _f. setMatrix (j * 2, j * 2 + 1, i + 1, i + 1, z_f);
for ( int
