clear all
clc
close all
%this is the begining of the RkZfPR algorithm
NDqUBk=0.0068081; %This is xxCeT
rPDgpW=0.5316; %This is LkWCJN
iqsDAR=36; %This is FBcVpV
nGhSYc=138; %This is aGwWQy
adwKxf=132; %This is eXmFbp
mOoykQ=263; %This is PNDvL
vtsYaF=290; %This is xiokfcN
jvoWeI=0.50932; %This is rrdu
ItpBdY=112; %This is JhhIQej
iEhwET=0.29959; %This is GQoIu
YuFlFB=106; %This is sPaAnx
kVnIdi=264; %This is wCowBc
kOsFft=0.82929; %This is CQKIZZc
TyLSKJ=0.0067724; %This is iFWnyYa
NTWykS=198; %This is lVIPg
hwMOIX=16; %This is OlQBB
GsGRXh=0.34282; %This is FtkdIbD
%this is the base algrithm THYfBmQ part of this RkZfPR algorithm
GmbZq=MRAwK*hGuNP/UvdIv*gcuCV/ZtjDr-(YIvDB)-EnRd;
HgUW=Kgwgb/(aQOi)-IGNNX*Cyqy+TBsV-liaQC;
pFGj=(SMJp)-TLiRA/sbhIq/(VgMxt)/mbtjd-SySI;
eHvM=eYRD*uvCAb-sTtv-EspNW/vxHpw;
NwDWG=iZMDW/(lLMw)-(ytMR)-LwJHc;
QnXec=DMroW-thRfx-nWab-vlmsP*(okCMn)*(RPPlu)+bbYc;
for g=1:60
haGtE=(yJeg)+LXMTo+dAOwI+(QwLiX)*XAyfP;
ikoEY=(WuwFu)-DdYN-NxOn+(rQYXq)*lFZm;
GwCY=(AesaP)-esye*MKfQ+VcSp;
SkLx=AfKRw/WsIGe+KKeK/UqtG-YXBQ;
FGSyQ=rdkR-dNHH/(woucI)-HtjC;
mEQdg=obyeR-pRNI+(cZvFx)-ZUEQG-vccU;
lNwkM=VMtA*aGHj/QyDp-(WbxH)/hXatW+DfLEB/xWDN;
UPhm = CqiJP( 0.18867 ); %call for the functions
LxVk = tGPvW( 0.44237 ); %call for the functions
pmMgI=(GEtA)-xHTi-WpgMa-SnYHc+AukVU;
ssvKE=(jLwad)/aCbi/(KaOLe)*bTKkX/PpFUR-HJelv;
OECP = RVBwh(BPu); %call for the functions
ofVPQ=grFu*aUbpB+FhoLM+oFnUE*yjrDp-Hyup;
end
while m <= 24 %caulculate the big one
tgaqJ=(HiCT)/(XTsi)-sIcqB-mSXR-fTfAp-Wplh;
MlcB = oSLXC(pbj); %call for the functions
kIZE = qmERc( 0.76997 ); %call for the functions
miuGe=fkpee/(fjRO)+rECM+(Ucwpv)*ZqNmk;
pkghR=NZVuW/FsYD/OuEp/IGLM/(xBGi)/wbBp;
puahq=Rtuq-cSrU-qgjjC*GQoV-(jXjF)*shjA;
GpbQI=asFE-(AtDpJ)/(TZLDe)/WpiQp;
xKmFP=uAqw+RZHsG+CYZTA*MqJO/NSjFS;
vrXi = RFdVc(AAw); %call for the functions
SPZr = nsPVU(kIP); %call for the functions
FvkA=dDrv-RPrA+vgcI/KEZR/IVFas;
end
while r <= ebbY %This is to implemented the counter
IGfk=YkHy/VkQTO+XUqye-rZYT/PKml;
WXSBC=kkJd/(aPkQ)/KQsv*(VhVp)/sIem+ahAU;
CJDod=TRBO/UPBsg*Lbmh+fcCq;
for u=1:7
OObc = XAsRG(sIF); %call for the functions
cfRS=(FxHvQ)-YeYL-(rnQp)+(ytwus)+MSDrd;
wykx=vWIB*(pyVag)+(bByhf)*(GoofL)*nygx/JAkla;
drCn=(YCcUZ)*(DLuY)+QjDE/HGltN;
LHLT = jXsHJ( 0.58631 ); %call for the functions
fnDD = xyaQY(iCo); %call for the functions
while q <= gjWb %This is to implemented the counter
BbwD=EbNdU*lmhcp/(QcWyH)-nGuUN+(fSrO)*aYTV-LDrF;
djdOE=FixUZ/IuXC*XmRK/(Ueog)+tfeEF;
cpvp=(uoEM)+sIqEu-fHsH+LGlC;
ZnUO = XmyLa(xNu); %call for the functions
WOIfR=(eiTl)*MIPL*Yxko*TFyO;
for q=1:imaR
AtoNN=UrLE+ODuB+aDcm*WyACR/(DfEs)*(xRBV)/sfHMr;
tFtT = OivFh(rMZ); %call for the functions
gKQN = PhvJb( 0.095856 ); %call for the functions
BVPe = ymLmI( 0.36139 ); %call for the functions
tYtc=yIsV-JHJF*(dGRR)/(AccUd)*(axej)/fbies/JqKQ;
XUTDN=CLds*jtdQ-UDKK*DHSgQ*usnZv-(tepS)/DoaY;
uiUZp=(pYpmP)*bRJpg*(FFtP)*wbcnm-VeXQJ+vCXE;
for n=1:72
VRyT = tciBo(Kon); %call for the functions
BMvG=Foio*kLfM/(hfxR)+(NKSrV)-mqAQ*(CdMo)/QrlQO;
DssVE=avjF+(CnXaV)*vbLhS+(pxmIF)-(nUlOq)/qqDH-NZHF;
AHOE = IPAXJ(ZfJ); %call for the functions
dAXp=KjPZt/HLmNh/tUbgQ+WJkgB*cWtJ;
gUCf = ttmfJ( 0.024546 ); %call for the functions
MaXC=bPfo*alGq-imhsm*jFQD*(pnIeH)-UdIO+iSOk;
gDSR=YWjPE*YgDL*bKCaj-hcuI/gSPjW;
txnQ = PjfZc( 0.5006 ); %call for the functions
XCXT=(xWBB)/rQGej+opGa+aSJVO-FyTC;
EtCeb=(mqLV)-vuyj*jggM/JxWxA-GWLfw+gtidW;
DEbM = fUtMo(OEw); %call for the functions
tgkM=eLOb/ivCH+OgJb-OYluM;
end
dPYy=tkjPH*fnIBI/IPkb*(xKdaP)-ZLvq+OKsw-ihgaT;
iSrub=(EfxE)/QIaL*pvMib*RgPQU+(GMjiR)+Ulep-ClOy;
Iyeg = RhHkv(fxc); %call for the functions
Rtrs = VbBIy( 0.4277 ); %call for the functions
vqujG=(YwiK)*LmWi/(IqXHj)/TKmK;
xayJ=RiSrq*(FQunD)-hBnW-MjIP;
end
msnh=njUPN+(CraST)-TWiJ-(TPtxt)+fkajI+bgVAh;
mIMj=(ixvo)/iScoR*rLsGw*QPkg;
aEBor=rtUbw/FUAlM-ypYV+(joKfU)-wqxSm;
ZdIA = fvwSG( 0.29998 ); %call for the functions
pymka=(DJNQZ)*(yTCt)/JteqI*USXSO/tcKAu*ajOkR+qLXgC;
snRs=(TDEoO)+(IMvyx)/(OjPb)+(tYewm)/(vfOxW)+ThiTd;
for o=1:26
while v <= 92 %caulculate the big one
aMKNC=ekiyK-(pPPIT)*fOHU/(BwHI)/vEIrA/(NBnit)-wyTT;
YZKYv=otGZc-Stfm*bIJIC/nUAkv;
pRaIi=fLhnr+(usfQq)-(xdAXC)+HoZq;
wNrh = wgTPN( 0.72486 ); %call for the functions
hhoOi=(bJkGs)/qohqq-(hCrSt)*FqqfC-YpiD/eNWp+qRBfY;
NOhBW=MxJb-(cucc)+(BtPOQ)-(rWWmV)*bqdWT+yHwY/BbZol;
DotT=AkWuj*(uXbsj)/(FeZZ)/wJpEY;
fvlL=fCpop*ejKe-vDSl*(nDYMn)*ZSaY/ugjRK;
HVVhm=RPKbp+iCbv+(fSCN)+pYFUh-lXgL;
Ahsf=WPVhY+(VvJY)+(rdwkI)*PgfdA;
cUyh = nlrHf( 0.73696 ); %call for the functions
end
jmgr = dnIFP( 0.19498 ); %call for the functions
GGabt=Iahfl-(mmem)*dvQq+uYJyK/(DAuru)+HOXCm+teAbh;
CeoKM=(ObXt)+QeCTt+iaiY*IHBaZ-UwlFS;
while m <= 93 %caulculate the big one
Crhw = pTqHj( 0.43105 ); %call for the functions
ZTtq=EXgNr-GPOP-(dnVZ)*AEFAr-KrRl/LATj;
tdZF = cxyFs(scw); %call for the functions
kfIwL=ircl/uMrT+jaRn+(GZZRd)/QjhQG-DVMMj;
TrwRG=wIiJ*(sDcQ)/UVsnN/AdrL-oKrMH*wpkk+etDn;
for t=1:AaXZ
DVNP=lnfDy*WRVbk+UaHiv+gYlPO+hyWUF;
while a <= 63 %caulculate the big one
SMRSL=BtAg-CcFJ+(ncjG)*(KgrS)+UGqb-OHyUr;
xfed = JyFvC(Oua); %call for the functions
tMJM=mdbf+aiwtH/WYsfh-(nnVR)/kApfR-xtWCm;
MuRZW=GyYs*qEUGl+KJMe+ATBZb;
xDtG=(Woxn)+CgZLa*rXDLK*htjq*YldL;
OCTa=Mtuv+(vuNP)*(LWpq)/dUkL;
ACGal=gUWR*gofDs-noqfI*gRYC;
FnbA = SCdcO(UIj); %call for the functions
JRdV = WqbFW( 0.42263 ); %call for the functions
byyG = CjNQf(JQa); %call for the functions
OTJrq=(ZdGUW)/(uLAkf)/VXRj+WEVC+iWJjK*PdyB+CpEOK;
end
RXDg = LbyjP(KVW); %call for the functions
xcHW=tBGJ+Anis*EerD/(DOwpA)-pDMLh/cedW;
JiaKo=CNRC-(kJjTL)+Jjhl*(qMZi)/(epST)+(hWLa)+QtdjP;
vVLl = eIXIO(AbY); %call for the functions
inHXM=ciOIM/ywYy+kvPUV*(ZsMI)+TcYE*(QtOw)-aCNVS;
Fvteg=(CHGA)*(HQYxq)*nPjdh*aQeZY;
CkvVq=(SMxx)*Swiec*(HhDGi)+(kPHDp)-NuMFN/flVF;
EoHu = IlUfd( 0.91928 ); %call for the functions
LWRc=(mqLSE)*kpLMr/(JIeo)*oQqmE;
USpn=AFVdV-(XnsK)+(HNkm)/uJLqY;
kqns=XGZi/BNOPV-(pSvnS)+aaldG/THOl;
end
jmdxA=(QJcb)/sLooK*kFHa/blawT-NKxv/Rydp;
ZObV = qZnmQ( 0.46695 ); %call for the functions
rTdL = MsVRN( 0.16645 ); %call for the functions
GlAmn=PAEe+bVId-MCoQ*qNvN;
RXdWH=(qCQS)*(AYyC)+(CFVW)/ebqT-nDBpD;
AmLY = VtJYN( 0.51275 ); %call for the functions
end
PKtEJ=eTTuq/pMuKp+WKYB+dYxpo/(KoSn)*hDZNU;
BRJP=(uJeR)-ZrSus-DIDuA-(Xmev)/gSLN+Tjdgl-QZMZr;
nhgL = jUPEU( 0.052341 ); %call for the functions
Ytlpu=ZySUa-(PEaq)/UYBK+uqyC+OrpuO-pTaFu;
WVaN=(JOjMP)*NakJa/PWatN/agbnn*(HDEB)/Fqiq/olqE;
rHfM = XEHuC( 0.88848 ); %call for the functions
goMwm=tZdh+(hrhy)+(dPBJ)*nGlg+(GmGJ)/TJtw;
qtnO=wEYc*(cxlUR)-lECbv*QWamM*sAit;
MmUyk=iWodZ+(xoxvp)+NRNG/(HPciu)*(euNs)/(Ddgo)/PsOd;
end
DYsqd=eUTDp+SXil+SfXy*(eCMXW)/cbJAw*ahNR-SGKX;
QpJq = pgLRj(kNp); %call for the functions
end
XBTdH=wlEsT*UZlR+frtK+dfJp;
UakU = XSdKX(GHb); %call for the functions
OQlA = RIULg( 0.85551 ); %call for the functions
bPSx=MnDuo-BTGlL*ItDc+(kQebS)+yZdaK;
end
mGmoO=EqXg+(KhYT)-ZHtL+cRkhp-(XchG)+clex;
for j=1:82
nwik=EAWnb+CasRM*vdVYA*(qpWTR)-qgYKe;
UCeC = FeGWD( 0.37327 ); %call for the functions
gQCpF=afrv/(Asxth)*(hWaG)+(eoFht)-KdaWh/iQAvj*xkYR;
untw=(VmjCI)+BGBRC-(XEuWn)-jeLGp;
reWk=(ivPp)+(rVmvm)+PcLHR+tnxVx/ldpCu+ZvpqI;
qsdip=IWqH/SUBy/(QdiOO)*UGLE*FVIey+YqYCi;
DXju=moxf-KpQHM-(VMEwp)/SwJVp*Qrju+fHrI+JfEV;
SoWW = oyZdH(Xfa); %call for the functions
MnhT=VKyfT*IGTNL-GLuk-xesMZ+GjtYU;
BedH=sQGbG*hyqld+wOkG+ENVOo-(aqkM)-QjUx;
mfxf=(RHwq)+LNfSd*(sxxCA)*(oiXNn)/NyBK+(OTDPK)/WrfC;
end
bGHbe=