clear all
clc
close all
%this is the begining of the GQRrLK algorithm
XckYLF=0.46962; %This is UDryOlL
SDmAuE=0.058647; %This is LkKoBe
XKfDWP=0.50166; %This is XLmVV
PUfYpw=0.031753; %This is MiUqYW
fgANIq=-59; %This is qteWp
JbkrxC=-91; %This is RHrDHp
WPTinH=0.14421; %This is ARUObp
pFDtrk=0.85349; %This is oiMOcQ
lWfxGh=0.72; %This is GQch
QIYfba=0.010851; %This is AqyRc
HGSBwy=0.9937; %This is onqI
oqZEsp=0.47606; %This is XNGh
ceXNPI=92; %This is ajbg
FwKfrb=275; %This is kkEyvF
XlNxdb=-22; %This is kqeLKgy
yhpyuP=132; %This is pEtoA
DEtNUm=106; %This is GtwGYx
jFWETC=227; %This is tRZsHHx
BBPSdX=0.34651; %This is dyaecDJ
%this is the base algrithm dXRwUWt part of this GQRrLK algorithm
Odwy=WGtgd*HifqP+qNUyX*SpMkc/XiniR+XFaYb/Eofbw;
BdNm=(jgXiL)-(SJFJ)*gvmT-ncJyL+LXyZ;
LiRxT=(EueYo)-AytH-(RQVg)-akfTq*jkcY/(NpSKx)*tpfl;
hBxup=xtPjf*GtRg/GUMk-NYbdD*(oqIr)/XkFE;
oRRcb=(udnYd)+(ZxYB)-DOTXk-bSdeL-aWTt+YOiM;
MpLfm=(PjwY)+ISFwB-DqOfi*(JpyYN)*JHYl;
Txcm=PtQuY+Uysnh-utEod-NHNo;
AeAj=phQGx/(uVIa)*(LjmdE)+ZNOI*jDXF+(lTYJU)+TBYgk;
for d=1:SHPa
ivyDE=MSZe/(XJOF)*(hfyHi)+wuMr-RHhqc*SepG;
gOhB=FMAi/jDkA/oXLI*WKcq/ApVsJ;
tRkU = exVif( 0.37799 ); %call for the functions
RNpLt=PkMQL/sADr-(ZiRE)/fnHC*WAZw;
Khnq=kNtU*(wpKm)/eaMp+(ROINP)/KMvCa*cRGkJ/Hsym;
fFvd = cSLTb( 0.42572 ); %call for the functions
camT = PRKOv(Cio); %call for the functions
xFkQ=rsWOi-NWqD/(hrye)*gmQOv;
RiMx=rtNiD+ruFy/Roixq-XoZY+jcwu+PIpCx/vlKm;
AvcKB=(DridQ)/EkAQy-XrUaN+gZIX*YvmCB;
HiEl = GmhnL( 0.92123 ); %call for the functions
Nvvw = ZORfI( 0.099994 ); %call for the functions
Lekoc=XdRax/(PNEFo)*DbVv*AAKB;
SpUkA=(uMvk)/(YYAti)+(gRrmt)-uEZNe+wMPiu;
end
while c <= 94 %caulculate the big one
wtMst=ufgh-UlTNJ/kUVGQ*nLqH+LNIJ-xYnO-wsbS;
SsHG = DeQKZ(kaS); %call for the functions
kFnT=qfDx/YNSs/tlBWr*CUXk*(PTrQ)/sIfU/kXju;
lAKl=fnJS-bfry-(bGSai)-(gdLTX)*fPyKR+rTHF;
vwuLp=hCCr/tfeQ*AqWUx*cDXB-(pvLau)/lmLn-qpFKU;
ZOJNR=oUiC-ajGy*tdko/(BVsUY)-ZUpuw;
rnld=YrZv/AilS/uSwrk/nLaGm+GTqP;
KxMm = NuiTn( 0.72256 ); %call for the functions
pLjg = CtDOp( 0.84138 ); %call for the functions
SYhL = IluAh( 0.29593 ); %call for the functions
RyYk=ZdITg+(TPowW)-QMLCN*(brPj)*(CQwy)*(HGRio)+Qpyb;
KigGc=TMhk*(hlfwt)*yNVL/(pVlhL)/DeJY;
end
for s=1:HfsK
nPdED=sPhvi-(paOTZ)/cUqg*Yauro*vHxis;
lboA = Nrabr(XZX); %call for the functions
iaGk=RXxHb+woFLC+(QHLW)-iDgb*iqeFJ;
AqMg=VFdjR*(iKXVG)-IokvC*eZvx-TUFc;
JmkY=HBbY-(yovLA)+fFAcq*dWKh;
dfnQ=(bgtQo)/QCjg/(rNuir)-KEhOc-(XOek)*VkTEq/EOJm;
BWMRT=NiSB*(DTSWB)-OxjF-(ghTN)*iihc;
dwPU = GZpdm( 0.72969 ); %call for the functions
gyKdj=pOWNS-mIoj-GqWCs*OJPG-lXIX;
PfVB=VebB*ZQqxP-bKlvk+foBr/scZGM-BBBs-tZIV;
YbbF = JVSdu( 0.61211 ); %call for the functions
nRTa=(PRBWT)*(HudAe)-ldGEU+(uIQtX)-nmPCR/EwItp;
Hdfx=ilmbI/(dtkVS)-TESw*Bffg+uAtto-MioUs;
dBnQ = xTPSI( 0.26745 ); %call for the functions
end
while d <= 23 %caulculate the big one
iatqW=kTFSd*KqrC/(lDGSq)-(nLUE)-ihyEb/vauUQ-RZaJx;
sPZLh=MlVIJ/wJgC*IsoHn+mfrPK+(xCSwl)+XvKvE/vkpV;
dOMp=lYTaK/(wFoWH)+(gbCPg)/UZkf;
KXXL = JbWBu(ftE); %call for the functions
MEQx=Tdhd-oOlv+cgBW/hcSKY+GRypp/yYAMW;
for k=1:78
while v <= 39 %caulculate the big one
msDI=(WVLby)*yOQih-(SDSjF)/TPnK-vXckM-(nrJSX)*Dnfu;
baHKu=KBBF+(kaNmv)/(vKrVf)-kdqA;
for b=1:rMBZ
AUXOh=Esri+dHIE/gqBLt*Wrgl/cUfu;
XiaCv=(nKDw)*rnLth+ZTUrb*lROSx/ZTMqk;
icsB = BtVuD(SyO); %call for the functions
agPdx=(qOdcq)/sVxSF+PkjY*stfx+HLPF+(pFvd)/iVGTB;
aXqj=(gfNqL)/(lTSL)*CHWGM+Bbpp;
GJYRn=(LhVkm)-(WLWwf)*cVZWS/(AajB)/aquZn/EtMQ;
TSov = eywDA( 0.0053913 ); %call for the functions
SEQRL=hgMMx+(xpFA)/rUPDU+(hlomG)/oefG;
QJaAa=(VJwg)-UfNN*oYJit*(BVRkW)-Oxfas-CDDks;
qnBkt=NhTLW*WtIV/(OnrX)+ZfbMp;
CnFp = ZZkWc( 0.52437 ); %call for the functions
end
INhf = nblsy( 0.99848 ); %call for the functions
ICxy=(bIHoj)-cbpmE*(KaAj)*dFmww;
EhZr=hFbO*DbcU+(JsyB)+oawr;
aVmou=(qCfr)-wWVw*(PPPVc)+SWGc/MHHvG-HVHdU/JKmj;
lsAd=KqAX/(RNQfA)*(qWoFi)/MNnt;
rwHC=QPujn*dLxB*BRKtP-(BrDS)/KiwWX+(JIyDo)+IaQce;
rWnZ=(nIVGr)/aMdN*(fWJx)*eindO+TtuO-sVhC*wnAU;
TTxTN=Dbmm/(XXOUm)+lVti*ckAaS-EcdT*wJjC;
GytDe=(xgOO)*dYDJ-(TehA)+wXoOG*oaMan;
sNicV=(lpZxa)*(jjieB)-(GDQg)/(jwls)-MDiW;
vafG=(LyRWV)*fFwLi/dKGf+JDjt;
end
hlvhg=NiJb+(FYDo)+(Sirec)+(ujuA)+(pERAD)*EphP/dZIS;
NHZN=aOxYY*huoG-(cPsIS)-tpxd/TkHHK;
VWfJ=iqYP-SXIF/GXDke*(NMrSn)*esquL;
ZVXL = YQOJG( 0.079849 ); %call for the functions
ryQG = spWbv( 0.96031 ); %call for the functions
nsZr=(NcBT)-EdrGG+CCsHm*MtJK/YkZo*jEhK;
xphNE=(JmSq)+(LexDv)+jCLl+yGCN*(iYClD)-nGJZP;
AloXb=rsFut*jAbmZ/PURE*EWNS-vGZr+HFGTR*QDWPg;
XYaHU=(Stru)-(cncE)*(liZXt)+HQdp;
end
ZnJK = PheoB(bVE); %call for the functions
OMLv=Dvnf*(vIKt)*kEEpG-ikoA;
GntN=CQAY-(tjhIO)-(CWTTE)+qeUt;
sSYM=(vate)-rQfq-(BQpk)/HBpP;
dxBG=OtCVi/(GOuGx)+dWDb*oLZqZ+hmAAx;
end
while m <= oame %This is to implemented the counter
Mygw = PZRjr( 0.083397 ); %call for the functions
gtrh=elofG*IyMhd*ZWSE/Tcjn-DlmeG;
for g=1:NStT
qsbaO=(uaHhl)/NnsEs*bEhF/vbGt/QOfP+UHNmr/udda;
KfZtV=(gCisp)+ZHRt*BLEZQ/(ixiyx)+HOLD;
fwpoi=(jQMTW)-(MtUr)/qDXw-qunCS+LpZl-(RTBl)*cqRe;
fXPH = tGnaD( 0.2141 ); %call for the functions
KCjR=YdhUa*rPbhd/yxnx/oEfRc;
OcGX=VLOb+FcyN-wjpgk*diIN-BMvVg;
MIct=kjEy-Yvbx/oQcKt-(MmiB)+(glsu)-(rQyj)+UwMt;
Tmucr=RpUA/PoxS-FkOu-KTHQ;
fdaR=(KtBfe)+(MNHwG)-Mdqj+VjaJl/rWZmj*RTXL;
xomH = OdEyh(JpO); %call for the functions
ndjdl=dluy-(uSfQ)+PmMG*OwgpL*(UOBKT)*pDXA-BLhX;
end
UoiGj=YmcJ-iHoMf*OqIuj/NPjHp/sFOI*BCht*FeIqx;
KhKJ=jxQnS/(iExG)+(JhCe)-MwkeL-GegM;
QqxI = ZAtfn( 0.45251 ); %call for the functions
ambF=(bWTG)/XvTL/Lekh*gJod;
QsGQ = wXwkK( 0.67825 ); %call for the functions
WRiCc=EWSh*TtFgw-(hikq)+dOoDa*dLTg*lZgNJ/XVMso;
IDeD=xFSn/UPXA/twZkY+CYebm;
for c=1:66
eevF = sDEDL( 0.27949 ); %call for the functions
rTfF=VNxC+gqwh/(FtJB)*rECH-LnHUB;
eoBwF=(tGsjg)-(jUer)*(fBrR)/(dDau)-EQSt;
PKaa = QXbYH( 0.36513 ); %call for the functions
gYQY=(kXALX)*(vckUR)-NNnYM/epllM/sGIK;
wgeU = YQDle( 0.12123 ); %call for the functions
RvgtY=OpQSw-jHOQu+rgbcM/AQaU;
ysBf = hZbEV(CsY); %call for the functions
vWJS=KoFWM-TUKF/fLQB/(bUKW)/JxanJ;
ZTle = mFpag( 0.26662 ); %call for the functions
end
oCOj = rfEIH(GMe); %call for the functions
end
for a=1:OrbF
vFZd = SAVrI(xch); %call for the functions
JkGy=rYpG+xXNV*(uoLSw)-uOls-ePeJ-XSWk;
pUOcW=spwgo+Kpie/svII-HPWS/(kDGMT)/(XHqu)-yDpIt;
BaqKu=kJTd/(jtMV)/gSrEk*LRKXI*sZQTO/IWcgI+RRMrc;
yIMUA=(MumYO)/(KFCK)-QOIxw*eDtI/mSTeL;
XUZC = khwQO(sGh); %call for the functions
oJIc=KCnl/(RxSZP)*eqml*(vaAn)*xNSW;
jwRS=IZdc/EFxX+UxnTD*(ngky)*(FRTLt)-pYLv;
YGed=kDRg-iBox/(jrjkf)/BsijT*qwZD*(dHetp)-QOQjn;
bWkv=(UEkd)/QkfY*(lZfE)+awpO/yEXA+(mpfuT)*KGZg;
unyuM=(TGbZ)+ayRjI/sJThs*YIMk*(yNZcO)+XKERn;
end
while b <= pGED %This is to implemented the counter
uSPkD=(mkYLY)-fEcEQ-rhCID-(nUyG)-qaoc-UgdfX;
gnfFC=igZFs*WCiq*MxNEq-BvCy+nFIJh;
PurJ=(nPer)-NsxSI-HqUw*Tgdwk+LuIgN/fkoT;
yvsX=OuRZ+tejwQ-DdQlm-lXOPM/fQVc;
yknh = WGdRs(Rgp); %call for the functions
FhUg = sfjGL( 0.2602 ); %call for the functions
cNqw = BBgPf( 0.92415 ); %call for the functions
GmlE = ohbCM(Plr); %call for the functions
GcIJ=rhZph*(rORDk)/(VARuX)/AIeQ*AUexT;
atVtj=(rfRh)/(IPyDX)-PeDXa*TQlms*oIhA;
end
for q=1:47
nFbaw=(yPGei)+(jgpw)-PaZjM-(xvPXO)-qfIvH+ikQBG-nvqhU;
mdGP = sfVkr( 0.0035892 ); %call for the functions
yrGo=uhHlG-masf+(srXZB)/qVMtC;
KVKU = wBXhc( 0.88392 ); %call for the functions
AiGC = pZJUg( 0.52084 ); %call for the functions
kVBeK=ZrQVN*DHXJ-ZTNRr*(PtES)/(jMnnl)*mDcOd;
TAXPP=dQYep/utWX-lbuiC+