How to make a 3d + Color Plot with multiple variables.

Hi there! The project i am currently working on is asking for the optimization of the mass fractions of a steam power plant and I think a 3d plot with color would be a good addition to our report. Basically, we end up with an equation of efficiency (n) as a function of variables x,y, & z. Where x, y, and z, represents fractions of the mass flow rate diverted by one of the multiple turbines. Ive attached our code with forumlas (sorry if its off im pretty new to matlab).

Ive tried surf, and slice but havent been able to make anything work!

The ranges of x,y & z are from 0 - 1 with and increment of 0.01.

*Also, is there an optimization function in matlab?

If anyone could help that would be awesome!

Code:

syms x y z n

%--------------

WT1 = 2824.47-2639.052*x-2756.124*y-2221.443*(1-x-y);

WT2 = 2947.5*(1-x-y)-2676.735*z-2294.919*(1-x-y-z);

WP1 = -0.10694*(1-x-y-z);

WP2 = -9.6937*(1-x-y);

WP3 = -3.1035*y;

WP4 = -6.5382*x;

Wnet = WT1 + WT2 + WP1 + WP2 + WP3 + WP4;

%------------------------------------------------------

% finding Qout as a function of h2, h16, h17

h2 = -(WP1-191.81); h16 = 417.44*(1-x-y)-h2*(1-x-y-z); h17 = (-WT2 - h16 + 3275*(1-x-y))/(1-x-y-z);

Qout = (191.81-h17)*(1-x-y-z);

Qin = Wnet - Qout;

n = Wnet/Qin;

%Governing Equation

n=-(7331870363841724791*x + 8112779943210907746*y + 2623084686859044591*z - 8561122989305755771)/(734886084214784*x^2 + 1469772168429568*x*y + 1469772168429568*x*z - 22485840445397007859*x + 734886084214784*y^2 + 1469772168429568*y*z - 23266750024766190814*y + 734886084214784*z^2 - 734886084214784*z + 23714358184776824055)