Get["QUADRUPOLE"];
  
quadrupoleSpin = 1.5;
larmorFrequencyMhz = 105.8731007;

  
quadrupoleOrder = 2;  QCCMHz = 8;  η = -1;

  
(Subscript[α, PR]) = 30;          (Subscript[β, PR]) = 
        80;                (Subscript[γ, PR]) = 120;

  
startOperator  =  
        0.4*Iz;
ωRFkHz = 
        100;
spinRatekHz = 
        15;
powderFile = "rep100_simp";
numberOfGammaAngles = 
        3;
t1 = 
        20;
Δt = 
        1;
np = t1/Δt;

  
detectelt = {{3,  2}};

fsimulation := (
acq0; 
For [p = 1,  
              p  <=  
                np,  p++, {
  pulse[Δt, ωRFkHz]; 
  acq[
                  p];
}];
);

  
run;
tabgraph["onePowderMAS"];