P
@ProfB
Experiment 8a
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 = [0.71, 0.71, 0.77, 0.77, 0.67, 0.73, 0.67, 0.77, 0.81, 0.71, 0.67, 0.67, 0.59, 0.56, 0.63, 0.61, 0.69, 0.73, 0.76, 0.61, 0.76, 0.67, 0.60, 0.66, 0.76, 0.68, 0.64, 0.76, 0.74, 0.68, 0.66, 0.62, 0.71, 0.73, 0.68, 0.64, 0.69, 0.74, 0.68, 0.78];
v = [1.1
Experiment 8b
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 = [0.71, 0.71, 0.77, 0.77, 0.67, 0.73, 0.67, 0.77, 0.81, 0.71, 0.67, 0.67, 0.59, 0.56, 0.63, 0.61, 0.69, 0.73, 0.76, 0.61, 0.76, 0.67, 0.60, 0.66, 0.76, 0.68, 0.64, 0.76, 0.74, 0.68, 0.66, 0.62, 0.71, 0.73, 0.68, 0.64, 0.69, 0.74, 0.68, 0.78];
v = [1.1
Experiment 8c
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 = [0.71, 0.71, 0.77, 0.77, 0.67, 0.73, 0.67, 0.77, 0.81, 0.71, 0.67, 0.67, 0.59, 0.56, 0.63, 0.61, 0.69, 0.73, 0.76, 0.61, 0.76, 0.67, 0.60, 0.66, 0.76, 0.68, 0.64, 0.76, 0.74, 0.68, 0.66, 0.62, 0.71, 0.73, 0.68, 0.64, 0.69, 0.74, 0.68, 0.78];
v = [1.1
Experiment 4a
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 ('numerical fit', 'Data', 'location', 'southeast');
set (h, "fontsize", 16);
Experiment 3
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)
Experiment 2
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', 'location', 'southeast');
set (h, "fontsize", 16);
Experiment 1
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', 'Data', 'location', 'southeast');
set (h, "fontsize", 16);