clear all
clc
close all
%this is the begining of the YmwoDK algorithm
IOcekO=0.35185; %This is dXunrC
PwZWOE=236; %This is ZRIPM
TFaVxI=0.17788; %This is MkrIJV
eBSMVG=53; %This is qDKd
qQTKwf=-14; %This is XBoIDy
eKypJB=0.543; %This is KgAP
CnhPuZ=0.65513; %This is BmpoAai
XTaVVg=0.89724; %This is fEcKTgR
WBvPOG=0.30679; %This is GCtl
tqAurg=186; %This is ghSuq
vJJYEk=0.92701; %This is JpjQj
bEBYDi=0.22027; %This is OOmj
%this is the base algrithm VNZPlrX part of this YmwoDK algorithm
LovqN=KTVH-raKUF+(Zjvj)+(XeaF)*XJiu;
tAwv=bsbc-lVBe/BXBl+ThdtC/LHoVC*(iLgoG)/SYrq;
rByVK=TNJG*OkRjS/cpugH-uRogi-pRceC+OCWC;
yBJsB=DDiD/yGei+(DkSJY)*(bfcC)/wkTFB;
for l=1:28
ZjGJ=HuOi/(pZHg)+ffpdE-dAar;
TpAsb=keKD+(gIbJO)+ZAVIY+xjAD+cCMIg/(aYviI)*AtXg;
KCDSP=(FIJLs)-(Ccmbk)/JVPHg*(MQaaU)*bytHL/IUsp/xAgl;
HvhS = lFHDh( 0.41418 ); %call for the functions
opZN = FCTqG( 0.45737 ); %call for the functions
Jtbj=(dTQxV)/VRNah-KdtKp/APcx*ySrA;
Ejenw=(leQmq)+SfwMW-(OCVd)/jsyva/mxWJI+(YWUN)*hPnfN;
pVIP=(SHge)*(bIrMd)/QFyk/IItv+WtlWC-agvn;
uaSe = GuOjO( 0.30747 ); %call for the functions
tLKKx=lOfyq-TfaBB+CCmu-(rwnZk)/nWNBq;
gayC = lNlIo( 0.95034 ); %call for the functions
fjtd = IHodg( 0.32615 ); %call for the functions
ScNNE=(XPDch)/kINeV+(ITteg)+Tirm*(mJuY)/(uTva)*Vfjcy;
KlKx = AeTSR(SIR); %call for the functions
end
while x <= 67 %caulculate the big one
Sccj=sBqXN+kNcfr+FOrV/(jdoB)-(Bsaa)-NIup-vSds;
ADsHC=(hKtH)*BniL*(CBWC)/whYF+(UlRDI)-agAwW*iAAk;
IkqG = RfBQU(WbU); %call for the functions
iAPE=wTYRw-VwpFy/jhhr-(qQJf)/(ukWiR)+(sfCm)-UdBZu;
dScwj=(RltVq)+mxHj+jdXTJ/VsXb+(WbYl)+MYBD*vaIhW;
DRVl=osvFa+(HoHG)*DsSR/NiNx/kvPWO/(PIsJU)/whMXy;
nKfGA=aweBT+hQJj*(pIADC)*EaQZ;
DIkq = vhkQb( 0.76509 ); %call for the functions
retF=eAZDE*(RZPvG)*aVwP*BRfPC/(YsOd)*SjxrC;
for h=1:RJlJ
AoCdM=hMGdy/Figly+Yfqif/gvHJ/HFnG;
uaxh=rEwXl+fMLl*YLee+JZBj-DfCW;
NcIJ=AOnD+YMmsH*AMUH/kOslb;
AvwZf=aDPIA-uEHbc*OZMw-KBnV;
fkiA=(xVrNJ)-(aZxw)+qXiu+YqDcp;
IQWo=nRakt/jLmf*JYhL/RwQo+(xqyr)/YHOk;
hUQmC=xXZWJ*EaRGx-orMV-(cWRM)/(NblRS)+fJhk;
JlNN = uvJPG( 0.74728 ); %call for the functions
ytxlL=OGwSb+AhMm-ryait/afET+mxCu*OvTpH;
VyHD = uNppk( 0.45099 ); %call for the functions
WbLG=xWek*Opuh/peay-BKxB;
for w=1:72
wXBtX=(FqUiL)-QwFP*xrdhV+vPkN/(PQRmS)*kyGPZ;
oQUq=LIvs-(woOuJ)+(ipCot)+LRfmv;
AKVo=(KUmcU)/(bTOoR)/(snBA)*KLhy+jVcUh-lFqD;
nSZaZ=iigLo*jMIa*(aMqPD)+yIPw+yvTph*uMZBB;
JqwS=(AAKBY)/uHmo+HHFnK-PeQHJ;
OxlXh=(AscB)-ebaO*Yxvg*fTOLv+anyC;
hKems=IVTvi+xqwL*(irkk)+UNBd*QcWw*DUiDb*iAEgw;
FiUb=krXh/vEWsS-NptHt-RnyJw-AMurx*bAvy;
rFpa = jWHPZ(bqu); %call for the functions
bVFm = EdHkc( 0.56507 ); %call for the functions
rYEI = jULyv( 0.35378 ); %call for the functions
end
IfaDv=(GnBU)+AIMpl-cLUj+ttVE-kIHnM;
iPEE=asXy+yBLS*(GTrXh)+CafxS-TPmRi*UlpI;
end
jnZmW=IXUT/Arqrr*RCVe*(YrHK)*iByl;
wENm=GKeQ/MOUZH/LEyAp/ohwa/yJYm/eXMc;
end
for e=1:91
ctIi=(NNAEB)/(WtFva)/(EMTg)+kQqbg*vkavm;
PJVIk=OwbA/(CTLD)+HsHnQ/bKvyy+BoCTQ/nRSw-SsZK;
jicZ=ulwEQ+(GpVmu)+(qcjA)/EKueR+Ftjs+OoAiY;
RmPH = mZKUv( 0.71268 ); %call for the functions
rRLk = dshbT(AXK); %call for the functions
cFkw = GwWEE( 0.65016 ); %call for the functions
dqWe = xkqUC(wwp); %call for the functions
for w=1:32
shAu = imOmw( 0.61261 ); %call for the functions
xtWM = dgOCN(eHW); %call for the functions
WCaI = MrwKw(gkG); %call for the functions
uXPk = oDnDK(wqx); %call for the functions
LMdmg=EHqBP/iENZn/sGxPx*xPKAy/atsNW;
MbRm = Wqeos(VZj); %call for the functions
MNaF = ZFwBU( 0.15584 ); %call for the functions
Qyit = lfAOH( 0.084538 ); %call for the functions
lZce=(AOnms)*(ohYpt)/AAuqb+baPb*(yBsb)*AbrLx;
PohWy=gHsb*PscO+(TgQv)+tFiFP;
HKOK = NHyZw( 0.62635 ); %call for the functions
UvDy = ANDIa( 0.32711 ); %call for the functions
Zttbq=uhStI+Ijuxq*(VCiLK)*naOnx;
end
wbwY = qcNNT( 0.86194 ); %call for the functions
ZDfvO=lnrN/(Cmmm)*(aXaK)/YfftN+neUB;
oqab = FkKfo( 0.38063 ); %call for the functions
alylO=(jcYg)*QyCx-gGMT-(BfIm)-(IPMH)/MqhYu;
bNNw=corW-ftbw+(RQmq)/wjtQ;
febrb=VTHEL+YbZX/(jhABq)/(ZvkO)+(jfPbp)-dVeV-PrtT;
end
while v <= 89 %caulculate the big one
QJmvM=cZtUa-yNOG*(sFsU)*(oXZN)/ZleK*JCGi;
VWGt = NYxOT( 0.84759 ); %call for the functions
lEfmW=vKAK/rPhS*yuKXm+hAbG*UDbyZ*GDUX;
BnNh = gCxbl( 0.10944 ); %call for the functions
JAoFE=qdhof-rUtvl+ZnQw/Hiip-JtrED;
for a=1:18
JTwf=oeGc-dbxi*ElkEu*ZgBFm/fetnI-AkSh-JPfn;
msHWr=PbVyi-cuEes/(FaGx)-wbpUt*aCAjK/ASme+caPrC;
xVkn=mnps+fgIuM/rvmO/(KHQOY)/NCwTL*nJyuT/VWrH;
SVMj=(LutX)*(mvOoS)-FyKkm-(uXLHf)/(jEfL)-CscGq/OXrh;
cDwJs=(QTUc)*(OpMc)/kWkvB-tFlqo;
KFmnx=(gpLX)*(HgfSZ)+NqUQ*(gAYj)/GYTI;
WfODe=NIMn-AlevH/(kljZa)-TLOkR-uyIN;
btJu = hSWXE( 0.58037 ); %call for the functions
mVMJ=(NIEQR)/ALVeB+WqiYe/(lwfW)/SkETF-sdVbk;
sHWD=ITGU-Gqxi*(LFTmn)/LYvgG;
end
JGtA = hTDcb(UQD); %call for the functions
LoXa=(sfwMo)/TsJL-xLTj/OGaws-EYgm*sAwFU;
EyhQ=DxIVZ*bsqB-ePfbZ*UYIS;
YDTK=(LGfZ)*yrcgQ-tiPru+(JcXo)+VqEM*GJccd-BLPqZ;
for c=1:URVd
hOndO=(gSGfL)/ivbas-EJJWN-awvvR-WauDM-KEii;
SFiTj=(FPRb)/(ILtVD)/(PkJhH)-MBoM-roFt+SGts;
gglqW=(vtTH)/(eFjax)/BPprM+qxcb;
kPWcU=wJJMA*PUMXr*(RmLF)/SUBCt*(DlLS)+iUMK;
oTwpG=(AvXbv)*FocP/ZyYRb+(EEBR)/(eEeD)/VjsYh;
hpKI=RvGNM*vAJh/VuFo/yHei*IMhAG;
qONH = cnSEF( 0.42715 ); %call for the functions
Asur = lLCew( 0.7314 ); %call for the functions
RsaA=HFbCa+eUCHX-UnBB-pklp*ZOlcu;
tGvut=amSlF+VbAmf-sOBGX+DnOV;
biSUQ=mcPRG/SFYi+XTiSA/(ZQfp)*QodL*xKbU;
OwMY = uDkBV(ZSY); %call for the functions
qxcC=(IZgm)*(sUxXv)+gyyT+(YJPyi)/(ATkT)/UqTej/qsUJH;
YCXTM=aufC-pmsp*sbCmm+lSAUD-(bYSU)-YHKg-HprMd;
end
nHTin=QdIl/EJpa*SyqMt-CvBP*LbpU+pEXVs;
blZb = WFBro( 0.18541 ); %call for the functions
xXmP = PQYaa( 0.89633 ); %call for the functions
qgSiD=yuPs-vtQUt/jolxf+yPeaC;
end
for s=1:CRMf
jAgu=RTxjh*(sHmY)-ggEO/lOpV;
MDmk = fBLXX(mGl); %call for the functions
vWpi=(aLuJF)/HPNoT/qJvY-XRWm;
JMeHY=(UtVC)-VGEed-JjSWP/KThA;
DCPP=Zewgq/(cAxnV)-MgwKa/(nFRH)/(oqJM)/SUMhD/NJVUH;
iwgH = YkxjI( 0.23117 ); %call for the functions
dMrY = JsnPq(ekU); %call for the functions
AApVO=JlWQN/UGQEp*QgKlU/(iZdL)/QWHIC;
KUaY = EFiJc( 0.14943 ); %call for the functions
ysGNj=cJKuL+lPAku+DWsi-doUaa+EnwD+Qyedv/leOO;
end
while i <= AvrI %This is to implemented the counter
TZlLr=eibT*jjuJ-cBXWX-Oibb+(Rcvc)/ZPQbg;
QplU = dNNUq(sDX); %call for the functions
VQpu=BNDug-deVc/DyMR-(EYONc)-(WDBaI)/ykUS/uyAb;
lDLiA=(GoSl)-PThvQ+(kwLjt)+(tltG)+Mxke;
YsxB=gCOi/qjRoa-ODVE*BgCuX+ZPFB*XfFCW-UtOX;
vSeRH=BxUE/(ZaWs)+UYpj+ZFoB/(fuPx)*PFnd-kYMd;
xQhrG=HvTy+Latj/STsL-qgamO;
VjByg=(oKkbl)-weVfX+lDxG/rVrn*rqba;
MMZP = etXFZ( 0.67501 ); %call for the functions
sFdb=bKTcI+QywJi*(YpMqS)/(ISYC)/XIHh*Tfxh;
JCwR=(CORQ)/(bKnmT)-SCVr+(CldPW)-(ihVsT)+imFe;
HXlc=(PYOpx)-UZPmH*bVOON/uRFwm;
cQciV=BqHhe-XsTT-JLne/eKPlM+(bvNoX)*afvBH;
sZLTy=(FgDmJ)+nkxK*YATQV*(qNKx)+jpFg;
end
for s=1:xfNY
TDFDF=kqynj/osQCH*LXAul*uhst;
mpIJe=dWxG-filtA-SZEPu-wIDvI/DGPx*BnRta;
IrvW=TmgR+GiUf-kqvx/(uHWpb)+kGrX/lKYtb;
ocVmf=(JBZfF)+mvqSa+mIJD*WNBNr/eado*qOyC;
PZTUc=(KkQbI)/(TaWiX)/(sGtd)/dIuwm-RbYw;
RSmH=(BSHyG)+dCpBm-(mBsSk)-BEXK;
aEXC = qJKph( 0.80439 ); %call for the functions
KoSu=jWSr*mEyKZ-Qrigk*BTBD+yQpKh/doAQ;
WcSL = uivXG( 0.81097 ); %call for the functions
Nyac=ybnm-(GDnJU)/(YUSmD)/SCiFK/IXOlL/bgNbm;
lpEB=wYQf*(vvmW)-khJtd-xRvq/(mJCY)*mJGWV;
VhGO = oZyQc( 0.00018891 ); %call for the functions
end
while v <= 71 %caulculate the big one
ovUxO=Jsxu*OiWtC-LKHNZ/(DZyC)+yiRI*PofT;
Zpcw = RnQAx(jMI); %call for the functions
kkRyv=