enderton/Diagnostics/ZcordA_FizhiZ30.m

nklev=30;
ZC=zeros(1,nklev);
ZF=zeros(1,nklev+1);
drC = [3500.*ones(1,28),1000.*ones(1,2)];
ZC(1)=100000.-(0.5*drC(1));
ZF(1)=100000.;
for klev=2:nklev;
 ZC(klev)=ZC(klev-1)-(0.5*(drC(klev)+drC(klev-1)));
 ZF(klev)=ZF(klev-1)-drC(klev-1);
end
ZF(nklev+1)=ZF(nklev)-drC(nklev);
drF = ZF(1:nklev)-ZF(2:nklev+1);
save('ZcordA_FizhiZ30.mat','ZC','ZF','drC','drF');
1	molod	1.1	nklev=30;
2			ZC=zeros(1,nklev);
3			ZF=zeros(1,nklev+1);
4			drC = [3500.ones(1,28),1000.ones(1,2)];
5			ZC(1)=100000.-(0.5*drC(1));
6			ZF(1)=100000.;
7			for klev=2:nklev;
8			ZC(klev)=ZC(klev-1)-(0.5*(drC(klev)+drC(klev-1)));
9			ZF(klev)=ZF(klev-1)-drC(klev-1);
10			end
11			ZF(nklev+1)=ZF(nklev)-drC(nklev);
12			drF = ZF(1:nklev)-ZF(2:nklev+1);
13			save('ZcordA_FizhiZ30.mat','ZC','ZF','drC','drF');