Using the diff command and solve command

can someone take a look at my code and tell my why its not outputing the paritial derivatives of fx and fy

the instructions were

Plot the function f (x, y) = 2x3y − y2 − 3xy.

Find fx and fy using the diff command. Setting these to zero, solve the resulting system of

equations using the solve command.

i used both diff command and solve command im not getting an error message for my lines of code but the code doesnt output an answer my code is below

f=@(x,y) 2*x^3*y - y^2 - 3*x*y;
syms x y;
G = [diff(f,x); diff(f, y)];
solve x y

Incorrect number or types of inputs or outputs for function 'solve'.

Answer 1

[Walter Roberson]

f=@(x,y) 2*x^3*y - y^2 - 3*x*y;
syms x y;
G = [diff(f,x); diff(f, y)];
sol = solve(G, [x y])

sol = struct with fields:
    x: [5×1 sym]
    y: [5×1 sym]

subs([x, y], sol)

ans =