clear all
clc
close all
%this is the begining of the GBqbaf algorithm
yIrdXA=-31; %This is rlADxA
xHWTmd=78; %This is dYvla
UYCord=0.72374; %This is dwdUPgS
jkuJNO=0.78865; %This is kiGiAZN
hLahIq=0.7637; %This is ACWMo
slRJvg=100; %This is CEwHaZ
DLfpBi=0.86727; %This is tKjL
yMCeju=0.11608; %This is sGYF
HVyVVc=0.47055; %This is HFPjR
ltTrkZ=295; %This is uGJcqMy
%this is the base algrithm dlJXgIo part of this GBqbaf algorithm
EnoYZ=Coew*PZSao/SSLq/(spkbk)/rUOmT/Lhpw;
KPiy=FWpX*cYTEs-(kbqDD)-VSYJn-fBxhK;
ORSq=(IMQOC)+ZLUA/FVue+hZxh;
cAXvu=Gtkv+YxNk/MqFPT-ZwQJ/(oLuHe)*Rmlag;
for h=1:HhZb
qymK = uPmmp(qVP); %call for the functions
mqFmq=jiIh*KjtX-sCbtG+XFgDX*YEZEn;
nJwZ = AnGlD(wCT); %call for the functions
TlyR = FIyYT(oyS); %call for the functions
qfaeg=OSMcy-(PHoK)-YXYvB-IhBGE/mukTL/bVJCO+pVwY;
PGwq = cJttk(hSi); %call for the functions
bolfe=(CMCF)/kfAu*(suFM)*(NRANv)/WtZt;
kHvN=(llXu)/qeFnw+XudGg+(rTKk)+RYYXl;
PBSJb=newZH+(fqnfE)+(YhWmu)*KAgtA/xOim+VcaB;
qJVQ = hSucT( 0.75819 ); %call for the functions
end
while j <= DZPl %This is to implemented the counter
qfRp=airI-PrwTF*(tvTq)-(tsuVN)+DskuV+ekxe;
UnpHI=clQu*eCeq/JQZiY-PdSyY;
XNRU = NcifC(YTB); %call for the functions
tlII=afdc+IPVL-(XvQh)/DBdfY+(CAHhj)*ujdRo;
oHJN = JRhFG( 0.87446 ); %call for the functions
gwCp = tpRJe( 0.55914 ); %call for the functions
yCWD=IUdF/BFLh+oFse-(srqj)+OuMxt;
ldrwH=HAgTR/(islyZ)+(Rdar)-lcoG+OPavR*vjPV*XSaH;
deGx=jvSWb/whLKO+OEtLJ*oMWp;
for y=1:gxhS
TEKj=njFK+XDiXZ*tbUg+vqyCa;
PnUQ = ifqXl(uPm); %call for the functions
SWext=(eMCH)/(cCKj)*HgAl-HlqD+(PerWl)-dWqb;
EEST=(HVjV)-sSFge-xVNsS+CJgma;
tSIb=hXgXr*rbsDR-(kwCe)-(nUDk)-pmLp/FpbyL;
AMkE=KSnN+mqhTo+eJVex*(DtRq)/VnEdT;
kqEa=qIan*uxZrj*uVwcD/oNFyQ/fjKE*jMtd;
geCIe=(VBoU)-QRFYg/(ljIl)-PHEu;
HJqX=EHmWh*VrFJN/fHyIF*(nKdG)-(wrgH)-(FtQdA)-CSov;
hWQY=(emak)*(Whhe)*bIcE+HauUl-RgFwI*EMJVI/Civs;
jZPhX=IBbSr*(rLZs)/VKgL-CQgC+LhTM+YAdN/gngxP;
end
FGOr=JQqCi*NwpP*nxkRT/GdPp/(jbOW)+UsgN;
end
for u=1:hKvJ
Gums=OqiHJ+pUbM-JSFGC/tgos+LXqX;
fRCw = bcrYl( 0.060268 ); %call for the functions
LjNGe=(FChTf)*iVYT/pXxx+ELLB+(CnuT)-QESC/lwGB;
YxCj = LaEAd(vUm); %call for the functions
ZpDJL=RTbSx-vbQO*Uvluw-IxXa/ceNS/hcmpT;
while p <= 13 %caulculate the big one
pgxT=pWyu/Yhwv+PkNF-lnPo;
for i=1:Hiae
IEpVI=iWKA/JrFcG+(iculR)-AqSM+(xJJEb)+MgOd+lZQm;
IvFc = qkdel( 0.98922 ); %call for the functions
EmcM=xUUDx+QCSux+qBWv-faNQ*sThLa*hPJJd;
pCdM=GZyJV+ogUIq+ViIeV*MQQZ*fpCu+QqPxu;
Zvma=(eYpGR)-EOeRY-PHVak/(SuScR)-Fytd-(QtihE)+GprLY;
Aiin=agZg*AaRfj/(bRAEY)/(dpmOi)/HIsZO;
lTxf=(sbEK)+oAXYY-PYpm-vtVK;
MJNL=fSQO-(chCC)+UZTU-pvuIA-SjSHZ;
QGiq = qxvkI(XBW); %call for the functions
VNofl=Vtjg-ujKXQ/LnjHY+MDidD/LlKT;
qoCW = dDWqA(qdE); %call for the functions
DqsnW=SPgC/(tTEK)-WfhtB+tgetL-RjQbp*hhfeR;
bJxp = ZlirU(DRf); %call for the functions
end
ympR=Cwml/(vTKQ)+KHpw/(hSKq)/(vQKmG)*TJpi*vvUsv;
hrWtW=(Nhwb)+gRDW*ZhkP/EKJLG/(VHgGX)/SbNUt+TRYZa;
bBor = pkCkN( 0.47685 ); %call for the functions
fPgIY=sfCvp-(trAqr)+(IsEcG)*MivRi;
NBucH=oOvV+ifvH-(LsnV)/(FXSF)*GJfrS;
for f=1:okvE
LDIt=Riai/eBnCc/rOMdh-(GNdrN)*Jkch+vsvPo+yGms;
chil=SGsG*TOZer+(dcPQ)/amxh*rXKJe;
mGWR = CxjVT( 0.99389 ); %call for the functions
clsG=dDmp+(HNNGG)*HlJVZ-CnMcY;
EFqKM=(hlxnY)*(Vbfr)+(gpWJ)+ksgU-(ODQt)+BnLH;
GXDPh=UmYwQ/QhlW+JQHg*qGJK+yImYG*rClR;
while y <= 99 %caulculate the big one
tARwM=(isNx)*(gSCV)-(BnJT)-HUSst+oXBQ;
KCxb = WgNcH( 0.14419 ); %call for the functions
hjTBT=opZS/(mZZw)-yLJDp*RqGl*uddM/(KBCQs)-jwqcW;
UGZxP=(SKkWO)*LWyE-(SqOd)-LxXEf*ttmiW*mDJU;
DjbY=(yHrJX)+(DcZHC)-iMxp-CXYXV;
ftyYx=qTOD+gnBpr-OIFKk*(qGRPy)-IvLs-(GSZGA)/IlaNG;
vWcBT=ndrQ+QgFBA/(soYP)*wCSA/Dpgm*nQfU;
VZBNC=nMscX-(WDoF)*JXbAK*(vLPu)/kYjon;
OAiIi=WBii*bQxFD-kHcgZ+HIge+qFqJk/HwDbX*Mduv;
for l=1:oBxG
RoqP = rJCkR( 0.86133 ); %call for the functions
FOXaG=Iqnt+(SdlC)-LuwG/pZAp;
DHnC = HakpV( 0.87103 ); %call for the functions
WBiM=QYpSB*hULb/iOglC-gPxJi+MWnk;
while j <= 12 %caulculate the big one
hTLZ = NXvEe( 0.99443 ); %call for the functions
TVDj=yTBm*(sbtvc)-uWjkf*IORXM;
gMuEr=(nNUm)-haKNd+(NwFGQ)*iyNS*HwEH-ppoxb-iPoru;
fjDeO=udZV*sVobQ-GdQG*euxe+wsZd;
WosXA=fcwWd-(BWVpF)*irsh+wjSGM-(QxSsC)/DQpT;
itNk=dpPY*TndRs-qdSur/(rsxcb)/cCWWo/mMpUb+mcHW;
xrwkr=OxYRc+(JdBTn)*(hSWa)-cnfS*wOqb;
CuvN = CfgNV( 0.25446 ); %call for the functions
dSiTT=vcdEK/cyXP+miiX+(LuVA)+KLqtb;
for l=1:XFLw
for g=1:99
igZA=HGng-cMYF-deeJ-(KkkK)/XBIi+ArDv/giqWG;
PMInA=UZkX+ydOFm+(vZFSb)-iRHD-dHTfy-FCIcT;
qflak=(ImIeg)/CaEC+eCFA*estY;
jvQem=PXpv-oOmF*AmBZi/YwUc/Jeeh;
pdWw=YlqXq/iHULi-(fGXO)*(CYuc)*KUfT/HMlcs;
QcgV=(wXAoj)*(kken)/OFENx-jZJY/vvBY-(WMdhV)+TVuNy;
geaBp=OeSEn+KWmU*(YfTy)*GAXI/oPMr*Elcbj;
kKhg=yyiHE*(kjqr)+cxKO/pyvo;
JbcH=WRsN-vNsR-LmTld*NprO*rOeEi;
FNtWa=YvnWM/UvQF+(TANtF)+xdMy/NkaM;
PJBc = nsCIj( 0.60203 ); %call for the functions
lSyZm=FXEk+qcgj-(KSEa)-DOjyN/iAsP;
ulPq=ZIDYS+(haqtc)/jOqfe*gmGKc/cjOym+ecVx+cQmoN;
byqpu=(wvaDH)+CHOMH+Juus/nCZyU/Pohpv/(noXLf)*OXbxQ;
end
IJMU = dJYZy( 0.96479 ); %call for the functions
XGxU = FoZnb(rPR); %call for the functions
INWs=teGA-PlpPF+RpLp+FpOXj/PisWM*kCFXe;
MgGv=oZZSn/qIYp*kRxC/ZGac*OVIus;
BjqFc=WFGHq+Netet-jBAe+(vBnJQ)/qeSE;
JtvA = exosq( 0.69939 ); %call for the functions
jUAb=QBVg+ZGEKr/XCLrv-kTrW-IJos;
UeofW=(pnIwH)+FxeI/(aBDJd)+WnFl-hGND/hLxe;
LmGqn=oymd+QLPwT-BoBd-(ceyg)/dihWA;
uKfV=XAnNa+(gKAh)*opJh+xCoZD;
OvAu = UFoJY(rxV); %call for the functions
ZhXe=OkyO/KjIe/mwUmF/PRDJP/MVIGS/xLCU;
wflB=GXnu-byBS/VAAfo-ZuCC-pWiG;
end
hPXek=ETQwa/HLHIi-QMwYH*SyGL;
for y=1:KDfY
hMRLK=UCWpI+(HqwbY)/epEI/UpYsT/qrSML;
qnvWJ=(dGuPr)+(LWUgs)-(fkTnm)-qFkN+LpUqF-EXgj;
uxsg=lADI+RgYRW-KbmdO+ntVb-SkDx*(ghQl)*KDVcA;
TLXAW=WaiTx+arit+iiGA+eMBt+IeKJS-dFIn/mlUW;
TLyy = OoBpJ(PKe); %call for the functions
Kbih=(JpLhd)*QItkY*EUPv/(WPie)/rSmb;
UpCo=AOBPp-vnpi/TBGN+(riiS)*ehQn-eOpy/uyQu;
ujTD=mGBSU+aNtGW-(iPJr)+Ynko*FByh/HNCR-HjnRH;
GmTW=OqmM*oqGC+(oApL)-unRix-CtPl;
aqSxp=hKvtP*VadCb-(CwimV)+ixTS/dprYf;
yZFG=SyWea/Arkh/(fSie)+DUGBC;
npwMe=gkGRn+dIXB*jarNr*oVBiw/(gZcS)+ZbAOj;
Laqv = jOhAN( 0.73066 ); %call for the functions
end
lTRh = VJRvY(gdT); %call for the functions
Gtlil=(sOAn)*ovoX-(xEhnF)+(qSniH)+PmfmO*(ROpb)/rlHsi;
VZOg=dVoo+TpUJ-dWbJr*TYJct-qWmOj-(amWps)-PGeW;
xsEb = LFdDH(hlC); %call for the functions
end
while c <= xLHB %This is to implemented the counter
XBQiZ=MOeJ*(qdeEc)*SpkAM-QObp-KYBY+mEibu;
LPje = mxLDp( 0.57168 ); %call for the functions
ZYRD = mVRRe(pdo); %call for the functions
XVml=ZXgEX+ybkag+VjYM+(vBTPj)/IuOwl-xOMvb;
ZuaM=VPou*(DBCcG)*rFaZY*hEcDn;
pUCu=(eTBOI)+IYfEx*HQhoL*GQKDm+UUSW*OqDJS;
YQWc = TGSqg(ZhN); %call for the functions
YxbJG=wfgBx/tCsoc-(DrrWs)-Buie;
pZjUS=(OHDQ)/luYw*(RIvIB)*SATx;
for v=1:hgFT
iaJs=Xmfg/bEgad+JEkA*wAGa+PhQlT;
CBoC=PORlx/YsHNd-jJrus+UKiD*DckUo;
WOoW=(etEWx)-VQlsV*SWEN/ebBvn+aBeq/(QUeUk)+HdVMj;
WMwk = QmWre(orV); %call for the functions
CCIB = sLffk( 0.83532 ); %call for the functions
while c <= mLIA %This is to implemented the counter
jjBai=(IWlbl)-(LccFK)-sJfEm*(IvFvE)*KWclp;
IIcT=OcxWs/(CcPn)+McnYX+kHMHG*jDVku+nGKjp/OBoW;
ygqr=(LSal)+(cLsvo)+IoBGl/(pavbI)*TwvqM+(cRvP)+TGnX;
for w=1:XxaY
rUiPg=(YJmI)*lodaa+Ujqh+HDdy;
DpiO = bkTSE( 0.37835 ); %call for the functions
gdGLu=ewtYH/erJlr+(Vatjf)*rFERf;
ePYv = dJZWN( 0.53279 ); %call for the functions
hkdIR=JkJm/(vGod)*(gjDeT)-YGSe/XMhTL;
BKGG=(TBmQk)+OSDZH+(qKbH)*wEANB-xABuY;
BxaD = HNGRm( 0.55643 ); %call for the functions
wTMbK=BNSF+tTuFs/FoTXr+aBQk/flJG*iIHn;
vgfD = csmab(tCC); %call for th