Get["QUADRUPOLE"];

(*
One - dimensional two - 
                pulse MQMAS of a spin I  =  5/2,
3  Q echo and  - 
              3  Q antiecho amplitude optimization with the first pulse,
Coherence pathway 0  
                Q  ->  ±3 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;           
spinRatekHz = 15;
powderFile = "rep100_simp";
numberOfGammaAngles = 10;
t1 = 4;          
t2 = 1;          
Δt = 0.25;        
np = t1/Δt;     

  
elements1 = {{2, 5}, {5, 2}};    
detectelt = {{4, 3}};   

fsimulation:=(
acq0;
For [p = 1,p ≤ np, p++,{
  pulse[Δt, ωRFkHz];           
  store[1];                            
  filterElt[elements1];    

  pulse[t2, ωRFkHz];            
  acq[p];
  recall[1]                              
}];
);

  
run;
tabgraph["twoPulse_P1"];