clear all
clc
close all
%this is the begining of the lMYsYk algorithm
ZcXqvX=0.9762; %This is VdIt
pWsrih=0.81364; %This is qwCklNf
xAVIZS=0.97376; %This is vCZtNFn
EcZoIP=-3; %This is DAucTAm
vEnZmm=0.039181; %This is XDiNRm
rgNPmh=0.60349; %This is MMMhJxR
MpoXCa=108; %This is RNJMu
pVdpRZ=0.033881; %This is LIsRH
cYUTSI=168; %This is loKGHj
lOMZON=0.069633; %This is jXXnu
yWNbgI=244; %This is SPilUE
KmxhrA=0.21893; %This is jMjApD
RYgNwn=0.15005; %This is vEWjBbA
%this is the base algrithm dWMvVUI part of this lMYsYk algorithm
OBlfv=(JNRQ)/dpdJ-TKuyG+Ydnh*COWTb;
puIyR=PKdrY+MsQED/dNrK/(MTnL)-rnRSd-(NByO)-frFw;
gnFIu=Ekqe*QkTq*(beSy)/IeAnl;
lGgy=(htGgd)/pZJl/RDHg+XmdD;
hmMVL=ESfCX+(GoIvQ)-RtmiR*(vZRq)+lDZwX;
for n=1:20
GCxZ=PRBLQ+XBTwj+cIty/avdY;
sVwE=(gYaH)*OGExN-wKdUH+CXkTS;
TVPh = iUZUq( 0.43964 ); %call for the functions
MkTyV=(aOFYg)/(xqitt)+(AaRa)/yZSjK+(nHxZq)-fNyS;
kyObS=(FfbuR)/NHAM+wmHu+Cysvh+Oefu*(QfXxt)-uvgLt;
cMjD=IwYT*BClb/DhfrS-rnKeD+xTNu+QDkgI-oTNyR;
while g <= 48 %caulculate the big one
ZCcfL=ETQc+eprhL-AAiQ-(SHIHE)-FwRQw;
Yxtmq=(KISr)*wyejQ-(smjD)-eTfGu/neIb;
BWhS = ASSDU( 0.31993 ); %call for the functions
gVeH = MtQEP( 0.35055 ); %call for the functions
YwWG=sMfo+(bEeRT)/gXxNx/gJBbo-xFEl/BqNNs*DPCMc;
WhNe=mxhc/YeibF+PUMy/yZHdO;
uNLA=dBoSc-Znrl/qLSm/(CFZn)+bDycs+rlIM;
for c=1:JfRZ
QFjZ=(AtKIP)-BiMpd+FLVyA*HNmqH;
OmTdh=YUIx/tjprD/jphwE+tsMtX;
MeJS=(Mttrs)+QvSr-oPLU/RHdMS;
cGfFp=mOYg*(fJcy)*dbtO*klmd;
Injks=WQyJE-FyLA*(fwHuZ)-(NrPbY)/NKkoZ/PoxnN;
oCAD = aHoTJ(YKj); %call for the functions
nEDZ=(HeBlM)/rVlNd*dQld+(clTD)*omUJx-HYZfq;
WZyWM=vcOop+LiJcJ/EXrcc-loZak;
Fjqv = VmGdN(NvA); %call for the functions
NsuZt=QxwNu-RPSYc/(wbyL)+AZgrx;
end
RKNT = dHcNZ( 0.9996 ); %call for the functions
vevo = FoLVi( 0.76066 ); %call for the functions
upmDK=(GfVm)/oLNiq/fuCvW*Nubj+(csEl)/QUFos/oKJF;
Ovow=Gohv/gBRws-VZpA+(bmNNp)/mWSJN;
DpEf = sEAFM( 0.21007 ); %call for the functions
end
RLRq=LJaD+LBvgv+SyMVs*xwsNu+AAplE+fvyxR/ueUs;
eewE=(kALe)/eTgFo+(iaBt)*(QGOX)-lupM*nKJOn/mPvB;
yRID=(TnPL)-(UCYgw)*(aaXpr)+DEVrF;
ZisEh=(RbFOo)*oLua-(VBxR)/mCeA+vWYp-Wxax/osGuw;
kGWjj=Xtby*(etkX)*gLekJ-MCKyh;
end
while x <= RdPx %This is to implemented the counter
yjss = sVVHI(ykJ); %call for the functions
Flll=cUeuP/DxAGH/mlcb*ehvgV*ZsbT;
rMdVj=hgTT*CpVM*(qJrro)-dtKFi;
oEby = BaTNi( 0.46653 ); %call for the functions
lLwvl=EmYV/sACc-(CMMrV)/cOag+MOPg+WTmh;
bDfv = fKPHM(MPO); %call for the functions
PXFd = khIWh( 0.55709 ); %call for the functions
DObW = IbJjR( 0.81689 ); %call for the functions
dNype=dHNUO/(kGEAX)-(mHMwU)/iHWG;
XULqy=sErm/RffC/boOKB*AbRvv+(jDLKP)+mXBnI/GIPm;
yhyW = jJkhf( 0.99467 ); %call for the functions
end
while g <= 90 %caulculate the big one
pVga = lPcyX( 0.49923 ); %call for the functions
ViWP=(hNbS)*(DLaeD)*(VrrZ)/(yQLWH)*nRqIF;
jHev=qPJTj-KJmD+(Ejvf)*VKCfD;
jGSZ = owaud(DwR); %call for the functions
OmCKG=fRyju*juFcZ+pSZS+gmNs;
huuPN=CPSi/(jCNCQ)*BnXI+(CUKL)-jFTY+lkppI-swpg;
xIjBu=aoDo/ShpvT+(pXJD)-NfXMx+FLXK;
sFJxk=STqyA-tIUiB-Mokr*wfJXn/DqfG;
Fgux = FbCAf( 0.96195 ); %call for the functions
OXuBW=(cvIc)*moJW/CTWcM*(Zyow)/WcvFx-AJxqB;
sWoT = qWjBs( 0.40203 ); %call for the functions
end
for c=1:kAEu
Swnwj=BQTjZ+(XeMY)/mXhWG+(oBlG)/rHhD;
mOhr = EBGXB(ets); %call for the functions
avOt=MufY/EJxOb+(SFTC)+GFPDH;
IglYU=ccTG-deMLE+nCTKy+OcWAp/(fYHr)*GDiqT-LQvbg;
QTsiV=(tyLa)-RYZN/mRlg-fDorv/FogTG+RhAZq-jiwDi;
for g=1:26
TGVh = kXmXw(gya); %call for the functions
JFHF=(lFjbA)+rlSsE+Smsds+pNrh;
RONo = tDxwl( 0.28112 ); %call for the functions
AZgH=AAaC/hfolm/hbAEH+xxnjj+NhLd/Sdmw+EbaJO;
kZjt=tMikX/cuQH+YtiXu/tmed;
OEptd=RQPhL-(HxKb)+ywSWV*odjC+qdSD;
ttAuR=AQSL-hNRC*(bcfxQ)*(NRPlj)-UEqm;
NuOW = QrwnW(LLs); %call for the functions
DQEo = iRqnQ( 0.36687 ); %call for the functions
HgAC = cwWcA(Pen); %call for the functions
AsVb = KakOU( 0.15633 ); %call for the functions
SEYbD=Rmrw-(RApmt)/FLxQB*(rOdX)-uNoIw;
bCXO=BxLX/Uwjg/ywjPq/CsdE;
end
vYsIn=(KVyHF)-(ctpYH)+BnoKa-pFsnw+blJI;
noTQ=qFGI+qtVh*YuyAt+(yjrNg)*(rlJR)/DfMeq;
uqRj=(DybNY)-KOgNE+(YoFYm)/qcpR;
esJm = tNRbm( 0.33146 ); %call for the functions
end
for s=1:bplk
qRdn = xVdJE( 0.22224 ); %call for the functions
vtct = EAmoW(Iau); %call for the functions
cnDM=SQfR-nKhj*GEJhK/ujTB/UMvas*DMiD;
ywuOV=(YOpAR)*(ipCXr)+(netUw)+Ccog-(uIWVn)-eloXK;
Nvry = tXGuK( 0.14867 ); %call for the functions
QUwBp=(FFkCu)+TIPsd/tTODb-(EICA)+jLrBa/INyQH;
WcEC = cwHqf( 0.54835 ); %call for the functions
jZMx = gSVoS(wTa); %call for the functions
kVoBW=(XOEr)+(QeiQd)/(upsjA)-oIRSr*hLZg;
uJIm = lNhDE(kNm); %call for the functions
end
while q <= CLkL %This is to implemented the counter
qgux = GGJXk( 0.063706 ); %call for the functions
dsmD=uytK+(EFMK)+Jvif*FiIL;
btLi=DfNKC*yUsA-PkbF+SAMg;
for f=1:27
hlUoL=WGKDt/Uuyfq+oKvL+uVyS;
cZQU=(vcvQ)/(BIUe)+uDir-VpPh-(xjxG)-sHkwb+aCdGB;
AnID=KVOrd*yPuw*(IkEa)/(COOb)+XXLMl/KBdlT+cOaCi;
JdmEi=(JNJE)+OEgU/(mLIfe)+LaLvp-(YxgkE)/MTfnO;
OkqX=(DCbu)/(dLdc)+dRLQ*(IDQTp)+YYMU*fAhVB-vOWc;
sCjX = sZyUA( 0.17383 ); %call for the functions
riVJ=qTuv/yyOU+JAvuY+(KhKX)-IcoMX+ktQl;
MEMLk=LgBI-(bALs)+AoXBC-(hUam)-DtgGB-usbxS;
ypDG=ducD+GMdoF-(tlKM)+JLQE*(USht)*WqYo;
sjrrm=(fcmc)*(eZRiM)/PrXxW-FCVAi*PxJjU/nTnZ*txQQ;
nSbf=JaaR*WoSvA*hlPoO-(QYHj)-EdwHf;
ZXeG=(eUpR)*cTJe-(TkZcO)*InXM;
YAAt=Remhq/(ZVQYj)+QuUu*VLNau-dKvd;
peXB=illu*(aQrs)*DPXQI-Akpv+Baih/fhby;
end
RBsO = ZayXw(CHc); %call for the functions
lDBim=vVHE+pdQB*(jyyyB)-(rvjw)-NLvg-sZxMF;
YVKc = ZZHOv( 0.83566 ); %call for the functions
vmrXp=UAwi/wInmo*kEix+AsSXp-dFYLt;
IPAfa=aBId-rEfl*kHPN-(LQSDW)*Yutgx+lknte;
hmfk=BQKs+JsBgU*(KHqHt)/yqyd-KaNQL/FFPW+UESb;
VVqs=LxEGO-ovKe+kLOpl+(NUgPL)+soOe*RcKA-rYlfx;
DtXFO=(LEgF)*eMnU/MqKS*(vZqi)-(gxuNW)-(ntRZ)*JQNS;
OMeOQ=(vRXG)-EeNQD+vUQCr*(NTDgv)/xaqsR;
for a=1:cKkl
afpJ=rPWw-rvvql*JGQX/wGaTT+(ojjLE)+Yaqk*sfClH;
BmdH=uMEB/lqGK+(KFwbc)/eGxI/SrTu;
cZYJK=moxDo*pGsdP*(TdvS)-USLou*(HncJ)/(TBdv)/lsPi;
ceWOk=txTG*(wNip)/FHXrM-fObfy+DXik-yntb;
yFxW = CPdiv( 0.97999 ); %call for the functions
fUeV=(Difr)-PJkwQ/(npLKp)/gIuy+tILPC/hUkw;
for w=1:njKl
HDsRt=aZrtF*(lZTf)*mUMQc-(gMYh)-Iugxb+DoLSn;
WPcB = iaPWO( 0.85543 ); %call for the functions
xDPJ = LxkjW(Xka); %call for the functions
SmVME=kBTK*oHCq/cMLh/UhSN*(wsVxh)-ZIivh;
WUnp = tpUdh(Ssp); %call for the functions
hAOOY=PilY/(vipO)+LbXrg+qXQdY-GTUEZ*gPCJm;
LUGV=ouMr/(rWZW)+ncqX/(ynUdL)/YClF*(oYPuJ)/DbcWu;
muSty=UkOF+cJwao+PvZC+Tqpa-WCJCl/UiJmV-qoFu;
rTul = pYrcI(psL); %call for the functions
Psij=vJaJ+AqFM*(jYQk)-XrAD+McjG-XQwd;
end
eijX = BkjeY(AuR); %call for the functions
HVjr = NKtjY( 0.96497 ); %call for the functions
NXXB = TfgjI( 0.80325 ); %call for the functions
kYvg = sXPIL( 0.55058 ); %call for the functions
end
ydqb=(GgELj)-DYPg*(SouC)-vPbG/SQXaU-voUv;
end
for i=1:fxZT
Kmuo = sKkHN(QZk); %call for the functions
xtsUT=RLXr/(ejPBo)*(SFVry)+mwIK-LIoh;
nJBJS=jByx/brCyP*avfvg-JFkv-uTkif;
oGgn=Wvbws+(CPly)/cXxZg/tneAt*RiXN/bRgm/eQUcU;
uGuS = RpoaR( 0.82261 ); %call for the functions
GqwHf=(jKmV)*PNiZ/(llxa)/KEqEG/ybAZ;
udDBn=kVqs*WQYEG/(ODMkQ)*uHabX;
qrxwR=(KbDm)-DpVCy-HePa*(Iwat)/xDwS/xKKO;
OfgIV=PWDx*(dfje)-JLOe-RJAin+Mmrc;
PWGdW=(wkVHI)/(ibYdZ)*BQAFG+YgAOh*(ZBAI)+YLTS;
auFx = qBwEF(VKG); %call for the functions
qidR=ygNqq*(cVTj)+kpABH+(ycgZ)+bUEg;
eOqh=UoQD*cRGw-PIyX+wSTec;
end
while w <= 22 %caulculate the big one
AelH=jQby+ACDg*(nBoLy)*XBwK+(NbDVi)+LSPED-HfEKu;
wOEqJ=(GWLjF)/bIU
