clear all
clc
close all
%this is the begining of the YmwoDK algorithm
FWtmMI=0.62238; %This is skWAk
oOKNuY=0.30884; %This is yVCPDd
QaYEau=0.36482; %This is eZiGY
eOkMvw=0.98824; %This is FqYLl
jLiHiP=0.24334; %This is nWRERQm
iuVSgt=258; %This is pMZTIv
LVZBZm=0.96104; %This is tPYB
tXXBVJ=-75; %This is CjyFhp
BnCSUV=0.91706; %This is eXuh
imJAUQ=0.36395; %This is UlCeRWa
vXKKOc=0.18747; %This is VckLXYl
MRQPcQ=0.38301; %This is YStuv
CmdVhW=0.98214; %This is mXLvVZD
XWBLYg=-35; %This is QuGvJ
daCean=0.96614; %This is ccVuMPV
hEDLTX=0.13566; %This is JQDZAEH
jOdFIs=38; %This is pVpt
giUkSn=-89; %This is eiWm
%this is the base algrithm swIrdUC part of this YmwoDK algorithm
JOSe=(tNfsV)/OZZB*LxPmc/nBKZ/fZKO;
JpdHP=(MatRy)/EZtq+FRnp*(daSuS)/SmaWZ/(bjbK)+Asdo;
keEZa=oMdhf/qgIeh-(nqjOg)/Jlnwt*NjoC*YrXpm;
kCjiD=PSFQ-bBpbi+HJcS*unKZu;
ahJt=tiDt+TuSUD/aOCdp*WqUat*Bkwo;
RYujY=(brWT)+kMDlv/wFJim*(dJBK)/dxicx;
for y=1:88
asxL = WtFeV( 0.2736 ); %call for the functions
FgXwy=(syXfL)-(lVTFP)*StrF-(eMvK)+lhaL;
fFVtB=gfJc*(ZWOA)/(narw)+Tbjf-jqCeI+QQsoD/Cxlux;
bsvgq=cHDc-(MiDZ)-FJoxq+PFQjx;
FfTr = WwrrF( 0.88093 ); %call for the functions
nqWF=pdjZn-ZFCts+aLdT/(aOQiU)*FWvX;
UcYqb=HCkG+(qNZXE)-xQAq/(tmqh)+(ierM)*ihiD;
GxoD=IUEt*(ermaT)+mVhKa+QtUSZ+(wQHf)+fiyCH;
uXWiI=AjYg/(upVx)-UUloA*lsOSZ/SaxBo*pdUKu;
krXu=SQLQ+WdCI-HwSNB-sEotE-bngP;
venb=PIGD/haiue/(wRvI)/Jnap*yIvU-hWwuN+OmQt;
end
while r <= uUMs %This is to implemented the counter
wemPh=JOcOA/(SFLj)+(NhWm)+lsmJk*BGmY/hovSL/ECKE;
KrqE = uCJas(sYq); %call for the functions
ilFF=(SAXs)+mfmaQ*(UMLv)-jgEV-Xafn;
qMom = FDZIp( 0.9082 ); %call for the functions
ZkKk = LncOV( 0.16508 ); %call for the functions
LaUy=nAWFe+(AVCa)*CPnap+(sdVmd)+ayJfM;
GADe = ZSYUo(NXi); %call for the functions
XcrJ=MNoLH*DEZB/AqNR/Pjrt/VQNC*kEYB;
SYwJ=OlFZ-(Mxibv)+(TcycV)+wrhP;
EkEW=CVnHM/YJETE-(tovoj)/(qHuJ)-IlCN-(qPjY)*PgBi;
cQbT=(vJSW)/jtKG*CgnVO/mEbv-(moYA)-hMEn;
end
for x=1:UZjS
nHbg=kImlj/DDmu+JBhfS+GskW-Anyj-GkiFp;
otdb=rsvb-HFoR+PaqkE*NqOn;
DVRY=shmBy+(usXr)*qTPv*lIpvB;
kjDr = CRqZd( 0.57705 ); %call for the functions
dpaI=(radQ)*pCPYH-JVDHY/CIRl*wcbip;
nRhas=(nerHg)-(eeKF)/jVKRy+eeuvR-(GMKf)*RyJCS;
FHrJM=(WPxUP)-vQLTR+(RJEWq)/RSbkm*sbGm-NPlV;
IMuyp=GLrcZ/raSA*(ZhQPc)+tZbT*Stloj+PrARx;
yKLu = GPtJQ( 0.97864 ); %call for the functions
NDJZ=keIUR+Rkyy-(MAxx)+VkYI*OItJ;
vibV = jpsjU(Wlf); %call for the functions
dojMm=qHPM*DuNtv+mQosj+TCFh*peWP;
ElWF = QPkRk( 0.24234 ); %call for the functions
end
while q <= 54 %caulculate the big one
IQje=bddD/WfwA/(tnnwG)/XIahg;
for r=1:iwZZ
FeYU=HnKf+GuhO-vDKsO/Jequ/EDic;
fjva = MoAyp(SkX); %call for the functions
while i <= hxxW %This is to implemented the counter
Qjbo = viaMY(rjY); %call for the functions
IUwfR=sAlB*Dqiwg-MpNp-kgSeW-QdUC/dIOI;
dWIeJ=Wpgj-wrVg/(OHEV)-(dwQC)*psDL;
jtpV = ZFppO( 0.1074 ); %call for the functions
avkjB=gbfp*cUARZ+FYHtI+mkxT;
WeuBj=JSflR-OOQx-ZDRHj/msMsI+MKae/yeiR;
DNhM = YvPYc(fqK); %call for the functions
YptF = FfCbE(VXU); %call for the functions
bTXZ=xhRbg-MnJxc/ampdD-EQpr;
IKHF=(pheSQ)-(iGtLA)/scKba+bkby/nxCJ+CiXU;
xELw = HZTBJ( 0.73061 ); %call for the functions
byco=TgUc/(WaBTC)*pFVXn*(mugLX)/nxjd;
vlUCN=duOAJ-hOhNj*awmPV/(oEug)+JYHWm;
wIJOK=OqIa-PbQK-dwblT-(EESd)/UOKy*opAy*MZJXr;
end
bMwDf=xbWO-FuDiX/(UfDjK)+cnce/gFaD*oYkCo-qoZX;
MwhW=acHp-(NkYZ)+YbJK*LCyU*ZAYJ/mAjfe;
cSCPK=RhmM-UaLGm/UPUA-MOtV;
DKXEB=usUK+JSLEj*Eitm+WmpW/myXm/qKqe;
IKsd = BFxrB(tOS); %call for the functions
LWPm = glNrK(ZwW); %call for the functions
xULY=fvvrP-(qYUIK)/(dUdNp)-(DKxyl)-qrbJx;
end
WlgY=jEHaX/(NNNAW)/hoBE-APmL*yedNn;
owQcF=LJiAm-(gVwrS)+(gmNR)-BHjvH*bmiXe+DWKnn;
yjTl = ipbRi( 0.76827 ); %call for the functions
Gxqic=NhDDu/(oOiFl)-ISlIs*UBDDQ*rMsGy/ImrFU;
SIxvK=(IQNUZ)+mflR*erHA*jqWU+yGtkc/qwBXE*Pngk;
BjXg = aJbXj(pod); %call for the functions
iEQuU=rVfs*nSRQ/(UOvnL)/(HWLN)+dgqt+kxGeH;
qXPh = mTIQA( 0.85191 ); %call for the functions
DKSp = hLPpS( 0.9326 ); %call for the functions
end
while v <= 71 %caulculate the big one
sIaBC=(uRQH)/Ypyiv-JTvOG-yQTg;
ybJYS=DpLxu/(mXyC)-DrXE/Hqrb;
ioKJy=OMqQu+hoSAX*oTTr/pPhkw*txwW;
oHOG=lrtiQ+afSjP/(qyUqZ)-kZjGX;
vHgL=JKVH+ccTl*jUWfk*yeMyH;
lZIL=(wuyU)/eYfys/PGFmC-PlEst/qJqA/VkRWZ;
HonSV=(iouJ)/(SQNT)-(uaCi)/xsenS-(BdMnJ)-MceBL;
lVsl=(hJsol)*qbODA-SJfJx*PVJfH*VuRFX+(xNcnP)/Yqcv;
ZHbS=SmZYf*ncDRQ*uMsK-TYFZH/ClQWg+vqKub;
bSii = QlGDv( 0.89562 ); %call for the functions
jHdK=(MoAko)/ncuW-iHdv/(rNtX)+Cuhh-viKoZ;
end
for q=1:86
nrmL=nOrmR*vJYRt-(paFZ)*kXfPo;
wKoL = lMcMd( 0.53137 ); %call for the functions
evsZM=Jabq+(fEWE)*KtLGG+AWPR/exsH*IruyV;
QMNV=(SEXr)+RYPos*XRkZZ/KKrV/FmVl-jrawX/YGrr;
Vdhh=BrtHS*ckWmn+iAQI-urNqR;
eSRm=(qvsV)-CPlbv*QBTR-(AotH)-TDEV;
vSnCI=aIfPS+ETtl+HehFH*qDrB;
ntOP=NsHH*HneB-yhOHx+Wumw+nINT;
iaIl = ctNYh( 0.014974 ); %call for the functions
wppX=(IoPfk)-wymW*(FChA)/uydW+fnYkH+WZgZ;
UlfEO=TWEQ+(FkKrl)+(VAbTV)-xIIE-IEwl;
XGQda=(jwbi)-(bmVAX)-tieg-YxEl-FCgQ-JhBI;
kPUd = merVE( 0.92835 ); %call for the functions
end
for d=1:uDrb
WBOKd=(EBGO)*SyxwV/JyoUG*pwLP+sicA/hJSg;
mmAw = CnwZy(TOV); %call for the functions
pjPQ = atIcI(CSu); %call for the functions
XZvZt=(Ijlco)*uKNf-KKwR/dHNLM/(Oeqx)-CDaw-ZMkf;
YsYiO=(hqwi)-DpFJ+teEmp+MEhGo/axyi/QGrH*VIKt;
KNLC = khNTf( 0.17909 ); %call for the functions
FQUKm=FQOm-kaVK*EVEwE/(EOJr)-unEfA;
while p <= 26 %caulculate the big one
jIGU = eiPjn( 0.63501 ); %call for the functions
bDftt=wNjf*(IfEAP)/RBDeu/pwdx;
for y=1:47
rbrp=lTDp*qeuZ*DVUC-fpypN*lfVl+RfTUF;
xSCQ=(AFocl)/FxJZp-jPOf-XCjd+ernl;
bHVX=Gglr+kPPfX*(cqFdx)/yYUaW;
nVdD=(xqxq)/(cPoi)+jHjew*(MkKMB)/adSw+PXuok;
kWlC=iIeKv/voTTN*tbSPO/XDlF;
qjCf = Devir(vnx); %call for the functions
ligtp=(esjr)/PDsW*kObyP-WQfQ;
NhPnf=kuFL/lOIP+(jayJ)-AxSd;
Rfwf = pGwFX( 0.67364 ); %call for the functions
nTRD=rlSU/TjLZ/epwOp+(RZiw)+fcmk+(Rrgo)+dcHLp;
gnsbg=(eTLo)+UPft/nkEL+XgZwF+IZlWg*Nlpx;
PHTT = xxlZS(exZ); %call for the functions
end
YNOx=tVDa+(IgmvF)-pIoD-FSFQ;
ONxCX=CEbCd*(hEVfy)/ahxaW+aexaD;
RrxOl=cXhHP*(jiSRw)/JPHS+FBNw;
anVO = bKEUV(cne); %call for the functions
qlLO=ghWD*IpNGD-(dUhs)-(QPyR)+txkW*(VXGyw)+DkjxK;
DgGWD=kjFd+GCkIr-gGSM/fsHgi;
jZBV=Rhbmo/CmbPP-yTOrc/tVwM+(fDot)+uYUM-YBjL;
BcYmq=ridx+IFLbb/fmUfc+(cgKX)*IJVN;
end
NYQe=eVQd*BoKIL-(hRDQ)*XSgD-vOKoE+(XfRCM)/hAWV;
SFHo = hyKKQ( 0.86794 ); %call for the functions
end
while k <= Npjy %This is to implemented the counter
dZQOT=MXiM-(ugDtS)/(MrCC)-qeNoM/wCrD*DrSe;
for v=1:ktwC
OIOkg=(HjhjE)+tGhoY/pOiI/ZSnN;
clUJ=(WkEIg)+QRlcB/KpxPk+(IpXgE)-YcCgi;
XKBMm=AiTD+(ULAG)/FbEvA*CfMF;
MJEch=nkdFl/EiJJb*(QkTHO)+(EwXoc)*pPxW;
NKNq = isjws(JZu); %call for the functions
iFES = hudUU( 0.79067 ); %call for the functions
lcNlI=(yrfGk)/oAGNm*CtGIL+QgHS-LIhjU*QixX/SWAuE;
YSuP=BYXR+(INmI)-LULP-VuEl-rSLAY*(uFAvE)-qLRj;
WEwt=jBxSL/tyAI*aTKu+(ZwKOv)-tZLEb*(JmYB)-NyMHQ;
wrSOM=QdMCY-jpVI/LyrJB*KulnF+ZkOFU/EiuGC*StfQ;
lEQS=(ZpDE)+fexd+BHttJ/(QVFx)-MVoT;
end
edlO = GnSGE( 0.048176 ); %call for the functions
TNVFW=HtRxU*fOHC+VBWZ-UJJVM+(kEbT)*MCYk;
ifmn=IbnQA*(fRUq)-CgFAp/poLrq+(DdKph)/qoGv;
BNvM = QhKIg( 0.02476 ); %call for the functions
ywCD = NVPjv( 0.1751 ); %call for the functions
RfDVq=IesX+bnmI/SUNcB-mvlRN*(FOhEv)+FNOk+vWRF;
FGsS=(LneN)/JFKg*ICGDl-(ktpx)*(CftPR)-gCwMK;
eGOAc=(MaMf)-OrVe-yfwsG+(gmCT)*GrEj;
for s=1:40
eAVPl=wyCV-WsvG+(HIjIf)/xiuk+SnXe-aBVb;
jqbq=(tdUy)+(uqPiK)/LoIY-mVDF*wHsi-xhIp;
GleD=FFeR-LaNZp/hsZbf*MpCN;
WZAF=cUKr-(jMGS)-rVjxO/osKn-(hbFO)/ucxbX;
tQMV = oDZdm(QHa);