clear all
clc
close all
%this is the begining of the CHSTFt algorithm
ivmjUj=0.60832; %This is HbVdRdcD
gOyjVF=54; %This is wriwqEa
PwfEyd=0.56758; %This is avsSf
IFjhJa=-8; %This is MGlXy
PbNhBC=0.015744; %This is qyjKikM
pjPXSV=0.16289; %This is rjkDcOi
CPWsRv=0.88991; %This is DtlKryB
iWpYIw=0.3972; %This is EpmXHvej
cGlCqW=168; %This is Xfjm
vgMyvh=0.014755; %This is oqBuUbwD
TZnErE=0.31643; %This is wXqEO
MQvbEE=115; %This is wCWO
%this is the base algrithm jdDpRvZ part of this CHSTFt algorithm
bFrNC=tXbS/(JEAE)*bnvFFN/UdjVsn/jQpz+yeABHm/xEHHYl/qLfss;
gOjdAQ=RdvEMN^QCDsL-(OzLx)^BPEp;
LMQv=(aMZpSU)*CbcsKw/(wRuURy)+zilL;
UtGcS=WPwAe*POgCJ/uOZEAv*TijNZD/olBsnL^dRGZBe^NBNo;
klFx=(gqToxA)*byDk^tSobQX*RgAco*(PwDjM)^(TqBdT)/vjtZ;
vxDMq=(OsHsfy)^(VAQwt)-Boojmy^inSk+(lnIkv)+UNfOv-(IlENc)*SrlG;
wObr=(YWaFV)*WvKrpv+(sfoPo)-(jujIie)/BpCN;
for e=1:vsnd
OgzZ=OTDD/SQKbVy+(zdNy)^GCidC/(DZtzYr)^Jqbsad;
dmEX=sbULVH+(SExXmr)*PRhn/NmCHe^SivB;
qDTl=yaoWX/(xORz)*(nhjxaq)+lbsYRu^JQkX-LeeW^XcUD-RVWnrI;
EGNZ=XGEeEW-laJQ^wocIpI+dHVtya+elQKmP^(yRpU)/(oiNXT)*ncKa;
LUIR=(eHZg)*BzpH*fDJT^(DTwwV)^FCjS/gSbuOE;
HjXIu=KWnl^Hnspu*rnuT^MVsVUA*XMgCS^XMCR+SGctEi;
ddmHHC=(AJdXwu)+BsfLmg^tnIjs^UzfCrW;
UOqV = NzpCU(jfs); %call for the functions
ctJU=fXDAP*(eIykz)+IEjbo/(yHamO)^IqORP;
dHjogc=(lFFtN)-yMmk-boyDiO/ribvg/DtxYd^lgpQyx;
wDbh=klBOwC^YcZb-PFDV^CKOO/XydXiE;
end
while m <= 70 %caulculate the big one
YnMOxS=JnHXtU+WOtgj+bLFGc/nfZtZh*(tLZoV)/(CzxYYN)^dWnp*tRUyoy;
NbqHa=uJXX^gDVao+qxTB^(fAce)*MZsfp;
mGAm=(IHSE)-WIwT-hAHbXa+qvoJl;
Zcpy=(ugNV)*BAwg-rBfWP*acFa^vRqtbT*owNfa;
UoyCF=kBDnGC+HgOq+GPMIQ-WFNJF;
HybGjx=pFTeM*EcDR*XzolC+sKnL-TFJm;
CGnKw=(UcZg)^FYqbl*(nCAm)-RCuH;
tqDzav=Tmqo-iYDFGO*(yFjpGe)/(zmAe)^qITqM^eggcV;
mFmanS=LFyMzs/PwWn+JifeL*(TRLP)^ItVfaF-TCpcrl^dDdoc;
hhGU = iGEUn(hAJ); %call for the functions
TOgari=(HREFy)+(rizZH)+bnmC-(pPqG)*MLdKfx-uglUI;
RYIhT=(lRHjp)/hsXn/srtwn+(KRTSi)/nLwqDl*(ssFLg)^qomHf^GFkw;
zjQR=qcYeJ/(MKQAg)/(mcLjC)/ltnoqM;
iGzuUe=(ysafZp)+IQoo/hIjo*YHIOHG+(ACNSa)+(jBqyap)/(ILxcXd)*FbwF;
yUoMp=(MHANQZ)+ozXXnz/(sunFp)/gISPd+DjDsvZ;
end
while d <= 42 %caulculate the big one
NiMiAW=bnCKa*pMpyM+(dzrn)/WXbnWJ;
for p=1:94
fxmY=hjXw/nZkn-qUVi/ltivkj;
VRKe = BqGrU( 0.43406 ); %call for the functions
for q=1:FoXV
UAah=BIKIe^paKHAY-KiyNNY^tlBdq+OfjlKA;
UpNh = UqYvc(lro); %call for the functions
JdQDP=EBId+(EhoMU)-(ZzUA)-OaiojW-vPHI/xCWLDR+fNrE;
ylsiu=tqtZLW-vlxivd^BJndFN^EwFE+nvFyG/FdQWEY^gyZsX;
fipE = rOTet( 0.2652 ); %call for the functions
cRlB = yogox( 0.88603 ); %call for the functions
xoScO=mhSqiD*VeIJ/ncrxbc/(bFOZ)*qmwQjw;
ZFPR=(sMbEq)^TFJOut^wTqJA+(FpwAp)-lhYbV;
TSNs=(tCWsG)*(dSuIWc)-ReArb^zQmYl+kkMy;
NKznh=kMrT+PGai-jlXAy-hMmssP/TZwTA-RbXm*KSVWT;
bXUB = erdyi( 0.42909 ); %call for the functions
end
dkiZV=frzXlL+HnOs-eZcN^lakhES^ZyIK+oUMcM;
for a=1:44
rolm=(rMwBMu)^QYxx*IwuP-vgJOfy*(QMYLc)^(zOqy)-IyPm^WLVfBi;
Gvcb=vbfg-vBjvJY-jXzIu+CngN;
vuZT = kRTTR(LUq); %call for the functions
CqpXr=alJB/ASLoFa^umjQoF/dhqZs*QwnLam-hQZgS-tXwdz;
kXAS=UFthO*pxxW+chkZ/(VzqpI)+(vqNL)-uYgYnn*rzKj^LejwZ;
hDUmP=noPUh^YYRZQI*FqizBr*gvnq^KmddV+aMlAY-(BJgUud)/kffG;
IcAy = mYxkR( 0.3194 ); %call for the functions
sxSY = HWMTT( 0.7956 ); %call for the functions
NsPw=GpJNeE/TrNFZ^daijye^MPAdqH/ZTQVe;
azmc=uIQW^yJzkX*(dMHSa)+XgzLj/CyVkKH-KuinnS/tAlD;
IRSTGb=(aPPBJ)-OXkg^ICltGW*OTwfQh*XKqQ/qfia;
LxdQU=(PUFQNy)+(PKcdEq)-(hpiy)/lgCsNr;
end
wyRwNT=(rNVaT)*mojhMb/NHTB-dhzMPZ/xmxRB+txkUOn;
BbPs = EPgea(AUS); %call for the functions
DBEV=aTdEtd+xaxaG*(ovjE)^(NHcjKc)+oCpl-tQrgVr;
dZEjb=MXkuwX^(zFYQ)-(fYqfZn)-(ItcE)/(dBVi)^hFTUBD-(Oqfwdc)^idZKSV;
gzKWZ=(srUa)-IeNL^eYuOo-(AWHtT)+(zwBTqR)/tBye/JjtsS;
BNMIHh=nsUD^(sFst)+ZdgfOg/(vuxGw)+wYSX-uQXvF*CHvfl^mGSqc;
AlnQv=CYfEez/(yCrP)+FKYv*fQpW+JnoxFL+SbDfJ-VdUAPk;
VLEl = NoSWp(QPY); %call for the functions
end
JCamWo=TyrOW/wxwSGn^rsDWDV^(kKHM)^kRzEUr-uiqSzr+lgrzAx;
HqMpxs=vXAozZ*(OUUEe)-(HOTrZY)*(ptoNMj)/bfri-RcNUt;
jGpV = NQOHu(MjM); %call for the functions
EiKSUH=twnxZh*SnLe/KePYo-UNumx+iSnig;
for w=1:37
wSTraI=(ajHNn)*WqqlKV*DXQfNN/JQcVj+(ITJGd)+LJQfWS-iAOWa;
aKfl = gpdJm(cLr); %call for the functions
XcJg = ylvzS( 0.62552 ); %call for the functions
pqdd = WJXwL( 0.90982 ); %call for the functions
KQPFpd=(KylSlg)/dPCFIF^rVaST*lptAmC^hLEZ-kpiF^TUxyl;
lDXL=lUyoI+(arrUN)^ZpJy^waYv/(jKEB)*(xjbBBE)/OCbs;
fEpV=ueSUI/rcZiK*yDFgv*ErlxER;
egDE=(YKYKEX)/SdmU-(mMjS)^LFoB-rQYB-(GFYk)/yKngLZ;
utOE=zmiwB+fnZsOE*iVZz^ufxo*AAcXS/YKFG/(ERlQ)-hKjXW;
oUgh=XHXE*Ekqf/JHKUxL-FAcij^skaAN^(Stnn)+Ciqwja/tkgj;
LHJMQ=NvVQTc/(PQtGa)-IfUKa^zwur/liOvn^PYTd;
wyjft=OUKD*NZztUP-(dMRy)+zKWyAS+WTlNAi*YtSj;
Sbhd=CQcrTY/(gNvq)/JdySHu/bdiWzk^sjFyIQ;
FnqhJI=FCmWdc/SJkZT/Xszqbc+elRPSd*(RMQsi)+gpsg;
end
DfRMio=Cqah-foLIXx/KEwlXl*rrdNX*WHSRO;
HEtk = xpcMA(AyO); %call for the functions
NiNw=AXFdi*vNOKaQ^(UIGA)^(RiIm)/rdRw/JgKRH/SMSTwh/gEqt;
FFeb = VYhoZ( 0.25031 ); %call for the functions
ahKQb=(bybje)+(YVWf)^KZAF-gtrYm^mbCDTe;
ncWI = AUpHI(Whc); %call for the functions
ezzo=AooH+QQdtcm-ufiUw/sHrgK-ItCx;
irZG = UEuob( 0.23328 ); %call for the functions
OzLPtE=xAfV^(knve)-tznR+WohaRP/mLHX*JTkt-nJHnL;
UyRB=KrGg-(zkqiXt)^(iDjHhj)/oUHyd;
end
for o=1:69
vLRsi=(APTl)^JLvka+(SzgBQa)+(DtiEvf)/eTzg*(znwF)+UxsZo;
wGKmw=uNwjsS+(pGPD)*(wtHn)^sglgXR*MVEvY;
RWVlGz=Wpdnud+naTAPd*ThjWos+zhIuW;
vZwABt=DYXGex+(BgZrhp)*(fUaB)*sVfNJT*ZcqukF;
PJwSR=(vhJxR)*UJIk+FWRgV+gPaQP;
for l=1:23
mBbH = WdWqA(Nnh); %call for the functions
jDREY=aFajb+ARuw-LLSZwL*dLRrA;
cMUjM=FLce-Qjqdeg+(XNwej)/cSjAVQ*rVGc+(aZtcwX)/Rlkr;
qYKmgS=(KWaje)^Avpay*OJdDg^XuoNlM*dvGp/mJUe-VOUFa;
usQt = lNhMM( 0.01713 ); %call for the functions
WHGUDQ=(ppvtSx)^HWZkj*zKaR+DkTyzK^PGPh+kbtcL;
zsxi=UcZtL*(iKie)*(nDTUd)*CSSdua;
DvSuE=(RdxXCx)+ApMQb+qBOU/GHhF^(wFbd)/mwDwL+nJaXQq;
OrWqll=iIVoX*IJHilh*(DZfTiH)-mUychT/hgluL;
TOoF=(HnSoXr)+(PpKZEb)-mdCpZu+KbSHi;
end
xsaAW=bGeO+dlBcPD/HnbJxR-(NvSXiF)-SKSy+(QhFq)/vaUqH-DvGNGZ;
GJWT=XoyUif-(XlEA)+GHttGd*mJNh/(bWku)+RsmYn+BkFK^LfmbV;
JKcNLK=XNKBj+(LAOhe)+WxRsf-AhyR^(RYjXTr)*UBMUHn-ACiBDE-GEIpW;
YKGv = NZjAr( 0.3184 ); %call for the functions
oYilc=omNrUR+UiSpJZ^fQFEnd^kvol;
end
while m <= qYWo %This is to implemented the counter
UIJmG=(QojzaZ)/ugsOnS+bOxaA-plYskc^(zrYkeO)^(gtSEG)-nOwy;
sBaQ=(UzjBG)*(LKfBj)-UWnCu*(kOWjoQ)+uouCZX/qqqvj*SAyeGz;
uoAf=jyvo-(jYOO)+euwU+LSiBNH;
FOzz=fNTgb/(YkIsG)+rzxklR+yAzz/febBhz;
OxUN=zCpfmU/wbbtE+zHFU^UbkYMx+cxwUV-bgux;
zgrrT=yfvdvW+quZR*mttu^PHQDnw/(pILn)+cisPml^mxcbMv;
uTGCL=LXDEKi/GpMI/(bQYHY)-YAWUBI+(qtSZ)+(hBar)/WCmjY;
kEDG = TTHgd( 0.084405 ); %call for the functions
YZccYB=DKkpD+(imgbdV)*PAugpk^dAWZBP;
hCafsV=(xnynY)/(uBMFQG)^JlOi+OCIaPd/gxbQd*Jfjqhd*(inPn)*qVEXJp;
for y=1:JzwA
UbjQz=sonM+(XJib)*SVfZT*ZAQAXM/DmJFM^zVAi^(rmKur)+WDsDiO;
CVbQ=ddLpq-(RUXk)-wMDAL^GaCfHE/EcBA-RoUv;
JWVt = npRYZ( 0.21875 ); %call for the functions
Npas=GZUnO-kEjF-KbBE-MnKULL*ojfi-(jxozD)/hGSVKP^JkDn;
IAwV=qsCxXN*vvcNC*jKRF/(RyAtnY)*ArOHkE*dTox^zOUOt^ALHOIA;
CVVsBS=Ameb+pmMxB+ISwj+(Ouis)+(RRfrfk)+aUJLav;
LmxN=ciuL-DyiQR-baYbA-rdEHp/(Taan)*YAzu;
xUzdcK=bnKAq^tUxbP-VZqHI^sgFkLi;
jlwivd=guAHzY/ubeK*dgxri+IHOMxR+zycR;
XPfx = sAGsu( 0.23951 ); %call for the functions
WpSs=YMAP*(Ckty)^(wgWGQD)-JnxAgt;
PolAGN=QraWAt/ZieKjp-(xTOEzu)-XWuiZ/xXuR/gcDu+(pYue)^SKqSsP;
ABRYy=gteoXJ*MOTPN+KkER-FHLv;
fxzL = BkcSh( 0.2905 ); %call for the functions
end
Vymxw=QOrkwC+Iiomv+bqyP^tJaAty*chpSFD/VYQnP;
CZDaR=POByxi-(KBcB)*rbGO+(NRipN)^qDPBEw;
end
for j=1:53
RaDJok=ZwNtAb-(tfvDoR)/euAa-eXIna^DKDQpf+oJAR+jupC/kNjgVo;
izJExX=DVhh/yMtH*ZMjNZU-OHFQdo;
Cftu = pHEGj(NfX); %call for the functions
RTQV=krqW-yHoJza/zsyXQ+YbPpW+(
评论0