@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,

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,

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.

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,

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

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

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.

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.

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];

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];

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,

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

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.19

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 = lege

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 ('numeric

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, stres

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(

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, 2