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

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);
[Atr,Btr] = freeprecess(TR,T1,T2,df);

% M1 = Atr * Rflip * M + Btr.
%
% But M1=M in steady state, so
%
% 	M = Atr*Rflip * M + Btr.
%	(I-Atr*Rflip)*M = Btr.

M = inv(eye(3)-Atr*Rflip)*Btr