.liquidFluid.getDensity(pl);
            rhog = this.gasFluid.getDensity(pg);
            <span class="comment">% Compute relative permeability coefficients</span>
            klr = this.liqRelativePermeability.calculate(Sl, this.porousMedia);
            kgr = this.gasRelativePermeability.calculate(Sl, this.porousMedia);
            <span class="comment">% Permeability matrices</span>
            Kll = - K * klr / mul;
            Klg = K * klr / mul;
            Kgl = zeros(2);
            Kgg = K * kgr /