Funcion solo para dos variables
function [val,ea]=error2(f,x1,ex,y1,ey)
syms x y
dy=diff(f,y);
dx=diff(f,x);
dx=inline(dx);
dy=inline(dy)
ea=ex*dx(x1,y1)+ey*dy(x1)
f=inline(f);
val=f(x1,y1)
ea/val
end
ejemplo en el matlab
error2('x^2+4*x*y',2,0.2,3,0.1)
pruebenlo