Get["QUADRUPOLE"];

(*
One - dimensional phase - modulated shifted - 
                antiecho MQMAS of a spin I  =  5/2,
Three - pulse sequence,
-3 
              Q echo amplitude optimization with the third pulse,
Coherence transfer pathway 0  
                Q  ->  -3 
                  Q  ->  1  Q  ->  -1 Q,
Wolfram Mathematica 5.0,
Author :  R .  HAJJAR
*)

  
quadrupoleSpin = 2.5;
larmorFrequencyMhz = 208.61889974;  

  
quadrupoleOrder = 2;
QCCMHz = 5;               η = 1;

  
(Subscript[α, PR]) = 0;           (Subscript[β, PR]) = 
        0;              (Subscript[γ, PR]) = 0;

  
startOperator = Iz; 
ωRFkHz = 100;      
ωRF3kHz = 10;      
spinRatekHz = 15;
powderFile = "rep100_simp";
numberOfGammaAngles = 10;
t1 = 4;         
t2 = 1.75;     
t3 = 15;       
Δt = 0.25;       
np = t3/Δt;    

  
coherence1 = {-3};        
coherence2 = {1};          
detectelt = {{4, 3}};   

fsimulation:=(
pulse[t1, ωRFkHz];              
       filterCoh[coherence1];     
pulse[t2, ωRFkHz];              
       filterCoh[coherence2];     

acq0;
For [p = 1,p ≤ np, p++,{
  pulse[Δt, ωRF3kHz];       
  acq[p];
}];
);

  
run;
tabgraph["shifted_antiecho_P3"];