function y = parpd()

y = [
   .085     % g0 [0.05-20]          1 
   1         % k [1-3]              2 
   0.4      % eta [0.1-0.3]         3    
   1         % Im                   4    
   5        % lambda                5    
   0.75       % T1                  6    
   0.3       % T2                   7    
   21.01      % G1 [20.9-21.1]      8    
   21.01      % G2 [20.9-21.1]      9    
   1         % Gv [2-5]             10   
   1         % Gs [1-40]            11   
   1         % Gf                   12   
   35.5       % beta                13  
   1         % gamma                14   
   1         % delta              15   
   .1       % phi                   16   
   1        % chi                   17    
   1        % rho                   18  
   1        % omega                 19   
   5       % tau                    20   
   0       % Fe                     21   
   0.99       % DA1                    22    
   0.8       % DA2                    23     
   0.7       % DA3                    24     
   0.8       % DA4                    25     
   0.8       % DA5                    26     
   0.8       % DA6                    27     
   0.8       % DA7                    28     
   0.8       % DA8                    29     
   1       % ae                     30
   20      % tau2                   31
   0.02    % Bu                     32
   0.05    % Bp                     33
   0.8        % DA9                   34     
   0.8        % DA10                  35     
   0.9      % DA11                  36       
];