function [dfx, dfy] = q3_b(func)
syms x y;
symf = feval(func,x, y);
deriv = diff(symf,x);
dfx = inline(char(deriv));
deriv = diff(symf,y);
dfy = inline(char(deriv));