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	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');