clear all;
clc;
syms x1 x2;

% Data
t = [1 2 4 5 8]';
y = [3.2939 4.2699 7.1749 9.3008 20.259]';

% Initial guess
x = [2.50; 0.25];

% Number of iterations
n = 1; % For example, run for 10 iterations
F = x1*exp(x2*t) - y;  % Residual function

for k = 1:n
    % Evaluate F(x)
    F_val = double(subs(F, [x1, x2], x'))

    % Evaluate Jacobian J

    J_val = double(subs(jacobian(F, [x1, x2]), [x1, x2], x'))

    % Compute the update step p
    p = -(J_val' * J_val) \ (J_val' * F_val)

    % Update the estimate
    x = x + p;

    % Display current iteration information
    disp(['Iteration ', num2str(k)]);
    disp(['x = ', num2str(x')]);
    disp('---');
end