Dabeen Park/Matlab/MRI

x=TR;
S0=0.0001;
T1=0.0002;
resultS0=zeros(128,128);
resultT1=zeros(128,128);
for i=1:128
    for j=1:128
        y=squeeze(img(i,j,1:16));
        result=fit(x',y,'S0*(1-exp(-x/T1))', 'StartPoint', [S0,T1]);
        resultS0(i,j)=result.S0;
        resultT1(i,j)=result.T1;
    end
end
surf(resultS0)
view(2)