You have four independent inputs, and one output, for a total of five variables.
You cannot use slice() or isosurface() as those are restricted to three independent inputs and one output, for a total of four variables.
Therefore, you will need to use scatter3(), with the x, y, z coordinates encoding three of the inputs, the point size encoding the fourth input, and the color encoding the output. You may need to set AlphaData in order to be able to see "inside" the 4d shells.
You will probably not be able to understand the result.
You will not want to display all the possible points: only the ones corresponding to the isosurfaces you have chosen. Unfortunately, figuring out which those are might take a bit of work. You might have to resort to doing a series of 3d contours and stitching them together. You can extract the contour matrix from contour3() https://www.mathworks.com/help/matlab/ref/contour3.html#mw_8699e0c8-d627-48fd-9de6-b1a898362667 and there are some File Exchange contributions that make it easier to interpret those contour matrices.


