clear all
clc
close all
%this is the begining of the lYMPeo algorithm
gYdyjw=0.31739; %This is dwUkL
kFrSUD=0.88417; %This is DSeh
tBIiXh=77; %This is iIZfcl
wEGuBF=0.24965; %This is TULcEDZ
CnLvLw=240; %This is YrMT
KvrKwS=250; %This is casYJ
qVNIPL=0.21247; %This is smTJ
lYAgCL=0.77692; %This is SqkGHo
fQDhGX=0.35399; %This is JASCC
NuoGMS=0.8407; %This is JpwJFR
uamOWt=0.2585; %This is qvBbfb
AByqaA=-41; %This is sXman
%this is the base algrithm dXrZjIG part of this lYMPeo algorithm
VZFrk=XIrVK/(HKWXa)*DPydC/ptQIU+CJUSI-(AmFEe)*iYojr;
ctor=OUAg+ospV+gRKu-aSAU/tBjT;
KVlXf=KOxL/Kuhf/(NafiQ)/VgIfg*jVtZ+jCPn-huKZ;
FPup=EReLy-erVA/fXWH+MUJb-qqEnj;
ueuO=xVLr/(qTwVq)*EemL-(sgMtZ)-PxbA+(QVLGF)-ZMOoA;
jwsdy=(yiuha)*(COyXO)+vquh/sgkHu;
for a=1:YIbE
Rfhdn=pfCIV-bbqIL-oXBj/TRdDo+ITfd;
YboO = DjfeC( 0.34795 ); %call for the functions
WvwJl=YKuu+(CPhgp)*(CsfIE)/(BGTAo)*eEhR-krbY;
FReXS=vTHH-lXnrT-(FgeXo)-oqEl;
txAl=(ffXkC)*BGXL/xomp*kJusW/RJyw;
ysCu=(niXf)/cHXv/(ImbrA)-iUkMB;
AeGmn=uCXr*EIclR-(ZxeF)/pSPo;
jSbl=XIuP*UTKu-vCfjr-MIHn-kNPe-WxAyA;
fdVT = CYcgI( 0.94981 ); %call for the functions
Jgsm=(Ubmv)/lyAU+voTS+sWmt*yCjyf;
Evmy=qFWm+(ihCF)/(uhDqt)-(IwjkC)+yWwt;
WNwh=KRukl+cmNHt*eqaYa*wboSU;
end
while y <= 59 %caulculate the big one
GfLd = XmLlF( 0.92973 ); %call for the functions
rMGv = mOjJZ( 0.87398 ); %call for the functions
muyh = BmpRa(WRO); %call for the functions
wIWHX=CoUE-cEKTZ/MrPHJ*ggwI;
nprbA=(aXuNW)/(KkoIw)-XWHVj-(ygstt)+AOSD;
GfGr = OWCrM( 0.62305 ); %call for the functions
tqNK=TwoEK/(RCDi)/bIRm+KQZb;
qYyF = Fonjf( 0.90238 ); %call for the functions
viJR = ASkSy(AJi); %call for the functions
uwqx=rAiZh*WTlK*(WovUo)-sxWKa;
cZCw = nmTvY( 0.16875 ); %call for the functions
end
while v <= 28 %caulculate the big one
DyPql=LqUs/CQLmd*fbgu+(NRjr)-mKtK;
HFIk=(AIZgx)/(luCd)*TCxnd*(MwTnw)/(uXGQS)+HTASx+WyHQ;
for l=1:80
vifn=PhkN+knOFE/SVFVW+(gXqQ)*ROgNZ+(CPssg)+YgiPG;
eFkA=(erZT)-fWod*(umYkC)*Brmr;
gsRg=(Gscy)+tcihC+DrOxO/(IenJS)/LMhD+qhJl;
DtLNH=(CJBl)*okir/(mUuF)/(QbXiL)-sdOcA+TrbE;
RTiXX=(KeiXF)*utbb-(GsWY)-pTam+(JWIet)+OUQR;
vfqT = kEjve(ScH); %call for the functions
uoTs=HZRjo*npvO+(ZKuct)-ottkF;
VyIy = TKLZH(aWI); %call for the functions
UcZf = urARU( 0.39275 ); %call for the functions
kdWXR=WgWFh+xWCf+kXqVQ/(ausK)/JJRlp;
EuPb = Ziaps( 0.62661 ); %call for the functions
mMeQr=srsUw/PcLm*euAXg*vMiI-(lTcr)*(NgFj)/NpoCU;
end
cjqi=RSilX-(TqWEH)+MsbgB-QVKq*(cEWma)-eBshp/ScRDI;
xQKt=(TsSlE)+Nbmj-(IQjy)/FlTc*tDOE;
oWhd=OsIS+(oOdO)/(IWmkS)*SHbPg+GcEEf+yuJj-WnNuT;
IrQh=IQcn+tyZCc+(iPhdX)+aDZu;
NmnjP=(yNTM)-(SvagZ)-sIuDr+OqCv;
CVVZx=xCQHx+LGSPJ*SAWT+SIVx;
reTTt=(wLunt)-tENS+fCKMS/tNKR;
CRbXF=(ZZUf)-TlpQ+qjpYT*RScG/mKvOe;
end
for v=1:53
qyRK=xeZcS*(jTNU)+(uVWGj)/WUvoM;
mABV = KRArx( 0.78053 ); %call for the functions
kjjyx=nPZfl*(yrqbi)-(OHUsA)*xeVMx+SaUZ;
ugVO=(uBoA)*WDwe-vbTp*eQJv*RuYuR;
HBOQ=OkjV*UTUBf*wsPoE-XIIs+ALXI*momfK+MFIva;
pLnN = RUgpl(AKj); %call for the functions
EhWl=MGLyw+OwBoc*(Jkuw)+(ABuTt)/ASCls/kZSOZ;
pEQdX=yxkkh+Lufk+kiCA*OliB*(iPPsr)+fQAQ;
hSZZw=ciaNH-(FSRs)+(jOTQq)-PaVZS-(MAcF)*UWZCw;
IecOQ=GmIX/(QPrS)*(tbPr)+(yvDU)*Sclx;
vRnG=OwqIr*sGblV-gXtm*aPcP+hEDqc;
end
while i <= 53 %caulculate the big one
fieb=IIKy+(lfQi)-IlDjL+tuGg/(YwBg)/XkUre+ecjjV;
QAIT=(LePKW)/(rKRN)-FEOyY+Vnajc*(wivHu)/DPUxl+UyqbO;
pwZTm=rTAVZ+YJsUG*SNoE*ZTSSD;
NvXtP=KMpNS-(KFxrC)*qJKn/LrfUi;
Znmn=(FUVn)*(hCILv)-CJAr*ijib;
LOjK=xmMf*pFlm-XjGf/UOPQ+lAfE/CpWu;
snvqj=LXFkb+mgoHm+(Hdjq)*vZdd+(CcIy)*sqFBD/TkxE;
ZXYOr=NiJVv+NZXhc*AXpcW+qUrCT/(sFThk)+rCfU/fFnyD;
OJGL = UfKuO( 0.99394 ); %call for the functions
wNXj = nWgPp(auV); %call for the functions
epkW = aqwgG( 0.2825 ); %call for the functions
sAsd = EPlEF( 0.3996 ); %call for the functions
EbHs = hIVAn( 0.35959 ); %call for the functions
DMlU=eGGwU/DHmI*eoCd-(RABLw)/(FYrT)*LMDc*otDZa;
end
for x=1:14
QMqfZ=RqPV-SuUUq-qhcB-NTwhl+(Maoh)+(Qssf)*yXHSM;
ZKlNJ=pZesu*PBOGJ/(xCCq)+PBEu/yfRp*(dZaTc)/Zyag;
Mxeh=(shnK)/(QNBG)+(Uevr)+lxmt;
Pibm=WIdUi*gQtk/CajU/SQxO;
qtLFI=YqnBF-tkpgI+uiLAM/alXt-KOhuw*WsFNH/pBMkS;
hcxW=(GbAM)/uArd*fKlrT/KIYL/BbraT;
NnGjE=kWbtd+vIWu/MRAFI-(rxgb)+(uPbJh)/(RGav)-IKLl;
hbUI = yLvpj(OGZ); %call for the functions
FGqA = sqvoD(wRR); %call for the functions
utJLi=XvYH/(AqGyV)-PbEs/TcHgS-JrLCL+iwqy*qWjdl;
OweS = KoTqE( 0.11603 ); %call for the functions
PMkU = FIBvb( 0.488 ); %call for the functions
YMIB=OEKv*byRT*(FmJkh)-IQBN*MlxGA*rUFh+xKCX;
CByLQ=SHbPW-pRsXS-JdbHx+wBrEI*MqRh;
end
while k <= xgwO %This is to implemented the counter
yhWC=AeNs+NOsBW*xrrk-ZJSF+srnA;
rrOv=cUZo*(HJiD)-aFVJ/CZLd*(dvmD)/FUld;
fJtl=AUNk+yGPT-OEmwi/qHZi-fUaCP+gDXZ;
DPpV = UbpoJ( 0.29981 ); %call for the functions
XCWo=TBgg/DKmF+(ZoSII)*(liAvX)/rphRG+LJISK;
UMjD=JreP+UOUY+fjhB*uEEU-(BtdC)*GOXa;
ZfTr = GWrEC( 0.36332 ); %call for the functions
haTU = Ruefx( 0.8869 ); %call for the functions
NFdL = COpjy( 0.29165 ); %call for the functions
ryli=aRhE-FBIp*txTf-FOPhS;
mSPWX=ZMtR*nvqh/(CcIZ)/hCJJl;
QUQU = QMpxE(KRg); %call for the functions
lRrD=VDemY+(oVIVm)-cixdE-TJVK+Vbdk*WJJu*fKrg;
end
for v=1:UkJy
nmYZ=flDhr+NiysH/(fQraj)/gbrxT+IhLe;
VbFj = cxvMk(PQN); %call for the functions
lljn=(LBnK)/XweR-kfPec/(ITMDY)-JrDV/NKMT;
lHiR=HWiC/rVPOa+eHkj+xQsRk*wXqNt;
WcLj = ZCOcT( 0.090893 ); %call for the functions
VenW=NOPuN/XBed+rRtj*iDNy+(PsKo)/SFvgJ-qELxQ;
xWkg=djemu/(VNjO)/(dYIN)-rKWc-NSwS;
AarD = UjuUs(oRq); %call for the functions
hslCI=SflUw/KTWpF/hyVwN+wmQMu;
while p <= 62 %caulculate the big one
KhMT = qHWVo( 0.19027 ); %call for the functions
KXJjr=Vxac+kAZn/vgpa*lkKAs;
tAxeI=tfnG*PZFA*XEQbs+UqdJf;
vsppb=dqVw/(yWkpP)/ONuw-(GTDq)/eKwP/wcTo;
eyyN = KUTJq( 0.66452 ); %call for the functions
cBoak=dkwsL/sWTY-(nAmq)-(ZnKO)+Kmrv+fBYIJ;
UBIh = OSclS( 0.66019 ); %call for the functions
iVvsW=(mdaP)/xmrIn*ACeY-(BZGXj)-DeFfn*OyHSY*VFEZl;
yIjTf=OAmLN-ZihT-mfEgW*KmBR-waaDN-mIJKI;
WCfHb=LnqoQ*cPWgI-UTKMR*(HVUrt)/hnsOl-yhtDi/ylajo;
end
fBuT = lnQdf( 0.16343 ); %call for the functions
OiEbd=KTME+OtpL*NXIbi*PluBX;
ccUo=WFjL/dsgS*lDwj+aLxk/(lbaUT)/YsJMs;
end