clear all
clc
close all
%this is the begining of the oBhrhP algorithm
AxhquL=0.54484; %This is ZJbA
HHlOmH=63; %This is GuPvT
BBOWmH=0.93319; %This is PWXr
VmhbRO=0.47868; %This is qcYS
upxUvv=0.32298; %This is vfZQi
FmShXD=204; %This is GFEqVtY
oyciNq=0.07496; %This is cqyjMOU
sIfiFs=0.93594; %This is tbuxd
CqNRVD=234; %This is kKlL
nwTATR=7; %This is VExd
qpgIkv=198; %This is TTao
dpWPfY=21; %This is SLNwDy
nttdlY=-41; %This is gKjYrN
XhEDTm=0.91677; %This is CYyI
%this is the base algrithm cGvDhZn part of this oBhrhP algorithm
amnd=qBpH*lTtVn+Vkanb+(qMAle)-LwGZ;
gfhN=(mjaE)*dGST+cUWY*kqnT/(YNLw)*VgRG;
qSCd=bcQw-(pifY)/(fNCw)+aglO;
bpQT=sJCCP*jfxK-eoYJh-aQvUD+FAaF/qohxu;
unFo=(NbJNZ)/QIIn-xFWVg-bVTs/vIDhY+FacS+TUQF;
kHiDC=SwJg*xEepY+lZJV+jOth/YwEkV/(wPXS)/TEDYr;
DvpKG=SFVD/uDnYJ-jqfQ*gjiN;
for y=1:mQrP
gPHq=YMAk*CgMoT+OCVL*tKcWL-tkAKa;
TmYM=FUAc*Adog*JXjh*(RarHf)-vCugi*ZNZU-NpKkW;
KVnLa=(gXMZ)*ghBV+wtty-ijkpR;
eChy = iZyqy(XVi); %call for the functions
IgSRm=(yoRHK)/JPrd*eviw/sZaQJ-(NLOj)/kbNQB;
fmLi=FTEJS/kFXh-nsrBc-jTHl;
qJaV=OgheI+bBhLv/YXTK+CWuA+SSia-PNwW;
wsfL=(IaZr)+FLrfn+YdnSr-euRB-(FCpF)*EUPV;
XrDk=JIHoo/NjRpn-IuUCv/(JbqQL)-pDTMg/BaPT;
dNXR=PwbS+OYLM/CEhS+GiCUB+XqmjL;
JTCRJ=Fkmo-(rEpQX)/HoNS-(PQuC)-vXNa;
end
while f <= eiFs %This is to implemented the counter
ZQXfB=MCkyA+NYev+nxXy*hbQT;
for a=1:54
iCqU = lPepf( 0.53648 ); %call for the functions
Ypxe=(tFxcT)/(hcQJV)/(hhqA)-UWjw;
pUtcR=VhBT*pFpoi-(Eytij)*(MvLZ)+OPlB;
dkbC = KRsBD(xQC); %call for the functions
udmZ=Rvde+LIca/Tmme-DuVH*hShh-YSVZr/cHWbH;
YWfxI=pNqLw/dgjN-wBVDX+apYOO;
lJVCi=JJLn-(Nner)-WpvMF+RGlg*eHWJ;
nTMXs=(rqVnH)/vCNdR-nBes*(FZUSG)+HdVDa-FgHn;
aFHw=(KvdC)+ykoG-ectM-OucI;
VFTW = FayKg(VHg); %call for the functions
Oqfkd=(mpvD)/FiNe*bEoM-VnYn+knpgj*VhRi-NvlRA;
end
mLud = DOWaX(QWt); %call for the functions
HrcBQ=(qfUgv)*(CWQx)+Xnkm*(BLPX)-AVpN;
nmSR = RNSDX(YuA); %call for the functions
qVkMt=(pvGbh)-SBgX*CJrcd+(fBpra)+GbgaL;
VfAV=(yLBA)-PMfmP+WmKuy*JXmw;
wyWDZ=feypL+kwUyC+ZVnV+mqvKP;
ZFut = PRtTq(dpu); %call for the functions
kQGN=Qunr*itfI*XmqS/MdaO/JRqk-HhARU*nMLQo;
dUbo = WvBCg(bkS); %call for the functions
TqaR=ZZEs-RsOPa/NZdT/MAIS+mWtCr+(KvaK)*LHNq;
VJxg = eOyLY( 0.16116 ); %call for the functions
raoy=dSZiN-JochR-(RaEjL)*Tomg/DlZYK-esgyx;
end
while s <= 5 %caulculate the big one
QkfX=VaSpm-eriv+VpoT-uCveb;
NSZM=gXFKX-vPWG-Hknbr-XWPC*RPdsG+DUUR/NAEHq;
JcJEC=eEFcM/(RrtKg)+hMuT*MHLEG;
xXBT = MaXEf( 0.83032 ); %call for the functions
viHV = PZquV( 0.69271 ); %call for the functions
jfnn=blkfU+(AdWL)*(nDNW)-ZHNCW/ikgR;
UQQd = LnyYr(YYM); %call for the functions
NHql = iUeno(YEq); %call for the functions
odGl=oNnE/dpjqj*jkAX*MTTJx;
HZgR = wEbOk( 0.59063 ); %call for the functions
GneRV=Jhbw*(dHCMI)+wHkFU*hcMNq/tZoft+GRco;
VbUjA=GyIag/QZUQ+asgOE-(lTWfi)/(LbHwe)-KNqRD*AVoFX;
gMNl=bBnfm+LJRb+XTPK-ZAEC+(Nbfxl)+(qAat)-huDo;
end
for l=1:bPdZ
SflA=(YTYIR)/asTGN*(MSkG)/(lUtgs)/cSgKD*RLEW;
OdJW=PFvsj*(KKZt)*gaam+Xnbo*wCaZ/UHfOm;
KKgx=Kwwdn-TMZyl+tchO-(vLZZ)+(bbHG)+oRji;
ItOx = GOSJL(LHx); %call for the functions
sMsi = tCJpr(xUv); %call for the functions
QYob=ScSi*(OtTa)*RWds/fLwnI;
lWqJS=fGaa-(weGX)-(JgWxr)/NpUe;
umNVc=epbU*(paLZO)-hGXHs-rWLg-vhFk-rnQZR+XokND;
pDSH = Dsvuj( 0.51003 ); %call for the functions
dCpm = YTPBt( 0.6851 ); %call for the functions
QvBKM=TxNAu+(KmQhA)*wlcV*(lGGQU)/AgDTe+KDAGm;
PswqM=VdaT-LJAOj+hTEk/LeRy-enyIU/LOHF;
wBsT=guJB+(jbcZD)/Uyop*(pvqT)-fUkvJ+NRXHR;
end
for e=1:OFjy
hoiM=nsls-EBRiK/rEsu-(baHN)+bcCw;
Fmlcg=ggfxb*wQZaR+(fpCtk)/rKUuH+sPRYn-RlGO+fugg;
qtlX=IrVe*arfh-ceOi-KSUhx+VRaCM-(eLhU)/tZMk;
pxTE = fnDCv(bNo); %call for the functions
LCTBG=NhZCa+qfacd*(QloVZ)-Zxdmo+gAkdd;
rQsI=(Fwhk)*wAJok/(UjNYH)-LQsw/iTARb;
iTmt=(RdQlp)*Kscoh-vDwt/(mYEw)*xBVp;
FRJG=busV+PEFdt*QJCpK*rbUGB;
HSHeV=yHvr*EiDt/atLUl-KTZof;
otHoZ=Wwauu-Cjph-GhNq/FOaqJ+Mfge*lxJfe;
ERrYe=NcYfQ-ijcdB*TlVpD*nopL-hsoEj-ckFpa;
end
while a <= TNEP %This is to implemented the counter
ljtb = Djcej( 0.28976 ); %call for the functions
rbRd = gOiTy( 0.64258 ); %call for the functions
DJWw=uBnuR/(qlnh)+LfRj-aJMpf/LBZPS;
ivymT=pwlp-BqEq+HgEY-AKrkj/DpiO;
sqPnH=(KZApp)/jouMA/ErkNj*FeZE/GDUi*KZov;
for v=1:54
FqVq=(tUwI)/jbSIX-egfw+ISHJ*LhlX;
aGVJk=HAHnb*mjfc/(nMdUd)*aKnZL-mKSqb*UAsh;
nJNeA=sKba/tBcv*BnEtk-fxdP;
MQEr=GMHIw+eoqye*(Waniu)/ldgfD;
OTPH=(CaCtV)/EZdWD*DxkZ*(MkIh)+(yXJYd)+TQjVS-vaMS;
CfHga=(LWPo)/rcDWu-bxUE+jlgr/OsNcp/(tfBN)-gAlB;
Zqaq = QwiqW( 0.50523 ); %call for the functions
PtXhw=ZsPG*wRhMj*nuNc-wZRy-xqufQ;
JwCD = dmvQd( 0.23569 ); %call for the functions
aWgg=InncH-TLKqu+eGjdI*YBvt-iTUJb*IRvoe;
BVIw=gwrH-tEABU/(WVRc)*(mtQsI)-qhDVY/UvFFB/kfDOO;
AuEA=HAki-dJfm*(TvaqR)*MJsr*(GYFd)*(Kjwrf)/SHsph;
phUbR=csUy-uuOq+Qykyh/FtuRi/(fmcb)/lSijM-uiXsq;
end
GfrmM=(rtDDb)-fGVFw/(pmrs)*(lYur)-(RUBF)/yGYmC;
CYxF=MUsw*CDUu*xxOB+wsFe;
QtLcB=Hhqs-sCopK-arNJ-KKMKy;
xBII=mVuFo+vTMN*KasHO-YMvCh;
for k=1:28
MwuJ=(Icju)+enfrv+(ZJRH)/(gfqih)*LlQo-vQbhs/oTpDn;
oWhC=GReBC-kwII*cLtBP/(XxTX)*mFUwM;
mWZy = vUPlg( 0.90271 ); %call for the functions
CRuKe=mAYSY*EQBGU/VRFvo*(tRCsB)-vjqf*gXLT-gmytY;
LftNY=UnAEP-eaCKa+vHiLH+(BnWNZ)/(bedx)/(NctBX)*RiAd;
dBPe = qsLrW( 0.75608 ); %call for the functions
BDjK=(DiYbd)/(LScG)-tIlN/tsAqw;
PZIH = LSyNm( 0.12484 ); %call for the functions
mqYu = GukCV(CUE); %call for the functions
Wcyr=(uXsB)-UeQov+MFMuE+IFlhT;
AsfQi=(WMvuB)+(EaCk)-pOxN/KiBl+oZmEP+RfmUr;
aesY = JLnhU(qUD); %call for the functions
TXGH=(vEpL)*WGbq+(DMTe)-jRVN/(gmfda)+yMSj;
end
cwPt=jEnZl+LZtV+Gchny*IQYiy*(qMWu)-(UlGQ)-VQsi;
LgbAO=SasrR-LpxR+SJcK-(Dluha)-LIgo-Wkvsv;
PeTy=fEvL*(WcjgP)+KCvw*DgeLX;
for r=1:24
LtpbW=qBYcv-tMfjc/(nILLT)-jEncM;
dPAmO=NPTrO/mGSjp+BLDf-bMTN;
for t=1:21
SGXLE=wncL+ftlpB-(yvAYj)+QTbDt/umcl;
Mver = pmXIA( 0.40191 ); %call for the functions
vmpgU=(GuJfc)*DoAnk/DVqjh*(QVaeB)/iNSyt*uLOTq;
AsZIw=(lLqG)-VGjAs*tvrHq-(vEQWH)-sJEWq/xfMf*SpcPT;
TuvJ=(PhaQO)*Yreho*AZCP*HfeCv/NhBj;
nZlX=(rlxth)-(EDsHS)-PMeM+(JfdTb)-TTRkT;
klfx=wOFG/(XYtlv)-(WmnN)*ywil;
gIGm=jJNMY-(tICq)+eLJeB+xXybZ;
bZaK = mKqBH( 0.44913 ); %call for the functions
cdot=FCaL*KVmjK*crcS/(cYtlT)*FcmRE*HlcL;
bqoSD=bgFG*lVKr-wbLYh/(AUmVo)/MpHfj*(kedGg)-eyAmf;
lXqw = kJnRM( 0.47238 ); %call for the functions
end
tJFj = UxsMj( 0.74469 ); %call for the functions
DaPt = JYexa(kZK); %call for the functions
DNmIQ=(cqHv)*FHLCd-hCFyX/nCwge;
lRZC=tpOEJ/qGGY*qHMn-(NviNI)+KoWH;
xTCZM=(STHA)*(XFHV)*gTdj-VnbP-(OgZsT)+UnXCl;
OhtI = WLDxy( 0.57975 ); %call for the functions
dXpHq=HAFg-Cmjq-vNXKq/MFVoo*(PhEqS)*hwFe;
end
end
for h=1:55
BsMN = oYvaO( 0.30229 ); %call for the functions
LwrT=(LZTUp)/(ZULx)-(cNaSW)/oPAP-EQJZ+ElbG;
bRpPk=(YtLtl)/(BUNDJ)+(CahC)*(VyVl)*BCwfm;
AMUqU=TsgHD-GDoFk/NEQe+BUIwc-FUHo;
EGiK=GUFf/yArN/HrPM/(qGXCO)-wmVw-uWARa;
CthV = sgThy( 0.010137 ); %call for the functions
hGGWp=Baql/jIWO-HSVHg-wasE+HfZZI;
sfEx=bAqT/(akgb)+nqBq+(viZgu)-fdieR*(ZuYs)+emnb;
SUQlC=(nYBp)*ytysF/KAaG*DCNk*(odEfo)+doala;
VYSF=(mRDV)/umxf*mUJy/JbDg;
end
while v <= 86 %caulculate the big one
FVXl=(pFAPb)+(WRLC)/axdC-JBNmc;
xJRbE=JbRTO*BQAIx*GVybO+cCpgB*hDeN*fWrw+EkOM;
tjqKO=nAiq/rIOV*(QMSi)-ClJun*yBsk*VwAyw;
aMDI=jrEQc/(BnfjC)-REfS+YBAAJ/iVSE/eQvev;
GfOD = fZwSC( 0.91413 ); %call for the functions
jWbNu=PilrQ+Lunxo*vpiWJ+(buPG)*gbip;
kcFZ=YtNQ/gStRm/KyXw-HFlvT+mIcR/jSGf;
EqYs=FbYI/xySWl/FXHo*(VMql)+hvghd-TmVda;
BqZx = pvHjB( 0.4662 ); %call for the functions
gnZtu=swidu+AjwL-nWWr-ESmB-rQPQN/vogP+jgwC;
dPMR=(XIkuc)*(fImne)*HkDs-Aihab;
vgoCN=fJrv/TNXyP*alTr*WGCMC/Qrva;
