# # # # # Duf = 'diff(f,x)'*a + 'diff(f,y)'*b; grad(f(x,y),[x,y]); Duf = innerprod( grad(f(x,y),[x,y]) , [a,b] ); #