ПЗ / 2 / Computing

clear;
%данные
me=9.10*(10^-31);
h=1.05*(10^-34);
l=0.5*(10^-9);
W1=0;
W2=1.282*10^(-18);
W3=-6.409*10^(-19);
%построение графического решения
%вычисления
f1=zeros();
f2=zeros();
z=-3*(10^-9):0.001*(10^-9):l;
b=1;
for i=1:size(z,2)
f1(i)=tan(l*sqrt((2*me*z(i)*W2)/(h^2)));
%f1(i)=(f1(i))^2;
z1(i)=sqrt(z(i)/(1-z(i)));
z2(i)=sqrt(b*z(i)/(1-b*z(i)));
f2(i)=-(z1(i)+z2(i))/(1-z1(i)*z2(i));
%f2(i)=(f2(i))^2;
if f1(i)==f2(i)
solut(i)=f1(i);
end
end
plot(z,f1,z,f2,'Linewidth',2);
legend('Левая часть','Правая часть','FontSize',14,'FontName','Times New Roman')
xlabel('z','FontSize',14,'FontName','Times New Roman');
ylabel('f_{1}(z), f_{2}(z)','FontSize',14,'FontName','Times New Roman');
%вычисление пси-функций
E=[];
delt=zeros();
A=zeros();
A1=zeros;
A2=zeros();
B=zeros();
C=zeros();
s=zeros();
k1=zeros();
k3=zeros();
x=-0.5*l:0.001*(10^-9):1.5*l;
psi=zeros();
for i=1:size(E,2)
s(i)=(sqrt(2*me*(E(i)-W2)))/(h);
k1(i)=(sqrt(2*me*(W1-E(i))))/(h);
k3(i)=(sqrt(2*me*(W3-E(i))))/(h);
delt(i)=atan(sqrt((E(i))/(W1-E(i))));
A1(i)=4*k1(i)*k3(i)*s*(sin(delt(i)))^2;
A2(i)=2*k3(i)*s*(sin(delt(i)))^2-k1(i)*k3(i)*sin(2*l*s+2*delt(i))+k1(i)*k3(i)*sin(2*delt(i))+2*k1(i)*k3(i)*l*s+2*k1(i)*s*(sin(l*s+delt(i)))^2;
A(i)=sqrt((A1(i))/(A2(i)));
B(i)=A(i)/sin(delt(i));
C(i)=(A(i)*sin(s*l+delt(i)))/(delt(i)*exp(-k3*l));
end
for j=1:size(E,2)
for i=1:size(x,2)
if x(i)<=0
psi(i)=A(j)*exp(x(i)*k1(j));
elseif (0<x(i))&& (x(i)<=l)
psi(i)=B(j)*sin(x(i)*s+delt(j));
else
psi(i)=C(j)*exp(-x(i)*k3(j));
end
end
end