% Percent, viscosity , T2 time (ms), numerical error (very questionable)
data = [1.00 1410 19.022840 2.487071;
        .98  939 19.964525 2.173574;
        .96  624 32.768710 5.898096;
        .94  437 36.278088 8.483888;
        .91  259 39.579193 7.836243;
        .84  109 90.326811 26.973864;
        .75  35.5 150.309497 20.761621;
        .65  15.2 316.628086 104.346099;
        .50  6.00 271.344190 87.715985;
        .10  1.31 354.612248 127.870084];
 
 %data(:,3:4) = data(:,3:4)/1000; % Convert to seconds
 data(:,2) = data(:,2)/100; % Convert to poises
 
 %data(:,2) = data(:,2)/293; % Poises per degree
 
 loglog(data(:,2),data(:,3),'o','markersize',5); hold on;
 errorbar(data(:,2),data(:,3),data(:,4),'o','markersize',5); hold on;
 
 T1data = [100 1410 0.0250512 0.000679172;
           100 1410 0.0238054 0.00278499;
           100 1410 0.0238191 0.00264131;
           98  939  0.0531649 0.00157548;
           91  259 0.0989349 0.00607012;
           80  60.1 0.308794  0.0434989;
           65  15.2 0.291683  0.039015;
           50  6.00 1.43547   1.00943;
           10, 1.31 2.13718   0.893303];
 T1data(:,3:4) = T1data(:,3:4)*1000; % Convert to milliseconds
 T1data(:,2) = T1data(:,2)/100;
 errorbar(T1data(:,2),T1data(:,3),T1data(:,4),'r^','markersize',5);
 xlim([1e-2 12]);
 fs = 14;
 ylabel('Relaxation times (msec)','fontsize',fs);
 xlabel('Viscosity (poises)','fontsize',fs);
 
 T1bloem = [.035 1000;
            .1 500;
            .25 200;
            .65 100;
            2  50;
            6.5  32];
 plot(T1bloem(:,1),T1bloem(:,2),'ks','markersize',4);
  legend('T_2','T_2','T_1','T_1 (Bloembergen)');
  grid on;
  set(gca,'XMinorGrid','off');
  set(gca,'YMinorGrid','off');
  title('Relaxation times T_1, T_2 as a function of viscosity','fontsize',fs);
  