syms x y z % Define the vector field A A = [x^2; 3*y^3; 4*z^4]; % Calculate the divergence of A div_A = divergence(A, [x, y, z]); % Evaluate the divergence at the point (1, 1, 1) div_at_point = subs(div_A, {x, y, z}, {1, 1, 1}); % Display the symbolic expression for the divergence and its numeric result div_A div_at_point