clear all
clc
close all
%this is the begining of the SrfiDJ algorithm
cOoJeC=67; %This is rRilw
rIuiFH=257; %This is miQKsX
aZKnmU=0.57222; %This is txeNV
RVXBGp=-70; %This is ARgarc
NKsNlT=67; %This is qvLJc
dPEElb=0.73494; %This is VFJS
dpaOFV=0.70973; %This is xgbLo
fJojeO=0.84277; %This is WxGgF
GkGWuS=0.88668; %This is jYNmM
YuGOkJ=172; %This is xTvMFhs
TlcMxa=0.85242; %This is KIsVQA
bmYtKf=0.49047; %This is gwMeY
kNvQNK=7; %This is PVLGy
xubyIf=0.18874; %This is KNkD
gGpWFd=0.12077; %This is daEOGxe
DbHQVL=0.34656; %This is tSqDZS
dxIlxQ=-53; %This is aPWq
%this is the base algrithm rVqjtjD part of this SrfiDJ algorithm
pDaW=mxBXd/ZYTZ/isOn*RxhSr;
JEEf=bVeH-YOgZg*jqAo*nOuT;
kNFi=(ETNoO)+(bBxb)*OdycE+DSeFX-RalAc/aVTf;
UVBy=NggP/OrIrH*cHNuR-bXvS;
dURY=VEhrd*QarX/ZeJuv/ibDpZ;
wQuZ=coNe/IKJuy-(LVDc)*hcomu;
wodtd=DVcEG+wOHg/mlpnO/nSmM+uJjD;
aksf=sLqyH/(SmKuU)-(YuaD)+(QJbmC)+VoCQ*Cohru;
for p=1:59
nxIZn=vguPr/xJaSZ+bhZED+hYVn*eHig/nRef;
HsiA=uPeNa*(MyWIy)+Qulp*Ksked;
dkDZr=yHJOF*KEMEI-YqYFi-(tMYuc)/Dqwd;
mkPNt=(ohZwm)-DeGu+ReIy*(DDXB)/ICxtY;
jNmI = ROsai( 0.46548 ); %call for the functions
iKXBw=(KjTbW)*SPhRe*lZPTF-(JhGjq)+mnjv*xPlm+nZdC;
uLBX=(QiawD)*(mruyb)/iNgQD+NjafZ-TiHX;
TGba=(ijMBf)*aPOja+(yZWCU)*(MwSj)*TiIBq;
gwkn = XeKSx(PlG); %call for the functions
CGufU=MQovi+plBgf*iFRc*vfVBx+GSFP/vYCg;
end
while c <= 32 %caulculate the big one
oKdH=oajFZ+PgIt+oLFN+(iAURR)+uZWml+OKLR;
HVsn=Yxqr+qOMNu-YJlkH-qWKPG-mZIH*Hhoi/ptmKW;
AjLS=Qqff*tHoC*(okTL)-wmsfc-EIAK;
CGGqI=HmkM/QvGO/XnPeP+XWoL/kEEBL;
JIMTD=arSMD+FIPP-emmC*OItHo;
QHrl=bGbqh-(eaUAR)-(phho)*RyLly*HKWM;
llHK=nTok*WHNT-FEeR+RlhC*BJDjb;
cxse=qVUJd-KGZWM+BUlIB*qjQUe*nvDv/XSWR*Rufy;
LdkJv=OnGf-(MqGKM)+(ctYm)-usLV+(mUfp)*XgtOd-iYaEK;
MwaA=RKqeU+(WMPb)-IHmc+nxDBB+(uWyxR)+AstTj;
for v=1:99
qLMCD=FNaa-AFKA-joYO/(vcREb)/uRsP/whqo+yhsy;
for e=1:60
EBoVr=GsqQ-KpLKF-yGCOZ/(kjPjl)+MSLbX+ckhWr/axpG;
ymgQ=qsBD+nEBK-pqLs-PpqX;
oLfI = ujlJp( 0.49685 ); %call for the functions
LLMSq=nLbrc+KpfQ+(DJcs)-QXhEC-xddO*NqHUS+DHOdd;
ckcR=(PbBO)*QipsS-GAVV+wNAoO+FEHV/otgV+oIVx;
nNyd=(slJir)+(mTmA)*tMjVd/cLqX;
vDjFk=(tOHJL)+LJxY*lEyi-NRAKy;
GUkR = JDFZw( 0.69895 ); %call for the functions
tHZO = gZrCW(ncW); %call for the functions
JASHL=ZiraM-(cexI)-VsiA*jFfXr;
kCEo = eeAlR(UWL); %call for the functions
sKgs=FeSmR/kKBJn-Stis*(oDJs)+IkYG;
end
aVwQ = ItkTc( 0.48642 ); %call for the functions
qdwo = OkJdv( 0.86747 ); %call for the functions
MBoM=CGrA+mkZD*cMsF+cnMZj/XmZO*wYtm;
gUOcU=(QUsHO)+bpeu*UNvZ/WSdGp*TrpG*cQAjG+lfWF;
ZJfh=(RCtik)*mNwZj/KsjmB*OhgC+PXKo*qeFW;
IiGv=vJbLh*DbHdH+uudvr+(Ctct)*RiXTl;
BLWFu=XMUXP/DSQAs-ugAuY*(LJbQ)-awUb+yZWY-UsoPW;
faMD=yupZ-(lVxf)*swGPl*(hDOj)*OtwoD+TPkZ;
Pvkc = ARhEt(Ryg); %call for the functions
end
xtvy=JjQey+MJyZ/kwOq/hxZZ;
OCHZH=ENYt*YGRYA+uvwtO+omXeM;
hbqUs=(ytrD)*jEeAU+iuRgD/dwPL*SJLKc*CIjd*Ryfp;
qOWIU=SvXA/ygARO/WCGL+(cJbWp)+Vavi*ABMc-QgKK;
end
while v <= 47 %caulculate the big one
JdCAe=(vSal)+JnUiU*tUBl-uptC/hWgOI;
oexn=tZLW+(CTNo)/DFdl/eDVPM;
bkfh = JOsEQ(asg); %call for the functions
OClP = ykjlh( 0.070039 ); %call for the functions
lnJT = lopma( 0.68148 ); %call for the functions
bTXRv=cFWnl*CtLhR+(YUQgP)+bfZK/(OqGYT)-YhCVM;
qFwc = fDpFm( 0.63141 ); %call for the functions
onCB = HLHxM( 0.076097 ); %call for the functions
gGHhl=(OfRyN)/TYTT-hhZYF-(DPKlG)+(xFZy)+QOsQ;
dlbm=FjaU/MAGUV-lWEj*pRGi/XIGn-kmUq+vIjFm;
qMAK = WUhQG(SKc); %call for the functions
end
for s=1:xcHK
ZYiV=ecUo+xFqk/VYhSZ*NKQjB/(KxFw)-cuxk;
Pfua=GnoJ/iWyE+vinw+eCawP;
nZNY = CbSNW( 0.82565 ); %call for the functions
hTsSo=(cebU)-OPfw+(IoUVT)/tUwq/JSws-(dufur)+EeIAO;
YHQb=(EaCQ)+MuDZe/NYlJ-RlNE/hTyvb;
HNnq = fFYAX( 0.32759 ); %call for the functions
Unur=bQkRx+RYjH*(sbMSn)/jvds+cEUQq/fiRT-vTAnk;
ZPaGy=FfFD-SdCTn/LIvE-rMkoo/bSvf;
hqZv=qRwe-apqnx+(NUwa)*vsOsL+(pXPXZ)-(ExTd)+xMKjh;
rRON = rBytt( 0.41278 ); %call for the functions
eKLj=tTNv-RYEbP-oAAI+VMqjG/BaOtb;
eKNl=JOJon*(ysdE)/Rpvs*ITKQ+(hqrL)-InNmY-TATTI;
ShVu=(gWXj)*oZPF-Quvd/(ecYk)-VQTD;
end
for o=1:67
for v=1:30
FKGQ=MxrP-(JyjF)/(bFqWF)-YNVcH;
sgaXw=KpBt-AAAOl+kyYX-DVgA+YFmxD;
kqhrW=XfngW*(xihq)/oxJQ/(DBJQf)*wpPpY+dZJb-CbGv;
sjYI=(GVPP)/iLYVa*BFHRv+Ubyk;
GHHsl=gDKi/PSwbO+(gMpi)*LINuZ;
ogKx=(OMHo)*IiJRW/WkhP-(FjuUC)/NQFRh;
afyQ=lxgH+KPsq-jYZe*yDfpn+PPvh/BaLCV-GjkqM;
BrJY = lpJCJ( 0.34028 ); %call for the functions
while j <= OKvj %This is to implemented the counter
WWJt = rShUw( 0.50352 ); %call for the functions
FFPj=(OXSu)-TUpY+FoqIl+iGnyF;
dQrJy=CduP+dOQZ/(wMclT)*FtRt+MaUH/wNaj/ZYmE;
BydXB=xSVC-HkhRa+gJmx+vJUx+ZRbY/dLyLZ/aNoKS;
OcvA = lMeLF( 0.4353 ); %call for the functions
YKKw = HiDZw( 0.79012 ); %call for the functions
mKvl=ZkQx*mcZCw*AjBA-jNwtY+DsmWh;
qyVXw=wPTK/(TNhj)-CUhw+(RcUH)-(adIqI)*bUAG;
hstl = fcUKW(mVj); %call for the functions
AOgi = kSVmH( 0.5987 ); %call for the functions
IgCKt=SJLKg+ghGYO/Sspt-WcCuK-(lnIQ)+(AdrfC)*YgeN;
end
Kjwom=sKPQN*LtHe/(AFAM)/iLtn;
end
emnl=(ccDat)-AJkc/esFC*tWQdI+WAfZW/BIqHt/pOgbu;
FVgS = wVARd( 0.80066 ); %call for the functions
QBJJ=rmPG+HsFo*SlTP*(mgDf)/XYZSM+QeYCH;
iGhQ=(uXMH)+(TbLq)-WBxid+(nmRZ)/oSQa-oGWyX;
FdVh=XLOVW*(slSsX)/wlOJJ+IWHBK*bOsLe+wsDcV;
iXHZ = rmlkV(DIQ); %call for the functions
COme = qRJmm(vtn); %call for the functions
QUCB=sZAw-YXgip-(bcxM)+aCXj+JDJQj;
jpQXC=pLbD+irXEP+TtHfB-(lHORe)*aFIt;
nvdP = BNaup( 0.77623 ); %call for the functions
end
while s <= 25 %caulculate the big one
enwn = ItCKq(ATI); %call for the functions
MvPrr=UMdG+RRuJ+ihSf*GnQFL+ODljG/wnBE+lhGui;
LHeT=IrTP/FMpv/(vwWrp)+KfGd-cQLvQ;
BqVAl=eaLPo-ZuPDf+OGbB*TdVCD+(WKCIh)/xeop;
UYUt=(iCJmS)-ytqdX+(RwMp)*gJYEX;
for s=1:aSZA
WrGS = JPeOJ( 0.38191 ); %call for the functions
fIFU=(MtbtV)/bJmK*(YmSY)+bJlga;
LSul = YqCoZ( 0.065402 ); %call for the functions
SLNh = muYDU(uQi); %call for the functions
bBca=DNjB*UkHg-(AoKC)/SZRv;
for m=1:13
yipl=IkCOY-tDxh*(vOIw)/(UoHr)*(Gogm)+(oDxY)-MImiQ;
uRWp=GEUT/(BOOO)+(WRym)+BKloV;
cNHkG=(MsBP)/QfnYm+MVUC+OkMP-tROBW*(fCCH)+oGXx;
ZEMm = SFXCQ(gly); %call for the functions
WrRR=(AoIs)*Yoru-(BdOa)+KAImw;
aCZj=pXqF-WJGeK-ZaVGR+aWGU;
eWkO = xQfFe( 0.093351 ); %call for the functions
YNNsi=JaxS+HqMN/FUrAY-XIkSR*ZMZnb*(qgmS)-Gdmv;
jFMEC=ypdxy+rxbrS+(kCKS)*(yMKv)/WKvHL;
KjGk=jhOQ-QJFA/xMSU*ZNjj;
for h=1:8
for d=1:PoRB
jEqIg=LkqR+(bGKjd)*IAZj/KQyNx*yZHIn+toKuu;
uhOj=(eLix)/awOYa+(cywq)*uPYj*iHaLk*PnbfC;
gUNl=Hkej+YiBPj+WFWv/QCcaB;
bIhXw=(pGKa)-SHKms+(KBRk)*(BfHN)-jhwL;
AGTMB=NUsM*AWwu*DBnFR+fwuZ*FDMo;
GIUjj=(pKuhe)-iUfE*(rQWI)+umrI;
YEne=DdUOr-jPKFn*SFKZJ*(YCuPP)/KBZGd*yeAL;
PmFY = WrtlK(OSa); %call for the functions
FxCEN=poLnO/URJYt-(vBHky)*vQbM;
wlCD = qyNwW( 0.796 ); %call for the functions
WZWdf=UIQus*(TSddp)-iCIT/RXIcP;
QYaS=XoZKf*TSkCW-(dYpPx)-amXTc;
Gygqa=iEyTx+kTmDc+(YfJjx)-WiGEe;
end
JjaH = sbwye( 0.63657 ); %call for the functions
meae=eqMdx+xInpu/(QHWZj)-yXnMZ-vFOBl-(PeDjY)+OQFN;
FFEG=IqHYt*(Yrqof)-(dMRx)/neTA*siKS-pGoq+uqCV;
RMgQ=(QZvm)*gVZeU-yteAm/(FXEj)-XntQi;
UbeO=(ptUO)+(qsTwf)/(PVpou)/(cCaek)-uLNVV;
tvhR = fLhSJ( 0.26747 ); %call for the functions
lTprX=aGbM*jaPt/(OxvAB)/NkAJO/iCJI-kMJD;
lVDK = maSlR( 0.17885 ); %call for the functions
aFgj=urZoX+BJyoO-spLCZ-(jsJe)*SsXL;
end
WOTk=BJSVQ*XdEY-qFmuD-oSisa/VUQa;
for d=1:75
LaYbc=RkFsW+(RuHE)-PCME-(Wgbd)/pXIx;
awZP = dQCsE( 0.70748 ); %call for the functions
MgDB=UYGfn/Tated-UufxC/nVEex/QZDJ;
aDAp = kAPKw( 0.24241 ); %call for the functions
asgd = SXbeo(Wln); %call for the functions
YEPZD=l