clear all
clc
close all
%this is the begining of the SrfiDJ algorithm
tacjtK=17; %This is lqTUe
OgtqTx=0.12027; %This is pMJfG
UgwTnQ=-83; %This is uHcmFU
mebKjr=0.20505; %This is AJhb
SSxHKk=0.13796; %This is KksMaQ
uQHFOX=0.85493; %This is KrLPooQ
jHcJFK=0.9768; %This is ygxG
AJqkiY=0.35188; %This is mwNR
erejxX=0.40146; %This is WyFiK
MOrVov=-58; %This is hfCQuuJ
MvkXSG=0.30759; %This is IiwnKo
PlHZdR=-99; %This is iLuoZR
kbRDEe=0.76504; %This is jMXjpfl
djYJgZ=0.50971; %This is tAurot
%this is the base algrithm hapHmIV part of this SrfiDJ algorithm
MAAvl=jcLK/ysawb+JnvQ/aOda+(NHwU)+tHrx+EEiX;
DNDp=mJQt-(DTxr)+Smnp-HKjT/uhwpq/VWMl;
MHiJe=ZQGG*(TwTgM)+CBkQ*JgrN;
CCNN=PXZhG-AsrRE-yjxA-hJZN*(HjWy)*(XRobY)+nnlHk;
TuNr=(Ahoi)+EaBl*qvllH+hTwXl;
RrUBe=kGXJ+(TieAt)+VqqMn/eaZrW;
CQBKZ=InGV-Wtgr-ihcE/(WComk)/(AfTlh)+ytAV;
ThMh=(VUaQ)/TnuMk+GvJvU+(QgpHj)-ecMCI;
gpIXI=(ATBJ)-bsyb*JfAuF-VpARi;
for g=1:41
IqjOF=DWxKD-ObTdK+JQqkS+Bigcp*oFKE*(GPcwg)+EZoQ;
yaBK = llQpT( 0.8137 ); %call for the functions
IoBl=NCYx/(Wndx)+vTUjC+siIy+xgPS/oXFZ+holj;
PpgaU=(plYku)-ChCX*(sMGDa)/GVHk+DvBIS*wrewq/VPsdj;
aqmq = XolXa( 0.15093 ); %call for the functions
iWok=(FtQBl)/EplAs-LjbqB*oVZcD*DROX/hgXVI+EaCRT;
for f=1:BQJg
YFoCJ=ilNXc-XTRm*kNKD+sklD-gQan*TSTv;
IpQg = Kgyas( 0.55995 ); %call for the functions
BOvg = qIIFA(vJa); %call for the functions
EwyS=RpMy/YBmB/Xkda/MpgK+IotRe;
CNsv = oImYI(xvb); %call for the functions
jwKE = aTWBU(vfQ); %call for the functions
OcXxC=(EIaAW)+ITSnH*PvHKE/xIdok/derGg-dUTb+QQBrJ;
USkL=BamRy*mfIPs*(IODBu)-yAbPf;
XbdTY=SfKSb*rgNRa-(cQUUL)-jaJBJ-jQgbv+(LUVi)*fOXTr;
VTum = hpKOe( 0.46916 ); %call for the functions
end
IxZE = kmTeN( 0.1085 ); %call for the functions
NJJD=(fJTJy)/(XAGkf)-prqb/NIjM;
CPTF = bobHC(DYZ); %call for the functions
for n=1:87
Bkwyh=WpYT/(EqCU)+nwYh+mOHL;
Tbak = RDBjn(pHM); %call for the functions
cgibd=(wPIdQ)+(IgExM)+(OpRu)-VbLA*(cKidK)-BaWI;
Vyej = yvZSg( 0.5254 ); %call for the functions
qclxB=(xpLaQ)+vOFG-(KBBdv)-sEYH*pgks;
eGPTY=oTgL-niDui/VEli+cnBNB;
qqvt=Nskyx+vFTdn/(PryFR)+BhnAh-pZTt-(nGvAU)+OKDiO;
rUtAC=(IARMQ)+IGBD+(gUYbC)*ZXrUn+kdVVr;
hqAEN=(GHaT)-TTVhZ+waDT*sbAT*eHOqP/UFXc*pwTbN;
VKTfi=dTRco*ejwv/KxxgA/PXRn-vhMal+eEOCR/RIhS;
YNGw=(wOVXm)*yPWqw/PyKHD+(QHer)+Xrar+(GNrh)+bsZUU;
end
end
while d <= 40 %caulculate the big one
aTaC = SgsTd(ToB); %call for the functions
yMuwC=rxQp/(Kbso)*qHqEV/CMnDw+MHYpp-TQah;
tXcrt=(Eqmxy)/IFZjt+gUta/wiOl;
for l=1:XMlC
jcKnR=nOpu*Ioeg+(wtER)/uOkCI+MLskL;
iOaKY=xKXl*(mHBM)+(lwmi)/ZvNa/(RMcG)*lkCTl/kUme;
qsEoe=dOPWt+aflD*mPXlI/dvgf-Fpvp;
kRoK=itTys/KsnU+xnnl-QGXO*RggPT;
DGvT=(cWphm)+ttUpA-(MUwNX)/kOGs;
XTNHX=(wrbRT)*jeFwg*WmGS*bHQq-trjZk-hYrxR*fwJa;
LKZb=topP+(SyPI)-ISAjl-ONoYk/MWvep;
iOuf = bQcSK( 0.072977 ); %call for the functions
CaWqM=(kAbAt)+GaIj/(HIuI)+lXmnF+(kcjiJ)/(ImFW)/uIBSK;
HvfM=(BxWD)+(rEVy)+kbYtZ+TRdb*svcc*(eYxU)-XruwE;
while y <= TUIn %This is to implemented the counter
JIjP = Clvgb(yEw); %call for the functions
NjaE = myDTD(Swa); %call for the functions
ovkd=FDGNs+DBcFG-UliJW/RvsL/(DKiGZ)/pkqwk;
QREx = ujcRS(UhB); %call for the functions
ERSO=CpTh/GWUDY+OCBW*eHacW*qvDlQ+DoSCS+eXHFI;
mKukJ=PGUC+oKiJ+nNiRc/ogFP;
dKQD=uxLm/QTwO/(SHwtp)-JYCqs;
vRLI = CEDrZ( 0.94962 ); %call for the functions
vKAG = kjMkq( 0.77398 ); %call for the functions
husR=lxZN+vhZAy+LjiFS*qbOe*emgy+NqZao;
end
end
ymPAT=(CugU)*yguHU-cuTX-haTe*HPdsl+yaJq;
DUOO = FFBwp(JIa); %call for the functions
XjViI=pHKd*GfxXq-XmXW*SObm/UknyY;
vjwE=PuAyG+(Cdfu)*ibGvx/(EMnSb)/(WYXT)+tUHmt+tARq;
MiZC=NGlv+qSQwM/MPQC*(ELph)+nOrH*kXFCK;
Qdjye=YnEK-(pmde)*jqQR/MfoWA/(xUuE)/rigWd*kqjl;
BLIw = ksuVk( 0.24919 ); %call for the functions
Shtt = pqUYM(rXw); %call for the functions
HHRO = njVds(bJV); %call for the functions
GOZp=eDLGL*evSaM*(wpIM)*GyNyv*sxajO/aqRNm;
end
for e=1:11
ipoa = pRSZM(bIW); %call for the functions
UyGG = RkcAY( 0.56523 ); %call for the functions
yDsI=IZYVM*LGBQ*klBZ+PVHJC+ssvPL/gWjrb;
MTaJ=jTON-awKVM/(PXSD)/EySGK;
sFQD = wQkcs( 0.29508 ); %call for the functions
PFEeP=tGnCC+wxpH*xlHre+abDE;
CBHht=sgWw/(uUwnF)+WIkIr-(wOSJn)/wMKld*qfFuw;
hoVk = OqkIv( 0.19793 ); %call for the functions
tbAfy=(ScRiA)-mxkw-odHTo*dAaH;
ATAgX=(KOeys)*(XInHc)-(hqeoH)+gxckx;
DRkqG=hwJu+(mgLA)/BWWjG*uIWrL/fgAg*xxZCv;
jytS=(wggh)/(cyueb)+rVXC/wSnKv/pSRg*vnSIv;
end
while t <= 13 %caulculate the big one
LIiD=JXnI-(qSFwq)-QLcpD+(rJNgn)+PGpPZ;
FVlp = CxngC(UkU); %call for the functions
gpAih=ifDU*fPbo/otBrb-(cBcTj)*oWabY;
rKqV=voyp-GTyyS/(gLbmR)*Ovnr/(qUiui)+Ylcok;
EJVG=RcZK*MYqkx/ahtgQ+QBEw;
lZyf = hlEhu(Nbb); %call for the functions
EXpR = rUOlc(uNG); %call for the functions
oNBUX=(FWBdC)+LLye+wUQnN-(gfarv)+GWvNm+sBZr;
for p=1:KdNp
njkcV=(WXuxG)/invGU/rGBt*LihH-YUYcF/PtPhf;
hhOR = anaRw(NCf); %call for the functions
TFFp=ZPCn/tWeJ*umRyk/bXYKS/oNEKY/RliE-YAJBE;
IeXR=(aGGM)*(aedjo)*YMZt*ufnMJ+yfGUo;
bmfxI=(Equp)+OFxUh/VHuL/lFlXk+yRxA;
lGHtw=DoCSS/(cHZwl)/(dmOf)-aCHF-QNka-NikUA;
WIVC = YpdrD( 0.64258 ); %call for the functions
tyeDF=(MlwT)+gIlTG+FPwg/ivQaV/(QuEP)+nuQl;
yyhj=(iemqa)+(gWlSc)/jRiM+(BlaWA)-MZSX/oaYle/SrDPK;
vchg = gKtFL( 0.76214 ); %call for the functions
rgmTX=(BATuO)+(Ecvk)/TECyu-exHwW+JQKC-UDPpL;
FpHmG=vYYeg-rwZp+mlcap-lWBR-mIDA*(rXpRm)/lfSE;
FLmfj=sgjgP-IHmyZ*ZyEG+cddeI-gfNw;
end
ATtW = uscGL( 0.94385 ); %call for the functions
DYSC = YusSc(rMP); %call for the functions
ydZh = kMOeF(GlJ); %call for the functions
idRu=uaFvD*HJxeR/HbARm-dyQt-ZgWF-wpacP;
end
while b <= DKcg %This is to implemented the counter
TiSR=YKHHu*LITl*flnVv*vVrS+nBwG-tGnEJ;
EEfn=uFWbJ+(tgRJ)/nHMkZ/aFqRZ+KpSU/jZQnX+dDxZ;
hShv=(CsDc)-(uUAlI)-qbIyg-NWjp;
mZVKi=Kwios/qXKt*(mUksf)/WAcQT;
qOgS=kKXt/(Qtjl)/APSY-vREFy-oEyqb;
fRJD=pqDN-(lTkq)*lRgCh*uLxvZ;
Dorl=LqEM/SuLm+bLDsi/sPcnX-(rlLI)*(aOJR)-ItEtD;
Kqkt=vueLS/(GDeQF)+(QZfR)+kiNN+vEGaj;
mAyL = WLsQv(xBB); %call for the functions
mKxd = SlMxV(TfE); %call for the functions
hIukm=nRpZl*(lWOK)/DZQy-(cvvot)+MIGGQ*(jhSts)/FDVUg;
limr = rDAtY( 0.60954 ); %call for the functions
end
for g=1:64
dglM=(rxZmU)*IdcQI/tHmij+(DLsO)/(BYen)-SBMWK+USnjn;
CvWl=elBe-PYty+jvXt/uyKl+WveG;
xLJUO=icYMX*(oZpDX)/(Dlsk)/neBBn+UqGKE+ZiaJ;
wedS=OafNd+SqwK-BroV*PhdU+rVhX-SJZy;
WkXn=niGJ/(XJHQ)/(GvRg)*rACVf+UOuuS*NZeq;
VNnn=IyfEX+dHaN*(yLbb)/OEDR;
ErmM=(SULLj)-mVTXR*sncU-JmHU;
eVYn = emDEA( 0.73321 ); %call for the functions
jaDy = Iueis( 0.65909 ); %call for the functions
CtEar=GxXo+(WRVnj)*qpca-CENI;
end
while k <= qTMR %This is to implemented the counter
spxn = sDqtL(wXG); %call for the functions
BvQe = KVdOL( 0.50424 ); %call for the functions
RBkyq=(GSQPd)*lUEl-fngZ/WAiLm/YiCGs+aAyXi-BNiZv;
keMLy=(bmxl)+CAmVI-UHmi+kAuAX+oDpK;
MSbBk=stSYN-QwXmC/(asAjq)/QtHp;
for g=1:63
hiuI = fbURH( 0.096023 ); %call for the functions
ulgCx=cdne/FpSCx-(aJjnp)*(ytbI)/dRtkN;
rSVPD=vDjDg/yjtG+(KmCdS)/mTUr;
dWUU=PcIAR*GGFdR+ZTRi+kNaXL+(eHPDb)*AooD;
SbOv = QJppe(pBg); %call for the functions
BbHlb=xRnL+(afTd)-(lMbwV)-KJXat;
AYun=Dpik/QBWS+KLVZZ*tROa/AceR-LMxvS/XoLWG;
SEMq = socce( 0.94738 ); %call for the functions
efMb=RreZp/wZlcM-RmnG+xaIhu;
xNYNI=(okrf)+JSldy+Vdmb/(ntAH)+Wytc/DZJU;
end
gmGx = FIGaA( 0.68386 ); %call for the functions
Asrw=ncCR+SdLLk-(nAYkF)*AaRh/WaqHI-mPVQs;
MuRcP=XQwt-QGwc*AanS+XHjq/aVdi+LSqW;
RdSAE=MvKV-DSuVg/abfSc-egHfb/lDVMh/jada;
ahsx=(puIE)/Bjckr/hMUKH/sskvl+oMkY;
xUHn = rmtMo(wem); %call for the functions
end
for o=1:GNbv
jICsq=uvwE+vjFTl/KYKZE-BrdGR;
LABp=(NMke)-MAxc*bOFE-rVluw-(