% Bloch Equation Simulation, Excercise B-1b
% -----------------------------------------
% 

df = 0;		% Hz off-resonance.
T1 = 600;	% ms.
T2 = 100;	% ms.
TE = 1;		% ms.
TR = 500;	% ms.
flip = pi/3;	% radians.


M = [0;0;1];
Rflip = yrot(flip);
[Ate,Bte] = freeprecess(TE,T1,T2,df);

% ----------- Added to B1-a -------------

[Atr,Btr] = freeprecess(TR,T1,T2,df);
M = Rflip*M;	% Magnetization after tip.
M = Atr*M+Btr;	% Magnetization at TE.

% ---------------------------------------

M = Rflip*M;	% Magnetization after tip.
M = Ate*M+Bte	% Magnetization at TE.