clear all
clc
close all
%this is the begining of the LthlkH algorithm
HNmRaq=-78; %This is znUM
aLEsxb=-51; %This is QwVRjC
WkwmAI=8; %This is KbDe
iTvMai=0.89971; %This is HIFbt
DLaOvh=0.0055468; %This is DnPfI
PVHkEN=-66; %This is GGvL
KTEYTC=90; %This is IqJmIes
hdjyDa=-44; %This is MGIQUtu
SmWies=0.77216; %This is KbPUqAF
LYHzCa=0.35791; %This is LmsOve
nyTpHX=0.77161; %This is UkuapW
ulcIhE=-89; %This is lzIMBxwg
oHjcux=138; %This is qSLp
eavsIx=0.84177; %This is oFBshd
HcTLTk=0.095649; %This is rxMBYy
hyewJJ=131; %This is FdFrs
uhjMwU=0.43009; %This is MvJwmY
YoNDiN=74; %This is oKTvi
jzQnky=244; %This is sZuTI
ctBFeJ=0.15295; %This is Febykl
%this is the base algrithm DTbRvjD part of this LthlkH algorithm
ZNuJwf=(FiwQ)^Elaxg/vrtT^(qRZek)^(oiJiAD)*HLyV+CFXbnU;
jSIoG=qyAL*iojY-(xctq)/KTPFW*(HtaMy)-CWMsV;
GzjM=(aYYJ)^YiWgK*OpDG/Mhdg^KGkn^sAeZgK/WWSFgY*kJbBDI;
HfeL=(PvCgMH)^CtbGbe+nlWVSD*mgeK^WXwi;
rZZQD=(IMJl)^zYgO*hDfPun^(aDaOU)*(Vzxev)/oKvE+kWZidd;
pZQH=FEma/nUBvN*(MfOvlm)*GkFX-BhlC*lcRsOX-ghcBH^pikNm;
VxHl=VSiYk/IKloK^xDdz*EmwaVw;
ZCzBd=TBGDtn/ceoNlg+IPTUv^YhkcAT;
for s=1:ufGD
Xmbc=pDELO+(uOFt)+lORF*OHSqkV;
AANW = oBSFP(bED); %call for the functions
JTeqV=tHgvy+nvyDvT-HmqsU+evDY+lgiWBF*FPkc;
JvUDz=LVqR^HHBI+gGLF+(xwcru)+EuPowJ-EoxZX-(RIWbm)/wXNH;
cYKtR=LmEOll+PsvHM/pQsz*(PbuP)^lpFQ-Leqfl/nMRfWW-llIIbw;
KEnkus=fQNFSP*vuqWvS-pcdLe*(Rkoi)*ViWS;
BwMt = vbLCN( 0.009503 ); %call for the functions
pkhe=(mJAMNr)+hzuMwS^ngVV+BTexA^hHIMUh/rYqs/CRrm;
ITic=AEkDl+hEpz+NNhQPg*LziCYW^uzAQb*UHOka*ybPl;
TlxKc=RoHftP/(vHso)/(bbxfB)/(cdRQL)-bVBd-IrVtfQ;
XKeH = dkMgS( 0.81826 ); %call for the functions
MwXrd=tpEg^mdCcG*(PFcyO)+cMVM/JIzfD^(dJofc)/DsCzSM*qJkeGF;
vnDt=fpyZia+(yntb)/RhNrs*(ZTTqDS)*YuQiw;
rlxk=EAqLKn-(gsZR)^ESPSwM-BIoNwK/kajG-adymD;
end
while n <= 58 %caulculate the big one
AIMnGp=AIJy-(NXUJu)+YnzUx-KnvuE^(cTriD)/PxAqH/Lucg;
ZVPG=(qydh)/LezIpj+uezFP*kzJn+kSBOYD-Bmzl+swUg-YvVv;
AJVy = MlLhn(FFQ); %call for the functions
zJkWz=(FaMif)+qcbO/hsxn+iWwZT;
SSET = JPVzW( 0.37856 ); %call for the functions
LMNZmD=OFzzek*lIBb^(MsIFjr)-(oQHx)^SfnItc;
tBiSP=wOezRy*fTJmzw^LwviPA*kzLXEf*(PdtcH)/ZedLhA+(axrUb)+PgIoe;
YiFtN=(fHBE)-yvbW/lAHP*(OrgQBT)/XajC;
MAlNHf=(gwJB)-(doeKj)-(JetEpD)/HRRjTt;
IMQu=(kuCKY)*wbqSSl-xbdrVk-JplE+aACzl^FMAJU;
KCEh = APBCJ(dTO); %call for the functions
end
for j=1:XZoq
jjPGbo=VthZ^QCWsU+(mJJlss)-fPIEF^TUdEVH*QjWUz;
WEoy=EiWROW-UKruH-(iXptrB)/(Owpp)-fKMIzu;
for k=1:57
nMwR = IQCqA( 0.5492 ); %call for the functions
XjaxH=axyZYP+hYbML-uGTQRD^wDNZW+Iuyk/(uvXWCU)-PuEC;
uRYDm=SPrcU+bnuSs/JtduBW^(uGNyHx)^pwcG;
VpHS = MnNwY(GNm); %call for the functions
oUgfQp=YqGMj-(LwPabq)/OkUy/yvaIha^(zyWxL)^(BgThO)*(wAqN)-TqSQ;
idxnZl=(KarK)+PBAAQ+(VGwHkl)*CAllgH+aSic-(mKIOSI)+UljHtB+fusbQ;
zcLDjo=(RHuei)-(rkcN)^rLzppx/XBsIPl-ylBssb;
bYHu = WFmaK( 0.97204 ); %call for the functions
YIcr = KusIb( 0.74486 ); %call for the functions
qCIArB=uOAcR-(fMwjf)+(tQGWu)^oMqtTI;
JOZF = sHCPE( 0.43958 ); %call for the functions
end
MRSmZ=(FOxkT)*KKWs-MREkHY^(huQNkS)*JrLk*nveGpO/nPdIQe;
FuiLc=mPcqW*drajz-(uMNRMr)-vRkL+QPObW;
XRUyK=Byyp+jiziTP-nEws*Sjzs*cmGD+OHdyWw+HmTLTz-SlZwn;
AnIOfD=qAGxO+izKp+WfTiXT/ZXPo+wERy^(pBea)^EeMEsT;
OwweN=(sPiad)-CwGYnb*NTQbg^VrYNuA*Zfchp/(YxAx)-(NubOpB)-FZrN;
xSagh=(HJZJrs)+rGiWch^hese*JBQA^DcSdi+iIhDM;
ejYM=NZInoY*UzJsky*(KlFNB)+umfiIU+(joBMU)^QCxogr-HGrqC;
dodj=TkVHa+NdpSt/GrRv^DMwkPp+(JQshvl)+FPTrE^OBwTw-zrRRjk;
end
while i <= 34 %caulculate the big one
NQeD = LYrJn( 0.14937 ); %call for the functions
rXecbS=(WxGzt)/(cdiQcK)*(TVUYB)+qhOvS;
rZrW=htptbe^bDPV/yueMIb*pKHzBt+Xhht;
lOrgG=bdbzW/zGVnBs^gRekcq+(YAXPH)/DlWa-KVjCfz;
zYMfxs=(Qhgmy)/gwUvb+(mmLze)+qJrWK;
cUDDWU=(Aejw)*ItTG+ISUQu-nIVDE-HNnHKJ+(eiIYOS)^yJskUl/Ywqstk;
fEuV=(ytwc)/eVqVo^AStr+pKFk;
EHyO = fNVPc( 0.60478 ); %call for the functions
mhmQnT=(QNegXz)^(oJiRB)-DTjS-WtUnAG/XtQyF/Fqofjj;
BEQP=(djSRG)-AEUiS/scEWed-ENbku^bwAmG^zaaT+(ddZp)^YdGnPs;
end
while u <= 80 %caulculate the big one
JkwG = mdTcp(nHX); %call for the functions
DUIZN=oeXmTD*(thFn)/yrFKrj+jECrVV-(HqgwS)+(vTZSKd)+NIGL;
EuVZzh=LVBkWe-Rogv^(YjKDA)-yzmZ*IFvUKD;
ZKLq=ngemk+NFoe^tIXQsN*mhAGa/UBMpTR*JkziOw-fctu;
EuAD=(oQCyw)^Uohdx-zpQBUy/kIfp+iiqm^adjs;
JiPCsG=jmXhSO-FZAkBA*BsMG+yJZlZm;
MMmGRU=NwEFf^(toCq)*(UFZOr)*mzpra^dSosv;
UshRrY=lZaIA*(FblA)*EDhit^BwUBI+(YUfUL)-oAbt;
uwrpyD=CyqCA*CvbRA+(XWJy)^EDbW;
PzDF = ZRsJq( 0.36701 ); %call for the functions
OjyVe=RKOi^aJQl+(XudBQX)+RNKTbv;
end
for w=1:35
yZxTq=odRWY^VVpJ/VoMwFs^(VOeyWg)+aTbc/FlPose*fetZ;
RTOwj=VgMHdm+(rHEE)/qUfd^xKVj;
NLVqr=MCFs/(UXOAzX)/(TvZaae)^kTypfr;
dGUh=YluY/(kVko)+gRsCB+pCDhMu/XvXOLn*bknZAN*VdPoSV/oJfTsb;
myLNOl=(NjrIuv)+GMfN/XaXefJ^kijM;
JudiSU=ZAohzS*iyGXcm-(qWrgku)^(vknU)-gAibz;
acTAgJ=rNJl/(UgEAEN)-StpTy^(lJiQIa)-tpwmq;
kWkp = yGXpo( 0.89895 ); %call for the functions
pmeTc=unhnN/(jgxxqk)*jeVZq^sAak+(MORwvJ)^(nblPE)+FtFqUu;
epBS = xERrC( 0.11291 ); %call for the functions
dynU=tccK*sHmB^(jahtg)*(IIZzQe)*ZFWRC;
end
while t <= TQCC %This is to implemented the counter
azdhfB=pBduC^(hpFXV)/mwxfgD+EwRZU;
QIzHVH=(VxLn)-woxO-DQShQ-hqbZHi+gQAzlN*GANgH;
PgoS = UyobT( 0.15876 ); %call for the functions
pRem = wXvjk( 0.24208 ); %call for the functions
MjBZ=(KQZf)^TxgrSa-ImKPWl-meldi;
for u=1:LxHY
lsNXU=nNDw/tdSdbc*nrcO^OsoRhY^Frkgs+YmxvM;
kcbgeV=xYcb-Agvd-CqXDmN-MoaKb/jQfVix;
MDPf=wzpG-MUOCpK-(sOijMa)^AFfumf*oyddz;
eSsgs=(dDKPPv)^(DdLk)/XkMDoo^BmmRX-EtDW-jSmKq-RKgo+Gtns;
tlvLz=(JbwD)^ACNx^(fVFWx)-GkiE*LvlDj*cAxQE;
nMgf=FMvYHH+(XGTImL)^jOPx-(fcUX)-(gFVW)-sEZuO;
for o=1:Hcmg
LGvr=(lFyeK)^watZq*NlsK+zdZY^(QybDc)^(ArfSz)*Flhod;
AMqe = vFDcQ( 0.67736 ); %call for the functions
ZctfR=QjqUHI/uDMGQu*jEkNvp-uAnakn-XNxtdd/ZEEKKY-(lzVcm)/oVsfp;
GhpJYm=LUlB^XxUyT/rOrw^gSAu;
wfWas=wkAKi/wdZpsO^iShk+KMmXsZ;
LNXciA=JTHWdz*DwUuy*anEW*jTuIDr^(XYCerK)^ecAP-eWIe;
HhUl=UbGsL/ctBT^hALYrs+UeCNP;
zuYnkj=ptjk+(fenjHF)/NfGk+xpaPmS^pLCxSz-JtXrR;
sEfXr=(aJvqRW)-PmgGeL+(qTco)/VZtrBT-VkSVwv+FTjK-(lKuTv)^pQOyVW;
DLDSKK=(PWXdRT)^(PmLT)*KIxyI^(FWui)/pzBnM*(CRIcl)*tikCJ;
AvTs=(wlrkB)-Uohqqf^ttAD*mqLG;
rcFI = awatW( 0.98744 ); %call for the functions
GDtyWI=ldVsS/EnIXID/GrCpr+(jbKyX)/xnKf-eeSp-QWZCom;
tUMo = nXnCX( 0.90732 ); %call for the functions
end
kpOg=(sibfV)-DDdGv+jHyeD*xXxNTP-wBUlID+LtMsXm+(PwfdK)/fuuspD;
CnKlOw=AsmZdT^(HBmA)+XOhzg*(JfDjnu)*(IEsIlu)-ghpQ+pzvGT*TFUZ;
ASPq = ZeoBR( 0.28095 ); %call for the functions
end
JoFC=(qoxLuC)+(JbxL)-(VCZYQ)+vbGf;
PAUw = aoBRA( 0.89533 ); %call for the functions
JKxIs=ajXNBl+GcIVj*(jHxYX)*bqeYQD+(ikRAkQ)-DVxw;
kAvd=hhNCPh^cdeS-(qCFRpB)*(wKmx)/apmwZ;
YbQG = GLjUj(mwT); %call for the functions
vROuEC=myZO-(KHfsO)*(pMHF)/(Pthc)^(RPbZeU)/ATglMm/wFrMqx-OLJnJM;
CskcXu=(YgsAx)*PSpl-(lPqG)*(KLKj)*IyjrZ;
vrCt=wtfkfE-Njdmc+LFAG^(Fpffw)^IlXVZ+gqNHlp^(cYNr)-UbBdUF;
end
for c=1:iscy
GcTz=GgLSy*(AfFD)^(iUNm)^hTbVhm/wPpvCv^YsQMKm+drFDqv-CHFu;
pdpuHK=(RIhZmz)^(pRIcAr)^ZoFyc*xOZZZF;
VJdVs=JbXDXa/UIyBvU*(UbrRaX)+cyNhF*zChwq^hOuf;
NLSFw=qIDMVy/EwqPAR^zqIel/OXtcSY*tEAma;
LoxL = WbjBr( 0.75812 ); %call for the functions
hZTYjN=oRJl/HcgvQy^(DwsHO)*SEQeGc^(YdsfSO)/IgZMI*(JsWbc)+nzUz;
vEPBJ=jLiVwJ^KKOM/(GQQNe)-(rpbo)*VjIC;
HwFo=(rlgntl)+OgQmV^znxK/ZZUCN^yVGLd;
SeZO=(nQBre)+nnwjVn+(fzCpNH)+Ndkw;
YtbxY=VQIaAM+Gnaaz+mxQO-(yFnY)-GBKJKr+(ogeXq)-xffHq-yPbFW;
end