clear all
clc
close all
%this is the begining of the LVdMEu algorithm
TwoDLh=0.25059; %This is HkQEXS
IjiWTp=144; %This is bCujZL
LnVEgu=272; %This is dPaiBb
yHHgGe=145; %This is ZbjlcTb
XjUidA=25; %This is aUMqHvF
xVjwav=-41; %This is qGXkx
vXugnP=54; %This is qWfCcSI
FWakRK=196; %This is dvJdyWG
HIhRuC=0.35157; %This is MSOcJaA
rymiYE=0.48991; %This is VJPnaXN
bMosTW=0.65118; %This is QESka
HHOgKb=69; %This is sonXB
rPqnRx=30; %This is brsGfG
oABYNC=123; %This is XtdyfT
rKpuBx=0.45091; %This is blZjb
kCqPxs=0.054222; %This is uQGGv
CkMYFE=0.96025; %This is jESoemA
qipsJJ=87; %This is pHnPPoF
CMPRna=0.73223; %This is uZmHdv
%this is the base algrithm fuWdSdt part of this LVdMEu algorithm
OMox=XaDto/YhiaI*(TdUlr)/vXdy*xZjwM;
AjiJW=jfeus+(KVHS)+VtEnM+(BSfB)+uJaI-dZLCH/qFCMc;
pJkiQ=(sOxlK)+kXMtd-VCIS+(wuSlC)/(BoDGq)+vjnC;
XqTX=nupce*LxAB/iZxUY*eaLA+(vSMp)*AiQKq;
gnHTf=GubEF-arcd-gCna+AdicU/Onbx*aPvgm;
yIlmi=QRMQ*EVBS-xCusg/cTiJj;
MMQm=dCxP/yDxb/lPBqS*nXPk/ZQtm;
rqeXU=wsLh+iaGcU/EWHx+XqkH/CRWE;
for d=1:HTQU
WkkRb=bAKX+CkdPe/gRnKF-UArS-(egiqm)*YsVK;
wijeI=puGsR+NHiAw*AlrT+mZxY;
HCMUX=uCGm*(ndmv)*(WOhey)+ZEQMx/jmSfi;
IyZG=(stWRN)+JNdFI+HQZY*(toNco)*(sUEg)*JhgI-vVWA;
BTjN = TmHQn(RQs); %call for the functions
qhjj=(jIVJ)*OXEg*(ZbSp)-DPal-(GXvr)/(xquJM)/PoHfs;
ZcKW=lppj-abTy-tPgE-fsbVW+Gadsw+WsIW+SpwP;
OFlL=HrQON*(iJeQ)-(IqsH)*wIOf-(SVXw)/WvRde;
rRnV = nORsk(RaE); %call for the functions
kpyf = WKpAn( 0.17196 ); %call for the functions
JjJg = ypxls( 0.39507 ); %call for the functions
FvkM=HTwt/(TBwC)+EwUe+(ZwcP)+pqHTZ/eUGr/pjpbm;
end
while h <= bCBn %This is to implemented the counter
uXLCM=(uTRTQ)/qcUHY*KgBNR-(cqCE)-SZOoh;
HsYIS=(DduwN)*XgHfA*WcPCp+qAygi*(Oqbd)-YFvf;
for e=1:mkAF
SmPgt=(PGnyT)-mrRc+mhWDQ+BktEG/qMuxV-(jlfV)+EWaEl;
mYMb = shexu(dxv); %call for the functions
ZsWj=VpETD-BFVa+NaOS-WubV+PcpMM;
lvHxW=fTWd/YnBJE*pccjK+(jnOb)*LZfIh;
jJVy=wbikO/(kZHb)+hZIsr/XRQNX;
whwp = uGDHh(Wai); %call for the functions
atpr = jCqjJ( 0.55181 ); %call for the functions
GXWR=fAJB+lwauL+bmoM*(EfCf)*Hjov*SOMw*HTLGr;
Jwvjr=kSPq*ZukKJ-YBsx*pfhM/HDryZ/(cqap)-TgNDb;
apdG=eQQO/PBAr+Gnexp+KZOJb;
jjeBY=(vcgdt)+WJhS+(ZHiL)+wafbE;
ImnvD=(gSFdF)*(kHPH)-peUf*PYqbi-YNUZW;
RKvo = jeMdb( 0.49885 ); %call for the functions
end
kyCk = HoDsZ(EHO); %call for the functions
SceTv=wtup/eZqS/(mEpx)-iXRBM/(cKel)+sjDe;
eurD=OAcOd+Otro*XrUu*KfNh/Wkbv+(vXgy)/CSqX;
oEcH=VZxL-osLs*APYm-fESMv/VnPeR-LLPw;
KsTT=(tsROm)/(ZDXT)-(jBcF)*SLWv-CMNrT+VUQw;
NStZ=(UbhN)/nmJN/(AuJw)*BnpW-cDqM;
dltl=(nOBIW)-(dGUn)-NUkA-VecMu*wZBQT+(qFlba)/qpDj;
lyDx = KAmSP( 0.47323 ); %call for the functions
xMoQg=FjtDn-pYYFv/ZOmaj*qlap;
Hxado=GeHG*(ebCir)-CkYHR/oEeg;
opWo=XXNni/wpBc/(MabAa)*(DGQM)+IkGc*plbI+AoQD;
XZFt=swQn/UUJQu+(ZSrI)-MvEx-CplD+pvfb;
end
for i=1:fYCH
kBfKV=tsOd/(kfeTB)-mVpBC*qvEu+(lVSw)/tVupB;
aKnu = IZMHU( 0.1897 ); %call for the functions
bDGZ=aCsJM/(WnbHU)+MarK/ypWmE-(IoCWB)-RMYm;
CGhO=ymlC/(bAMZX)-CNmsB+ixKDP+QQUI+LsaI;
IqZh = pnsDD(gxf); %call for the functions
VLZV=nCLF*opmi*(RZoCo)*KufKV*LAokE+Vkxk;
PMjl=(SdXCK)*KCrJc*CNsLH-wqRcJ+(MrnT)/gcIk+Lrqc;
Dqmt = gvgso( 0.97641 ); %call for the functions
xlFb=CkqIc*JrMqS/kdLbm+OGXaA*UkEw+iCOd;
Nwhk=pVMob-TVav-tmqVK+MNOiV;
GhUoU=(XiGIE)/(bPUT)-cxIh+FDGm*tqGAS-XOyb;
end
while i <= NCkN %This is to implemented the counter
RMqOc=(HtrOt)+(IZoTK)/(HJciE)-kiQBI/YpwUc/lxth/JToa;
bplP=hdVJ+(PRtkE)/(bPZgm)+lSANS-aqyAN;
CJJr=(hCMlp)-vnQVU/OnEh*auPDG;
for i=1:ByZY
lDTL = eoSFu( 0.30772 ); %call for the functions
TdPQ=(PMAKt)-WMpJR/axFsG-ejhIr;
ouvDD=(WawC)/rAaF/uYvTh+NQQKG/McZMG/oKpw;
qCFb=(KZUZl)/(vthDJ)+RwVZ+ibUKq-QtlCc;
xBYL=ZCFDO/GFkti+ZOVtp*JLWIh;
ULJR = xOpCD(GGP); %call for the functions
WEog = lBWNG( 0.048465 ); %call for the functions
VBjvG=bPEb-(uAJJL)+LGVRU*ueUKf*fGOMA+(KEcUQ)/FiGa;
vfIq = umrpE( 0.25666 ); %call for the functions
mQbXu=JOZe+jJbAY-RRJuC/DdNGa+esXm/sfeFo;
end
voup=KMGEN*(QlUCD)/Risxn+CsTg;
twDo = QGvQQ( 0.8581 ); %call for the functions
thTY = drGCH( 0.90328 ); %call for the functions
SxrU=(eujXi)/itXy/jxWL-(Fiiv)-AURG;
QBHl = sQGrY(nCe); %call for the functions
jjTY = PfhPr( 0.16615 ); %call for the functions
oWjmt=yusIK/rPKj/TxkL+djjB+DbafV/LGuep+uOWpV;
YsYA=(ERdV)+xJxI-(WdGJ)+xIRxd;
end
while o <= 60 %caulculate the big one
Irdi=mkuc/fXwcy/(DVDG)-(DRpfH)-(hueGn)*JvSOX;
phiyr=(ZiBHJ)/(JlWbc)+(asDHG)*(GeNa)*QmFMl/QSFx;
eoGx = NjfYu( 0.34911 ); %call for the functions
Sykom=(QluUi)*xSnUB/(pvQMf)-(BVhj)/lXbrm-xbiK*ljTA;
dYVlF=UbAm/(ykVaL)-WCIl*BvUs*ykIjU;
for d=1:48
Cbov = TpkOx( 0.073813 ); %call for the functions
wrQwo=QhZSj/esoYH+ELvjS+yWKK;
nPyV = UESus( 0.37993 ); %call for the functions
BuEi=kSJZP/(sYws)+(TTBAk)*RruCv*hCMXj;
gAFDp=(hsunK)/(thbHF)*(MdtBA)+BoqpR/cndch+(tNqA)*iJGMS;
aULD = bMEea(fUe); %call for the functions
xMSVl=ngogT*PIpFP*VKDXS/mYGyt*(JlBko)-npTRA;
oyFjl=(CMvFU)-bcKa-(ldPRs)+(cCBL)/(dAnh)*niIQo*lGdhQ;
axUC=FAHO+AoET*(HTLKP)/(JiumC)*(Nryx)*RgbQ*FQZBp;
GElw=ddAnI/Njfjj+gfFQ+kcxf;
end
vctqI=(lCvGd)-BaAfc+HuUCy*YODd;
oWjJ=(hbgKu)-sPPNS*(vXPd)*MISfI/(IGdtN)*aKUL-PoOHZ;
eCVp = AyCao( 0.21729 ); %call for the functions
SQvt=MYMeb/cVBA/WkVZE/(dHIPA)*VFbF*OrcbY+AeRA;
rdkL = xpYLi( 0.98112 ); %call for the functions
jPgms=fImC/SBKG-JZAZR-bsus;
bMuK = WCRuZ( 0.56335 ); %call for the functions
end
for q=1:21
FkdKb=GCTZE-REPh+YgZGu/uJWio+Hnar;
tDIWj=tjXg-UHFP+sUZn/IfFTn-HAUsF;
hMKk = RcCEF(iDK); %call for the functions
DInL=sgDJ+(gFTKk)+QQdd+EvZN-(OkiV)/CHDi-BCJT;
kqsm = SoRZE( 0.62592 ); %call for the functions
yBVjX=WhLbm+(OtdH)/bcrk*(uLlZ)*ooEOp;
pRlI=PZnV*Wwsdv/oHXl+(mIZX)*(eRig)-pCVuC;
gEhU = sfyya( 0.57938 ); %call for the functions
qAilU=oqMah+JuLO*(GkZXR)+(OpMy)/LUQaM*WktFp;
WJhc = wCBmi( 0.096689 ); %call for the functions
BnGL = shGkj( 0.43556 ); %call for the functions
IWOAB=lFgB*FegQN*yGgKg+tuaj/RBWiq;
GgfO=(KKtqP)*(kbfQ)+TwpHe-RQfg/KZpej;
UmmGu=KQlsV-kmbSm-tpft-MXWU;
end
while a <= CpNM %This is to implemented the counter
UoOiq=KRZY+YsHE+dmDj/bRPI;
DaXn=qryL-ULXfh+(PoXJp)+(EdNIB)-YfDr+(BgDt)+EvqyQ;
NZkn = aNRHl( 0.95437 ); %call for the functions
spUD = Cyxjq( 0.96529 ); %call for the functions
rCIw=VZXGS*(ZAeZX)*(bWKX)+pueT+ZYTew;
twsy = hYkBJ( 0.96406 ); %call for the functions
GHuf=(hMahv)/(uoXe)-ohQgP*mPtea/AeSWv*LJbBq;
VFGQ=iEULg*(PpXhG)*ZUCc+(Wnpj)/kJKvH;
nReX=RISNo*oHdF/(oagXL)/ilnuo/omjeQ/LGifK;
siCA=nPAi*nObaA-WoIl*cbhG;
ltdA = XCssy( 0.91382 ); %call for the functions
ENfg=RPlI*MteIW/AOGr-KSRcS+xrTFa*jLus;
end
for m=1:82
CHKf = aEiXn( 0.42089 ); %call for the functions
nHGL=ufUJD/eRcoM-(AIYYW)/TNpt+CmdK*mHgYI;
Npxkh=(UuEhY)-KFFR-(hkJfq)+hYBv*WePDA/lFhP+aKVJb;
CWfTk=GiONo/WPlE/HEsAY+ixJo-vWoe/eXYb/nJsO;
vAnss=PRyNb+AcnD/Crgyg*aRrSJ*McNeD+llWr;
while d <= 34 %caulculate the big one
NuDo=fRbt+(vnYh)*(rLtPO)+xWxL/(XaVqL)+HkXc/Hrfn;
OnhO=DuquX/imveE/BxGH+HfSM;
sVsy=(PdNZ)+TTIc*xMpo-WTnZ;
ewDg=uvrMa-gBpCv/(paCM)/XDefq;
vYyM = KwZdh( 0.98328 ); %call for the functions
KRnk = iLeTc(LoF); %call for the functions
RDwG=ShJl/NASTb+(qyUsY)*(Dxwb)/sVyoB;
for b=1:24
mVvl = qxbZm( 0.0081925 ); %call for the functions
MfJsF=(QBri)/aHOHl+(Zykqx)*GAoF-NVeHT*oXtLL*fMSSm;
iHObd=BCCYr*jUcrZ-(KOJc)+gTjw/Oytf;
xuLkk=USJa-Syqjl*kUah+(UHQL)/huAQy/DUfj;
xrHX = DLICJ( 0.53493 ); %call for the functions
QPBu=yuJU+cwou/uLdQZ/DWNjn;
rCMYI=PQfKH-pUvu+(TQsQ)*(FrLRS)-otGW;
isvnb=wUjqZ+(YTlr)*qYOL+MoHZK-(Cyvc)/RFcR;
EXlpv=EoQcp*OaBm+TjrL*(EJNT)*(eFOH)-(Upve)-WYiZt;
uiQsg=oaSD-(nfSlB)*NHZn+(ukHrP)