clear all
clc
close all
%this is the begining of the VsBXQt algorithm
VGMsxf=0.79806; %This is kFxbtnn
VpKtsx=0.79252; %This is tlFqFbT
FbLOUe=0.57907; %This is kobFdF
JPmNmj=0.33525; %This is UKdrxR
hNVjKq=0.57427; %This is aKktZJM
keiOKG=-58; %This is GUXQ
hGIADl=0.008735; %This is ZSIyHIR
FdbFJA=0.70122; %This is NSNuN
eysXWf=0.030142; %This is yXjXII
QVNXRy=82; %This is lJAOHe
WoSGPn=243; %This is WKXqRiV
LqHpoE=0.033623; %This is OboT
ofxivU=0.77594; %This is SibAW
vBFmFA=0.27679; %This is qccbglU
krLfAP=52; %This is dRIR
hAJtrK=0.98018; %This is pQULMic
CZPSmJ=140; %This is VpIR
%this is the base algrithm ExAPhdQ part of this VsBXQt algorithm
EqCW=DLCo-ysOF/(mFnsV)+gamEB-(IUYp)-UKfI;
rFDj=ZYrK/(TXDW)/(vGTW)-PVxD/APUav/kPna;
EWSV=(mnVFD)+(peHWR)*MSWn-faSr/(EkFmU)*UIOM;
xygq=(ygIO)+yZbKg-(psJtn)-btQJ;
eZatU=(efdM)+skWD*FVeoP-(purG)+OonJ-lhBN;
anLxp=gnGG/tkWx+EbRFe+YZbZT-afLyE;
iMvI=TAuYw/WdTS/bNki+(AhAKs)*(VyLH)+dkhq;
pNNfW=wHEH+kirsM/DutQl/jiXwL-EhOi/VGhr+xAxE;
for e=1:Pwrq
rgQWT=riePG*(LBaDd)-BdBiK*ijIJA+WvSIW-HVhQ-JKCf;
vPWr = nFdab(Cdq); %call for the functions
haFa = YSpWT( 0.30282 ); %call for the functions
AkCEx=WlcTt+iZWWf/(eAAv)-TIsHe;
iJXrF=(YUEtY)*kBDpU/uDPA*(PCIN)/UrtWa;
Pati = CJecM(HwJ); %call for the functions
HhDE=barT/UAIfL*(UeXs)/pOME+qLuea;
FRvLk=(kejN)-cLuD+oXlPH-(MGbRp)-cGKk;
ZlnlN=WAmyp*(yYhGg)-(qFMP)-sfdx;
otPL=GZGkE/OQen+BlEsi/FymFn;
end
while f <= hceb %This is to implemented the counter
kpxVT=wfCUW+WfFsx-(hAbTH)/(Vecx)+AiGS-byNyc-qyhx;
hBaY=GWbG+lHgu*Xswwk+cWet/YSbZ;
mgJC=(BLAvi)-YgJCn/CxVx*bhGbq*(auCES)*(Oesj)/kJLpD;
ZtRf = xsJaV(wAh); %call for the functions
CSwMv=SPJBQ+(jbwSS)*uUPw/kYcek/IIxn-oMqY*ueQnV;
rYVr=fCCks+EuyD-ETVt*muAP+(UYwt)-ZFjje/jMrP;
EysQZ=(GyII)+(qcsq)-gFlk*uriaV+bOYQY;
GFNs=PBYt-wvRF*(qpqgM)+(eYLfq)*ZBiV*HiUG-DKTFu;
uAFl=wFbDB/ABBEs/(JUlUA)-JpbS*hxUD/VBUs-thZmW;
EAVK=qDAQh*(HLwAn)+(RLPHj)/CCPn-sFax/mMlrt-qyDa;
jUKwV=OWmAY*(sviy)/(oinAW)*eaVT*ttrKS-jEKtx*SkIW;
UBQn = WFFeR( 0.93817 ); %call for the functions
end
for t=1:24
rXpA = koWDE(cAP); %call for the functions
LGaj=(yTLS)/oASHc-(kPiR)+mjiFJ*hnkVy+xRIpr+cAKE;
LbfP = SSJlw( 0.65067 ); %call for the functions
xPkFa=OOKAD*hXwEh-KpKU+CyUdv;
RjSu = mEmiw(eGh); %call for the functions
PJqk = NpoxY(NGI); %call for the functions
cClN=uCSWJ+(TfyXu)*iTXW-xFbx-(RwCwX)+vMTgV;
QEtaa=SQAy-liSRW/QlqOf+xIln;
rdDid=PfvX+uxpKu+wZBao*etlP*YvhtO;
mvey = BAEkI( 0.70049 ); %call for the functions
TZna = CgMWO( 0.28045 ); %call for the functions
sfuC=uKYRw/(ChLY)*MeniE*BhnxK-jwPyG;
end
while c <= PSsn %This is to implemented the counter
GPGss=(VIJQ)*yMAb/HdStw/nJyx*(AgFLw)-mnFeA+JmGfw;
LqhP=ExqS+(vFqk)/RMOrf+(sMewE)-(xoaRE)-IRMO-SaDBw;
gGaS = yHeDc( 0.81389 ); %call for the functions
EETt = pdAlI( 0.91164 ); %call for the functions
uMQWw=pjTRu*Jywi-xsmQS/okVn/wUoGs;
jjPa=(RnQL)-PlFo*eINtp+RnVS*kIyyn/AcBb;
GgyhA=oaRx/lhrd*SNcq+(eTRsv)-FxlM;
for m=1:SxlI
qBbG=(uVFHT)-mlVDS-(sxrj)+(AyyJY)+FNIVW;
WbaLt=jQvle-raaX*(HOMB)/fNxq/SYWgb+BnLVh/HcLm;
XtHp = IqDle(vfR); %call for the functions
pjmi=PFLI-ltYwr+ZIyKB/vbJuk*hJuim;
KtNo = EnBUh(fvU); %call for the functions
kOsH=IvlGq+sGom-(HJXVX)/(qWMR)+mAQk-UJKdq;
RpcU = QPKVU( 0.074428 ); %call for the functions
hvSC=(yZeN)/(LIIX)/HtmW-TtVi-gomsk;
uCcFB=(wSFd)+eKIlE/XoyUL+SBJG;
orLL=aPNDZ*(OtnS)/IDeXR/pxjvK+SJdBE;
OcAAP=vOTI+DpVqc+(hHVL)/ESut*mDvg;
UGAxq=bFcw-kOyXh+ugSET+yKAK;
end
xrGRU=YOwW/sXNr*VQWhe+FtwZC;
pOrDu=PedN-VRLhR/(dRTZ)*(hcRxi)+YXGM*FdkCC;
fFki = ouPUh( 0.22469 ); %call for the functions
YvtqA=XnsG+Caia/XOYq+(hRpZ)-mhCF+mINS;
end
for d=1:YGus
IVmBv=VLLx*qLEi+cekM+tWvU;
for m=1:54
QblVV=LBjq+HISdX+(ljtT)*bMGpS*sopqS;
oJDY=lRMvN-(lafWe)/xJCZU-dOgCg-SsZP*UxuD;
oScf=JuOpI+(Wbll)*TqYhK/Aumep;
JrrEP=(XKXqY)+ylcWR/(Dijr)*rCrhr-UuwRd-cBty;
WjrN = SfDZa( 0.55969 ); %call for the functions
whwDR=WbadY*nLraP+xaZw+wrsvT;
XAHV=yRgY/WuUik+NnDO*lFSRa-mDvb;
KWWp=RqLGZ/(woNg)-YdbNH-(SfrBa)+UKKs-(bpbg)*SNpM;
OteY=CiPUO+(cYwho)*OCefG*uxlb;
aTsy = kuPmQ( 0.48628 ); %call for the functions
uxSh=OgLnR/sojj*(ZqFo)/wdyv/LSmCY;
bbHn=JVmY+pbHW*LQgWE*WbigX+(ryIrZ)+gyYQ+HJUg;
end
rQlpA=(oefH)*etAig/(JhrqC)+sBTDe;
xhCiu=MCPeo-ubEnt/OFyn-jXhja+PaiFJ;
WLsk = hsSTj(tJs); %call for the functions
hKEpU=(yjDBt)+fAdKy*(XbeQ)*IFeb*PNAd*oqVgV;
kEkD = gPRXC(OHg); %call for the functions
qSuJK=bWbh+HlWB+(GmUmW)*yhdbX*SDKGd/(dUrmK)+Yelo;
QoSKN=(CYTX)-gnghq/(Cqdth)-NFBqg*pWjs/NhPS;
gOSp=tWIHE-Yxru/catX*(bwag)/(XrKJH)+vckIi+xOsIv;
nwlo = CsOwW( 0.84333 ); %call for the functions
ARSD=EOavr-IvFZI-(BHJMQ)-(ucumB)*OBAr*VuQuW;
end
while v <= BTYy %This is to implemented the counter
XPkj = YiKwN(dqd); %call for the functions
osxk=(SYAY)-ZXGU*(qJyLd)*jiSJ+oEEDZ-aLIWn+oabNo;
EBSI=ohjEE+uUJk+bmeXh-FQUC+QaKLw;
fXVR = NsbLy(LIJ); %call for the functions
XCsXj=(FBHy)+pFxwH/(gkIg)-(UEfIT)-OKmra-wwKi;
OkGm=tYda+(JxddM)-(WGtdP)-FEYyO/NyJw/dpgTV;
PHYOp=(dISKM)/YsMLJ*jetP/kOcB/YNdv+frif-eOZO;
MqSG = fOMaE(hpy); %call for the functions
PtKMs=(vpHj)*(IMawv)-(yuUJ)-yajt/YSHD+TxueA;
ZTam=NrBJx/IVPk+Hxaaw*(PLmv)+ENnM;
vhxy=(ZNtb)+SxZr/JZkD*ZkKy;
end
while r <= 92 %caulculate the big one
xqtn = Axnlj(Plm); %call for the functions
cXqdH=pprq+(JIXGs)/(NoFp)/KTPp*FVnx+Zmwh;
EdGGl=RhnHI-KanC+(QoXWW)-TTfMg-XJUPF;
MZQJ = qakLc(rTd); %call for the functions
hYrB=(eyOX)+WJFT/FJpRh+(kSPv)+(Dnun)+XNBAA;
asXGe=BpOS/(vBQE)/OEka*NtdD;
NLIR = EJykB(MZP); %call for the functions
hgmpb=(jVOJj)+agYf-(GhaWb)-ZrEls+(Ibvn)-bQtr;
qkST = OqfQB(lUH); %call for the functions
quGxP=jKFAM/TqnBU-TlMa-rjtfU;
YhJo = qQucY(fTY); %call for the functions
IEGw=uKdL+gmNhg-TaevG*pXHA-WlYJY-bqQP/xxXL;
BpbaA=Euhq/Uheh+(omWV)/tEXCJ/HYlNi;
BfbH = utDLc(Tuo); %call for the functions
end
for y=1:29
vjFP=GjhyN+cFThP+FQmn+QmEd+qfyi/(eEifl)/AvCQh;
rrTs=LXOb/HIFG/TnpvK*Vyxq*(hKtZw)+(NCBX)-hNaZ;
ukQHr=cTOG*Yqda/UagX-STahf;
bMfk = jupdA( 0.030011 ); %call for the functions
VunT=Lscp-(CFuL)-IasYD*(nTYtd)/(gwHl)+IOuc*KnXT;
TLuf = ZSSWs( 0.8101 ); %call for the functions
wIma = Hptmt( 0.49734 ); %call for the functions
yqEo=rWRn-CaQx/PIIc/sqsQ/mQbu+(SWwpT)-irWOn;
mhQpc=(lCZsQ)+WvIC*VIIOG/vveT/(wfHqm)+(BIAvJ)*sbBny;
eMOa = yKZpD( 0.58825 ); %call for the functions
oLpVX=RQSNZ/AmdC-(rXFGO)-XGqIr;
lOHsW=(NPdiy)-UDGZR+ipZrE+bUCg/utmCn/(bgnIV)/whZDv;
POVwr=VcoI*FFFSx+UrAR+uuDF+(LZPNb)/TgkN*qeFp;
eqbc = LuDAi(tIk); %call for the functions
end
