Given are the functions
$$f(x,y) = x+y \\ g(x,y) = xy$$
and the set $$M := \{ (x,y) \in \mathbb{R^2}: x^2+y^2 \leq 9 \}$$
The question was to visualize the level set and to mark the points
$$\max_{(x,y)\in M} f(x,y) \\ \max_{(x,y)\in M} g(x,y)$$
I know that the set is a circle with the center in its origin and a radius of $3$, but I don't know how one can calculate the maximum here.
Apparently the red dots here mark the corresponding solution. Can someone please elaborate how we get to those?

The lines that are plotted in the background all have one thing in common: The objective function is constant on those lines (you see the value of the objective function by the color).
Moreover, if you were to plot every such "level line", you would get exactly R^2. So what you do is to find the line of the brightest color that intersects with the admissible set. In the first picture, of course, you can't see such line, but that does not mean it is not existent. It is simply not plotted. What you also notice is that the diagonal lines are getting brighter and brighter the more you increase $x$ and $y$. This is very obvious since your objective function increases as $x$ and $y$ increases. So just imagine where the brightest line would be and then find its intersection with the admissible set. This is where the objective function attains its highest value.
Maybe it would be helpful if you were to command Maxima to plot more lines, but I do not know this command by heart unfortunately.