файлы матлаб / grid_functions

x(1)=0;
N(1)=Nef(1);
Mun(1)=MUn(1);
Mup(1)=MUp(1);
for i=1:length(W)
x=[x max(x)+W(i)/Sx:W(i)/Sx:max(x)+W(i)];
N=[N ones(1,Sx).*Nef(i)];
Mun=[Mun ones(1,Sx).*MUn(i)];
Mup=[Mup ones(1,Sx).*MUp(i)];
end
I=length(x);
dx=diff(x);

x=x';
dx=dx';
N=N';
Mun=Mun';
Mup=Mup';