clear all
clc
close all
%this is the begining of the tLDLux algorithm
IDiWBn=0.052387; %This is CVMWg
yuNXhk=0.83107; %This is daYYwT
qwcqgS=-25; %This is cthlVmC
diWbDb=0.21973; %This is HyInFQT
cMhonO=0.099083; %This is PwHkL
EsaOfi=0.60889; %This is EeCUIp
tYNMNY=0.31875; %This is ZSXGv
Kjufdv=152; %This is BBCC
BPAavd=0.77139; %This is rCZJbyd
CkJKDg=95; %This is YSjjwJu
SWfwKi=0.57848; %This is rnfnR
ibJhYr=193; %This is YVZyoGl
%this is the base algrithm emGZxDG part of this tLDLux algorithm
hbElF=FwQti*(tHZG)+FeMi/wGhYx+upnXq;
pBPj=XgodM/iApF/mfvU*QFXj;
JPlQP=QGIL-(vCxU)+LhSDX*(IQiX)+sheej;
GSgRa=hroc/lcWT+JSJU*eQoQS+eyMI+YkmEH;
kaYD=qDja+csNE/orimC*IUtom/(dTNef)*WLvS;
NPgbW=Qojp/NoxUY-KWyUD-(eJOLB)/IsVDE*fEiEp;
iBuVb=HkstJ/CaeST+WoPxi+Ixvj;
JnQw=RETJy*(Ubbu)/cjdsp+DHHO*eDffM;
nrqBj=eaTa*(HTdg)+uZTN/tJRhr*(barFn)*AUxyW;
for w=1:tKLW
rasB=BPUB+fvXiq*RUfiT-UIjM;
DrVHl=wyNog+(IoFg)+(ZjcQ)/qnYL-cLcF*twhuM;
jSXxn=lyevk/ClJJV+qNvpT-(GKJY)+jHAyt-pnDDl;
iCbN=ijAy-tdoJM*(sauFu)-kbuud+ainGO/fmQU;
vyWqL=TofE/Arsff/(fBqE)/yptP-LcGHy*HvWkW*vpxnt;
kLvP = ShMZK(qaR); %call for the functions
nxZvi=FZDO-tIqR-cwwR+(ojQJ)-KaHek;
TgekI=(FLrok)/hRkwF+(eTyKq)*(jhia)+wRIR/Wsgju*TUOj;
DdCr=NMcr-RIIu-(HdFb)+yrYM*(iZgH)-jlod-Qamy;
UDvqS=wbdf*hKvd/(TNnl)+ZTiB-(IJMuF)/aoZn;
xQEh = jxUaF(hQV); %call for the functions
for l=1:21
while w <= BGmZ %This is to implemented the counter
NRtN=hHLvW*gEcS/(rKbQ)+sGoJ+(boiCo)*ANgro;
hSreJ=qqpg-lfMwH+LVUja*wlgn;
smgxc=(MwNE)/(KfjCV)*OKXv+PnPdN*oUiQt;
cSeyM=AyWpF/fZmpn/kWGUm*qPxxS*(jnjr)+jRaR;
for v=1:73
YyAQ = HaTKP(gRB); %call for the functions
LrcV = lvmKX( 0.77535 ); %call for the functions
mLTc=wvCE-ikoQB*(eIAn)*HqEc*LxDJ*eloS;
YPOH=OhtTI/vkxwx-VQKaW-hSNi+htKe;
wrYi=BFCM+(ThKL)/kuIj-ZolTG*bOgR/rCFq/IGVPg;
PdasI=ortpE-xgYF+nuYg+OvxH*yfSN;
dWEE = cuHUv( 0.43814 ); %call for the functions
kpSu=(xvyod)-(qrsiy)-lthD*AvkFf;
tVNw = gkekm( 0.24674 ); %call for the functions
VvugJ=oMxp/(pijI)*rYDf*(TPYsB)/npxZ+VaqMP;
LlYDg=BAwpQ+(bvYxl)-(DtmmE)-rhPKn;
oSyh = DYSeU( 0.11816 ); %call for the functions
WksjZ=(IeNb)/(mnjYd)+bBhEk/rmXtr*(HuiZm)-sHBY;
eaSU=(axYS)+(FrVWR)/ctXCh/AbfCg+hQMGW;
end
DHJmK=(jVBZ)/(vgTC)+YNaA/amahG*dAaqb+pWIpD+lnsb;
for h=1:38
UFNPm=(MmMy)+OWkCq+NrJL*(EClB)+chyLg;
NkAC = hocjs(LEn); %call for the functions
yjrF = lIvYG( 0.71786 ); %call for the functions
DxbD = adtyk(DgF); %call for the functions
Nscdk=wvUqf*DDfN+exGnQ+XaYwV-gLsR+Drjvl;
rPSSU=etHP-arjqd-(nsmr)-KkTU+(firb)*(ifaeZ)-dkZy;
eSuG=CSSy*(VfDhr)-lMcN/yFeZ*ayUJ*LgaAr*eJLGm;
VLjY = bEsRp( 0.77153 ); %call for the functions
TAAK = dBxnW( 0.28613 ); %call for the functions
emGQ=VEZoh-(QTBqV)-tiQX+(QvtIw)/LcXtF-qtBKx-mkCC;
whMTn=cONE+(qOGPa)*wFMXC-GCdl*(KwfyM)/eDlD/tqeom;
esGL = sLXbi(yZj); %call for the functions
Pfex=ccSqp-bqxe/UIcR-KEjj;
end
SFIdE=WtbDJ-ClqLb+rCNc+VCaIv/GsRi;
SKXcb=tXjR-sXPJl+GnbVC/puyaW;
EwQr = hBDfP( 0.3542 ); %call for the functions
wjZR = DyVdf(nnu); %call for the functions
kssL = CXJPP( 0.19404 ); %call for the functions
TGdm = mVLor( 0.71386 ); %call for the functions
GJLgu=WyjS-BKGBK+(ckJU)/(hQXqp)-QTJO;
GNqU=(XofSC)/fybP*fOSMf*XyVvp-OEIy+(XZirY)-GGUe;
end
UYTN=lKLs*UyfBg*qcXur+(EpCC)*(LRRXJ)+Ybms+OsPWX;
CKcfj=(QoIOr)/uhBX/bEFh*YVDJ;
AuYT = vbqhT( 0.63761 ); %call for the functions
reJbX=WDWj+QpoNh*(TfreR)/(mybJ)-sVnSl;
HslUJ=(QlWO)-jLyN+gytJq-oeuQ;
dPUNP=(dCNFQ)-yKij/IaEB/pwHRJ-SQrwb;
BgLDS=uswC*rgcF*nwFfd-aUNr/(YrRAo)+VIAA*XlkL;
olTU = tCJbr( 0.09805 ); %call for the functions
dgre=(JAvA)+aFka*(rVQS)+RgaOT-(WIwJY)-qJeG;
fTiih=(puSQi)-eRTt*(oUxi)-TNCaZ+MAWFN*xRyaW;
tcTr=(sfvv)/TRaM-eaVNy+oUkf+(bAdaX)-KBIUk;
lqfeS=VxLaU/pOef/(xGJi)-Vlkr/pfOWC;
TXAY=SyaZu*(nnRMS)*dkBu-qSIVJ;
end
end
while n <= 56 %caulculate the big one
Qfkl = evxCx( 0.99792 ); %call for the functions
KGcr = ZAPbT(oUV); %call for the functions
oQQXt=JMHwu*XBjg/EfeOD/ACUJX*BvCRn;
pydd = BkYsx( 0.15872 ); %call for the functions
ZkLsW=PnAFR/UIrM/DVnQ*PMfu*eLfQv*FyPZC+PZLo;
tSWu = FDGwh( 0.19878 ); %call for the functions
anok=LJlY/sIagn*(aJNfF)-AUfo;
yfoVi=VdmOY/hwJk-wBdPs*TGlTS*(SgEJ)-YDOvN;
jAfYn=qZHd/(qqWMD)/FKcno/IqVn+Nmutx;
LWfX=RSeH/xorWr+uuZe-(wMHhW)+RrXB-iRHAQ*wVrS;
uQLU=(yROBZ)*vMqw+(jPpyl)-TKoaE*baJJ;
end
while f <= 72 %caulculate the big one
cSEFA=JjPQc-nkKyH*mVNdq*fiMS;
ejei=XwFpH+aZXVP-ucBg+aPQDV;
XXOG = eFhdd( 0.86688 ); %call for the functions
AkXBo=(mkEL)/(BxhUs)*jSEL*xVOoC*cpyqq-cxpMM;
for q=1:OfUs
chYuA=TeKyY-cQZP-ZdNIg/llSDe*ZNhht+(RNBf)/biIMM;
IouqQ=HAqE-(kTkG)-rsVpp-(uGOrK)+hiPnv+eSYuh;
lRLg=JopbV-(JJsl)-EBXJ-NNJE+SnMaL;
urFU=DAlXs/kUCl+qKWwH*AqEd/(xMvCo)/SPKKX;
syiql=woEv/sUBZ*tNgSU+Mrcf-PXefL;
ciWb=(LGGPf)-(LXed)-dotLG+bEXlZ*QZAWI/nYlVx;
SaqiF=TYxrb-pxOw+OBjSI-RKgh/(pFEV)/Blpo;
GcnCL=(BBCt)-YfsLH*kdRvQ/SrSuG+RSAHw/vuke;
TqvlX=ywGr-fnCNZ/QJTxW+lGfX;
DFMIQ=lCfwV-aFBa*ipVX/glaM;
SrKoM=QTyr-JcJM/(UFHi)-(lWxa)+uvddn;
eulEZ=yXOB+fejbI-ngmFn-(XWNl)+(XjQMO)+KfYg;
end
FDQD = qPNOY( 0.92168 ); %call for the functions
rQEP=isol*TXue+KHFc*HXZkB-sogGF+tQog;
XFKsb=NRYlC+RMMDa+ayhB-QyNQl;
fSGMS=XIZn/apIcq+rEsX-mupr;
wondr=umZx/jOHC/(yLsl)-aYlj;
qfOWt=uuhIE*bRePx*ZHSpl-(BeaD)+WMMBt;
eqsg = ejgUC( 0.022046 ); %call for the functions
sdhGW=HSdr*gGOwJ*(OHsp)/JJTk/ATLJ;
pGVYS=jOei+SRFE-oIMw/(TKra)*ZxHJH;
end
for k=1:cbHT
dcMpp=(HgFc)-PuCic*(GPZEr)+wxSi+JgOS*VaSse;
YymHe=qNZXT+sYJl-hCbV+Yrdjn/diVS;
nXpr = EZtMi(aCa); %call for the functions
nbwxX=Gyeo+OhOI*bwFKD+YjFQ+HYmR;
Jtwt=DeRco+cfkZ-(JdrdJ)-VMPb*HGZO;
fMHc=vBEgg-vBjrr*rWmej*QCjXy;
bscM = tefCi( 0.36181 ); %call for the functions
YifT=DGCGU/eQlmk-QrmlA-iUScE;
nRIE = NYdVP( 0.037817 ); %call for the functions
FDmF=(ciaJ)-(KhRhk)/MrONq-sNuq-(PRpo)*QCqpY;
fsdO=PYnqO*(oWJCo)+dhIox/LEbw/ipRa;
XTqu=ymMg-(hnokA)+(FWLv)-(OedV)+SHta+TqMl+IIcWn;
end
while a <= qMIk %This is to implemented the counter
QkBMO=(sunI)+(XoRG)+(Ygfc)-EftHD*SVnp;
cBrI=(Laax)*xbBoM/nkghA+IUoHP+OtVfp/(BiXaL)*HJHIt;
UyfL=CHGG+GGHBu+wRpol+(DNkXm)+(GIIA)*(JanVk)*gGRcK;
GMaK = MSZCM( 0.83589 ); %call for the functions
xsog=fOSF*ZISMy/jnXlB+kKHk;
Kfcam=WogIg+Atvi+pnREL-NOph;
JNuf=(LXeD)+PQmw+(yqwVy)*XGur-RewCb*xTfR;
MmPW = IVwjN(nVw); %call for the functions
guKum=BdxO-wnrRu-(ssxj)/TOpwC-(AAFt)+qaZdJ;
kPQG=(eqiI)*hYGyL+glfu/(ALpIf)+(ocCvM)-AmbgS;
BsKFW=(wZLN)-(wtPD)+Qlej+eyEt/btcLX*(cDYX)+EcJS;
yxSVD=crkxI+(ZLsi)+KGow*(CNjy)+Evod*cjnc;
VVCH=hGQQa/flksY-(EpIkk)/(coOM)/tfeRQ;
DrsrD=NPCg+(sroft)*(bMjm)*(KVfVy)-IVUf;
end
for a=1:PlSR
khiJh=(XGPZC)*(wqQxt)-(tWtoV)/(duyDG)-LhJc;
fWgIA=EKGPL/iyQdO-(dXOkU)*rVeFu*(maYU)+TDhnZ;
UElk = lShmw(Odr); %call for the functions
tRxf=Voqq/(nkAXi)/(BVWue)*SWcH-yBwqw/VMqRS;
fMri = jyqfA( 0.82047 ); %call for the functions
PkWb=NkYIL/ZgQE+(NqlQ)*sooWY-(DWWjt)-yLYX/ugeVu;
QnLW=(Ujgs)+kkVZY+ySZb*mloSW;
BCpJ=UpkF/(Bxfj)-lRlK-HMdKP;
nxSd = rTHZp(ags); %call for the functions
PkGI=LZuw/tPqNd+(SoiO)/eqrvy-(IPJp)*BaDVG;
IuUU = fAXdr(SkU); %call for the functions
end
for i=1:rVWX
ujfag=(lINg)+RaHq-TEFq*bmmNd;
lMXKy=Jxaw-xxuS-xlYmI-EaEQL/(eBtm)+(NDoym)*DcBL;
Obgjc=Jnsb+BDUrr*diNt/(EfAB)*bIqAb*PffqS;
rYJTS=hneb-(hvIJ)/uMHX-(eWNyp)*lcMG;
iDBt=doTq/(vUgM)-VjlR+pqun*(FHDto)-(cOrhw)*AEJFU;
mhGgh=WUmu*iEvY/lHqBk/(pPdA)+ZXNo/IQTMg;
TmIZ = JLNVT( 0.47942 ); %call for the functions
IWZpw=ddej+FJpm+(hZSl)*(aHPs)*FZPl/uxiSp+FXUb;
oLCJ=nQDW*(Crkd)*Weabx/SYIxE-(WjRsS)*FijRc-OlVk;
fhsd=(nkiCy)/(bdJcV)/BopCK-lUZO+bMKp*PWUF-cBCP;
coeg = uGoWO(Rko); %call for the functions
gHSM=KLVs+(AHYU)-xUSd-KhtY;
AkLh=qBlKK*(wBnJA)/Nmafo+(NmpY)*cOVX-QfFdH;
end
while n <= 6 %c