% Vm (V)    Im (A)  Thetadot_m (V)

data = [1.64    0.044   1.5;

        3.66    0.047   4.0;

        4.86    0.051   5.53;

        6.87    0.058   8.07;

        9.29    0.067   11.02];

    

Ktach = 0.0286;

thetadot = 1/Ktach * data(:,3); % Convert to real ang. velocity



Rm = 8.6;

e_backemf = data(:,1) - Rm*data(:,2);



result = thetadot'*e_backemf/(thetadot'*thetadot);

