kill(all);

load(mnewton)$

van=facebook.com/LearningCalculus$

f(x,y):=x^4 - 2*x^3 - 11*y^2 + 12*y + 41 $
g(x,y):=y^4 - 2*y^3 - 11*x^2 + 12*x + 31 $

s:[solve(f(3, y),y),
   solve(g(3, y),y),
   solve(f(x,-2),x),
   solve(g(x,-2),x)]$

[s[1][1],s[1][2]];
[s[2][3],s[2][4]];
[s[3][3],s[3][4]];
[s[4][1],s[4][2]];

X:[rhs(s[3][3]),rhs(s[3][4]),rhs(s[4][2])]$
Y:[rhs(s[1][1]),rhs(s[1][2]),rhs(s[2][3])]$

for i: 1 thru 3 do
    for j: 1 thru 3 do
    (
        z: f(X[i],Y[j])-g(X[i],Y[j]),
        if z = 0 then
        print("Answer ", X[i], Y[j])
    );