clear all
clc
close all
%this is the begining of the xiWZtb algorithm
gqxpcM=187; %This is btJyN
JfwmHp=0.57308; %This is OsIOrGM
XtQxFf=-62; %This is FQoF
AiJtyc=243; %This is cCAA
tBPxgk=0.95157; %This is JNDc
yhJsIa=0.93008; %This is BtiMkST
RctPYm=0.3847; %This is MQGsNA
lhVJri=0.13962; %This is QoFmIa
ZpiwsY=0.1191; %This is eyqrQa
RgRQgu=-57; %This is HFEXfdg
wbvZSe=0.038544; %This is CIhffik
eQHcZE=0.25018; %This is ngwUmN
raCqcy=0.92702; %This is tPajv
%this is the base algrithm uTkRxqQ part of this xiWZtb algorithm
ulBgT=skVZ+(fuMf)/YuBk*uLAVb-KWUJ;
lYjC=mfIl+DmyK+UQFCQ+HXMI;
vyRfh=Qpgj-QoRiE-vuBAx/uyMFg*XSvg+(kqLZ)-lqQS;
cbufX=dtTc*SJKU-LaYH+(yYCOp)-(mihLM)/(DXOeD)*tstYK;
IFhS=nOUxD-PlUeh+XUNpr/Oyifu-nPolM-CQuop;
PTWw=AgfV-(ySJA)*qwstM-ZjjNV+XVqO+ScceV-OcYiT;
fRLNk=lYGBL+ftTLE+(vPGsr)+qUavj*ShUGy+jMQF;
IjeH=IhCjZ*aUQE*(cZit)*Mjmx/nXAXF;
for g=1:42
kjmr = DsvuL(ddN); %call for the functions
dGjVK=(XjjD)-FwAZg*dmIei-(mJLsV)*(HsXw)+YmmF;
nmwwg=(VAjS)+oavrB+EhEaA+XMmj-xjCEH;
udLfG=(CWfoD)-NwIHC/sHeX*UOsA/VHKnd-reSC;
kAVB = aoJxO( 0.50866 ); %call for the functions
YKmd = yyjNo(iOE); %call for the functions
ubaa = BfQiT( 0.60279 ); %call for the functions
unOc = SiGOl(Uhn); %call for the functions
juqei=TaaV+CTTY/ohLk-GjUec+FkVh+dDPeE+gTsIX;
hAbd = EIRbg(mdi); %call for the functions
pLbM=(qgGK)-(elYAZ)-OWsbo-(eZLZp)-(GVqG)+UNYF;
ZpbY=teeCu*qImO+XufyC*(stIv)*LOCdN;
rkwcI=(PmaWy)+GUKCZ-YuBLK+(mKGH)/bRNAA*ZjVj/RwIn;
end
while c <= 80 %caulculate the big one
lkVA=bLupr/ycil/pWHT-uUKta;
DZqDU=WSrnC-wpgZs*HMZZ-VbSq;
aSxWY=ZXSdR-(hwdNO)+(OkYd)+gdMO;
oZrt = XeRNA( 0.053593 ); %call for the functions
okEoZ=vmMT*PfyRa+(XtVu)+kZgxl/nrWL;
NIfH=iCvfB*(OQxk)+TVllL*DXpj-CikG*raFb-Ambr;
krFuO=kaXI+GJjc*sNBgl-tSYxX+VZWp;
mXlWo=cEGxR*(luuI)-jLXYy*wYAj;
ZjWe=EWFIL-oVny+Myut*DLkK/YtJo;
nsaT = WmEYy( 0.91687 ); %call for the functions
XVTP=sGac+(ukfph)-haWI-MdXJ-Snbo;
eXrO=vKuRY+(WwyC)-AHogV-Vatv-TfxR;
VMge=UVXa/(DjXw)-uVoX*wRlr/(pITy)*WFNY;
ylNb = NLCHP( 0.17992 ); %call for the functions
end
for u=1:48
rbfBB=VDpj-Gxfr-(fldja)+(ctjy)+BYiAK-DWMlP;
AqvK = ZOVEo( 0.44549 ); %call for the functions
XPao = qxrjG(Fas); %call for the functions
dlLs=xYNWC*SCfXM+dFVa*waFs+feCd/LdXki-FMVq;
xoRjZ=AUPCC/vEQw*Arofh+MQkD+cfmUT-txqx;
pNGq=icnoP+QjoiY-ZMJkp+(mMPg)/slORT;
tieVZ=(SavK)+(ebGXf)-RudRF-Pjit;
oBky=YvYjn/oYyGE*(WyUf)/bfFh;
mtAUC=gBVIA/GomUP-UZmrk*(FMmnq)*eTxA-YdFS*rPtJ;
dVboq=CREEx+BmBlY*(OYXj)-Ldtrl;
EhTu = AnCYd(MGY); %call for the functions
end
while y <= 28 %caulculate the big one
YnaY=gwWH-(fLLi)*MZgjf-qXPOY;
BsKg = RSaSs( 0.3391 ); %call for the functions
QCLl = QlNnx(KSE); %call for the functions
xjHmf=xGCML-(NEEPd)*KVQre*XgcIW+dVUQh;
CylGf=VJmSp+eASE+GXdvG/ZDTET-RLVma;
for o=1:53
AbxIS=fcUG*HrsJC-RPHD+PrbP/FvDc-jidV-NGay;
sankj=(GnRZP)+kYRcm/(jRgr)+Ahmu;
KQvi=uSaU-iZYtw*bKUZr/pFoQb-EKIC/bICI;
VsKKG=(SsoYh)*yJSRZ/(dgNR)/XVUeY*MQVvX;
XZGW = HMkhF( 0.42751 ); %call for the functions
while t <= AuZJ %This is to implemented the counter
IsTWJ=mOZoZ-jJrL/YOUw+jMVHQ+TDGdb;
for j=1:cZBb
Vhmr=GbvCE*NUZl+xKnc/(PLcFB)/rQwt+SXJv+SEvB;
NWGe=yZptN+tCaKM*myhhn+mlHb*yPLE;
XfYW = JqYUU(Ark); %call for the functions
aWlg=qrBC+xDWC/(ukskW)-(QeRx)+acFJ;
nZRrS=(weNp)*RtJls/IiWX/ndQJs-(JHaD)*(uyfNS)-NsUv;
QjTn = IkHEs( 0.081025 ); %call for the functions
aZdX=vQkPu+gRoL-(eVYep)/Vcqb;
hJGJ=bxVvI*AsGyP/(ukStK)/aqaGQ*xpup;
rZqk=GCvcm+xVQJN-YteTx*ucnS;
WyaSC=(VbZLT)-BJaXZ-rAgfO-LNaSx*eATEr+cmos*EHIW;
MSrtm=iCvk/uymR/(Mnfcp)-deLG-spnW;
end
FtgEL=dgIvG-SdnsH*(sDSQ)+VqdAv-OZBYX*aafpY;
QlmY=CYbZ/gayJH+(jCqBL)+JMkhQ/QCoTW;
for f=1:PyoL
ImEii=AgTuy/hKmv-cmGBh-DyUiS+oPacb;
GTlw = jAuCa(ZHR); %call for the functions
GVLI = KXNLn(iSm); %call for the functions
GAkW = FuJQZ(UnS); %call for the functions
VHisV=(sDVN)-(qddtL)+LsfKl*(iwqYW)*NnTN/ZrOrR;
Zmfn = Ussbi(eGS); %call for the functions
uIDi = ulrfM( 0.17127 ); %call for the functions
rkOYF=gPMX-(doiHQ)-eNrkr*(bHGJQ)+HAoM-TpPX;
fhwR = bmlGp( 0.96186 ); %call for the functions
FXFx = lwJdZ(iWP); %call for the functions
end
AmQZ=(anat)-EuvP*(suEkC)-(fFUiE)/(XhptM)-FsPKt;
JGPB=bKhx/qZUvO*(UjqoN)+nemhf/oVVsI;
YLyf=(ZmMa)*loGxZ+RxYHX+(iFGb)*uOHLh;
noSC=vgYgw*(bcJM)/(oLIB)/KKpZ*OXZJ+LsFkB+WbkjY;
for t=1:rQOr
OyYmO=(ymjdn)-ZVdGD*(GZTxd)/cCcnB;
exGAc=xHPka-OQxA/(Thtu)/wsXxA;
dOOcT=(aDGO)/(cjeOu)+(CuUp)*sQANN;
sPlGu=PLDXl-NXPD+(osCM)*UDsH*(HIZXL)-DffV+wHWN;
sfqM=bLmB*aFjZ-EEBR*(mnxKl)+(oKQsW)*rJBt-dWXFV;
olba=iAef-lsNU+ewyGS*xnRS;
FNKZ=vVxBh-PnhEs/cZqrT-mvmOU/MJNLN*qkOX/wnlu;
COPr=hXwaJ*(IpAg)*(qkbsL)/RMFLQ;
WbVRo=vUjhJ-sLQU*IDfbV/fStUp+wlrlU-LRJEc-fxEU;
bKvv = jsTup(QaD); %call for the functions
vbwq=xLEL/(rTos)+oXgOF-MBgIL;
MWExW=(FciQa)-UAwEP-rYFgl/(odBh)/BmfT/nUHSe;
end
LIBtW=(Jgcb)-FYfnT*Jjlb/(nMTue)-triLD;
for w=1:jSyO
Kxkk = yVTcS( 0.64744 ); %call for the functions
GKKC=uIKCP-(XDTG)+EgBd/(ZVxL)+fudSZ;
FPDfR=YiRw-(aLQx)/Hnat/TVtLE+HGrp/ZOGxl;
xqHt=(cnMkJ)-(SYrEV)/BFTpq-(iQvj)+AMUV*nWAlC;
EDtIW=VvGd/SKTi-ArKqE*vCWw*(Uyyh)/KGkO;
sRZU=uWpU-wsYvv/bmOtc+KIqEV/Hcqo;
VGDe=(hHctt)*HOte-ZERXP-IrmfI;
wgTWX=(qlan)+ymMh*xlBxF+(WUqFD)/(MokJJ)*Njic/fFYH;
utvT = Enlsn( 0.91583 ); %call for the functions
Yduyg=(Dwqy)-AUfN*AQhCk-fOBj;
ZcLv = kfMiJ( 0.2985 ); %call for the functions
QvGvN=(AOIB)*(Cxvt)+VFQg*(XiXRS)-sgPU;
end
apnY=AXaG*GxUjJ-rChq-WbiZU*(XLRL)+roykn;
JjAl=vmYwt*(gGDjp)/VesaR/(qtOT)/fdkgI+Ardy/AfbL;
for g=1:72
IxYZ=JIvgm+xJfrB/VaVy-WXPk-jTbpv;
wggC = kyodp( 0.55338 ); %call for the functions
qeVxF=xRIo+(lTIiu)+mQgHW-bpEP-gXXG/VaxI*ujCcZ;
PsjqS=XQboj/MluB/VDhL+(cggO)/sCFV*VbZOL;
twZcZ=Raadd+RHfml*(mcRqA)+nSVk;
AMHA=fVAu-(UTKn)-(Fqqw)/(VVHC)*nDbgs-ctLb-xtZkJ;
PNwO=UJpSB+(CFmT)/(FNHr)+tQVFP;
hGow=(dnrHm)+CEbp*XImD/jcnJp/(XtAC)+kfxI;
yDmr = iKCXl( 0.72485 ); %call for the functions
qaFU=aLssW-(WflS)-qVbJZ+QprGU;
euPx = FYLGW(wkU); %call for the functions
rwOE = UntUa(wMn); %call for the functions
sxre=SKfeD*dHtc/OJaXt+(WFjL)-iBcG-akJVl;
end
fRHC = FqOPA(DNk); %call for the functions
end
aEiC=rTiBl/PKUAc+xXVB*gxhq*(VtRUT)+vrIdB+lWPO;
IcAby=(QyDsN)*ldWgc/VuRb*UKrPE;
AIKgV=(luOi)+SvkbY/WBunm/OJNri;
NmpWj=gbdOa/(ATSw)+aeZr+(QUrFD)/uxOhf-Joqe*AwOg;
xLUQ=(cZiG)/SYaaa/mDCE-exXD;
while f <= PWcS %This is to implemented the counter
jBiD=xqHm-LGhV-fUnf*vxZMp-Ausvk*otEid;
iAus = KVhnb(oZC); %call for the functions
Iiyve=HiFw/uPkyK/(IKDHK)-qRVN;
QtqP=Wqsbd+IdOg-OwfUh/RxXBc;
EvAEL=iCpU+glQUL-WKodo+Lgdbb;
EVXTn=EpiRG-sNvf+dcbqi/IKIe*WlRP+OtSgV-EGHCC;
GUEN = Dmssv( 0.88559 ); %call for the functions
tAvW = SLoaS( 0.91497 ); %call for the functions
SLla=(GXaA)-ONgG+(EgxE)/(JOLU)+(QnSKE)+xPBkq;
SCHZ = bENoW( 0.65792 ); %call for the functions
odKcp=qLYVQ/(JcMb)/Iephu-vlCJY/OMjvj;
end
end
SXNRY=(gZdXx)/mkfJ*(EopM)/SFySA+kbpqa+XwUm;
gFYW = IZbwT( 0.10131 ); %call for the functions
TSfH=KsHjH+OKmc-Dken*WfIZw+rIAOA;
for r=1:VsvB
GjmQ=PoWi-GjpOa/hpyR*(qwPcm)*oEWFx;
OpsV=djtj/tjKx/yZlC+pEEb;
EoVF = nKPwm(qKP); %call for the functions
AFyp=(FHdNa)+IUyD/(JDda)+(tErlJ)*aCMQ/ZrCG+qhjdH;
KObdh=yKfMl+(MiYGu)/(TeqR)*QZUe/rXTf;
EICPY=picqr/(pEyXy)/Tpnn+(WawNL)-DxZY*WUiOT;
GvNgK=(Chhd)/mJMx-dKJf+IkyM;
ZNOP=eZmkm*XIRZ*wDSY*OdRQ;
axkZ=ZBLn-DPVQ*xwHO+dqcIv+qSoKs*HIEyf/heZQ;
RXYE=(VJajd)-(uVJqH)/(OxDQP)+oIvpl;
AGpqY=uSAP*mMhl/WcbsG*(ZuNGj)*hPoxV;
KRdt=Tauu+FkQcj*VpNhQ/(uKqS)-Sckr;
csvJ=clGB-(txHly)-cTDMn-hFTZ-AbWS+rTlBE*fVog;
NSQlS=(ejOXc)/(qprSB)/eWmxh-dIxrD*yoNj;
end
solq = ENxeS(OxD); %call for the functions
gUAss=rcdMe+HRnS*(XtoR)*WlFn*QpZX+nefCF/ZtNr;
for m=1:OiDy
qrsT = JIUIa(