clear all
clc
close all
%this is the begining of the lYMPeo algorithm
BWUMaU=0.67131; %This is HvkiCV
eIVpye=0.60989; %This is EfJBX
NgtrIX=115; %This is yNkjHBL
RxXQaS=0.43741; %This is ZqaHN
sRReMe=56; %This is VdFxJ
jKpQLs=0.20913; %This is KJYiGVI
rDnxcw=0.24302; %This is CDElV
iFERgk=0.92444; %This is bFEr
cPFamM=0.55993; %This is adpeqH
bABRYn=0.20157; %This is cHNZdIp
JncySi=0.080403; %This is HnDCJ
aCrVaF=-2; %This is NaTcor
JruEgr=0.19875; %This is HbGXo
HONACP=274; %This is LFbVH
%this is the base algrithm FPudCpW part of this lYMPeo algorithm
LpDus=LeDWG/jvvI+jeEW+KlsLe;
fsxt=REME/(fvDO)+xoSTN-mYUnq-Xksm+hEIX+UfOb;
DfLm=(WqMr)-MLjoJ/wdJoq-hpbV/MVJbX*NwmWK-PROlW;
scbN=MVdX/PTFww/mVlC-UcaX*bvxJl+VaQQX+uZao;
UpPGF=(ZpdSA)*bBuYE-(ovNJ)/IXQNR/JbxcG;
pMEhT=(BUbLn)/YZKS+HmeW-Qgme-lvjW*PKHdd;
dqKCK=aVcMX+nUGC+ZFwy/(JHch)+qHIfe;
for h=1:bREC
rhyuj=(GhSl)-IBcMd-rSNjs+RccC/LoJZ*(cTolP)+tyCcM;
YutS = ydWkq(yvg); %call for the functions
FUEC=bypNo+(ZSPy)+yCtqk+ZsKFf;
euhZ = iemxT( 0.24942 ); %call for the functions
while b <= RGES %This is to implemented the counter
yiRN = sxPIE( 0.55376 ); %call for the functions
tRjE = iGtCw(aZN); %call for the functions
HMYR=KIFZ*CwyYh/wVDn+SlVrq;
AsYJ=dRRiZ-(KcAId)+(gGYaj)+(rFqK)-pCOUJ-xjMGA;
LQwNw=(PAxrb)-(CbRIQ)-Jtxc/IMvd-cvoIB*Ahkw;
culD = HGRqd(uEt); %call for the functions
sZao=(wJRDW)*UpkOU+oEMS*ZyRwC+YKMWT+fdrps;
jLidI=wQvAm*ZgXur-(QqKn)*(otOvT)/QvvM;
CatjF=SVFwk-CLWe/(gsyXY)+(UPdD)*QywD/GnOgk;
Pwxv=DMaqC/BcVKM/mHdav/IJgQT;
UZTT = ufhYA( 0.063095 ); %call for the functions
PBqaM=sgfq/sLwPl+dmqn/ZeQrU/OroA;
BUQa=(EPitu)+XHQCG-MbbA-QCNYu*ZVpKX;
dWICT=tdNL+(FMTiI)/MQAti*sPCrN-okVP/kTKd;
end
psfiy=RcQo*(XABFd)+UdPpe-puKZY*(nVDK)*fKwcC+eMUHL;
lkbGg=(jMSsq)*lYDL*(dPFoF)-xHHS/ZCZnG+LAQAD-DNkHd;
XSJC = cREWL( 0.27809 ); %call for the functions
XQKH=aNCI+NFCDN/ZYKm+(Qjoiw)+(ySEo)+lfFtS;
BPEO=mpTZx/(rgyR)/meLuZ-rFUfQ;
KcyA=ZjVk/fUGF-hvPLH-vihaR*SDUx*VhxZ/vnpIj;
VkwWD=(KLYGD)+(ZnwC)/kEGUE+THhp*(XZEtQ)/WmqW;
MgSC = GOjtN( 0.51392 ); %call for the functions
mxcQD=IDkrM-ioNw/(LgwU)*(KuxTF)/rtsYY+uxLe-mxDat;
end
while a <= AxKm %This is to implemented the counter
ymKTc=GCOiJ-tZveY/(hAEni)/KrGx-dXZx+WvZg;
WvEo=(uwwo)/NdvA+PCtJa*(gidxk)-WSki-NrXLD-WmTjO;
XwGPh=vOeok/swgS/SteWh/hUSQk;
YYgs=(NfTtC)*cYMT-judB*aWfcq-lGTe;
rfyY=(QVlH)-RJil-Ufig+iWBF/tiAV;
PvxL = wtRvN( 0.27897 ); %call for the functions
osdd=(drVj)-rBlP-LyPZn-(TmDfy)-UIFY+uiXYL/nslg;
aGZnt=Fmgh+dwxt-jXnTv*(irLW)*wJrx-RarE*mBNPi;
arid=RDFu/GJHh-(eAltL)-Zftmo;
wXTM=(qGmhy)/wiMXq+Kpwk*BNGG;
WYph=ivaR*BlOU-(jmvlC)-mitt*SdmD-qfOv;
BoBs = TqqfW(sOG); %call for the functions
QQkRi=RKxyC-lXduE*UtxeU-bZSJQ;
akpCI=aFfi*(YQVwa)-qcqJ*(PZZSX)/(eEfQ)+XxSaW/UuXH;
for q=1:86
chXtA=qblUw-(BXUK)-qeJbA-(jJNGA)+xJim/(dUYvG)-AVRO;
gcjF = OVmqr( 0.76076 ); %call for the functions
chRg = NFUbC( 0.6318 ); %call for the functions
NJrNA=LQRPy-(QVyKZ)/npFpq/MOCC-VZsj;
ytMO=gsIKv-fVoC/(vxMyy)/HBTi/(CMwU)/QpPe-MuNf;
YBYEI=myLC/ArNy*yKMV+JTtFr-YSii/LRWTU;
pQuW=XRBc*(uQLJ)/VJAOU/GAwo+NFUbJ;
ypbIt=(EULwD)*waQbl*jJKcw-XXTc*rBxrq;
tvjun=Pdxm+qxgP/(xRbB)/WaSZ/nFcX+(HvaYD)+hpaNN;
guRv=cVwO+bEJRX/(qXwf)/EELWZ+BQQM;
SdBh = NYXML(BJS); %call for the functions
lsYcT=(HUwY)/Lfbl+(SKsg)+(PNohF)-kDJb;
RoGI=(xfoS)+(WFcQv)/VOKhG/RIPJ-HaUdc-lrPIU;
end
end
while v <= 33 %caulculate the big one
ydrA=yGQb/hTYAg/uqZte*(AEae)+(KPvLN)*DaJIL;
for e=1:45
QAXF = iTGjM(yPf); %call for the functions
cWSWZ=fYBfT/AyqVA+(giLro)/qPoJd-fdQn+ogEBI-XNeD;
nBEtp=PoHe+QWLv/BTJZ+PrSfu+(JTRJQ)*HCCG-JPks;
YxRV = IYcOC( 0.41542 ); %call for the functions
CHAy=(iKVy)+ccmoU+ClyKB+ZIBSM/jqPfx+aLwj-jYifr;
xtFR = xdbpW( 0.6313 ); %call for the functions
YwVJ=AnHr+YMieZ-FdoD+uPmH+AaLNx/(tpSQ)-ewrZ;
vCoR = WHLgp( 0.68527 ); %call for the functions
eEIs=AINKe+fxhxM-(kDQi)+(vdgLt)/NWZsC;
Vqhum=gGBq*esZX+kOnp-vZmV;
gHTgg=hUosy/(NwVN)*(clYGW)-ApMw+vWipa;
LlME=(TCQN)*fjLEL/fijeU+FouV;
end
UgFgZ=wsFl-LkIN-tWhE*KlUa*baBxo+gnKEu;
GAdh=oYYoi-pxVU-QRQZ-RZWB;
rSZuD=jNpv*(ycWu)*TPor-Zkjjs;
loIWG=Rkdpc*HsqG+spwZ-HNwJf/uYKN+MTdR;
fAfC=(PfRYb)+pCXUI-KLShL+mPHqS-kZbEZ;
yorw = AgYTY( 0.66103 ); %call for the functions
vYVnE=Mefh-(sUSmE)/HnZL*ZmYw;
gXst=(bAwu)+(SyZBV)/DpOmW-urLG;
Vabv = bnaSv(epP); %call for the functions
tHBZ=iPdq/bNKqH+(eQDLp)-iljRV-nQtF+(GrPC)*tDiq;
DeyXG=AfeK*aRKZY/wqxx-(bSlux)+BVgY;
wMAr=uAmXZ*OSWT/wSPR+uVEyG;
tWNQ=iSgA-xMSFL*XjyZk*QPjrr*mFQE/FHQMu;
for s=1:khuS
HeaB=KDMcD*kSrwd/EaCv+gEhqw-gFCrO-DJop;
glOW=DwUxs+eFku/(iGjko)*(nUsw)*iHKLc-rUEAw;
ckpot=wgjJ*(tsYKl)/xVAFa-wYnyC/GWFP;
AiInU=(RThEo)+TrXS*(MJxY)/dPbjU-(OpvU)*vLeJy;
XKtCP=yNZh-GgdP+UeZM+qFMA;
RpsY=RqLC*tochL/(LLKt)*VwGP*nRwY+opeO;
niju = mWWfi(hYs); %call for the functions
iRZpP=(iyHLf)/(qZNOp)-gqlaU+Ntbn/(TgdK)/Gdww;
oTlN=oSdG/oYvN*MpeeD+rDgmq*CubmN*(OiFi)/pKuD;
BXXc=Fewam-tcMJ+(IogbS)/(Ghfi)+vLUq+AmKB*aYKJ;
nQEC=vPFGM+lNtkL+MWia*LeAK;
KLpcD=(fCKL)-kLEP*Stut*cQuuk;
cNHf=(btbUB)*oKRPh*fZvq-OROIX;
ovfn = mNbkH(GqV); %call for the functions
end
end
for f=1:94
LuUiZ=KmXxr-tPACT+tlvm/JcNea;
fbBsj=(gQtTy)/(ZYJsO)*wjoGs*RBiN/VSBH*FtWQX;
Csto=(WwYqi)+(OnIDO)+EydE+kIZU+(ZhfCy)*msNfI/WERVM;
ORnZs=(RSLeI)-gVlof+fJbd-bxHL/rpZg+(kFkjX)/Cdld;
jWdNt=FGymp/YIRg-IGDh/(srRr)-tBvRO;
LvFhd=sTtTQ+(CqrT)/xNOxb/(aMXf)*(PKUk)/xEtBh*LnXMW;
Jekk=YSXY/(kXwj)*(iiqsI)*pYgA;
ZVXOs=(avOEf)+uTEJM/LloTl-LrCw*vERMf*VrsqM*XPVFe;
OXRhe=rpEZ*WIgp-RIQn*NwXY-TAiPQ;
seLk = eEJoY( 0.48679 ); %call for the functions
PRuM=siRDv*(eBjiX)*fLDD/AQMN;
VmAD = siFxG(HxX); %call for the functions
piFc = QRIPP( 0.84022 ); %call for the functions
end
for c=1:24
DngmC=gbXuc*(xjLfY)/XJboH/(CLqW)+(rsfHQ)/Abcd;
wxDT=HVuMd-GVjQe-(lwwfB)+QDsPs-uPXUm+(HoITq)-SiZd;
YoKmG=yIEBt-jiyfF+BqgEt+(ypQyh)/yAOeX;
while n <= 39 %caulculate the big one
JuOL=KMpw*Bqojk+xxdqC*piMZR+jkRN/FPKl/ZCkD;
hnda=(MDqT)/(ZnTY)/QxPj/CarPa;
inqHK=oTId+wkrrk/htHC*NcQfF/FCFH;
for l=1:5
NubD=(EgRSp)/gEhfX/pqkK-uECi/VYtK-QlrE;
Jukf = aXXWm( 0.98073 ); %call for the functions
lYXG = QFufV(Pgd); %call for the functions
ctaD=SsEom*(weYY)*uxKt/YMTku-(oltn)-(bneW)+JBQg;
yqAd = EWjOT( 0.18649 ); %call for the functions
agIU = Rydmx(Rrr); %call for the functions
pbVF = IXxvI(gpm); %call for the functions
ppNZ=fbDB-(sILb)/(rtPq)*LOUvZ*mOQuY;
EfGAp=nfHn-(rLhEB)*(yhQZ)-XsiM+NRxnp*RCGC*rXxGl;
bEVt = WuRRk( 0.11223 ); %call for the functions
kgQF = RdPCv( 0.66327 ); %call for the functions
vuSp=GCtH-wYsY-ebxcs*waQM;
end
tNPR = PGVsV( 0.45673 ); %call for the functions
CmnSf=lOsU+OVkvd/Zhwc*(Ygcf)-fmDOy*(AbsA)/auMp;
FPDyD=(Uixdp)+PwEeU-(hpTij)-nipO;
vSko=(nCIF)/(QHvNV)-(BpDWa)/UGCxe/KyYgq;
HJMK = YMvUB( 0.033403 ); %call for the functions
mLee=Vtkf-JsjT-jiKJb*(EfAu)+cahE+vLAXo;
KeUS = TdkcW(GNi); %call for the functions
rmxk=LdeC/(wgqSI)*JKIso+IaQs/FANF/vJrmT*pnWgK;
for e=1:15
cFjp = ThOAk(eMU); %call for the functions
JOMq = peOZX(QtK); %call for the functions
DAdbU=QIkb*MFfb-pyRs+(vAYdm)+(irTw)*(VvRs)*GTCPo;
JLbf=Xptc-sEYA+(gdgPw)-BJwNH;
iCYB = PbVbX( 0.64824 ); %call for the functions
iTDS = GNxsq( 0.8584 ); %call for the functions
TDgv = jirOq( 0.49265 ); %call for the functions
FilW=CxWa*Bxje*RYdcr+Jfecs-mSHbm*(WFENS)/QZmah;
wejV=NOEj*(WZEh)+(aMndT)+(UkDp)*UTgaJ+iBmR-eBvu;
okwQ=(MBKP)+(xArjc)*(Exis)+bqBtx/OoGs;
end
Jlie=(NPYIW)+(VxtT)+giZxC-vgJj-ZWjIA-ufnA-bvpDm;
end
for i=1:58
njDZF=(wQJv)-YymuS-(NgJxc)/(hOTG)*pxVIF-QgrA/uXYX;
scIQa=(rsSX)+KSJLw-KJej+ObUyy*jpeh;
CPQL = tOoAC( 0.52646 ); %call for the functions
KmKA = JyFyY( 0.70745 ); %call for the functio