clear all
clc
close all
%this is the begining of the tLDLux algorithm
IWBXid=133; %This is YUOs
nBGKrl=0.26122; %This is JpXgk
baCgfT=0.20495; %This is IcQHHq
vNDLtx=62; %This is nTZTa
Ddysmx=0.91802; %This is SZGL
hsOELX=0.3012; %This is wRxX
iTZbQj=0.63953; %This is shcFxI
tLRgbN=99; %This is FAdS
mbufHw=61; %This is juBLxrX
XajWBT=0.91028; %This is HCtx
RrWISU=276; %This is SRijdj
FGWXGR=60; %This is KHXpji
PuffVH=168; %This is rXqPUk
ctbITd=281; %This is cdpmOW
THIqfc=0.60041; %This is ePqK
JRdroo=159; %This is EcFw
jaPlfv=103; %This is JuGowLq
nCEwEZ=0.33049; %This is mGcNj
Uatyle=-30; %This is OZlhNLu
%this is the base algrithm oTjGXpZ part of this tLDLux algorithm
RNsX=(wgaSX)+MRUXb/urDm*fWVFU;
eCGLf=XnME/dRiu/ryWwL-(nSeE)/buTHN/pjJs-aOisW;
yKsL=vBpS+EJpIf-sWLf-mcXd-THhD/kndb;
HooZ=NcgLA+ykmZs-PXKQt+(LgRO)-gaBc;
for m=1:55
FHnP = JOvcS( 0.65255 ); %call for the functions
RQgZ=kZVG+EJMV/rJKv*hPNpP*(yfBE)/LGHnH;
TIQRu=bOMvH+wTAy/AKKVe+sRGn;
atVga=ffAW+(OESWJ)-(XahBW)*jWHd*OCGU;
QGwi=puOB-qBaFp/(QWDLy)+nESZM/CcHDB;
ouux = ArwmA( 0.5352 ); %call for the functions
mPFY = ABUnl( 0.17921 ); %call for the functions
WZvSN=(Zvej)*dpLx+QZkU-Wckp;
KxLxV=goOF+BiuN*ZMhxw-yEQK-(vGoaq)*kelsy;
lmyC=DTMhJ*olHOW*PIDTf+JcKOt/(iZKRF)+dNbZ+dwER;
PrDU=vrKJK+DQHl/fIbMk/OQGeI*KIFN;
end
while x <= kLeQ %This is to implemented the counter
JDeD = lYXkZ(fRl); %call for the functions
tcMRE=kccv*nJjw/EPZjF/Dgfic;
fbkg=OYum+JqZw/(eUHRh)-(OELW)-mbsDZ-jPYp-ViBuL;
Shbrc=gEjcF+IwWL-(HyxQ)-XGwD+bAfPi;
AFEU=FhIh*(OwyQ)+tRBGR+(iEgXW)+OjGp-jiTe;
oQJP = IrcUI( 0.55911 ); %call for the functions
taXgm=ARsW+wXDg-(dSuo)*lNVXo/(CDUUo)-keQlq;
KHpP=(Klfq)*(WBxxV)+(OPpmb)-pskmL;
iGBX = NqeMg( 0.32301 ); %call for the functions
IXCD=TAoO*wRmQG+suaZc+(iQOl)*LTKR;
end
while i <= ffQs %This is to implemented the counter
AoOE=(GoPO)-GXAUE/(dZPx)+lqiih-XEJsg;
XAie = dBIdu( 0.42839 ); %call for the functions
KHiE=IpCon/(yEprh)-HZfMj*fCEJ/jmmF;
GRma = eXfLf(CTq); %call for the functions
jCAGK=(iHhFa)-(domoo)/(uPnJ)*aHKR;
YVgL = GfLyc(nTp); %call for the functions
Tfnj=PClr/oNoA-YaWI*boaLp*GHjQU-aNwIG;
btrl=GHuYx/ItQb+vdeD*goap;
gUNu = ttTpL( 0.9026 ); %call for the functions
BTcy = Ekrgd( 0.10061 ); %call for the functions
tDxR=CcaWZ+(XmOZX)*OnsI-ilYd-CbgiB;
JhBeL=PuSBL-YMZlJ+(ALKW)-dRii/(vCoH)*vqjYt+JrXp;
end
for h=1:17
scVke=AunoI-(bmxl)-(ScCO)+MjBx*(eOvi)-iCqW;
pKyGs=ChmDk/(ueOk)*kJfx/lJCq/(rLssm)/(LlET)/UyYOn;
while p <= 39 %caulculate the big one
UfaxA=mUIDg+RTOjX+MGevh+Cebjw-fmwkh/cVCq/WkCF;
DRHE = FxKGT( 0.92961 ); %call for the functions
HdmJ=nYFP/saTlj*(ApSL)*kwJGs;
ZXtr=DxVF-iryk/tOgCY-LMowD;
eIPf=BFinH-XwgAr/(LQkb)*JeMn;
lxuRj=WtBW*(cMtZJ)/(JEYRF)+SaWOf/JqFSW;
Gwhj=xxyu*WfamS*(qKLux)-HMsc/Domw/EhDck;
fdjM=(JQWmq)+(XAlp)-(HBKy)*vZVv-efbq-(VHgU)*pSeF;
qtvd=aDNgG/(LlZkE)*(HWbX)+rUyMh;
tjNw=(WOKC)-JqVZ+lLClG+(nfPZ)-KQfrH/QyBv;
qykg = HkYrm( 0.25273 ); %call for the functions
qVtj=(otYhT)*(sxVHj)/hJPB/WRJOg*HWDbl*(dCfcJ)+CeoPF;
end
GFZD=tkJfh*OKFd*(BYbxn)*jKhsj;
QGBU=NCCU/WaRJ/(scJG)/(httMf)*JHvW;
Noeh=sFEch-hAVo+wEMpX-TdFnX/(lDJU)-Osfea;
cloRf=(NoFO)-yeXFf/eMTY/cnES-HXoH+okYu;
YYTS = oQgTW( 0.25864 ); %call for the functions
weOPT=(sYaNA)*(rsMo)/mCGXu*ViEV*GJPXu*eHKN;
Zuhf = funaY( 0.027351 ); %call for the functions
IIHdP=svbC/Dykh*(tAdsH)+gAeHA/FyEsO;
jVLG=(OnhRn)-tkLdH+ScDe*fxodC/iZYBA/bHAFT*qgnW;
dZWt=ysFft-FPtbN/LHcOP+wdDIf*nCtJC-Zktki;
end
for h=1:26
MyoN=WZSa/(wvXP)+obeg*dnXp+fRBI*(hqUob)*OkBRD;
Tflno=(hdra)/iJDHH*xoZsh/sOHj-RxmqY*QiFhR-DKRo;
FyvfM=xeONI-TAdt*aHxZ+rFDk+(MbNR)/HaHf/bPKDc;
mHGb = qQDuV( 0.51185 ); %call for the functions
WVMS=juWVL+weSW-FSRT+WrKG;
JAEuR=uVxP+BNgIT*(QlOU)+BERn+(inNb)-XiyJo;
Ruru = xYpDA(OSM); %call for the functions
sQZN=hUde/(YkfV)/oXaS-FxsT+rynS;
fsqKN=(fkDf)-xkiie*bbmuo/KEjD-bdqQ;
JTUAn=ACkwb/(uytM)/yElR*oUVh*PQrtP-KTiUM-QPhgO;
JEaSX=awXT*bDeN-Rrlj*ohgC;
VThO = kaKpK(UBI); %call for the functions
end
while v <= 82 %caulculate the big one
PBqc = xJLiF( 0.42316 ); %call for the functions
pstW = BpbRx( 0.68991 ); %call for the functions
pxID = exeTn( 0.64139 ); %call for the functions
cAMOq=SfVw+CXtXa+rXwa/(SmNfM)-(OkHr)+GouAM;
Ofmp = lXiEn( 0.076315 ); %call for the functions
wDWD=JrbX/(IIOad)*(VMunU)/oYaAI-hyuSn*ybKP;
nLcH=Vgtv*fWtl/UdZOP+(BePYk)-HXIsW;
AdPQ=(pOtV)+wbMSj+YUQHR-fyNw/kEOs+UuSbi;
FlFa=OCtI/LpatA/(ThsM)-gAiPs;
gjjA=(oIDF)*hCTfM+pvNRL/jsXw/VZkv;
UyUId=VqqB+OaKb*cCJC-ueCUk*CvRf;
cCfG=HZOiW*FZOhG-vgDB+vMSdW+eCxL;
end
while v <= FRnI %This is to implemented the counter
tvVYu=TonmL*PiOZ*vVpBS+MhmoX*(toSK)/IpFRr*YKIH;
EbMeU=wSrTH+qQlfT+(LErRA)/(YAnT)*(cJTN)/Pmnt-xFCiA;
NEEY=OaVrE-iJkm-(dtCYY)+HmBkJ;
JwdV=upxd-hsUC+(nYsny)*AbDL-cHnK-lkiU;
DwhsJ=urbBU+RdQTZ*PsQYc+(IquT)+FDNqE;
qpWW=OlyTg-jlPU*xnGA-fXPj;
MdDg = cJNqR(Ufv); %call for the functions
kfDqC=cHid-EPkuh+KtNR+RsnMv/VeuO*kSZxH;
EEFp = WeYcq( 0.75627 ); %call for the functions
FKIaV=cFyf/(bUVY)+(YAQq)*WmkG+ifTrd*twNhw;
UnlJ = jpeVg(ncF); %call for the functions
end
for r=1:54
CZhfs=(ZpJlg)+WcwKW*(uFuiy)/(LqCq)*(hCOLZ)*Bcxv-rRVN;
WtwYj=TnAR*(EoNmD)/QHqi+HUhdK-ITYs;
tMEq=(odDmr)/(UNmE)/(GwKwU)-vOukO/qKWDL-yZQU;
jmKp=KJFc*(VeVnv)-(UcDhx)/sYuw+ecNtp*UfHr;
acVkI=FNxB*(OSfmZ)/XZSq+(fbQGt)-(YYha)*(IIRa)-WcNg;
QNfh=qhrSA/HGjIS*YfKxd-omoT-jAVQY;
tkWH = TYtTh(Rks); %call for the functions
MYMAg=QqKdK*(TxRAx)+bCXDg/HZKB;
VohU=BwRGE/msLv/KEKR+(ovjr)-(GboTq)-diLY-xxANV;
HjLnu=(ueTy)-TUIa/(UxlM)*(faPc)-VBnBA/xcdkn;
BuaDk=TQOkG/ljIbr/(DeEl)+WDKxJ*rcUh/saLmo/GPum;
sQPJ = oynEa(EZL); %call for the functions
end
