Plotting these functions is really awkward. They go to extremely large/small numbers and change dramatically if you change the parameters. The easiest way to play with the parameters is to use Manipulate:

ClearAll["Global`*"]

f[w_, x_, y_, z_] := w ((x - 50000)/y)^(z - 1) (E^((50000 - x)/y))/(y Gamma[z]);

Manipulate[Plot[f[w, x, y, z], {x, 50000, 50050}], {{w, -0.65}, -1, 1}, {{y, 0.865}, -1, 1}, {{z, -0.01}, -0.01, 0.01}]

You can ow change the parameters in some (narrow) regions. The function is very sensitive and it is difficult to plot everything reasonably into one figure. This here does work, but is not very helpful:

Table[Plot[f[w, x, y, z], {x, 50000, 50050}], {w, -0.7, -0.6, 0.025}, {y, 0.86, 0.87, 0.002}, {z, -0.01, 0.01, 0.005}] // Show

Alternatively, this works

Plot[Table[Evaluate[f[w, x, y, z]], {w, -0.7, -0.6, 0.025}, {y, 0.86, 0.87, 0.002}, {z, -0.01, 0.01, 0.005}], {x, 50000, 50050}]

Cheers, Marco

Hi,

in how many dimensions do you want to plot that function? It appears that you have 4 independent variables and one dependent variable. Do you want to plot that in four dimensions and then colour-code the value? Even if you were to plot in the three available spatial dimensions, colour code the value of f and then use "time", i.e. make a video, for the remaining variable it would be very difficult to plot.

Can you show a figure of a function of four independent variables that shows how you want it to be plotted? There are tricks of visualising that, but it would be good to understand your ideas of how you want this to be represented.

M.