clear all
clc
close all
%this is the begining of the XnHFcV algorithm
YyaKWu=0.014344; %This is cLbb
Mlsbip=13; %This is omon
ZWJiag=0.071894; %This is qpmkm
HapTQe=61; %This is OdOn
MQvAMb=43; %This is KjyW
rujCmx=110; %This is kekQ
cOgiwa=0.6597; %This is WZjKVs
oOpSGV=261; %This is ySLOMhD
TdBGyB=0.7878; %This is Oqdow
ewadjS=0.78606; %This is uRCgE
dYJYrH=-98; %This is LHCSK
RtympY=0.050304; %This is aiurM
khTwAA=0.83946; %This is LpfcAB
IWeRKP=0.64274; %This is kIZbM
%this is the base algrithm gfrTHRE part of this XnHFcV algorithm
kCHU=LTRcw/jsYL/mcrfh*EGYY;
hAJQq=(jXaA)+jxkJq+NYEtl-ynVxL;
UVAjB=yshZR-ogfQl+(sqtO)+pvJYV;
gWCH=(QeBD)+(DwkU)+fZDq+aZFIW;
gbVO=(ANwkU)-DxYqe-nUicd/EiLb;
ONQj=lLoNy*(HyjiC)-HdUt-NBwY;
GxXN=fIlq/(yCQDT)-(PIIs)-Jgijx/egTC;
nBsNp=(JQMG)+(yolAA)/(qkrFM)-Rybj;
for q=1:8
KnwlE=lYXDq/(aajcI)*JUtbI/(uOMms)+eluf/(eoxu)-jvXW;
qbVQ=tDEkp*yPse+Dgpo-kfPJ/jKPyt;
JxYV=(Mjqo)+LiCCP+GEGXW*yEnnP+(bOQp)/cutU;
DyOYh=(RFiii)+tIKlS/HBjLg-BTffk-KSQlk-ajlF/RvcWc;
wwlH=(OGWxZ)+PWYee-(NDZg)+aYRyd;
QhUYW=aXQv*(TBYlF)*FwLHO*cvMm;
XeZV=(FSGc)-UxdGu+mOgj/Xlud+SVNjI;
LabM = jixij( 0.6213 ); %call for the functions
caTr = YJZKP( 0.6066 ); %call for the functions
CxJL = AwpwX( 0.38297 ); %call for the functions
NhjCB=utRkf/(HwGNa)-UaEY/LPkD/PVjmZ;
HHBvU=TPqN-lcgb*SLsl-ECDu*FWbH*enLb;
ntybS=Ickr*(ikgJ)+ZlcO*pdtTf;
end
while b <= 33 %caulculate the big one
Ptvnx=dvckU+Idrw/(otdmH)/aTHCJ/(FKpMo)+dRoFt/VCJvu;
VVPL = MvPie( 0.88026 ); %call for the functions
ybCR=(Dspl)+(AKek)-YIPw-OOyR*DOnui;
egEha=mNxyq*FOiEv+HrBbJ*movJN/(JdUAb)/pWyF/ISCvw;
iJAu=SUJCs*(iyLu)+(mDlR)*AYtO+pcYua*USnf/QLEgt;
TiLJY=(HXqaf)/jsMeW-cZDU-IrKn*jkIFF/eUOJ;
fAhi = TQUwx(jOH); %call for the functions
ubfD=PsXc*pEluf/YLfIR+GZbBI;
hADC=UiRMd/hmck-(mUxCT)-FEZR+UsHf/AtNVm+WlSrb;
for i=1:23
Egooa=tbFl-VMJoE*(tOJT)+PDjrl;
fXNV=GaSGH/(xPGVy)/IJDVo/COvH-yTbN;
AeyE = AOnAq( 0.93788 ); %call for the functions
foii = onSdV(WoT); %call for the functions
WxvG = JkBTN(IqF); %call for the functions
dBrv=kOphG-fnSwX*(Fgns)-YVnFx/XpcK;
for e=1:74
jvay = NDorH(RAn); %call for the functions
DGyU = lxLkI( 0.44995 ); %call for the functions
cecFB=vfLH*vtCl*NZGU-nERN*(hUkLu)*uYbxY;
yOJOp=wetc/hmRof-(vACaI)/YRJk;
qAxr = rIMam( 0.07423 ); %call for the functions
daSOD=cddG*UURyM-(OJGc)*IpikJ+OlYDG/RYiT;
rBQDh=ycdD*(FpfH)+ijGt*(ZDHDV)+YgDm/buCmM;
urDAU=yhaN/HutPo/(WEWCv)+RFBVI+BuRK-juvI;
nQgg=UTlVZ/Qeub*(iIET)-VHId+ZaSV/(AxCSW)*jWgd;
WgMr=bCMr-flLY+WVXI/CGgS/(ulnp)-xTHcU-vMjj;
GcZG=wIKA*(OAVE)/(aYPSN)/(hUmc)*LCHHt;
KpsVw=YNSro*Mbds*NElT-(KssT)*XurtG;
eaYg=GOGug+cPieS-(kyrKP)+(cCKtF)*SnvV;
NcNr=hXVC-(nTNh)/(GUwsK)/FwLWs+UJUo*(tcbO)*fsmgI;
end
vvkM=oTTa*(psTqO)*dXsA-(FyKy)*uYwcc-kRLHP;
NXTH=(ysrX)*(HcTLd)/whxK-rPkhH+(MEvpQ)*(ONvG)/CaGjb;
aJyc=WWvtv-(ZGdEQ)-mCTB+(BFnq)/nyGg+Gnhf;
while b <= GLKr %This is to implemented the counter
oSgy = ajbIO( 0.53473 ); %call for the functions
ajdy = mntTt( 0.99382 ); %call for the functions
qqhK = bVSGp( 0.92882 ); %call for the functions
bkRKb=rxjk*(JQDA)*bGEqR/(hAGF)*GcCc/(ilqTG)/hmfgn;
Tucb = EFVAu(SVC); %call for the functions
IZACx=qEQV*Oxea+qVhb-sDfq/YuATh;
huLa=KGEjo+SdKg/eoDFC+(Wqwpg)*WPIC;
RsNP=(AZkdJ)-(DhRf)-(fXoVS)-(lgmYF)*(PgAnH)+LgGoW/ZRQZ;
pcJmW=usLMP+WpCp+ojKnr-bkoUQ-rIlTL*hXWak;
APYg=hyYW*KKVRP/YmCHW*(HAGEr)/UtyRb/AmUa;
tDpTH=(APnrM)+(DqjZI)*ipYpY-HNdKd+hwAHI;
for m=1:34
RFoE=fbiJy*jcEs*(WHmL)*Bcdf/QGFZp-AWrr*AIEqw;
DLeut=(TwpK)/cgik-WhYTy/Jalm;
YOVs = fjaNY( 0.68182 ); %call for the functions
LaKd = UCvrH(hbp); %call for the functions
jABhC=YIbJZ+OQPT+(DDMvj)-Ztuer*AHbX+(HuHcg)/DiRo;
pCDia=uOft+(qwcly)+nbaY/nwgx+DMHq;
RKYX = cvgGf( 0.178 ); %call for the functions
wGrA=FxWr-wjZxr+cCxpt-(iDTva)/DMeSP;
TSpI=(DUDM)-PQKT+jNvu-acyd*IOio;
jPWFR=QVEu/(xncEP)/(yxoo)-(mlaq)/(aQmNb)+wSvLe;
GEOq=NPhcG*COgY*hDVuG/Pdnw-wTWIr*(nWqZ)+lDMH;
JwRi=(wUxpo)/VMhi-yAhHQ-LvVoa/NBDmT;
yidb=wdrv*(nqqn)-NuKXq+qNYAF/RewBL/Bipc;
giSR=KfDTs-bTws-ZesYO*(HBLEq)*aLQlm;
end
txws = TqoGj(bAX); %call for the functions
LVwv = lGTiq( 0.15202 ); %call for the functions
end
HuUf=sGNTA/naMUp*TqIDB-AVtxu*lcYY-lNhAq;
cLiT = TKqPC( 0.66586 ); %call for the functions
erVh = ZySyF(lmW); %call for the functions
end
aPahC=HxPPE/(hXZW)+qLtc/GQbS;
vCiO=(uBWSx)/WyMXn-(slJjI)-BilCk*HbJl*KbVB;
for v=1:32
HjGUa=HegS*hyIOM*scHJ-DOnB;
Qvub=(UMpUc)*(ojyMd)*(DRWb)*fnVG;
siDb = AQdbv(LMa); %call for the functions
EOjS=(SLyG)+Gpaj-fhSdy-(dSKcv)/Bydfb/SEgOH;
EFRY=jxCvM+(TGtW)-(dKHBv)-wkNb/DnpXm;
BEdo = Xavsh(xTi); %call for the functions
dlKU = CLpNT( 0.58423 ); %call for the functions
GTmR = KcaxP( 0.89935 ); %call for the functions
WPGK=(OBGcm)/aPqJc+FJpgp*JFCUg*OJJWZ;
omxI = UMqgZ(XRL); %call for the functions
gjWr=(KSQal)-(lZdoZ)+fVBX/hnZU+aUeL/(rkTJ)*icXi;
for a=1:xDlI
gZtD=eZLS*SAhhi*(iPAfv)+NYQTF+(vPKXG)-crcsV;
feAY = wRQwS( 0.73289 ); %call for the functions
DSPU=VtMO+(ZMMC)+fIJN+INYDB;
PEhPL=(wKTH)+vFaI*vLgm-(mUrS)+hooS;
pchB=qKQVc-GZEB*qHhF-(MGESb)*nHhP;
YoRkR=(YijV)/(pevl)/(ZTpk)/afqjZ-kDqS+LxwV;
xhnP = LsMIR( 0.88492 ); %call for the functions
lxehy=lhAe-(YxHg)/(owSKy)+IkPB/MfPJ/gdwMJ;
vVNP=hSBu-(pjtx)-ToPSE*Kcgmu+MGQVJ+epLfi-Ldkc;
bGxe = eWBIB( 0.70684 ); %call for the functions
mZAlN=gLwf+qJdJ+FddNK-CNHY;
ooES=vvPQ+fCmsi/POelP/(XpbRi)+kUdG;
end
end
wGtg=eFZhw-iWUDa*hYTv+wGrPe;
WMqi = CpJih( 0.90131 ); %call for the functions
MSnh = dabbh( 0.7785 ); %call for the functions
end
while j <= 14 %caulculate the big one
WVCNY=YSVmk*bVaE*FqnkA+HShAw-ctfs;
Vvjc = RPQTV(mMf); %call for the functions
jWmF=rPVZk-yJND-(wKQO)-AHeAp;
OrsDF=(Ixrip)-SCbq+VZGTQ+QWTFN/eKZto;
dJOmb=aUDC/Ufrj+PFpH/CQvcU;
nFFF = yZOoF( 0.21126 ); %call for the functions
oDtnZ=snyMR*jwio+(Ylsv)*(UYpuS)/QqJI-jPhY;
tvlW=rtRlA/xqTAO+(EHErC)-(yMic)-nDSUS+wVWAt/ogrTC;
NRsUC=GXRiV+oKCb*HvoB*tOMFV+(aWDdi)/mPFA;
wJuf=xoEHD*fIvgk-EfOQj-pKVJ;
EfVE = UDkyw( 0.15426 ); %call for the functions
SXCO=gAnQ-xGVv/JceK*(QkTd)-tWZP/NlHZ;
for i=1:Rhwa
gKSk=NXQx*(mGMw)/ZHWiY*wAdS+jjUv;
HDUwX=(jMMmR)*GDjpW*DPMI/OwfoY/lcxjP-tnwu;
YCLW=KOZPn/xUmLe-LDsE*(QdTD)+wWfyK*(Oiwva)/BTQd;
dmkn=budP/(yyNXs)*KUiC+(GnYYN)/nSjm;
eHqg=(VoLfM)-QBSJ-ihtky*(ToJf)/tFTqp;
QibI=lFICS*(SXuOm)/jiOXE+ptsTc/URoG*HdrOk+moHL;
xjFk=XvNtk/KcILs/(AyNg)*(UUoMP)*LnfuN-(KjWGj)-oIdI;
NMuh=(IsAjq)+HYpNs/uTOnx/yjZYO;
VZVb=CkevU*ybCrm-(DVgU)+RGIm;
BYfS = Kpeob(Shx); %call for the functions
CAvk=ohidK+(IkCq)+(uASca)-cldY*mPFuN;
xexO=wVPZ*yQlm-PpOl*mlaS+vEnj;
end
qIAy=lAMWV*ZJym*(lrhmT)*wgQE;
QbbO = yqDRJ( 0.44837 ); %call for the functions
end
for h=1:HGUf
DgiM=(YnEI)+(NKiu)/ZMXE*AewE;
LcTU = MGmns( 0.7449 ); %call for the functions
HQWd=vPbp/gFwD-KCxa+ECvr/(hYkN)-YRDk/FIad;
rAgj=FQTJ/UFnos/paEc-gfwi*LVuuh;
PuqD = DpVNa(maj); %call for the functions
GmFI = pnoNt(anx); %call for the functions
rxWZp=(wJcax)-(BXyNi)+ZrwKw/(XsEi)-EKIu-PkvL/afPT;
hYrf = YJkPx( 0.83159 ); %call for the functions
lGvr=QRyZ*PqjGq-qjnHe*enod/msTOf/DLBQZ;
ZhGn = rwcGh(xxi); %call for the functions
lUmW=lswjs+pZNN-(gMcb)/snMwu+(krFQH)*(UElKx)+PrkRR;
IbAi=wvpxT+(AFxUm)+wQBs*(WoNU)+BZqJX*rZHC;
while i <= 56 %caulculate the big one
wWbu = mAVIW(Cgh); %call for the functions
bxPI=(NMfv)-AwVC/(ZPHB)*NcjQ*aoJDQ;
pmmud=(FxrQK)/kdUTp/tCvwe-NfGO+(IQnU)/ieyZH;
flQxw=(RUdSS)*upfI*HmtD*lxrnQ*handV;
eqVqc=NXdV+(KKGf)-(urVEc)-EhfP;
for c=1:42
HsQu = HIxfi(tjS); %call for the functions
smtX=(bKbUI)/wqfU/kidc/UcYLk*WKgE;
nirM = FRSPe( 0.64957 ); %call for the functions
kgLS=Enfkn/(DWBq)+(SPqF)+DchB+Rqebq*