kill(all)$

van=facebook.com/LearningCalculus$

z(x,y) := x + y*%i;

f(x,y):=z(x,y)^3 - (2+4*%i)*z(x,y)^2 - 3*(1-3*%i)*z(x,y) + 14 - 2*%i;

R(x,y):=''(realpart(f(x,y)))$ R(x,y);

I(x,y):=''(imagpart(f(x,y)))$ I(x,y);

sol:solve([f(x,y)],[x,y])$

for i : 1 thru length (%rnum_list) do
sol : subst (t[i], %rnum_list[i], sol)$

first(sol);

second(sol);

third(sol);

M:matrix([a,b],[c,d]);

jacobian([R(x,y),I(x,y)],[x,y]);

invert(%);

plot3d([R(x,y),[x,-12,12],[y,-12,12]],
       [grid, 50, 50],[mesh_lines_color,false],
       [palette, [gray, 0.1, 0.8]] )$

plot3d([I(x,y),[x,-12,12],[y,-12,12]],
       [grid, 50, 50],[mesh_lines_color,false],
       [palette, [gray, 0.1, 0.8]] )$