@ProfB

L = [0.93, 0.99, 1.06, 1.03, 1, 0.9, 0.95, 0.96, 1.05, 0.95, 0.9922, 1.034, 0.98, 0.94, 0.97, 0.93, 0.94, 0.79, 0.94, 0.89, 0.91, 0.91, 0.84, 0.84, 0.91, 0.91, 0.95, 0.99, 0.96, 0.91, 0.91, 0.86, 0.96, 0.94, 0.90, 0.96, 0.89, 1.04, 0.99, 0.99]; d =

L = [0.93, 0.99, 1.06, 1.03, 1, 0.9, 0.95, 0.96, 1.05, 0.95, 0.9922, 1.034, 0.98, 0.94, 0.97, 0.93, 0.94, 0.79, 0.94, 0.89, 0.91, 0.91, 0.84, 0.84, 0.91, 0.91, 0.95, 0.99, 0.96, 0.91, 0.91, 0.86, 0.96, 0.94, 0.90, 0.96, 0.89, 1.04, 0.99, 0.99]; d =

L = [0.93, 0.99, 1.06, 1.03, 1, 0.9, 0.95, 0.96, 1.05, 0.95, 0.9922, 1.034, 0.98, 0.94, 0.97, 0.93, 0.94, 0.79, 0.94, 0.89, 0.91, 0.91, 0.84, 0.84, 0.91, 0.91, 0.95, 0.99, 0.96, 0.91, 0.91, 0.86, 0.96, 0.94, 0.90, 0.96, 0.89, 1.04, 0.99, 0.99]; d =

L=[0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65]; T=[0.903, 1.011, 1.095, 1.193, 1.262, 1.338, 1.412, 1.488, 1.562, 1.632]; lL = log(L); lT = log(T); p = polyfit(lL,lT,1) plot(lL,polyval(p,lL), 'r'); hold on; plot(lL,lT,'k*'); h = legend ('

y = [0.228 0.189 0.127 8.28e-2 4.93e-2 -9.757e-3 -6.161e-2 -0.108 -0.143 -0.211]; t = [30.705 30.739 30.775 30.808 30.841 30.875 30.908 30.941 30.975 31.008]; p = polyfit(t, y, 1)

t=[0.040, 0.080, 0.120, 0.160, 0.200, 0.240, 0.280]; y=[0.009, 0.036, 0.085, 0.158, 0.251, 0.365, 0.501]; lt = log10(t); ly = log10(y); p = polyfit(lt,ly,1) plot(lt,polyval(p,lt), 'r'); hold on; plot(lt,ly,'k*'); h = legend ('numerical fit', 'Data',

t=[1 2 3 4 5 6 7 8 9 10 11]; C=[0.28 0.399 0.49 0.569 0.636 0.699 0.755 0.809 0.859 0.906 0.951]; logt = log10(t); logC = log10(C); p = polyfit(logt,logC,1) plot(logt,polyval(p,logt), 'r'); hold on; plot(logt, logC,'k*'); h = legend ('numerical fit'