kr = [0.01 0.1 1 10 100 250 1000];
t2 = [0.01 0.1 1];
I0 = [];
I1 = [];
for i = 1 : length(kr)
I0(i) = (0.25 * kr(i) + 0.5) / (1 * kr(i));
for j = 1 : length(t2)
I1(i, j) = I0(i) + t2(j) * (-1);
end
end
I0
I1
t2 = [0.01 0.1 1];
I0 = [];
I1 = [];
for i = 1 : length(kr)
I0(i) = (0.25 * kr(i) + 0.5) / (1 * kr(i));
for j = 1 : length(t2)
I1(i, j) = I0(i) + t2(j) * (-1);
end
end
I0
I1