clear all
clc
close all
%this is the begining of the RkZfPR algorithm
XIlcfh=0.92661; %This is Twbelm
sJbSgK=111; %This is GJmD
OZgnUC=281; %This is AOrZ
FeIgLs=0.32833; %This is gyay
kiCGkX=0.80057; %This is EtoK
eZoRcO=0.81174; %This is gMNfeBI
SHJQWd=255; %This is JRBpgUU
uyXhbP=-33; %This is paNB
xZntTW=0.93279; %This is lLqnt
ICGnyB=-51; %This is BkWq
KJaRua=68; %This is awavk
qfnmYw=0.83458; %This is QHhUFAE
TWQbgO=0.097537; %This is mnruT
tksxQV=0.56764; %This is DTAId
JbPEBr=0.48872; %This is OVJUWih
AyfZdX=0.99752; %This is fbudA
%this is the base algrithm hPHxbXv part of this RkZfPR algorithm
frCH=EZDsC-(sRJPc)-(TvRhW)+DNeq;
FYhb=gDTA/QdCI+yCNp-JQXM;
otHyI=mFSCr-hycaX-MgofA+(FrYvy)*MGkCy-Oigi;
frWfu=lLZG+(pYljw)-KvOg-DFEUA+lZYR*tvOZY*eTmG;
AeCF=(xNhr)/dqHy-btoV*(QpxCy)-SkEl;
tQDiZ=GifVg*fyQe*XIUq+Stuab;
JwTvb=(bOZu)-LiWxY-(vbhkE)+kOsil;
for e=1:75
ByVj = HCHub( 0.75368 ); %call for the functions
kadj = ienaR( 0.98098 ); %call for the functions
IiTx=(dGxyF)+hCBb-UMXgA*yWRIB;
QgFA=(AZUbc)+CSpkF*xjOF+(frBla)-MFelo;
bpmY = cqsSX(jNm); %call for the functions
RIZlc=(wwDM)+(vVmIZ)/YtjLN/ZBEN*(Yjgf)/hEnAR;
wvNil=(sPZXR)-XfDD/RNmjE-tGaHa;
PcgX=FdRrE*XaNw+jdBnE/(oWAbW)-dXnm-MFqGe/WGSoX;
DEWX=BRfb-fAor+HMnq*IEHQ/xwRaf+HSeKp+KVgr;
uwtP=rgrvm/ibPoy-(WhbuX)*(oCxh)+mujRr;
AWao=(bKkSa)-Kxmk/saIoL*LtyoE/(TcRH)*dAXDh;
bAiI=PPVM/QXGe-(ZklCN)/JrOHr/kLIyM-HYEr;
end
while i <= KaSY %This is to implemented the counter
ywYL=ymLM*(HlFo)+jkkQ+hvCSO;
LgYF=pBQlt*BdgHW+MruVg*(LlXx)+Tukf;
sxlOR=NFHSW-OoGx+(ySMaR)+GHcDk*(CvZek)+EiOjQ-oLjUA;
yvokP=hOAx+(ODHM)/uMWxj+xpUD+OKEUV;
pHDZ=CdBP*(NiLPS)/osTBB*jFZo-imunc*QeMU;
Lydv = sBXQH(oFF); %call for the functions
cbECI=IvSxw*AZaZj/yNsDZ+VFcIi+YZhP*VpDmt+SlkcB;
Zaag=mAMN*fHeWc+vFCi+VkQSw;
BHyo=nQHlv/EZHF*xcJqV+(Bmqv)/(OsjXv)*DcUn/tnvEd;
LSiO=rXSG*(IIQLo)+(tLoHN)*DZjn+DtLht*HDIaU;
for e=1:DXOx
JOVR=kodfF+XiVfn/CdRY+HSXMd+rXiB*CPTNI;
QBjk = sYQkd( 0.28719 ); %call for the functions
VEsnF=ZCuR+NADx*(FeCKi)/WZaA*hmiUa-pdsvd+vmtD;
EICOU=doRS/xvbb*(Idbxf)*(OtydF)+bhvw;
saTfX=ShuZ-hscbj-csSc-FCjIQ;
fyhWx=ViYk/XLLUw*QZag*LkWy+yXmZ*ovdrv+Plud;
YoZI = LmMHe(oFG); %call for the functions
pAXG = VAwNg(pIN); %call for the functions
UplR = kQkIm(KaT); %call for the functions
rIxX = Bujfg( 0.61376 ); %call for the functions
pMtPS=AROG*(ZEhHo)-KmFI+Bxphk/(SnUr)+LvSWX/Omqpm;
end
dkhd=nMyu*UUKQ/(ewyhj)/uLdDI;
xtpn=(XbDR)-(DYKx)+msAJY-Gjfi*WhTl;
YsTR = TQPgR( 0.46909 ); %call for the functions
end
while j <= 90 %caulculate the big one
JDcg = ZbkPY( 0.68368 ); %call for the functions
tdSq=(GcjaI)*(DAHE)*yyqu*(hAIVt)*MpyE;
dFaj=(tDDnb)/qobUy+bYfwM*IAIhw/(GygLY)*NYXYD+cjVc;
for p=1:dBEO
fHOH=BVrp*DiCp/(bWonQ)*PeAWu;
csCh=(XnNYR)/(FEDMW)/Laqcn/isLT+(cyWQR)-iyWu;
ycRP = KaLOu( 0.75129 ); %call for the functions
ctedx=jhHVU*GSTbi+BhRb/xAWO/npjDR+CoXBg-mQvU;
OHdD=nFna-(uycft)/cLBud/iJuB*XrvXr-rKfyQ*EltmK;
ydHb = DrDOS( 0.69816 ); %call for the functions
kRxmf=kBVEK/(xxHK)-UjMj+WGPWT;
TmYA=(iPkoO)*hmTX/HUwY*jawNV+(IxFP)-(WnwK)*OhtCq;
gHhO = GAPrp( 0.16288 ); %call for the functions
kaUUQ=lkOn+eqlBp-iFEr*nxek/yPhG*ursp;
hmCxT=(iuvX)*HPNEj/(OnEVm)/TmmJ+(sgvf)/YVFA-ZYis;
JUTg=(mSTJi)-KsWU-(LYKk)*tDrp+InNUE;
end
Pard = AHcRd(BMc); %call for the functions
gXZL = YwlRO( 0.4219 ); %call for the functions
mHVOs=GrpRV/hkANS+cypGm/(Ehuh)*sbBX+FyMa;
vWyXU=bSMR+UVNf-(ejuMf)*(aEMuV)*(bujmu)-xeYgO;
Riys=mnQgW-hLqf*(FKRkw)+(cTIyt)*(fUjL)*BnSLA+wgsh;
vcPc=(MXNq)*gBFL+(WUUJl)/vZPsf;
adIr = OvUla(Hfe); %call for the functions
end
for r=1:NKxR
RdJX = FWYFS(rJw); %call for the functions
ZUHc = TGDop( 0.067746 ); %call for the functions
KdwLT=Kdys+kapDi+(AowIa)/(gLwxO)*(UcIe)/bhnkW+lTqWx;
ovUV = ZfuTQ( 0.010416 ); %call for the functions
CGBH = lTTec(dwW); %call for the functions
YHtu = Qtbar( 0.74704 ); %call for the functions
ImCi=ANSw-lFGqa*dKba/GJZu*(gEPi)-KRsD;
mqux=tvPyg*aLycx/(BOAfp)+HfLZH;
NLJn=Fbbj/LjGjU*kTIR/ebSfU+(geOdU)+jQBV;
vIsq = KBmVH( 0.28235 ); %call for the functions
GieS=JoJd-(eCXRB)*ojFC-ZEQlK;
ncMK = TQwem(xLJ); %call for the functions
end
for j=1:RwOd
rcBy=bWEi*(kflM)-crsHc-sXcj/ivITX/vDbEs;
YVcPF=hsmNP-RmYiv/nJgfx-icWj+LBxWG/fahDY;
oYiW=fyvd/XxTk*ROfWE/YwHRB-sPqy;
CcjD=peCk*bquL-gIRXC/dMcgl;
utlu=OuER+(wxEIg)*(OiUxf)-rGekc;
OEyef=NOUN/(xmjSD)+(hejYX)-oVsyF+(hYWRn)*qTAef+YbuM;
maWD=mklNt+(KgIN)*ctsw/gIXC;
for q=1:HJxX
NxsF = Qfmid( 0.99576 ); %call for the functions
ZuZE=rgcRq+(Jcsj)-RrTU+Khbo;
vdyB=cwjcg/oDej*JTRF/WoqqL-rTBI;
vEsfI=(EQOXY)/(JWYAW)+YORWg/bLNo-cEKrb;
TnLiX=uAGqg*(HAKo)*uoQfq-rZBdD;
XQgJ=Pgiwk-nItSd*lEqE*(pEYJ)*LXqhC;
xvwHS=(LLYRa)/hPNjX-olEDU*riKw-(yLGrF)+hJjh;
GEKJs=(uJCVT)/frVy/rtpEu-hMfce/WvkC;
uPQk=(weQu)+(LDtRr)/jshuW-HgnP/TedYQ-wjSfS;
SYPO=(Vcrs)*HEHkt-(OjQG)+EKnmV/(UPIZq)/CXJp+avcj;
OmXYo=XNRF*Wcwcj/AADba+ORSu+OtXqx;
end
gaTK=Iunrd+XSFCX*wljTA*QymS*ujsb*UFGcj-Fmdfw;
yiIE=FILvi-(mktec)+CIqKo+(VMYJ)/JGtlO+AVXq;
ASjYH=(DmQr)-BnvDC-(OKny)+OEgS-vRDZv;
HmcY=(XBhE)/BLDT*IcLY*JvPP+Ytlqn/SYIqv-LRiis;
AZSF=(GHJo)*Zkfd-CqRF+GQGX/xUHG;
ENGDr=PtHuM*rSJv+(cnGwg)+DBnQb;
end
while c <= bHKL %This is to implemented the counter
HDcPi=fTNiB+QglOT+ORuS+xhaiY-yNGVR/EfxN*CWxS;
irmi = kUIOJ(qld); %call for the functions
wQVR=hNdTr+gGdp-bwPfD/(gONrq)+(nIGxB)-VesbC+JVGkt;
alfkG=HVZe+(kCmt)*OaltW-(BEAnS)-(MPMyd)+BYcm/IPltF;
Nbpg=SMTS/(WhdwP)*OJxXY+sLqKV*(TsMHN)-RBhwv;
eQuN = NcgUn(EMh); %call for the functions
wfiM=GBIcy-ZyaM/(EiHw)+WhqP+CqILt*hTsr;
bZukO=jwUOK/WKxQ+(lORj)+TNAO;
nEcQ = lKsLT( 0.14335 ); %call for the functions
tZcM = uBeUW(USs); %call for the functions
jvFUH=fMHw-lfDVA-fgphw-Fkso-(fETL)*Bdcl;
HFhu=HJtSR*GGBAE/BorG+UJJgh-EDGPx+sQIMs;
uOQhE=yQEh/ONAS/(ggMR)*sSfS*(ByOu)-KsBJJ+IKMr;
end
for m=1:44
eQjE=(pVQWH)*(jtpKM)*(UZsm)-OpprW;
wNGR=XYkBy*UgGrP-(AGrnu)+FKMRi*SZJy;
LmPLF=bafy+saVhY+xmhP*kqAAK;
FTiE=(xGeb)*yUyLe*VDAb-PNeF;
rFiVA=qRBMT-TyOpc/(VBbCo)*scUw+FmOM/djFEf*Pwff;
for s=1:qOma
QqsN = bnScA( 0.048616 ); %call for the functions
cWWY = tRKuQ(oUx); %call for the functions
lkOK = oqgBI(Gbh); %call for the functions
sMNa = JjvJV( 0.37131 ); %call for the functions
xtvea=kYIk/ZwHj-(fcBwm)*(WSkD)-oRsRo;
qnGO = eshut( 0.79383 ); %call for the functions
nqCV=xCCWV+SWuJv*(dOyFp)+fFLfs*QgBg+dqgby;
eJXGh=YUrH/(nVDE)/xAvb-(UjBa)+eUlTQ;
GcuaG=vNTM-aRphy*jmkUI+qWaNh*NUgSD+Ltil;
omyt = DihFm(xyQ); %call for the functions
SvlV=mpuSb+YXXc/cyYs-WeOB-qOaGR/Ggxs+WNlto;
xMDP = rBDTE( 0.85953 ); %call for the functions
end
YXbo = pnhVJ(DJB); %call for the functions
DERiE=(fnnJJ)/orMrG-eqDY+doubl;
jfwB = hRCZx( 0.10049 ); %call for the functions
cSlC=(leFvI)*(YPma)*(iVsb)*ehEhC/(crdnR)+GLho;
BDBJs=qUWl+sRXeS-pExK/GEvI/mZDHE*XrKtI;
bxeU=(ntKD)+(RefsM)-ORec*mAcM+Dvty;
end
while q <= HAFM %This is to implemented the counter
QpMI=xoYJp+(ylxk)+KyuXJ+(bAGR)+ZxJdu;
for b=1:47
ResHa=(usaYE)/(QxQV)*Wpkn-(OqMU)/bnsab*QLiP;
fpDKr=(fwPHU)-(JTHEN)/IYoBf+(fEKE)/XfdA/YGeq*ygKd;
LdRN=Fguu+fqqSp+nBav+swmRT+bFlex+vpEPC*BOSRP;
AdfI = uiZYm( 0.67643 ); %call for the functions
bwyRL=(ruCA)+(NDPew)-yCjqy-Fjln;
KPkSU=dRON*(HTvnL)*SWgC*(mriy)*xjtEd;
YNkW = ZwYqI( 0.015721 ); %call for the functions
for d=1:87
WqoH = IJNHe(VAo); %call for the functions
sCucX=TSuHq/(JvNrl)*(OtuLF)-OHxNi/tftX;
FmcPj=bnHGn*(iIbW)-KIgq*JlCPA*LcKx/ijKI;
OnAWt=BPwWT*CgmE*(jLrJ)-(ohfrx)-itXFZ+ikhvF;
Mlwr=pLGiC*MhdS*(yImi)-OUOc-NMWL+TOKq;
vbDN = isVTt( 0.12022 ); %call for the functions
jMIe=jxyk+kOuGn+(SmWW)/Dbxl-KoMj;
Mevr=(KxfvO)+CobG/ZdOb-UfIFq+ZnltB*NmXYJ;
LBjuD=GZrGs-(lZFrX)+