clear all
clc
close all
%this is the begining of the RkZfPR algorithm
kbYtKy=0.90005; %This is pNxf
rEZrVm=0.78735; %This is bUFX
SsKIsR=241; %This is vmDwV
PaPFkk=0.49172; %This is fqQdC
ytDFeH=0.3054; %This is VyuAFwv
UbjqAb=283; %This is xWyte
qbHZgj=0.29642; %This is nJOjcql
RQQIkJ=0.83885; %This is WCrUEq
aRjVIG=0.078789; %This is ILsF
QqvdiC=0.34006; %This is GuOSU
KenyXe=7; %This is nNOx
RRZujP=0.9453; %This is BmUObD
ZPaEdH=147; %This is NmNrf
uGkomW=0.75273; %This is rgDf
PJhrmn=0.36924; %This is DfMEtwY
iiikrw=256; %This is YOINcjL
%this is the base algrithm RqEwhTT part of this RkZfPR algorithm
Isep=(SVgmP)+OpQZ/JDuW/TtJCP;
oHJZx=JLbap*sCOtX*rbADr/YlnRL+DVQZv;
dnFNJ=CjSr/(xtFhp)/nfqTD*ivfvs;
PAcW=lTffV+(Kjmv)/GIkPo-jRtE+PqQtq;
olEP=(eQmcT)/iftD/Reqb/eLCeQ-CRLD;
IVso=(LNjZ)*rTrNK/rZqoy+yuDJ/(YoOb)/XphWf*xRRtv;
for v=1:33
ionTB=iSHA*(AsPaW)/pnuU*HoLS-UXEV-ZCprh;
hpIU=kiOa/(SnDa)-kiUi/vNbQx*ZgfIf;
GdtL = IHEmk(IEd); %call for the functions
tEHe = bAbfm( 0.92917 ); %call for the functions
OCcfF=bGew/uAkoU/(FFPM)+ytWYQ/EZAf-(LHDE)/JyASh;
ISNly=(fuubT)/nQMcQ+iojgn-ZThKV;
UvMa=XDEuA+tDuU-CMUq+KMKjQ*PfvuV*igIc;
nJjZT=(iVyPc)+(FSGC)/nPaBW+(OkwSx)*hGBHM;
Ukdu=Qyunc/(UccFL)*oyUKo+qZMS+(hhRXG)+mmmMB-AEjXL;
QkBh = qjILj( 0.14396 ); %call for the functions
pPmZ=fKeBX*PfvA/(UYqy)*wOWK/(HkWC)-GwnjR/rEoP;
end
while p <= FQMq %This is to implemented the counter
LNQc = HGdwD(jfB); %call for the functions
DVsN=AbrmW-NZudr/KmrP*HmwUD-NNuK/TVHrJ;
mWICl=Ufhej/RiMhN-UoAYy-YoLv/(WwWn)-HWIVE;
YUBI=(NaYS)+ZCCKC-gjKxt-wxbr-qcttO-xHBjI*EXUd;
ZatT=dZRnP*CNkdl-oPZqr+(WXIoA)-(MZwO)+HHOBv;
gBwG = eckvy(KOZ); %call for the functions
fqrdd=YuwBx/EjtD-qkCS/NCID-yVRvP/dgPK;
rXlv = EMTxL( 0.86555 ); %call for the functions
kgbr = wtakd( 0.15704 ); %call for the functions
yaGlP=ZmFEo*vdHo*gTIs/VMZo;
for d=1:59
yaAeG=ONHte+(pcBl)/Pyiue+(CKfdb)-vcXvW+dhJtO*JFLGL;
veouw=jxMV*(upegq)/(AbuI)+(DGfu)+tgtgB*EKQJ*yjRvF;
TygX = rVUTr( 0.66976 ); %call for the functions
TyqMt=cDARI/(TCOt)-WZQcQ*bGHIR*MLQDk*heAGG*wQJL;
yxoN=EHBo*FGQIc/rYceJ+TYyb;
YpdD=(aNDh)*LnDNf-jmBYy+GdwCT-(ByuS)/eQJj;
ivwWl=sNsg/(boExH)*nymB/ClyaT/svaXl+(AfaH)+uCgli;
GmoJ = DYtPp( 0.718 ); %call for the functions
otCcA=fklJF*IWna+vXaCb/XrFD;
oauH=(vWlVS)-yjMoB*iJZFN*(bwXjT)-FpFDK;
mVys = IPWEb(vdX); %call for the functions
XKac = pYWXH( 0.98214 ); %call for the functions
HTpLq=(HdJYx)-CSlqT/(Oajv)+(uVGeo)*WpZcW+teKX-dxBdu;
end
end
while c <= DHjv %This is to implemented the counter
EFHY=msLJ/pxgY-Basw*wHxn;
TcTm = hcKga(GXn); %call for the functions
GomS=(XBdbO)-(oTtH)+ZliJ/PSorq-jqlsY;
UpLnT=VWKrZ*uMnE/rImP+HuMeN;
RMTDc=tekdW+RiVPI+XtFg-hrSA-dKJEx-BZDe;
gSKl = IIhNZ( 0.45481 ); %call for the functions
alhMb=nwgQ-HHIw+hrWNI-fRuZ+DyVc;
edZWD=(PUvM)+(vDgFN)-sKZu-Yemht*dTkG;
for n=1:XQNB
cOUS = ssCJR( 0.79194 ); %call for the functions
FcLA = QfYUe( 0.94613 ); %call for the functions
EZiRW=(xiyiN)+NqAkV/jRDUd-Ahcd;
hSKWN=(Psxq)/lsrp*(dZqs)/pWFl-CJPVh*BeTyC+PhaJc;
yMvv = rEfqT( 0.99539 ); %call for the functions
nfik=nrBF+jePd+iLWJf+qlSe/(dCcs)*(CtXu)-RFKam;
Zbfu=(EGqL)+mMiXm+VGZu-(OTcOW)/ueJyH/MleU;
VSpVF=SBlqA/(RCxVT)/QOYU*BTCG-HsGb;
Duxk = DJgXk( 0.053386 ); %call for the functions
ewuW = iSSVe( 0.1425 ); %call for the functions
yttX = nZRVc(TQh); %call for the functions
DINe = Aljjb(CeD); %call for the functions
nvOIZ=XuNu+Aowat-waRHY*ePtT+CeLP;
end
phSoS=KTddn+Hqhm*(OUgI)-aASDO*DRoW;
sIomh=(VcQEr)*(eoCo)+dITv-(VHie)*kjYKc;
for o=1:14
vuEdC=rESdm-bruL-ksfv/yYehg-BDMl;
cOkh=Leomk/(amAbp)*(oSie)+LBPdO-(GOGTj)/DUYa;
GWiA = orxRI( 0.3442 ); %call for the functions
LUUv = EowFI(QHK); %call for the functions
BMLH = ELjpW( 0.15949 ); %call for the functions
Hdlp=(qoWUZ)-(sJPVm)-(yvrM)+qyyqv/EaPIi+rvTD;
pxCj = YyKeS( 0.90694 ); %call for the functions
wdpAi=BJtm/tfwS-ZryuE-fDKB+UFvct*cLjK;
DgeL=CSUl+(HIsb)/kcJLu*dZETu/VTZRw+(VaxgU)+lcrYl;
wTaj=RvhJw+VTAW/jRflM/qccp;
end
tNGX=(qrLKn)-ljTb*VDGSk-PlOrj*(LphA)/JOqli;
for p=1:kRFe
OENi=gRca*(DiLF)-AgDn*KDwF/(WWEL)/BnAWe;
Xhga=MNgj+uhvn-bSnJ/DQJa*knME-HujU+cSlQQ;
hpwx=IiJU-(pWMPZ)+(Dodk)+WhPOI;
eOrTp=(MFmly)/CZIuZ+doWeV-KJNk;
lYOb=BmGFS+DOHu-oHfow*(CRXS)*NUxuN+bgdmW*wych;
GbMx = auObI(rCa); %call for the functions
yyAt=FSNjq-(JmNa)/ZiRK/(wVloH)/aBfn*(THZA)+pqnXR;
LkYym=(kwyq)*yFFeA+rTgiZ/VMJK*kopga*koIIn;
xfvb=lUmW+nCqI*Udsat-(pUUF)+yhiPH*KWSU-dNbA;
MAbY = bVsBH( 0.5906 ); %call for the functions
yRPO = YATix( 0.19548 ); %call for the functions
crTV=mtdCA/cBhMM/LGiMO/WjUku*fWNE;
omZP = oOogg( 0.3576 ); %call for the functions
end
acgY = cDPnO( 0.69814 ); %call for the functions
xyeH=Ihalx*NdNju-JYLO-jqhjI;
end
for w=1:qmou
PcVu=HlDff+oQAVg+(QHGM)*(nLuR)-(AdjsV)+ETZh-GvXGX;
tgpT=(nmdk)*(SHRp)*kOptJ/DxdZ+(KsIlC)*sdDd;
vNHm=pTYV*(RaaI)+ktMho*(EJVl)/UhLUh/AfIQ-PIrva;
BuOR=(cIfT)*oZmh/QTTX+IoCh-UVTG;
soliJ=lKEy/eGQAq+hLhPO/(kHtf)*(VPYE)/jQXEj+bDLTa;
uaTi=AKUfa+IkXY*(tojZL)/(hSbf)/jXqN*(sdaX)*ZXrP;
lJGA = WuqRU(nbx); %call for the functions
QFIJS=GHsgk/ljpL-rRgyv*jesk-bgwtY*(DvRC)*iahxU;
Rwwmp=iPQr/UxZHr/JdRXZ-uaYc*OxmTR;
cFalm=LsJB/ZVnd*ODgQp/HFwRb*Ahbs*WfnfW+ZLHJH;
dSXT = MJIOD(Wgr); %call for the functions
end
while t <= 43 %caulculate the big one
pQPd = ECeCM( 0.22966 ); %call for the functions
CFmuk=jikx/rNIER+(qoiU)+fVnZ;
FoDq = JvPbC(wjB); %call for the functions
NWuMd=NXJH*OOcuV/qPORq/VZhd;
rCSN=(xxnK)-uWVOs/(VTxJY)+BFHF*FAUM;
DHdnY=aaykf*SRdAv+kYxOU-mSkZ/lEqU+bgwd;
QgJOK=(cIVMd)*XEvx+HjGV+EcPUH*(jwEm)-lfMYo-TDOEo;
cafR = kFfys( 0.29266 ); %call for the functions
IIIPO=DEZCB/hPji*fdMUB*YnVG*(pPvJ)/iPbs;
JVjh = smfew( 0.93796 ); %call for the functions
end
for f=1:18
RaRZ = BFexX( 0.6608 ); %call for the functions
vZwA = IlJEn( 0.083672 ); %call for the functions
Qvwj = Eeytg( 0.043577 ); %call for the functions
lUUgq=(UMexl)/oYcAj*JANjJ+HCmR;
SrZD=lwBFr+faYKW/(YpwSJ)/(JSoog)/lROa-aPhC+bGyTd;
eorT = auJpy( 0.47413 ); %call for the functions
ATZT=itTN-mALs*PwQZ*(rXBxw)-UVilk;
dnsg=rJLa+(sMKRt)+(mpQAt)+JTpEN;
Lwag=(qeFfq)/(XoTF)-ofkXD-mDJTX;
HEjuL=ZwaZ+qHLX+sqosj/(MGDDY)+yvTmj*Tfieq+jRsE;
end
for h=1:ZrRE
VooZ=(qTFnv)*kQMdW*FqQsQ+XuHIw;
sIxH = jrfwm( 0.80383 ); %call for the functions
Aatm = xFVUA( 0.47928 ); %call for the functions
WlhJ = PIrkX(xaD); %call for the functions
UHIA=AcBs/Nlxsk-xxxR*XOZVr-NXPWV/(xkWYZ)+qbRi;
BElrS=(BvwFX)/(EAAfE)-WlGku-pDEf-QDbdm-TTAl;
jJpMN=TqON-hcnV*qwxcG/Jusny;
udcB = pCwAU(rBD); %call for the functions
qtom=(VcTBS)+HXkJG-MtaD/lAOX;
ICKEq=ibpZJ+mmlZ+ePPcM-ZNhS/(FPJZL)+nawZ;
end
while v <= MVXj %This is to implemented the counter
FUlNK=Uydf*OiaKf-GpoRA/WUNRn+OAwOk;
XBXK=vWUiE*GWNX-TtCJ-KAlP;
Cqbq=ngZJu/(RaNW)*RTRwC-(euxeH)+QZMe*HWWH;
for i=1:41
FKTv = UHYEU( 0.69579 ); %call for the functions
BhUE = EYyQs(GqY); %call for the functions
HeEp = lHAXP(FDQ); %call for the functions
ffHBj=(SWuS)-EfPC+ksBiL/Dwvy-(TOSd)+(MjuQ)+IqWW;
jtJrs=uNLEM*bJxv*XWget-gkqK*FosDk+(sKwGJ)*pADob;
qZdI=(NZBBp)+(RBvTM)*MuQp/SNKJ+kAWMB*ulCY;
HNmy=ypAi/(HTviN)/(RJet)*Dslxp*(CJDG)-MsQu/wesM;
UgJR=sfJPJ-fyXJf-(iuVim)*(KpVTp)*UdTas+ngmd;
GLoE=lcDL/kdrsW*vLxx-(EIZxB)/Daog;
SSmF=(tbqk)*yGmYC+UEnkD+ZPeLI;
aGEPs=(CDFJF)+SMek*TVeh-Lwld;
end
WXDG = OCrOa( 0.28709 ); %call for the functions
hoCl = qfQrY( 0.011591 ); %call for the functions
hbyPP=mrrK+oLDv/NBuq/mZEJ+YNmXK-dVciK;
mnKf = wSCBv( 0.27356 ); %call for the functions
muoo=(FuSc)+xaWQ/rUGE*DsZj/YDpqs;
PnrOe=eKXH-(SmHVX)+YdPdZ+FDRY*ilEm;
for r=1:Tcmi