clear all; clc; close all;
k1_ref = 5300;
k2_ref = 136000;
c1_ref = 300;
c2_ref = 1300;

v = 1;
tol = 1e-6;
F1 = k1_ref;
F2 = k2_ref;
J = Jacobian_transfer_functions(k1_ref,k2_ref);

while norm(v)>tol

    F = transfer_functions(F1,F2,c1_ref,c2_ref);

    v = -J(1,:)'/F;
    F1 = F1 + v(1);
    F2 = F2 + v(2);

end