this is the error i am receiving while using the jacobian function,
Error using sym/subs>normalizeInconsistency between sizes ofsecond and third arguments.
here is my code
————————————————————————
function [func, j] = sysNonlinearEqs(X) % X here is a vector
clc;x = X(1)y = X(2)z = X(3)%Define F(x)
syms x y z func(1,1)= x + y - exp(-z); % Row 1 Column 1 of Vector F(x)
func(2,1)= z + x - exp(-y); % Row 2 Column 1 of Vector F(x)
func(3,1)= y + z - exp(-x); % Row 3 Column 1 of Vector F(x)
% Define Jacobian
J = jacobian(func, [x y z])j = double(subs(J,{x, y, z}, X))
———————————————————
i am trying to call my function in this script and the error message pops up
syms X%global x y z
% Initial Conditions
X0 = [1; 1; 1];maxiter = 50;tolx = 1e-6;%computation using the Newton Raphson Method
X= X0;xold = X0;for i = 1:maxiter [func,j] = sysNonlinearEqs(X); X = X- inv(j)*func; err(:,i) = abs(X-xold); xold = X; if ( err(:,i) < tolx) break; endend
Best Answer