clear;
f = @(x) x.^2;
d = @(x,u) f(x) + u;
disp(d(x,u));