last lab / lab

function Task03
clc
clear

%-------------------------------------------------------------
% PARAMETERS OF TASK
m = 9.1e-31; %mass of electron (kg)
h = 1.05e-34; %Plank's constant (J*s)
a = 3.25; %width of potential well (A, Angstroem)
n_max = 3;

% Интервалы
x1 = 0:0.01:a/3;
x2 = a/3:0.01:2*a/3;
x3 = 2*a/3:0.01:a;

for n = 1:n_max
% Функции
y1 = (sqrt(2/a)*sin(n*pi*x1/a)).^2;
y2 = (sqrt(2/a)*sin(n*pi*x2/a)).^2;
y3 = (sqrt(2/a)*sin(n*pi*x3/a)).^2;

% Интегрируем каждый интервал
I1 = trapz(x1, y1);
I2 = trapz(x2, y2);
I3 = trapz(x3, y3);

% SСуммируем вероятности
sum_n = I1 + I2 + I3;

% Выводим результаты
fprintf('For n = %d:\n', n);
fprintf(' Probability in region 1: %.4f\n', I1);
fprintf(' Probability in region 2: %.4f\n', I2);
fprintf(' Probability in region 3: %.4f\n', I3);
fprintf(' Sum for n=%d: %.4f\n\n', n, sum_n);
end

end