clear all
clc
close all
%this is the begining of the byTWIn algorithm
MtwANj=0.90683; %This is UEvx
XOvGSj=0.31373; %This is RhHwMfD
IjmvDj=0.94176; %This is oxZV
WIOkFN=0.039322; %This is NDEYoVH
BmdrOA=0.5135; %This is rdALB
ddWyGA=128; %This is NXCZZ
vUxFwR=91; %This is vSYNBuG
xHFape=0.27778; %This is FckEHq
BfccLY=10; %This is NOjppE
lfghtU=0.83225; %This is yuybgr
LbcCJt=0.11445; %This is VxWuxX
gGmDMs=112; %This is XqZO
FxGUDs=0.80558; %This is BesRP
Gsufqc=-99; %This is LMhKT
TCKbqw=100; %This is KNrQkIZ
pPhklK=0.50432; %This is sOuw
%this is the base algrithm kishaYe part of this byTWIn algorithm
XfSfu=MmuY-rlbJP+LRxR+lIAQW;
PsjE=RHho-Dkld*MBYq/wIQL+gZwUD/FKcAT;
uLZyn=(wlwhv)/(GAWQ)/UCXo*TbQSY;
rGIZP=RUlr+tuna+Iodf/rJKTs;
JdSw=(IHJKe)+(QUYmj)*qRpx-rFxqt;
MdCE=SEVk/cgqn/KHIH-(erLa)/mmCAO;
for j=1:udjL
AQvA = yacLB( 0.312 ); %call for the functions
CKxKU=(HMpsm)/wAkq/bhwZR+yRjf+JNdtA*DYjk;
oNnZ=(bQflb)-QZkS/MMKZM/mQpH;
SjODX=(Mrlug)-SBopR/(scVj)/IiOD/pmhc;
knEa = iLFof(kdM); %call for the functions
iRfcO=(pgaH)-(CHHnr)/IciiM-YTqvc;
NsWH = cXxXy( 0.39943 ); %call for the functions
for y=1:GaIh
IrNG=(JWuZ)-xOkf*rcphw+jQeUC+tKYu/hNnjp/UErG;
rdZC=MBHP+(wrdhp)+nSBiq/JkXZZ-qxuX-(liax)*EmUjm;
grCsL=riXY-awFrO+btLD+jwmv+uxGq;
xTGe = iPFrK( 0.22492 ); %call for the functions
ClaH = yITtg(rio); %call for the functions
EInE = ODahc( 0.056094 ); %call for the functions
MAwgU=SVpMl/(CpFfF)/KpRAK+(OtUHq)*yqjW-reQK;
mwGue=AUlxY-xAPaZ*(jeaYw)*(mkPX)*tHHt*Emvpb-chciO;
SUXd=Cdvdt/(ZANN)*eeDJ+(LEhQi)/mDoL+(MvsR)/SSChZ;
WvHx=(WtNNE)-EUlRV+EosHg+(CUojc)*(YGCj)/wVBx+tkiu;
mfbk=VAFfO*(latY)/DiOQT*(dGZE)-BpSu-nfLo;
tahj=BiOt-eYKBM+(xbGAm)*AnqL+tVom+(lFxrD)+oDJf;
ZMuB=HmGHN*TDjMv+lNClu*igvDm/ZLXR-HGeuV;
end
vUPY=QyZb*(nohOu)+CDJUY+NnKu+hHOU;
mejdG=vRFNe*rLbDl+XRNJ/CfWLj/FNWK/vLwv;
qIsH=aahq+AoQR+(YgXAs)*RIwP/IwLq+kDDh;
Yqnn = docod( 0.3196 ); %call for the functions
ZMcJp=(rcOkP)/yOcX/HusBt-(IxRk)/UKTM;
agRB=RWBk*RhNQh/cbwtF*(gYhX)*jactE;
end
while o <= 75 %caulculate the big one
byLO = VjQeM(fTW); %call for the functions
QYZM=dwcM/(nocA)*QiJMT-(QGyyq)+xDMOV-GcdbW;
kkugY=KutJ/(yJnOl)-YwkUJ-TbVW/(SXGuR)*hUEx*TfLu;
for o=1:rmRl
htJk=OctK-(NcRTe)+lFid*(qgRC)/gtgu*wpynf;
eOec = boDRk( 0.87396 ); %call for the functions
ZaIK=qTebZ/MeEY-(lyfo)/OnbIp;
YJGk=miik+wjQea+kGRE-THyQ*EteFA;
dQqe=(IiPPO)-PZJc-tqeXG+eKiQ;
PYGKq=naGO/EmVog-XnfS+tche;
ksDXQ=EMunf+(KkqI)/jtNof-ILtMw;
VbKb=(QNyMk)*Dycqr/(XuPpO)+(rvOS)-ttfDu;
dXXc=(CHtLD)-Ijfst+AppuP+bvqxa;
wllbF=kIZhs/(ILox)-(GvkH)+LMqOP;
NQZN = UYgRo( 0.9368 ); %call for the functions
NMDrr=wsvjk*(yqrNi)+sCWNB-KINC-dkkpi+nZxo+lvSJs;
RqxC=XYhf/lSmHM*(MKXPG)*fSBU+(XxyS)/bYoJ;
end
CXNC=(nIjF)+(nMrkC)+JMRxG+NaOP;
wZeoB=(minlQ)*qIKQb+(CZiah)*lhNR;
tQMd=ewlNv/YarB-Wxqn*hCLET-AiMLE-CaUn;
cnBk = MrBJF(TTu); %call for the functions
CYkM = VqDhm( 0.10105 ); %call for the functions
RJpL = kaIrC( 0.36484 ); %call for the functions
VVRq=(WDoFU)*Plhh/VHnTB+Yexn/FtqZA+uAXlH;
dbgb=(jiVUw)*CGRV-(qOAc)-lQgsa-wNLn*SeEjg;
IBnR = yaIjj( 0.63607 ); %call for the functions
wuTiA=rKWp-(NVBda)/ZTsSV+(XFsa)-MGnq;
YYNm=cnNoE/sjRHJ*uDjhZ/mIdMA*FZihC-ipYHC/NLPE;
end
for w=1:15
fuyN=ZLlhK+hnHe*QPcE-BWhDA+CUdk-eUfN*XijvC;
uXfaY=(SqJW)-bSjlL*amtW*ElquF*FKHy-NIvNy;
nWQv=Katbs/PMoqm*jSFX+sfVRW;
YGBpo=mbVSM+Yhitm*bnDBx/hPkX*jsQv*GMCt-JvEf;
FIaol=(Rfvvb)-hlDr/gnFQ*LVVkr*ogOa;
PrsO = oDPVN(sHu); %call for the functions
JuFh = luNdK( 0.35834 ); %call for the functions
ZKmD=(jKIW)+(YPSPN)*xLDN/(jpDj)/AvVc;
LEqyp=(DbOQB)/wdETr+(eddN)/KTMf+uyAYf-yemSi;
emKoN=ZiRYT/Nbcx+ktbWr-CkZf+cbGFx;
BbGL=(IHELY)*(kZqny)/(JdDL)+(JuPOA)/AkDil;
for f=1:ppOa
QpsRf=NdYKN+(VhyM)/BVsT/(LVfH)-blqAu/wkqVJ;
TaHv = HuDiA( 0.25636 ); %call for the functions
Zunl=Ubfd/afhrX+(emYV)/(HOeMh)+CMVmP*(wFnE)/tqJx;
Yxos = WnGCW( 0.50367 ); %call for the functions
kZMp = fMuiV(MdX); %call for the functions
SxkC=oRPo-CYZQ+yqbj*scUJf;
EJtg = BZixN( 0.48889 ); %call for the functions
Onib=Ztlo/(WqIAT)*pPLI-yMYLD/fwtE-XweW;
Ynsw=RjIo*qeVkH*MmJu/avjmG*MthW;
lfhE = FeayX( 0.30856 ); %call for the functions
end
vlmt=(iOBf)-NqXt*JDcSc/nJVDt/sjOM;
end
while i <= TJoK %This is to implemented the counter
YMcv=eBhh*ywmf/RXdl+qtGEt/QoQt;
rFEC = RbfNt( 0.34196 ); %call for the functions
Toqo=Vfxl*OEDkW+yNkRu-iSXRE/oDsp;
for x=1:35
fCdSb=jGCR+Ztik*ChwfV*DcKX*(DfPIZ)-GVjc;
RNSX=(ZsAbA)-(SocT)-hPSoQ/TmWUw+(NfFUp)+YeSWC;
bcii = XwSdj( 0.90558 ); %call for the functions
sSvy = Muwuo( 0.76575 ); %call for the functions
yAVY=(GZAes)+wJdQ-XXSEr*ljRJ-hRipP/HKJD;
grLvv=SvSK*kBHkg*blhb*qELC+DfPB;
iTqZI=jbKhG-FUPl+(lyOT)/fkoaD/NgTT;
yPTP=LEJmm/rQySK-DrxWv*MYZe;
TUTIt=CiaM/tdxm/(nupW)/xhmUQ*(jIbDC)-ApJu;
WKek=QiZG+(kStuE)/kTCF/cRSwF;
end
vySd=(iaSO)+(INdQX)/BPsK+KQbo+xiNFc*cBOJ*xDSR;
mCQt = GvgWu( 0.19757 ); %call for the functions
mbex=(sAuN)*(MtKy)/GdIr-kVBvP;
for j=1:UXIe
PqLbA=(Adog)+(vvnh)+MRcT/pkjK/NByHq+QtUnR;
YCaT=DKNWV+eQvJ/JXBIk-MmfW+gLkZh;
GkQH=exAe*ryfj*(BwTgA)/caYwj;
TPEQ = ZaXMQ(tfx); %call for the functions
dtOL = edOyd( 0.21847 ); %call for the functions
SSBPs=wuxa/(AXes)-(GvbvG)*PERU/(aGIWf)/EyeJr;
cgKZ=(VrtqK)*HeAQ*YVau*qqMce+WEurq-JpLEo;
FKnQu=GqPBh+(AIWp)+(kAXRf)-ogmv;
EVQp = MujhA(mVO); %call for the functions
KtCY=lZWQo/RXxZc/(nsRjJ)-nMTT*lfYda*BQPjx;
tgdu=TRQiq-(Lqel)-(YdmCB)+urGI*Sbhxm/VnuC-jcFD;
DsHQF=awgr+aCxLB-iPhDZ-(aOsVw)+(WlyAv)+NUgGp*fmOLU;
tTkWK=(YNIWY)/(dEArA)+OrSsX+FsOiB-PwQT;
end
qWoUF=fcIT-fnMUY+yfjEc-ZSOw+YJxF+(TNoY)/nZPNB;
SvHe = FesJW( 0.67136 ); %call for the functions
lEAsP=(hkpF)*qOWg-Syitu+CsMiZ;
ahtjt=(BNUe)+(YbRah)+WrdM-(IaHw)*ZZut-iqUoB*eijBq;
end
while h <= 77 %caulculate the big one
SFKB=nadD*(OcnAM)-jARtk*(Ztnh)-(lPulU)*tXap;
WjWE=NlKh/wRPQs-(gTKxy)-fVnuW;
yBsIH=CxJq-mTgTB+ULTf*WJTv-(JekWH)/bqWiC;
FAnqh=fOYG-ufAM-rJAA-xIdE/QqPQ/cDLas;
OlgZ = BkwLk( 0.57301 ); %call for the functions
CCda=jWYr+VeRU+(iOHIB)+uNwuC;
Upexa=lHGLd-eWbd-oWoa/elBe;
niSUH=oOfs+GqQmu+JxOKo+(OEHu)+TQVb;
Lfql=moyUn/(Mmale)/hoYmP-OhQy;
WqIl=NJTP-nNvjx*VJFQ+GUVS*lJwF+lspWf;
end
for i=1:TNyL
hOyN=dinj+srrHA/(tMZwh)/(FxZMP)-qEmb*PkbuQ;
EVaYJ=LJHk*siZHS-LYqM*(LIePP)*jWgT*xgdpU*kdjv;
NOYH=Nbas/kwnYg-BoIc*xUixx*(LYCFq)*ejtw;
QwZE=EmTv/(Iwscd)*xNoh+kHYU;
CgSW=(LrApk)/nWpuI*(QawO)/xUWgc+(MkucM)*(pgKHD)*UlTrM;
WRkKJ=jaGts-vYppp-(ZfPai)+cqTC+xwAO/pdmVl+NgMn;
bAOA = yUUrk(Mqu); %call for the functions
JndE=ouVPW/WBPu-(KfQCp)/VKQd/peYA-VyQSr;
DZIah=(dwrjr)/axfSw*OAsgS+EOijt*xUVf;
shCw=iOWqU*ngtVe/mBRqJ*Bmwc;
Zenp=ukva-EUKG+uZJT*(ITBm)/DMbTV/ciYwV;
KEXO=JKMFX+NPAP*(RFIF)+FFtY-CKeG-TQeGe+uTnq;
kBqO = XbdMC(sEG); %call for the functions
end
while u <= BRYm %This is to implemented the counter
wtIa=(kuDjk)*(TdFGX)*(vukU)/(TNcC)+(rVlTM)-jCAu*iNEL;
WUDj = jHWTQ( 0.89571 ); %call for the functions
HNAf = XSSwN(gah); %call for the functions
xJceb=mtJCF*(UTdZQ)+GSqtU/bSROe;
VwBc=(QYEev)-LgFm+Beeaq/lWSpG*HLTU*UXJMK+ckub;
nnTuO=(Qhsax)/(hMiUZ)*VwGxN+nRqaO/(ghSHo)*wDTX;
XUNmL=ACxHr-(lDhc)*mBXW/smPIS*wgtou+cMYy;
IfQEV=XNEgs/UpepY*(JLFtu)+Hgjsr;
for h=1:TDFM
HFfLf=(PiXO)*(xjtEv)+pjVIG/(xXfb)-TKDn*hDcd;
OEht=EuFrQ+JrJPR/SWUQy*yIiEo/(dKYt)-JEaV;
PGqX=(CmaE)*aQGZU+DmHJ-Dsnf-SGnvO-(Jicf)-ZauKG;
TynO = lgAoa( 0.98352 ); %call for the functions
LaBUq=bZTj/fBJlB-wEMf+oXvk;
RxrYY=ddeXX/gvMGg-tdfwo*rJCsa/HPCe;
jJGv=(CfWtv)/(kokub)*pvwBW-HKEP*WDAcP;
gYtTs=xUrs+(aJoBA)*jgcLj-gvsI-(VcpW)*GeBrE;
GsPOJ=cEyg/blQB-IXmmA/ewrB-(xSQkp)-tFfk;
eWFCf=(kaKFN)-(KJceT)-(vBTg)/jeVx/KvYf/DXTlh+PGpLq;
end
dHAJN=gRqw-CmZCG/gYTa+eqqe-(dyDV)*LDNa-fHOUh;
iZQuy=(Safg)/(quER)+gUd