% Im (A)    deltaT (s)  delta_thetam (V)



% for Kt_over_Jm

data = [0.35    0.5     25;

        0.2     0.5     16.9;

        0.18    0.5     11.7;

        0.3     0.5     20.5;

        0.25    0.5     18.4];



% for Kt_over_(Jf+Jm)

data = [.2      0.496   2;

        .490    0.490   4.9;

        .6      0.5     6.3;

        .3      .5      3.3;

        .4      0.5     4];

    

    

Ktach = 0.0286;

dtheta_dt = 1/Ktach * data(:,3) ./ data(:,2);



Kt_over_J = data(:,1)'*dtheta_dt./(data(:,1)'*data(:,1));