%Uppgift 3, Anton Hagelberg och Filip Johansson, fjäderkonstanter
clear all; clc; close all;
k1_ref = 5300;
k2_ref = 136000;
ck_ref = 300;
cs_ref = 1300;

v = 1.1;
v_gammal = 1;
tol = 1e-6;
T_bar_k = 0.7;
T_bar_s = 1;

ck = 0.3;
cs = 0;


while norm(v_gammal-v) > tol

    F = transfer_functions(k1_ref,k2_ref,ck,cs);
    J = Jacobian_transfer_functions(k1_ref,k2_ref);

    v = -J\F;
    k1_ref = k1_ref + v(1);
    k2_ref = k2_ref + v(2);
    v_gammal = v;


    k1_ref
    k2_ref
end