Extremwertsuche bei Funktionen mit zwei Variablen (mit Nebenbedingung)

Die Nebenbedingung lautet:

restart:
readlib(extrema):
f:=             ;# (* Funktion eingeben *)
g:=             ;#
a:=             ;#
extrema(f,g=a ,{x,y},'lsg');
# (* lsg enthält die Menge mit den Gleichungen der Lösungspunkte *)

#(* Beispiel *)
restart:
readlib(extrema):
f:= x*y;
g:= x+y     ;
a:=1 ;
extrema(f,g=a ,{x,y},'lsg');
lsg;