function [ sDot ] = SpinEOM(t,B,S)
%function [ sDot ] = SpinEOM(t,B,S)
%
%ODE function to solve
%input vector components are rows
%

%Preallocate
omega = 13.534 * 1000; %Hz/G (13.55 MHz/KG)

%Equation
sDot = omega * [ ...
    S(2)*B(3) - S(3)*B(2) , ...
    S(3)*B(1) - S(1)*B(3) , ...
    S(1)*B(2) - S(2)*B(1) ];
%Make column for ode solver
sDot = sDot';
end

