%% HW01 Problem #2A
% Define some constants
dt=1e-10; % Time steps in seconds
gamma=42.57e6; % Gyromagnetic ratio, [Hz/T]
B0=1.5; % External magnetic field strength [T]
N_Steps=1000; % Number of time steps
M0=[1 0 0]'; % Initial condition [A.U.]
t_max=dt*N_Steps; % Maximum time of simulation [s]
t=linspace(0,t_max,N_Steps); % Vector of time points [s]
M_H20=zeros(4,N_Steps); % Initialize the magnetization vectors
M_H20(:,1)=[M0; 0]; % Define the initial condition
dB0_H20=PAM_B0_op(gamma,B0,dt); % Define the incremental precession
% Simulate precession
for n=2:N_Steps
M_H20(:,n) = dB0_H20*M_H20(:,n-1);
end
% Plot the results - Note F/W differences are hard to see on short time scales...
figure; hold;
plot(t,M_H20(1,:),'linewidth',2);
plot(t,M_H20(2,:),'linewidth',2);
plot(t,M_H20(3,:),'linewidth',2);
%plot(t,M_H20(3,:),'linewidth',2);
xlabel('Time [s]');
ylabel('Magnetization [a.u.]');
legend('Mx','My','Mz');
title('Free Precession in the B0 Field')
set(gcf,'Color','w');
set(gca,'Color','w','XColor','k','YColor','k');
set(gca,'Color','w','XColor','k','YColor','k','LineWidth',1.25,'Fontsize',11,'Fontweight','bold');
set(get(gca,'Title'),'Color','k','FontSize',16);
set(get(gca,'Xlabel'),'FontSize',14,'fontweight','bold');
set(get(gca,'Ylabel'),'FontSize',14,'fontweight','bold');
grid('on')
print2desktop('/Users/pmagrath/desktop/','HW1_P3_D')