mathematica_1.0

xzak = {-10, 10};
yzak = {-10, 10};
fun := x^2 + y^2;
Graphics3D[
 {
  PointSize[0.5]
   Solve[fun == 0 && x <= xzak[[2]] && x >= xzak[[1]] &&
     y >= yzak[[1]] && y <= yzak[[2]], {x, y}];
  N[%]
   stacjonarne =
    Solve[D[fun, x] == 0 && D[fun, y] == 0 && x <= xzak[[2]] &&
      x >= xzak[[1]] && y >= yzak[[1]] && y <= yzak[[2]], {x, y}];
  wyznacznik = D[D[fun, x], x]*D[D[fun, y], y] - D[D[fun, x], y]^2;
  For[i = 1, i <= (Length[stacjonarne]), i++,
     If[(wyznacznik /. stacjonarne[[i]]) < 0,
      Print[stacjonarne[[i]], "-punkt siodlowy"],
      If[(wyznacznik /. stacjonarne[[i]]) < 0,
       Print["Nie mozna roztrzygnac"],
       If[((D[D[fun, x], x]) /. stacjonarne[[i]]) > 0,
        {Print[stacjonarne[[i]], "-minimum"], Point[stacjonarne[[i]]]},
        {Print[stacjonarne[[i]], "-maksimum"], Point[stacjonarne[[i]]]}
        ]
       ]
      ]
     ]
    wykres =
   Plot3D[fun, {x, xzak[[1]], xzak[[2]]}, {y, yzak[[1]], yzak[[2]]},
    ColorFunction -> Function[{x, y, z}, Hue[z]]];

  wykres[[1]],

  },
 BoxRatios -> Automatic
 ]