@ProfB

Y = 2021:2100; C = zeros(1, 80); M = zeros(1, 80); N = zeros(1, 80); T = zeros(1, 80); for i = 1:length(Y) C(i) = 424 + (Y(i) - 2021)*2.4; M(i) = 1911.8 + (Y(i) - 2021)*12;

in=[0, 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85]; r1 = [0, 4, 7, 11, 15, 18, 22, 25, 29, 32, 35, 37, 41, 43, 45, 46, 48, 48]; r2 = [0, 4, 7, 11, 15, 18, 22, 25, 30, 32, 35, 38, 41, 43, 45, 46, 48, 49]; r3 = [0, 4, 7, 11, 15, 18, 22, 25, 28, 32, 35, 38, 41, 43, 45, 46, 47, 48]; dr1 = [0.5, 0.5, 0.5, 0.5, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2]; dr2 = [0.5, 0.5, 0.5, 0.5, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2]; dr3 = [0.5, 0.5, 0.5, 0.5, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,

D = [0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0]; y = [0.022, 0.022, 0.030, 0.032, 0.037, 0.044, 0.059, 0.034, 0.052, 0.058, 0.061, 0.059, 0.062, 0.081, 0.085, 0.075]; lambdaD = 650e-9 * D; p = polyfit(y,lambdaD,1); plot(y,polyval(p,y), 'r-'); hold on; plot(y,lambdaD,'ko'); h = legend ('numerical fit', 'Data', 'location', 'southeast'); set (h, "fontsize", 16);

time1 = [0.53, 0.35, 0.79, 0.82, 0.51, 0.48, 0.68, 0.72, 0.57, 0.51, 0.40, 0.66, 0.51, 0.50, 0.73, 0.51, 0.58, 0.51, 0.62, 0.60]; time2 = [0.87, 1.37, 1.24, 0.98, 1.18, 0.82, 1.46, 1.08, 1.22, 0.95, 1.16, 1.22, 1.19, 1.44, 1.00, 1.31, 0.98, 1.49, 1.24, 1.16]; hist (time1, 5, "facecolor", "r", "edgecolor", "k"); hold on; hist (time2, 5, "facecolor", "b", "edgecolor", "k"); xlabel('time (s)', 'FontSize', 16) ylabel('frequency', 'FontSize', 16) axis('tight') set(gca,'FontSize',16)

V=[3.1, 3.2, 3.3, 3.4, 3.5, 3.6, 3.7, 3.8, 3.9, 4, 4.1, 4.2, 4.3, 4.4, 4.5, 4.6, 4.7, 4.8, 4.9, 5, 5.1, 5.2, 5.3, 5.4, 5.5, 5.6, 5.7, 5.8, 5.9, 6]; I=[0.079, 0.085, 0.085, 0.094, 0.094, 0.097, 0.095, 0.097, 0.098, 0.109, 0.106, 0.109, 0.112, 0.115, 0.119, 0.123, 0.121, 0.122, 0.132, 0.132, 0.131, 0.131, 0.136, 0.140, 0.141, 0.149, 0.151, 0.147, 0.152, 0.157]; p = polyfit(V,I,1); plot(V,polyval(p,V), 'r-'); hold on; plot(V,I,'ko'); h = legend ('numerical fit', 'Data', 'location', 'southeast');

v = [3.35, 3.34, 3.29, 3.21, 3.10, 3, 2.89, 2.78, 2.70, 2.65, 2.64; 3.52, 3.54, 3.46, 3.32, 3.16, 3, 2.83, 2.67, 2.53, 2.45, 2.47; 3.84, 3.94, 3.77, 3.48, 3.22, 3, 2.77, 2.51, 2.22, 2.05, 2.15; 4.49, 5.18, 4.41, 3.70, 3.30, 3, 2.69, 2.29, 1.58, 0.81, 1.50; 5.22, 6, 5.14, 3.83, 3.33, 3, 2.66, 2.16, 0.85, 0, 0.77; 4.49, 5.18, 4.41, 3.70, 3.30, 3, 2.69, 2.29, 1.58, 0.81, 1.50; 3.84, 3.94, 3.77, 3.48, 3.22, 3, 2.77, 2.51, 2.22, 2.05, 2.15; 3.52, 3.54, 3.46, 3.32, 3.16, 3, 2.83, 2.67, 2.53, 2.45, 2.

v = [3.35, 3.34, 3.29, 3.21, 3.10, 3, 2.89, 2.78, 2.70, 2.65, 2.64; 3.52, 3.54, 3.46, 3.32, 3.16, 3, 2.83, 2.67, 2.53, 2.45, 2.47; 3.84, 3.94, 3.77, 3.48, 3.22, 3, 2.77, 2.51, 2.22, 2.05, 2.15; 4.49, 5.18, 4.41, 3.70, 3.30, 3, 2.69, 2.29, 1.58, 0.81, 1.50; 5.22, 6, 5.14, 3.83, 3.33, 3, 2.66, 2.16, 0.85, 0, 0.77; 4.49, 5.18, 4.41, 3.70, 3.30, 3, 2.69, 2.29, 1.58, 0.81, 1.50; 3.84, 3.94, 3.77, 3.48, 3.22, 3, 2.77, 2.51, 2.22, 2.05, 2.15; 3.52, 3.54, 3.46, 3.32, 3.16, 3, 2.83, 2.67, 2.53, 2.45, 2.

data1 = [0.66, 0.73, 0.79, 0.80, 0.81, 0.89, 0.94, 0.95, 0.96, 0.99, 1.03, 1.09, 1.09, 1.10, 1.10, 1.14, 1.15, 1.22, 1.27, 1.29]; data2 = [1.66, 1.67, 1.78, 1.79, 1.91, 1.91, 1.98, 1.99, 2.00, 2.02, 2.05, 2.11, 2.12, 2.18, 2.20, 2.23, 2.26, 2.28, 2.33, 2.35]; hist (data1, 5, "facecolor", "r", "edgecolor", "k"); hold on hist (data2, 5, "facecolor", "b", "edgecolor", "k"); xlabel('Data (units)', 'FontSize', 16) ylabel('Frequency', 'FontSize', 16) axis('tight') set(gca,'FontSize',16)

data1 = [0.66, 0.73, 0.79, 0.80, 0.81, 0.89, 0.94, 0.95, 0.96, 0.99, 1.03, 1.09, 1.09, 1.10, 1.10, 1.14, 1.15, 1.22, 1.27, 1.29]; data2 = [1.66, 1.67, 1.78, 1.79, 1.91, 1.91, 1.98, 1.99, 2.00, 2.02, 2.05, 2.11, 2.12, 2.18, 2.20, 2.23, 2.26, 2.28, 2.33, 2.35]; sorted1 = sort(data1); sorted2 = sort(data2); for i = 1:1:20 cdf(i) = i/21; end plot(sorted1,cdf,'ko-'); hold on plot(sorted2,cdf,'ro-');

data = [3.2, 5.3, 1, 2.5, 6.6, 4.5, 3.4, 2.2, 9.1, 3.7, 1.3, 4.5, 4.3, 6.5, 3.5, 4.2, 4.8, 5.8, 7.5, 2.4, 3.9, 8.9, 6.2, 7.1, 6, 3.5, 4, 5.7, 2.4, 1.2, 3.9, 5.6, 5.6, 3.5, 7.6, 4.9, 7.5, 8.6, 7.2, 3.8, 4.4, 6.4, 5.5, 3.8, 6.4, 4.5, 5.1, 5, 3.7, 2]; sorted = sort(data); for i = 1:1:50 cdf(i) = i/51; end plot(sorted,cdf,'ko-'); xlabel('Data (units)', 'FontSize', 16) ylabel('Cumulative distribution function', 'FontSize', 16) axis('tight') set(gca,'FontSize',16)

data = [3.2, 5.3, 1, 2.5, 6.6, 4.5, 3.4, 2.2, 9.1, 3.7, 1.3, 4.5, 4.3, 6.5, 3.5, 4.2, 4.8, 5.8, 7.5, 2.4, 3.9, 8.9, 6.2, 7.1, 6, 3.5, 4, 5.7, 2.4, 1.2, 3.9, 5.6, 5.6, 3.5, 7.6, 4.9, 7.5, 8.6, 7.2, 3.8, 4.4, 6.4, 5.5, 3.8, 6.4, 4.5, 5.1, 5, 3.7, 2]; hist (data, 10, "facecolor", "r", "edgecolor", "b"); xlabel('Data (units)', 'FontSize', 16) ylabel('Frequency', 'FontSize', 16) axis('tight') set(gca,'FontSize',16)

t=[30, 60, 90, 120, 150, 180, 210, 240, 270, 300, 330, 360, 390, 420, 450, 480, 510, 540, 570, 600, 630, 660, 690, 720, 750, 780, 810, 840, 870, 900, 930, 960, 990, 1020, 1050, 1080, 1110, 1140, 1170, 1200, 1230, 1260, 1290, 1320, 1350, 1380, 1410, 1440, 1470, 1500, 1530, 1560, 1590, 1620, 1650, 1680, 1710, 1740, 1770, 1800]; T=[26.62, 29.09, 31.41, 33.60, 35.66, 37.60, 39.43, 41.15, 42.78, 44.30, 45.74, 47.10, 48.37, 49.57, 50.70, 51.77, 52.77, 53.72, 54.61, 55.45, 56.24, 56.98, 57.68, 58.34,

T=[45.5, 49.6, 51.2, 46.0, 45.4] mean(T) std(T)

M=[0.02, 0.04, 0.06, 0.08, 0.1, 0.12, 0.14, 0.16, 0.18, 0.2]; T=[1.411, 1.424, 1.431, 1.428, 1.419, 1.430, 1.417, 1.422, 1.430, 1.432]; lM = log(M); lT = log(T); p = polyfit(lM,lT,1) plot(lM,polyval(p,lM), 'r'); hold on; plot(lM,lT,'k*'); h = legend ('numerical fit', 'Data', 'location', 'southeast'); set (h, "fontsize", 16);

L=[0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5]; M=[0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.02, 0.04, 0.06, 0.08, 0.1, 0.12, 0.14, 0.16, 0.18, 0.2]; T=[0.903, 1.011, 1.095, 1.193, 1.262, 1.338, 1.412, 1.488, 1.562, 1.632, 1.411, 1.424, 1.431, 1.428, 1.419, 1.430, 1.417, 1.422, 1.430, 1.432]; m = 0.49777; n = 0.0034532; x = L; for l = 1:length(L) x(l) = (L(l)^m)*(M(l)^n); end p = polyfit(x,T,1)

f = [0.20, 0.39, 0.59, 0.78, 0.98, 1.18, 1.37, 1.57, 1.77, 1.96]; x = [0.017, 0.034, 0.051, 0.068, 0.085, 0.103, 0.118, 0.135, 0.154, 0.171]; lf = log(f); lx = log(x); p = polyfit(lx, lf,1) plot(lx,polyval(p,lx)); hold on; plot(lx,lf,'k*'); h = legend ('numerical fit', 'Data', 'location', 'southeast'); set (h, "fontsize", 16);

m=[0.02, 0.04, 0.06, 0.08, 0.1, 0.12, 0.14, 0.16, 0.18, 0.2]; T=[0.25, 0.37, 0.46, 0.52, 0.57, 0.65, 0.69, 0.78, 0.76, 0.80]; lm = log10(m); lT = log10(T); p = polyfit(lm,lT,1) plot(lm,polyval(p,lm), 'r'); hold on; plot(lm,lT,'k*'); legend ('numerical fit', 'Data', 'location', 'southeast'); xlabel('log(m)', 'FontSize', 16)

strain= [0.000, 0.014, 0.028, 0.042, 0.057, 0.071, 0.085, 0.099, 0.113, 0.127, 0.141, 0.156, 0.170, 0.184, 0.198]; stress = [0.000, 0.214, 0.432, 0.656, 0.863, 1.024, 1.282, 1.571, 1.659, 1.957, 2.229, 2.371, 2.686, 2.664, 2.984]; plot(strain, stress,'-k*'); xlabel('strain', 'FontSize', 16) ylabel('stress (MPa)', 'FontSize', 16) axis('tight') set(gca,'FontSize',16)

strain= [0.000, 0.014, 0.028, 0.042, 0.057, 0.071, 0.085, 0.099, 0.113, 0.127, 0.141, 0.156, 0.170, 0.184, 0.198]; stress = [0.000, 0.214, 0.432, 0.656, 0.863, 1.024, 1.282, 1.571, 1.659, 1.957, 2.229, 2.371, 2.686, 2.664, 2.984]; toughness = trapz(strain, stress)

h = [0.5, 0.48, 0.45, 0.44, 0.42, 0.39, 0.37, 0.36, 0.34, 0.32, 0.31, 0.29, 0.27, 0.26, 0.24, 0.23, 0.21, 0.2, 0.18, 0.17, 0.16, 0.15, 0.13, 0.12, 0.11]; t = [0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 110, 120, 130, 140, 150, 160, 170, 180, 190, 200, 210, 220, 230, 240]; p = polyfit(t, h,2) plot(t,polyval(p,t)); hold on; plot(t, h,'k*'); h = legend ('numerical fit', 'Data', 'location', 'northeast'); set (h, "fontsize", 16); xlabel('t (s)', 'FontSize', 16) ylabel('h (m)', 'FontSize', 16)