Surface plot with four independent variabels

Dear Matlab experts,

I have a function f(w,x,y,z) that has four independent variables, w [1xm], x [1xn], y[1xp] and z[1xq].

In order to plot f, I use w as a parameter i.e. multiple plots with a single chosen value from w vector. Now, there are three variables left. The "surf" function allows only two independent variable matrices of order nxp created by [mesh_x, mesh_y] = meshgrid (x,y). The resultant function evaluation is also of order nxp.

But I have a third independent variable z. So, I tried to use [mesh_x, mesh_y, mesh_z] = meshgrid (x,y,z) and wanted to use the third variable matrix to encode the surface plot color C (in surf(x,y,f,C) syntax) matrix. However, this meshgrid produces 3-dimensional matrices which are not accepted by surf function. The surf function requires all matrices to be of order nxp (n and p defined as in this case).

Can some valued member please guide me how to encode the surface color in this case?? Or suggest me some other solution?? I have already checked following threads:

https://de.mathworks.com/matlabcentral/answers/114036-plotting-3-variables-of-a-single-function

https://de.mathworks.com/matlabcentral/answers/317838-4d-plot-representing-3-variables-function

https://de.mathworks.com/matlabcentral/newsreader/view_thread/154664

https://de.mathworks.com/matlabcentral/answers/104616-plot-4d-surface-from-x-y-z-c-data

but they don't apparently solve my problem.

Any suggestions and help is highly appreciated.

Thank you ..