function [pd, t1, t2] = get_nmr_phantom(n)

mri_phan_pdens = [    
         0.9  .69   .92    0     0     0 
         0.06 .60 .80   0  -.0184   0
];

mri_phan_t1 = [    
         500  .69   .92    0     0     0 
         -200 .60 .80   0  -.0184   0
         210 .1100 .3100  .22    0    -18
         210 .1600 .4100 -.22    0     18
         300 .2100 .2500   0    .35    0
];

mri_phan_t2 = [    
         40  .69   .92    0     0     0 
         10 .60 .80   0  -.0184   0
         60 .1100 .3100  .22    0    -18
         40 .1600 .4100 -.22    0     18
         -10 .2100 .2500   0    .35    0
];

n = 32;
pd = phantom(mri_phan_pdens, n)+0.1;
t1 = phantom(mri_phan_t1, n);
t2 = phantom(mri_phan_t2, n);

end