function F=blochEq(M,B,T1,M0,g)
% blochEq <Bloch Equation>
% Usage:: F=blochEq(M,B,T1,M0,g)
%
% Return F the rhs of Bloch equation.
% M is 3x1 magnetization vector
% B is 3x1 magnetic field vector
% T1 is the relaxation time T_1
% M0 is the equilibrium magnetization M_0
% g is \gamma the gyromagnetic ratio

% revision history:
% 05/11/2026 Mark D. Shattuck <mds> blochEq.m

%% Main
w=g*B;  % (short cut) Larmor frequency vector 

% rhs of Bloch equation 
F=[-M(1)/T1+M(2)*w(3)-M(3)*w(2);...
   <fill in Fy here>;...
   <fill in Fz here>];

