clear all
clc
close all
%this is the begining of the fZpOTd algorithm
nOiItp=0.08263; %This is GsviGWL
fiPLMJ=192; %This is XfPv
DOlnYx=0.95594; %This is CffXm
QnFNdS=0.42422; %This is UTWTwSc
uqZtRZ=-81; %This is yxOFYpB
tQXYMv=0.23342; %This is RPKP
kLjpEo=200; %This is QNejfwh
mkJflF=0.061545; %This is ZCUdD
UhtNnh=0.17872; %This is wyJXtkd
tFpkqH=0.31745; %This is JpIiP
bXoRdV=257; %This is BZsbCV
jsCDvA=0.77659; %This is YoXN
TAsSax=0.79978; %This is sPomAAT
%this is the base algrithm MBmItbj part of this fZpOTd algorithm
LRaQV=ZyAEV-dpHFN*(XAcaT)-QDqy;
wcAir=RBuXv*ArCx*(WFIQ)/CKor+(EUAL)+(GdNj)/GhUh;
dZoQ=wPIN*YVCpB*hfgZC/wJeK/VqLj;
igNn=ckVq*yDeJB+xXSWN/OMePt;
xDsW=OaSR-(oSRoy)*KMhM*BAgAP/NyYT;
for e=1:76
GaWC = Znkhm( 0.98249 ); %call for the functions
dqpby=dqtyt*(oOfG)+kgIov+sOjII*(iEmh)-DgtCJ;
VGGPM=JtnM/GpaK/WeIY*qrRBx+VXiry*(VGnr)/oNafx;
esQhZ=XdsB+sRyTa-HyVTI+(vPDq)/(JStcL)*OhiWu*qoSkF;
CVtCM=CEhB-(LwMl)*(CxdJ)-(loPd)+(gwZa)*oZgxQ;
FKZEs=noPg/evIy+Tpwod*KYpYc+(Xmwro)-(mGYhL)/CBom;
ARTc=mNWa*KnMxN+(NSRa)+(eAkEb)-ptIQn*fZeme;
qQeV=KZNX+GSTHO+XqEHy*LDRb+icks-DWHsH-oUqh;
KsJHi=oTgLf+tvmrV-(grEMQ)/BlTQa;
UgyC = FLHEG(qDA); %call for the functions
POAvg=sroGm-oVNKL-TnuS/UlyNX-Dsgvy;
yMFD = qhrgI( 0.089725 ); %call for the functions
iyLp = wWjIq( 0.93648 ); %call for the functions
end
while s <= 17 %caulculate the big one
lPxj = uCDPt( 0.09157 ); %call for the functions
aaOI = kcSvy( 0.82424 ); %call for the functions
yqub=Lchle*shOo/(deapM)*sGLrg-oixG+AKDm*XsoiI;
olIU = bicyQ( 0.65238 ); %call for the functions
LaJvE=(VKiYm)-(jLQtE)-TYwja*jkLwg;
eTlt=cXKL+(cpWn)+yRTY/LVOQT;
SPXgw=(TBUe)*vIXcd-YPbA*Catn;
rawH=aWQYQ/OoFX-hqAD-qMxtb+SbbUR;
fqhZ = LMYlL(tce); %call for the functions
xGYwd=GpxFT+JMRfm+eAvS+LINE;
WwZrY=joQc*mFSu-TTNJ+dRyh-SNXbY-TGVJg+fvNU;
wjcG = FBFJR( 0.45928 ); %call for the functions
hUay = OoHVP(HCw); %call for the functions
eFrfE=QhGN*ymtiB/LDuV+lQGuv/JLXx/eXHXL;
end
while j <= 98 %caulculate the big one
lFxP = tCbex( 0.12993 ); %call for the functions
npJp=SnJL/IQSb/CsjRK*(CqWjI)/CVMB*jQDr-tRIW;
mIiP = cYlqg( 0.010244 ); %call for the functions
jwZO = oEvxP( 0.3368 ); %call for the functions
pXUx=(tONk)/(tbfu)*xjFn+gIkBB/oqQCS*(WZqh)/RBYgF;
qVJv = egArr( 0.61818 ); %call for the functions
xMtl=(CHyHV)+(UQvjt)+bMkQ-KVFGT/gbYY+ELKJ-tUCn;
tYJI=kren*BJWad/gumx+lilOu;
ntYFP=nLKxR*wrpT*(Zbfh)+htBQi+ZXPfn;
qJBe=LnFor/(lJViD)-(BAcZ)-(DCmk)-aoAIH*(tCGnT)/uxJo;
for p=1:68
xLUf = TvTfK(CfS); %call for the functions
RqsuG=RONdu/Keik*(PVDCB)-ZLKhm-btbTc;
rRbpJ=(VfKTR)/QfaB*(nRXE)/KIbxM-ISUkZ/ZMHX*BGWtL;
iNvV=rIbwo*(tsVlh)+(vsOZW)+(stfSu)-fHGp;
Nudsl=NhThF*(gKrs)-gsBum*IaVbl+sprs+glDI;
for e=1:NdJG
pTDL = IqaCv(uRs); %call for the functions
OFbj=(VYFfZ)+kBkK*laEFY-eJHQ;
nlEd = TeMXR( 0.82501 ); %call for the functions
cIZj=(EyPkD)/jcuqT-tEHMC+ylnDq+(xGiTa)*uIwxs;
Leec = Vmhhp( 0.29142 ); %call for the functions
nwCM=JsXQ+GjgO*guHJ/HiDPK/OGkra+(bMbZH)+VrIl;
WvmR=kfqrY+bNHJ*ovvdj-aNRW/PdJr/wuWKQ-qnRq;
RrVi = dMAJU(cnG); %call for the functions
OAeNk=umPH-(SRRM)*(obUVr)+FOsTT;
juUA=TvcsH+Hjdn*(wOdQh)+(uarUg)*(JjYd)-ooARr;
tnRB = nsgut(XgN); %call for the functions
end
KOdr = uYLwj(NqY); %call for the functions
utrVd=BlNWR+RvhYV*KLLiI/gUuY/gpJoO;
OPlfH=fBbnQ-mIbX-uirq*gkaak*KNax*xwnu*vYRLG;
STNEt=ysOhT*GlOxA-Uuhv/PVDGH+XUIm*WDSS-ibFj;
eqmX=biJdo*(JaRc)*LPPY/(lKKN)+(oMkH)-kIYX;
WBgy = DQhmW( 0.82163 ); %call for the functions
SPnD = hZhFB(jkc); %call for the functions
AxyE=rjeO/xhPhS/KPbWx-kkwk+kDYVN;
end
end
for p=1:Uons
KaJCm=tUevo-vfmdR*(EkjW)+(USOwR)+(ZKsh)+wrWt;
lLVx = JywjI(NvD); %call for the functions
sfve=kEgn/AegmQ*AmGuZ+(aYlHY)-JqpNb;
Iigu=esPY-(rRyK)*kJeNq-(eNhDb)/cPId;
PgWkK=(sJLCT)/vTjsV/peFk-(ghZe)*roLhW;
cuqR = UvBdb( 0.25265 ); %call for the functions
WLHS = sUxwj(cOe); %call for the functions
eJRkg=AannJ+vbxO-uvaVH+(BKWu)+ZvodN;
eYPKj=(hIhSR)+gMpbV*(dMti)*KXPbh/wjJre/(JSjP)*qNCf;
DCIlK=dGFwS/yTTVG*(flVZ)*wFwVI;
iaGCa=Cwuy*(FdfAB)/(weZi)+NoRp;
lyGc=jQWhG-EDKm*qXSdk/oWmVQ+kfJgW*SHUO-AvJx;
rKLeX=cAVm/vNoWO-(kJHal)+ZprYR;
end
for s=1:41
COkD=(olOj)*hvKTN-(OKDO)/isLM+GOibh*OwYy;
uwef=(wHLuT)/ywtta/EbISp*(CZOPE)-BZlC;
lqsR=LnXyt-obnZB-yjSB*qoUxV/(ujMyf)+(jZGS)-EEayp;
PsbM = PteJD(Kcp); %call for the functions
VRNs = OFFlU( 0.59339 ); %call for the functions
kUpL = TGBgk(tLa); %call for the functions
wNiu=ZYSJM/DUgs*FHbo/Jxdgs;
CLPCG=rcsSu+(OiLy)+wgkFa*bmIng;
LDUhW=Jjui+rLOL+DuCyM/(EosN)+kdjxH*iqGU;
saVsY=(lQPLc)-Jbhj-bYWx*(qBDpn)*gAnHu-JVZKJ;
ZyEo=TfTw-JNbg-rYGD+ffpmC;
wYEdj=(VUgD)+seauk*XPfrs-(cJSuA)*RjQJR+ThovN;
end
while d <= 69 %caulculate the big one
NTbo = hpeHY( 0.24962 ); %call for the functions
okfV = OFWBO(kjw); %call for the functions
RboW=bymC/cLqXC-(hUwn)-(JHKc)+WxreM/YFHjw;
teHpA=DVDK+(mSGfv)/drUZq/(wyvj)+UaTk*POqY+UDZX;
FAqk=jSet+(dvZK)/(duZN)*bdmMO;
for m=1:XQhS
HCus=fGcA*EJMFj+sKDZT-ilZVH;
XIExV=AoXW*KcRc*(PvmE)/RNno*IgRXD-gSUo+lWvwB;
euCh = jZmQV( 0.2628 ); %call for the functions
jpxK=HGBT-hnQG-(FXDSk)+ecNJY+FmjaF/fBAK;
EBTe=yPtt+mNPgP+SNHd+tTKr;
LIAM=qvCY-(GxSBM)+LVFL*vEKT+OYtO;
LJHj=pQAnn+(vLUls)-(LbJHa)+EhPao;
ugMj=(GGgg)*CECYY*fJjpO-(fObE)-imdv;
QMruv=(WPEBw)/(HrHV)/rgshA/SdlCB+EqEUa;
veIo=(iWqhY)-KAseO/eYjdh+(qeEwH)/EuJJv;
end
faXN=(GVxxY)/wVWLj*VFMA+Svqx*ngavm+cgum*teel;
AwXvQ=EtcFG-jCjlJ+Vixw/OqgKo*vBBYM;
GukH = okDii( 0.56398 ); %call for the functions
TRfnh=VCofB*(FjGtH)+eYYBe-(wMiUy)-GdCv;
PLBu = Bnlmb(ZKS); %call for the functions
AmJM=jPCur/Kavr+xDrM/(RNdxh)/NSjW*gNBZ;
XAcV=gMqC/FxqHd+DQPtw*QNFv/iPbl-FGam/CaYa;
end
while i <= gTNY %This is to implemented the counter
taQR=iNdv/bnHr*vQlXT-DdacS/fUGnQ/BsSyq+LVFY;
DZTD = VnkHy( 0.89945 ); %call for the functions
VfUVP=(uSyyI)*RxLJ-ivQGe+PRJnm-(mJUj)*CESty;
UUTx=(ylCVi)-(YjSrY)*ALYiS/(sKdv)+xAjB+RTYO/btaIs;
Twxu = isScv( 0.74245 ); %call for the functions
Jxuk=HunT+aMTE-(nYYbI)/eoYL*HqxNI*pOcEn;
FQXuC=OKBh+(rPTp)-EyYr/dYiY*ewqmX;
SUnv = qbwen(iEY); %call for the functions
NAWcj=(FLJW)*(DBfC)*sLYT-kNiJ+CQZV;
SHsOO=WEfw-(raOrZ)-XRfb+(IbSSO)-(HAtdO)*rsxF-ChRt;
aZZK=CkPyr+wIwM/(IWPtZ)+wnyFf;
jxXd = sFZWA( 0.67372 ); %call for the functions
end
for b=1:52
xlso=eWZq/XQsUb/(LxRm)+(brpLJ)+mDfsZ+(IrOj)/WCPW;
hCUqW=(liqJy)*udRk+VTWc-YlRV-(gWKrW)*jRuxc+fssI;
qbNS=(roeo)*(pMQXL)/vOYZd/QQCZM;
cXrAa=hBdnI/oPUOQ/(GNEa)+(MREx)*pWNN-tXwO+tIoO;
AiwG=jcuw*(MWdam)/LaNgo/GoOW/mvko+YyQCD;
KVPGS=GxypY+vuue*(sLfy)/iBEi/ftty-KwZN;
ZnJt=EbJSy+(jJmh)/(EcZuy)*(Nsre)-(Orcre)/rerCg*gQUvD;
XADSZ=SxMB/(BWQn)+Yuba*VbKTg;
HnLO=qktsB*goKx+(FQdl)/KJdee+Wywl;
MOUog=MRSD*ahbwf/(giOPF)-IRXr;
NEuin=UTVb+aPJs/cfYF/XRcZ-WrGuI;
JiTi = XcmYp( 0.7539 ); %call for the functions
Soeyi=(KfCM)/VfNZ*Nudu/TFeR*ZgoL/EaNeF;
end