FA1=90:-10:20;% degree
FA=FA1*0.75; %0.77 r1=3.28
theta=FA*pi./180;
%% 
[fitresult, gof] = T1cal_t1fit(theta, M1);
t1(1,1)=fitresult.t1;
[fitresult, gof] = T1cal_t1fit(theta, M2);
t1(1,2)=fitresult.t1;
[fitresult, gof] = T1cal_t1fit(theta, M3);
t1(1,3)=fitresult.t1;
[fitresult, gof] = T1cal_t1fit(theta, M4);
t1(1,4)=fitresult.t1;
[fitresult, gof] = T1cal_t1fit(theta, M5);
t1(1,5)=fitresult.t1;
[fitresult, gof] = T1cal_t1fit(theta, M6);
t1(1,6)=fitresult.t1;
[fitresult, gof] = T1cal_t1fit(theta, M7);
t1(1,7)=fitresult.t1;
[fitresult, gof] = T1cal_t1fit(theta, M8);
t1(1,8)=fitresult.t1;
t1=t1/1000%ms to sec
R1=1./t1;
P=polyfit(c,R1,1)
x=0:1;
f=polyval(P,x);
figure
plot(c,R1,'o',x,f,'--')
title('r1')
xlabel('[Gd^3^+] mM')
ylabel('T1^-^1 (sec^-^1)')
%%
% sin(x)*((1-exp(-200/t1))./(1-exp(-200/t1)*cos(x)))