
function endor_mims_echo_nitroxide()

% Isotopes
sys.isotopes={'E','14N'};

% Magnet field
sys.magnet=3.35;

% Interactions
inter.zeeman.matrix=cell(1,2);
inter.zeeman.matrix{1}=[2.01045 0.00000 0.00000
                        0.00000 2.00641 0.00000
                        0.00000 0.00000 2.00211];
inter.coupling.matrix=cell(2,2);
inter.coupling.matrix{1,2}=[1.2356 0.0000 0.6322
                            0.0000 1.1266 0.0000
                            0.6322 0.0000 8.2230]*1e7;

% Basis set
bas.formalism='sphten-liouv';
bas.approximation='none';

% Spinach housekeeping
spin_system=create(sys,inter);
spin_system=basis(spin_system,bas);

% Sequence parameters
parameters.spins={'E'};
parameters.electrons=1;
parameters.tau=200e-9;
parameters.n_dur=50e-6;
parameters.grid='rep_2ang_400pts_sph';
parameters.nsteps=200;

% Simulation
answer=powder(spin_system,@endor_mims_echo,parameters,'esr');

% Plotting
time_axis=linspace(0,2*parameters.tau,parameters.nsteps+1);
figure(); plot(1e9*time_axis,real(answer)); xlim tight;
kgrid; kxlabel('time, ns'); kylabel('intensity, a.u.');

end

