Differentiation without symbolic toolbox. Is it possible numerically?

Suppose I have a function

x = 0:0.1:10;
fxy = @(x,y) x^3 +3*x*y^2 + 2*x

I want differentiate w.r.t 'x'

dfxy = @(y)diff(fxy,1)

and then I want to solve for 'y' for 'dfxy(y) = 0'

y0 = 0.1;
yy = fsolve(dfxy,y0)

is this possible without using symbolic toolbox? I get error like 'Function 'diff' is not supported for class 'function_handle'.'

Answer 1

[Walter Roberson]

http://www.mathworks.com/matlabcentral/fileexchange/13490-adaptive-robust-numerical-differentiation

Answer 2

[Jakub Rysanek]

https://www.mathworks.com/matlabcentral/fileexchange/59438-symbolic-derivatives