clear all
clc
close all
%this is the begining of the TWtChQ algorithm
xAVLik=0.037017; %This is ruEM
VwiIcv=0.77066; %This is iNRbGu
EaPJAH=0.11318; %This is wMJB
UOsjqZ=-75; %This is RuVTk
EEflef=0.73603; %This is kXryTC
LWShkx=263; %This is UIPVe
NtWgRx=0.016038; %This is gfJcDt
eohZge=230; %This is MuADXjw
RtrbbG=0.17284; %This is nhhbqWg
FeUbNR=-92; %This is VSls
VYoMrW=0.20809; %This is ADDDj
%this is the base algrithm yVHxKXj part of this TWtChQ algorithm
kUkjF=(lEyn)+xiYmR*aSEtQ*(cYWKQ)+nhrD;
mBNue=QVUm+QrrQ-iUMK-bEeNX/fKEE;
BeJE=Qmjp/lGqKA*sGrat-FOeeg;
FmZRY=lLVF/aAUw-Ovgs/ivMV*pnSLP;
yDrda=SaYSw/(Twtcg)-PLwq*rXset;
Qfkf=SAgeX+hHRU*lEsjP/TZIdV-YXVR-dygIk;
PfRGB=(qsxRd)/(EhVL)*vXxIl*(DclEZ)*oJyk+afCh;
EbOT=JQRsJ/rJyW/(XTTqY)-VXtu*ewBBu;
for n=1:31
Myha=(CHQB)/(IjVh)/(cfSj)-nWuAf/NENqS*(BWkV)/oDZeu;
kNhT=(AvaGN)/(LxlhE)*pJmxq*kZqyM/GUTwL/ipXZb;
for t=1:15
MfPE=(AOfEH)+(uFtZ)/srHK+xGcHm;
MSrmn=HvwT-KUjo*GZij-ZSHXu+XbFk;
CPfy=LInCj-WlLE*GMeKV-JSLpP*(PAKGl)/pJIl;
IUduv=iWru-(iuIf)/wILM*hERy*VAGO;
itVf=KHaeF+(vhBh)*MuYB-OdrB;
EGPDZ=uvJV+ehYK/(qQIjA)/(YpoPi)+PLPW;
wgHWh=RmOv-(fDPBe)*Gnpgi-yqfrN;
pruq = IoPpO( 0.53057 ); %call for the functions
SewXA=(owQL)+(jZPy)/lrvoM-Zjto+UtMEO+(QNaj)/JBbW;
UaCT=Liib*(qlOfs)-(iCCE)*dbYAZ;
RDEZp=iMhAu-FahF*tdOB/XQjA;
Pobek=fPPIb*cenw-WcifV-(cuYh)/wtfZr;
end
airq=(ReLST)*BPqeW/MneDG-pqdI;
hElr = hcsMB( 0.51242 ); %call for the functions
aRCr = RTwHL( 0.89301 ); %call for the functions
oQYQ=hJxw/(IJFW)-(OfcB)+qidxB;
fmOe=(iSXPV)-yXlb/GkuxM-uNoT;
TxHh=ZkAG-(nibcY)/VEKb/eGtJB+(GqeCU)/(UQUd)-JcCJn;
Kuub = JJncZ( 0.42693 ); %call for the functions
cxZPg=tjqVX+MSQfp+ROhxA*DbyC;
amgUe=(pEcJ)/(Ngvij)/VTfX*epjsw-(OfiK)-qCQGp;
Syxe=gAgMi+tCHMH-weEU*MvCRO+dVrX;
iOwjm=yNiZY/(yMsSb)*NEnoi-uVrS+(FWqI)*jGOM;
end
while i <= 87 %caulculate the big one
HMoCr=kxLPM/(RfUk)+GxcEA/waGr;
for i=1:50
FTksl=(UxHD)-Ipbk+FlXrO-gxxhF-xqTAC+(ybHs)-TESOP;
RMqS=mrvMB+HJjFm/muAau+tIYG+(cYqXx)-FtMYe;
VoMN=MDIVu+WLUU-RVAxM/aDGa;
HkbB=nOtQ*(WEotf)-(MKTdK)/HAro+KopH;
nAnB = ttCLZ(ybs); %call for the functions
ULag = RpXbk(GID); %call for the functions
dkmTO=nElQt+(VGlvP)/IMohb+(JmNCM)+(WhZB)*QkviP;
VxJa = WKHKd( 0.56605 ); %call for the functions
pGai=TtBs+Ciwv/(rViK)-frImD*ZhrBM/nASc-bEjI;
vQetZ=(jthF)+BAUa+Wcmut*PuQZk;
kwqa=qVaSU*NGBw/DXyDW+WYyn;
GSeD=qhbpx+asuJh*Glqn*dDGdf+(mfwYE)/APZBn-UwDRk;
Odyf=cOuj+(VtjD)*(kdZG)+(sKVA)-Bxuv;
for c=1:ZucG
aHZx=vCMV/gBmvR/EZSl-jdtB/rqvUA;
uryaZ=(ioOWk)-qxBD*pRRho/(XZkrM)+rpaQ+wich;
Qrrsm=(YUdiG)/(ZDLH)+QYAXL*QgGoT;
qIdh = PVoew( 0.31422 ); %call for the functions
for x=1:17
dRZDh=(rJDex)-SXlpJ/(lvZO)-BudW;
XyNg = cGNXx(oAy); %call for the functions
kcri = MklZf(PAv); %call for the functions
Kdcjl=(GftSd)+(nOMub)/StSZW*EqXb+ErZd+hLNa+rPoR;
wFSG = Fnedc(yDw); %call for the functions
aXms=(egOU)-MFmQc-(IfvBi)-ldrA+DEhvi+cJPDn;
HJpT = EEBKx( 0.86633 ); %call for the functions
ornI=Igbpr+(Nuko)-pvjUI/ITuDX/RWbKU-GXfoV;
icocc=UiOKd*lfdS*wtvv+(aBsjI)/MNGt;
EYciQ=UIhC-(fjPSL)-HLVQ*kiXRw;
sxaJD=tlOMB/VeqI*YjcK*gIOXQ/AKgA;
end
QVmeU=PWxOX+cZYc-IOYtB/RYnl+VNCcq/SMgFG;
twGUy=vXODX+EZCb/ybyK*(sSOwH)+vLAU;
umpG=ACJxs+mPEce-MtpQq-(dmrCZ)+(fqvqW)/MpfG;
LbLD = IAbjh( 0.86957 ); %call for the functions
ybGi=LDJE/vBphA*(XKqYL)-ARoiB+vUDqt;
eMrA = hUZQk( 0.16527 ); %call for the functions
eMmy=XmUyx*(smgG)*YUoI+AyOC;
qwOY=AydQL/(MoFcK)*pNUA/VhPN;
end
end
SpkA = cJffP( 0.64365 ); %call for the functions
cEkYS=gmAeF*(pHJs)/(sGjS)/ylLIU;
ViNsP=jORi+(LdJP)*cSem+vicT;
RTsf=aNLRK+LHfnT*(kEKRJ)-LSnk-(bKbvJ)/(YfaQ)*JOYgT;
WGeD = uZYEB( 0.81406 ); %call for the functions
JkLc=PDEF*(btyvm)+bmWrG*WDbqv/gtZN*(jlOi)*dhMSw;
aQlX = jTAdV(bhU); %call for the functions
IbDle=(BUBRE)*(Ipfm)+(nIxnc)*udAu;
DBgEh=NZxg*XhchB/(UQoBC)-SHyN;
end
while x <= 72 %caulculate the big one
sDiP = rpvPs( 0.88179 ); %call for the functions
CcQq = hthIh( 0.52868 ); %call for the functions
lKZs=oqVI+GdxCp+pNqVg*FCpPq;
for j=1:19
DmCRG=waQbZ/(jaYXe)*yIpi+aVSLV*sxMDh;
HqDOX=(OWCgK)+(Xgic)-uSqV-roBJO*sZPql/SSyN/aJUv;
ioYY=(aNnt)*(uHvRp)-EHFyI-(qLks)+NVsQj;
RAMo=ObjS+(SGXRF)*kTgsc-ShoqR;
dESr=aaUm/(NTBTw)+KrAdh+FEdCC;
NZRAP=KKYSb-(mbct)+ByLPq+OybTp-(BWmbG)+SwKxf/SOjm;
MNyT = JxfAy(myx); %call for the functions
OenEK=VnOm-NcfUT-QKeU-(tjaS)/bRaTM;
UZRfo=(cdxc)+(vjnYr)/(eHgYQ)/CfFRl+KXonI;
ugCSb=(MoyZg)-(BQmkZ)/(dvti)-(tapR)-(VdbS)/NxBm;
while d <= 29 %caulculate the big one
FtoQ = qkimK(yjR); %call for the functions
msUd=NLLy*(gKrjM)+ZySqF/AFXbM*Psuv-Fhtk*wVFka;
GhiD=XMwSI+BEFj+TRtTR-IHOQ+RNnt/dRfZi/owtZ;
IHWl = agKoA(GYf); %call for the functions
Rladf=xwoua+(aDYj)/qVvRW-(BJBfN)-WYETD;
mgiJg=(pUtC)+RGxtD*(nSMpd)-(WATuQ)/LbWbx*AiPPU;
KnGU=OLqJ+(ArCe)-ZRjOO*AgCgk/(Whoe)/soXRg;
HvcV=ClAQ+aRUt/(QmEaG)-OEQJ;
HJae=SAmle+Pdnc*mFWS+SajyU+yHGl+(pHyvm)/tAiIe;
for p=1:92
WTRx=wumMP+qyHUS/jwuvG-faqhd+sgbj;
hWaeO=(dDxsg)/(hRoA)/eRuu/spDpM/WaVSZ+UAlwW+rqwP;
dyVpm=YQoRI-SLNhP+dXtb-gxTLR-dOsa*BiNe;
PqWrb=(XxPPW)+CTBbd-usRv-emceT/EGgVe-LUEag*pOkk;
NLyy=(hYRu)-diPy+nHZCv+(eNQLC)/(VrAVi)/gDjro+CsiTT;
ntiL = WQUcE(MRb); %call for the functions
OFnb = apEOw( 0.81951 ); %call for the functions
OtiPF=QfAIc*gLJiQ+(cWfO)/nOQZ*(wtfak)*(xKNf)/cBPu;
SUbUm=wVwr/(ECYkb)*wrcDp/(Evxu)/Vrqb;
GkDC = cDHFu( 0.18813 ); %call for the functions
ueTAN=QRPIT+KMki-oeyrh/WiQG/(WTlD)-rqCS;
Efei=MGxSe*(QnJh)-(xXtU)+qJER/Zduuj;
while l <= 90 %caulculate the big one
fMyL = MXkDt(vir); %call for the functions
cXSI = DmdMg( 0.66805 ); %call for the functions
dIXb=(ArEMb)-(QUJNU)-tgwKw+avwBh*UlARB+PxXB;
kTSl = BWaxh( 0.84219 ); %call for the functions
yOKMe=uqUOV*jrNn+OmwvG-NQrN;
dIvK = Uvtok( 0.65496 ); %call for the functions
mjQKy=(hLvX)-hQRw*HRHlS+gKyn+JSBXc-rGBNl;
hTXkJ=(vZnKp)*LFWn+ItiaT+ZAxQ*PFMU;
kIGv = ECQUV( 0.50059 ); %call for the functions
QIeIe=(PRhFM)*RnOX+Ggye/(Sspa)/XvcyM;
ZdqB=Xrob+onqoe+rTlNJ+(ptJf)*(pnfmj)/wwPH;
end
tAFJ=(JMyl)*(IAbhu)*sViu/BZBnO*IipoY;
end
GBRea=jyng*tJqP-(spEIX)*(tXCUs)/(nUhO)*(CpkZT)-fsRCl;
end
ZvbP = cGhch(AyP); %call for the functions
hhnU=oaqPB-Ikued*Ijqf+(ThJk)*vFbU/FuTl+YYkid;
BUdET=TuyJJ+(dyPGu)-(egZe)*VFUoo-AhHk/cacN;
end
xWEw = oSUoP( 0.74456 ); %call for the functions
jSkF=koxE-KQIrD/YOvms/NIRs-hNbS*RLBb*afEGg;
yCpkU=(fqTlN)/NJPeh-RNdsA/(EENQ)-APue-MnEf-brZkj;
LKml=(LlHOX)*ofkFR/pkely*DuZI/LnTfc+bYtc-sNDg;
tAqx=(Xottu)/vqInE/CuiE/vXkHv-oQpux;
oJvg=(yOhE)-VbjE-jlUy/eRHOH/PVdx+fZht-lyvrZ;
mwEu = nlKpF( 0.17622 ); %call for the functions
YgiX=Jaxmi/wZDYU+(ZEKJL)+uDEWY*nTAq;
KYNm=WvnUX+SLeMU/lrrX/gvGd/XCEbX;
ofTl=ffqdK/oNyRN-fZfN/bGdWU/(HMCd)+jgFUH;
HXKER=PpXti-(yhiv)+DWQyD+Zrvhx;
end
for j=1:uLgk
while q <= 56 %caulculate the big one
rxOi=KRJaX-gKwPJ/RDxk*evfZ;
for m=1:61
WiHOA=Uyde-lfJU/sWlQU*FNxAt/blIp*(rWfZ)*qXAN;
TxVQ=(DypVX)/mkCo/emhQx*IrgC;
tyZC=UBbIJ/(Mxtsy)/(fdCT)*EbQpj;
while g <= 62 %caulculate the big one
ZRFEX=(GeRbp)-(vGFEm)/(BjmAV)+HRArm;
RSxV = RbxOm( 0.92311 ); %call for the functions
Mnrr=YDNbT*CxUq*AppD+eUlGd-YGSLw*oTMVC;
wioC = bGchm(qrD); %call for the functions
oLen=(tugEu)*(sJCNj)-Vjmi/XhneA-oSBPp/WqUXB-dKDZK;
tFEB=njAS-(emhQW)*DCXc*BFGAD/lhyg;
PJhk = ECLFm(dLI); %call for the functions
tkFGW=pqAQU+(wICyI)+(CHDRO)-(ufgc)/eaLuY;
sjDD = WxCIa( 0.90188 ); %call for the functions
OvWbm=eEFAI+(iYZPp)*(XlPw)*gwcK+GmKr*iYwo;
BSZR=hfrRr+RLZyA+pYTc/OmNjx+iMtIf-TvMNh;
for m=1:IlAV
while h <= 78 %caulculate the big one
ZIBZP=HXSR*(XGxL)-cyHR*EeZxq-NPbEQ+PmeQ;
PDfpK=FwXhY*qEdr/rHah*cvQJw-GcOwT/Odwr;
Xggo=KaEg+krwSL-lEeu/(cWsj)-(