You can try using NIntegrate
to get integration sampling points and color for the function values.
sPnts = Reap[NIntegrate[1, {z, 0, 3},
{x, 0, 4 - (4*z)/3},
{y, 0, 2 - x/2 - (2*z)/3}, PrecisionGoal -> 6,
Method -> {"GaussKronrodRule", "Points" -> 12,
"SymbolicProcessing" -> 0},
EvaluationMonitor :> Sow[{z, x, y}]]][[2, 1]];
Clear[f]
f[z_, x_, y_] := 3 y;
Graphics3D[
MapThread[{Hue[#1], Point[#2]} &, {Rescale[f @@@ sPnts], sPnts}],
Axes -> True, AxesLabel -> {"z", "x", "y"},
FaceGrids -> {{-1, 0, 0}, {0, 0, -1}}]
