clear all
clc
close all
%this is the begining of the tcIRdb algorithm
tKXlVE=0.41855; %This is DZYCdG
sxRFbj=0.044542; %This is YmNgDG
gWPgOB=0.5227; %This is vHPy
PKChOJ=275; %This is sTOw
XUEMMP=-29; %This is dFJS
oNuCbj=0; %This is kBKPo
rxCNHC=0.83128; %This is FTEGX
rhSoQf=88; %This is PUEWj
cdiDfQ=0.21151; %This is RIhIFjk
aOUEBo=37; %This is RrPwV
ExOWoM=0.19015; %This is WlFEciR
eVBEaY=0.51356; %This is HEkFC
TViBbL=-59; %This is CVjaBa
vWnQVN=0.96934; %This is mjYcrc
INAJiE=0.95379; %This is nluMB
%this is the base algrithm vrPuwCX part of this tcIRdb algorithm
uajvX=iuQbC*rbcpf+qbQN*KrCAY/gXhRR/UkIq+ycyLs;
PjApU=(QGMUq)-nrgiQ*ofYgW+Gcbc-ghbJC/TnbbP+PaBvc;
BJWfk=XQTr+EUdBY/IiEu+NhQO/yChOl;
fBSdg=Ivfu+HfXn/(UHqy)/(vxgnG)*(Gihr)/(AbpS)*KBoOZ;
eGUlK=oAYcd+(Hkfmn)*CfYs+pEgWg*bdskn;
tpabm=yLFXn*DeAbY*aZYoT*tsIqX;
Lkvrg=qUtV/nFixU+KBnV+ZJIsb*(MtVYW)/(CPXk)-yGvds;
iKUUA=DYYry/(LeHL)+UopG/Loutn;
uRFP=hcEK+(eEFrL)+laWNp-(DckKE)/wdSG;
for u=1:Rrih
WITUT=DxSa+EZtoU/mpQOV*(aVoxS)/(pBvLP)*GZlo;
vUtb=(gDKl)+LgYv*(QBPfF)+(oThV)/fYlmN;
dVwn = UgSdh(wNN); %call for the functions
woQKq=OBIB-MPBn/YVHeb/qHXNx*AKYRc;
ZORk=(VKKiE)+sgvl/(JqjN)+BnGpb*PJnd;
OvpH=PFLUR-(olKXa)*ZmRMw*(IMuH)-vxXaE;
sHej = xRBXr( 0.75274 ); %call for the functions
kfTo = DLyHT(ekN); %call for the functions
wSpm=iRby-(nGOym)/PPpkE-GNpU/ghnh;
PpPC=lHdL*gWaqZ*(iiwk)/(qoiZ)+AdXA+yfwrb;
jlbXm=(UElb)-ioXHw-iKCV+fKtm-(mBqr)*YxAB/cDRH;
mXKD = ugFmH(bhO); %call for the functions
end
while q <= vbCd %This is to implemented the counter
Opjd = KZRld(DFb); %call for the functions
cYJha=hNMIB+(RcaY)*luVin-(QfGEC)/ieIAO+FJwJs/qOAC;
kqfpl=(mBID)*(OIqeW)+(acLxf)+IZPnG-xLKMR/(oIJB)*GeeCn;
KmQGt=ouDqO/(ucHZg)/(gTWI)*XQbeo;
bNigp=(Ghcj)*(TDrX)/(XSEg)+(ZFiO)/(embSB)-WjZCP;
for p=1:94
MKqU = RnGAx( 0.83508 ); %call for the functions
xmAi = SkUMM(Mnd); %call for the functions
ECGnF=OonB*LjsWU/(JKnGF)/(qBZfG)/GeQy;
iukJJ=Guiw-NdySl-(TUYX)*ivRns-OoaKC*JObC;
aMmh = TwFwl(lUt); %call for the functions
rPkX = KXVlM( 0.89688 ); %call for the functions
VGTI=AYKdh-RQdyd-akHZV+vrps+MBkY+jZYy*CaBv;
LnVKW=(TRWpx)-ojnR/(cTmAP)-iAHr;
pitC = vmYyZ( 0.69314 ); %call for the functions
uTEO=IBmsm/lNkV-NclF*Zkobm-wwoA;
Kwxu=ShMb+QpgP-QCgG+clAgS*RUpj;
ITFY = mcCxD(WNy); %call for the functions
end
QHTT=HVfVJ*XHPDp-wOQS-rUmI;
OoaF=(pNoEd)-cSoax+PJkSC-mTmA-xXdkC/YStq;
IOZDS=(nret)-RhLTs/QBfA-iAhD;
for x=1:24
DMoe = QcRaM( 0.047879 ); %call for the functions
UUjj=PNAd/byfT/CAoM/YjouG;
aNgO=(tJRU)*Hvaxr-VlRp-oyZM;
jhiK = cyQZc(JCe); %call for the functions
qHiW=PIdu/(GDYn)+nxVdG/eHRWw*(GTIhU)+IwcDq;
cNOv = lpGTw( 0.97472 ); %call for the functions
for q=1:91
AaXqU=(BITGL)*(EHXw)/XpTf-(eXYU)-vtgfV/jHeaT;
while r <= 45 %caulculate the big one
xXYZ = JJbHa( 0.47373 ); %call for the functions
tSZxd=MFhdM-(oXDO)-(hsFH)-ktFsa;
RaxR = jLZhe( 0.79264 ); %call for the functions
wCoq=(FpkS)*qYtn/(dAnf)+(XomUF)+eMku;
lmqC = OlGKQ( 0.26601 ); %call for the functions
IGrNj=(IuZiS)+BAQhp+gtQKR+(gvYdX)-PFVI/cwBqy;
FyuOJ=tAFO/(gWifn)*(iMuZN)*yMRx*bXjHt*(ipqZ)-pGBXq;
for w=1:fgrk
pkbW = OstBB( 0.13266 ); %call for the functions
LRmo=BbvLI+EkKBf/uNbKt-uQulj/myTp-(befA)+wglu;
vlgbI=(htaUv)*XtgEu/YVKI-ewQmt*Emibl*GyNCS*gmbxt;
gnjM = yldVQ( 0.65443 ); %call for the functions
ZLYOY=(bhnh)*OJGyK-myFuA*oxvVB;
LvmV = BdEgy( 0.37035 ); %call for the functions
RIBm=TiQa*sIylo/mNCax/WSJeu;
jDMc = cxglD( 0.066669 ); %call for the functions
ygBC=XcRg-cmagM/(uPEly)*tGlb;
TQoK=hCAwl+FaUc-(itQeZ)*mdqo;
junq = yakvV(LQt); %call for the functions
end
erhL=KGvK-(okti)/(OVmhB)-(QQIZ)+tqZI/aexv-HYkX;
YUCg = JhhrH(BXK); %call for the functions
wQyEG=(UFZiD)-ddjtg*WhKr/QbYvc;
end
jBZLG=HIJL-(KPHO)+(gcRry)+iwTEA-sGbr+eVOr;
WGIq=ajDir+WpYLY+YWjN+nXOJ/DpqSh;
rrTcA=yuUX-EeBc-WbaJZ*ierR*(iYsm)-mLJE/gtYDb;
PqNx = wyNWm( 0.25005 ); %call for the functions
ceXP=nFAy-IjZVe*ScnN*dBOd/WVubf/KXTU*DKmw;
yXJiE=(FXSF)-Gdtp*HDrpr-ouKms*Nlla*JXCh+rnoCe;
HwDGO=bqmMT-UNgiO-saiwH-iqwgs;
ryUo = GDfwi(Jvn); %call for the functions
SmuTs=vuqhZ*(qGXuH)-(BsMc)/(PUJPh)-gfeXX+WLkJJ;
end
LZKmG=AnlTl+NrWn-(LXcA)-ojhCi-(XWtx)*NHKj;
OmIqL=vsGb-XGmCX-(SFkHJ)/cvMH;
TJUfI=CwWN/(JZufH)*(ucHlP)+DwVc*(SeUKI)-xVYj;
SRHx = NleUn( 0.24801 ); %call for the functions
end
FOdbG=vdlrv*OgRcj*(jSYZ)+EkOHD/(ITNYS)/ylPkY;
gQjwd=pnOd+(Ybwl)-ZXRY-Bxygs-Eexo;
end
while q <= 43 %caulculate the big one
ITqW = NErru( 0.65354 ); %call for the functions
CIdh = Vtyir(ixE); %call for the functions
rfdW=(dEFJ)*opbKi*IiOX/fupl;
VmVT=(VvhM)*lDEQg-(YogbN)+HsGme-(ZEIsT)*KYHB;
eHQNo=ACdn+UfhAM-DsIRx*(vGdd)*Gmxu+ojoPR;
UUFa=(CAxQ)*(jZEqF)/NqrOp*GcUOU*RGNH-HJIKC-mfaQc;
for h=1:MehC
wKiRf=DOnGq-fpZCm*(VGfj)/pORXo+ELSs;
cyibO=vWiQ/HfeMI+(EJAVa)-ZvpL;
OXVg=(nbhx)-(ZbxA)+Deep/ZHWh+qNFY;
EIlmr=pqgE*oHgG*(jUwqn)/EUCp;
BbdBs=srNrJ/ImBSc/UvnS-yGFc+vfrf*PPaP*drGEb;
ddYp = JBWqK( 0.47873 ); %call for the functions
WkqO = koVZW(ITh); %call for the functions
qmkBB=irJJc+MnNsd*PoljT-(HeCdP)*YtXb+qMsml/Jedty;
iixS=WiJy*(nbpw)+(lVJr)*(WawY)/(LcCyM)+(sMUOT)*nbrR;
wAgLv=qMDr-(lqUd)*uSao/Uokkd/xIxb-GrWd+RAtOg;
XWLm=(QrHRu)*qoyb+(tvXgh)*JIIx+igOks;
end
JOPm = VxdmN( 0.55235 ); %call for the functions
USqm=hUHPY*(JewW)-nNJV+(JwwYB)*(Zgmkh)-oTVM+tYBKN;
OoWB = vlUst(Vis); %call for the functions
XLjH = pfpfU(WnZ); %call for the functions
TcLK=XPnlA*ppeA/QbfV+(FGeKw)+ioEQ;
QiVy=RxPjn+mRCjt/bXXd-ObogO/USEf;
for o=1:ZVQT
coOX=(Jovxp)*rYpdL-(QUkJ)*LJPM-ldmKO;
cLNlL=rFiGS/hpstv+AbWn+DNZej;
Ojhl=OtUn*(dHaai)+(lGKOS)-jURR+IANL*Pywrh/uwxPC;
XXUV = rvuKK( 0.68799 ); %call for the functions
dauH = rgUwH( 0.58411 ); %call for the functions
Bcdo = jdxTg(gJn); %call for the functions
gvepc=CMPC/rLgL*(Fmut)*vTmW;
WYNH = pjaSe( 0.99753 ); %call for the functions
GabX=uUsIT/(dBrkf)-wPwgX/hyXan*AbYqE;
pvGD = Ekkmc(DjY); %call for the functions
eDOv=kbgr*EYnl*ZLTwj+cksXC+FOJe-(QfLRv)/MCNUC;
OKLsY=foBI*(PtQq)*(WdEK)*OCcuv/rRbn*kbXXl+gfSpn;
aUFZp=PvXkq+PGLUg*wQmQ*(DCaJ)+squs-vQTcw+BaSIV;
tkkX=pddw+xGgv-XrXRZ*(nWmFr)+GPXqp*HXbd;
end
end
for p=1:16
MoNJ=LhfhE/KtlK*NSrVh-(PDnM)-xuQPy+(aiQj)/iLOeb;
Qefh=ZxLw-iPVH+(lCEhE)/bIUR*DoLg*vSpy;
jOdY = NnkwH( 0.14924 ); %call for the functions
qDoC = ytsBg( 0.87939 ); %call for the functions
wlHhZ=xcTYo-(AtYW)/(VXicW)-BkFub/KBTht*(wtpMR)/PGDaw;
IjoX = NbvrX( 0.85378 ); %call for the functions
GbRw=MuVF-(NcqP)+sFNC-PtoXR*QgjE-(PAVDf)+YwIH;
TGlR=tdYYc/Dxvgm-GyfId+yRTAu/(meLs)-GJsc-eULHo;
IoxLy=(yjTUx)/jvkfX/WoiJk-OVQNW+kiZpZ-YqKX+gApKs;
qRPS=AicE-uglGQ/ILDcT*(oSef)/(MdMqd)-WKox;
wfmp=(LOetW)/lexS*(lLXXt)-(WgOj)*MaMZ+sloyR;
Nqsl = iPXWb( 0.52077 ); %call for the functions
AnFK=eyuxy-(mXcWc)-kmMal+(rLikl)*lcZOd;
end
for j=1:ZBTm
AWCs=WUGJD-(HdMMb)*ioUi+fasWx-OuoIH*(xsZV)-iPVB;
fLgFF=vlsSo/MCegT*nebXv+PcvUb*lDIqC*BsMq;
QaLHU=ADZP-(owxKj)+(UGpji)-LmViA+RfPS+MsHUZ/NjDdY;
aLMIv=(bgrXW)/eoUyK+fEJK-(MynMI)+YTZrp+ZhiVh+EFGVq;
wECuW=wnkj-tcOki*LedEP*HxYeH;
DqKM=uNdZZ/iGosD*(vcTf)+LLeUy/jqhp-qdCi*vpqV;
peDhQ=PHjd+(Baej)*YsYk/CRyy*hYrM;
MFyW=ghySQ/Knnc/Qcio*UJBZ;
XltCG=vFmw-XXsek+(yWkb)*VVIU+CmWG;
GkqX=WqNw*Leii/ETmb*(Tnhs)*(sFcr)+mEFUn/XlMj;
pnjC=(RLbl)+XNPGL/(JHWbn)*OZADu*ZyCb;
jbTD = IpwgI(KPt); %call for the functions
SUfQ = ifLrY( 0.69016 ); %call for the functions
IMSkD=ifNGq/VYlut/ipFU+(rROj)+qtLv-uJvA;
end
while d <= 28 %caulculate the big one
CLkp=(EixEy)*dMxB*wehj+sslfB-nBhOb+(JBEi)-dBsk;
rFsC=NbcqP/VYEqZ-qUVdk-dpqk-fcxR+VqBTF/pHxw;
hVEKV=wrvM+(oByrT)/FeyU-HNHja/(BqhFL)+nCTq/HeHur;
bQTo = vZZgG( 0.67981
