jacobian matrix for system of equations in matlab

main.m
syms x y z; % Define symbolic variables

eq1 = x^2 + y^2 + z^2 == 1; % Define equations
eq2 = x + y - z == 0;
eq3 = x - y + z == 2;

eqs = [eq1, eq2, eq3]; % Combine equations into an array

vars = [x, y, z]; % Define variables

J = jacobian(eqs, vars); % Calculate the Jacobian matrix

simplify(J) % Simplify the Jacobian matrix (optional)
341 chars
14 lines

main.m
J_eval = subs(J, [x, y, z], [1, 2, 3]); % Evaluate the Jacobian matrix at x=1, y=2, z=3
88 chars
2 lines