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