r3 Solutiondriverm syms x fx x2 expx3 gx cosx a065 ans p

Solution

driver.m

syms x
f(x) = x^2 + exp(-(x/3));
g(x) = cos(x);
a=0.65;
ans = partial_eval(f(x),g(x),a);
vpa(subs(ans,x,a))

partial_eval.m

function ans = partial_eval(f,g,a)
d_f = diff(f);
d_g = diff(g);
ans = d_f + d_g;
end

Output:

0.42641415292460