M=10

Just learned that ParametrizedSurface3D is a surface object, not a plot. So we first create the surfaces and then we render the plot, which will be the sum of all plots of the previously defined surfaces. I guess you can replace exp(-(x^2+y^2)/n by any other function, but that is just a guess.