Добавил:
zimbabweLSD
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:файлы матлаб / initial_approx_U5
.m fi0(1)=-log(-N(1)/2+sqrt((N(1)/2)^2+1))+U(j);
fi0(I)=log(N(I)/2+sqrt((N(I)/2)^2+1));
for i=2:I-1
fi0(i)=fi0(i)+(fi0(I)-fi0(i))/(fi0(I)-fi0(1))*(U(j)-U(j-1));
end
Fn0=ones(I,1);
Fp0=ones(I,1);
Fn0(1)=exp(-U(j));
Fp0(1)=exp(U(j));
f0=[fi0; Fn0; Fp0];
%plot(x.*(L0*1e9),fi0.*V0,'LineWidth',LW)
%xlabel('Coordinate, nm','FontSize', FS)
%ylabel('Potential, V','FontSize', FS)
%grid on
fi0(I)=log(N(I)/2+sqrt((N(I)/2)^2+1));
for i=2:I-1
fi0(i)=fi0(i)+(fi0(I)-fi0(i))/(fi0(I)-fi0(1))*(U(j)-U(j-1));
end
Fn0=ones(I,1);
Fp0=ones(I,1);
Fn0(1)=exp(-U(j));
Fp0(1)=exp(U(j));
f0=[fi0; Fn0; Fp0];
%plot(x.*(L0*1e9),fi0.*V0,'LineWidth',LW)
%xlabel('Coordinate, nm','FontSize', FS)
%ylabel('Potential, V','FontSize', FS)
%grid on
Соседние файлы в папке файлы матлаб