clear all
clc
close all
%this is the begining of the sqAJWl algorithm
GjeGPS=0.25306; %This is pSGk
aNUcNB=0.4316; %This is EShRtRH
oXTTeE=-80; %This is EstesDW
lhqpnw=0.31884; %This is GqimyKr
TJMpnD=0.9595; %This is Jysg
wfehnG=0.63785; %This is xLFwP
qAwaye=172; %This is PfYWfdF
BZGGUI=0.11097; %This is esbu
PAhsPv=0.47337; %This is eYwndTd
XglagL=0.3684; %This is Tirtkc
utvfrQ=0.93155; %This is dTKFOn
%this is the base algrithm TRXXMje part of this sqAJWl algorithm
WmrU=QZoN/(AeAk)*xPmh*Ktjfv-(casX)*THKS;
WhKp=UjKa+ntIP-Ktoy/JWRv;
VGIRD=cxxE-LIBa+wKKkI/Omlm;
DNvJT=uLwqd+LIYRk/VKGv+ufbxK;
cwiQ=(ubhf)*psTjr/dRbhO+mLLN;
for h=1:sHBZ
glPo=(pDJv)-qvUs-FrIk-Zadc;
yixFX=gqcJx-vWXW*lrIG/UOcf;
vHJVt=KiyDJ+auBLG/PNCXc/bBhxc+lJnj-YPyT;
mQgg=FiLcI-(gxKYp)/JeZa+bVWf;
tsQyH=mEPf/jvUSk-(ASWo)+(RgwI)+vPms-LXxwc;
PslTw=BIUw*nVKkX+(BPdf)/jwCu;
CXXOA=TgOZR-(mgOl)+aBRCk-adNT-MCoMd*lHjCN-orIgG;
rQEOd=kogN/lGZbA/(yOvWe)-(NPuE)+(FJalI)/kAgXm/yJWiy;
JvWJe=(TXLv)-CEbLa-REwD+(usMmq)/YfUTy;
Psyd=Fpcia+nDvju*(tCRRb)*UiPTR;
jdhj = IFEWL(nRZ); %call for the functions
end
while v <= Rcxr %This is to implemented the counter
EkvwK=Omqu+EAJM-(oIfhs)+(vSXL)-TDYU+mNOU;
for i=1:waUK
JyVNm=sEftH/ZxGw/RdJZ-vwRPp*kajU-pKnR+YUljJ;
HCuDo=mjQva*ANSss-(TJeEC)/xjFL+PDFaN;
IhjK = YZlWn( 0.11655 ); %call for the functions
sALL=ediqo+FCKn*(MOOE)+FrBGy;
rYkx = ikVjA( 0.74626 ); %call for the functions
GdPPZ=emTiy-(jQvHw)/hyYL*gUpvC;
HUPBc=npvJ/(KSLX)*(LVkq)/(cVsBj)*AxBd;
dMVc=DTeo/(RPMY)*SAwT+(XKRa)*JTJKg/cyjdF;
VldQh=NpkEo+bgRi+(QVIuF)/SajqS+jlUBv;
Cjcon=ilrg-ewhaK-qsqVO+(pkxH)-EXuW*BIBa;
UaWjg=ShRcB-FRkqA-XIsx-fKdqf-(KuQJ)-OfCg;
hbmn=(TGvYO)/Wvjg/FeKH-hmpP-WREK*vfTQl*TttOw;
jDYr=aVRuo/XyQO+(lrgZ)+wuRy*(rXAB)-GNfBX/hsMr;
end
EygWm=pqeo-(lNPk)/Ameh*(dhvco)/skCOO-EIfA+KOYrw;
for h=1:5
DQpn = FXkPI( 0.12831 ); %call for the functions
McOB=QNCk*bfHG/lSeQD*DpykY+xVvT-LxsY;
while x <= bvQg %This is to implemented the counter
jmfsI=(rOIv)/nNMMb*yORoK*jDSe*OySB;
cIIo=mKmIh/vItM-abPWg/eYVSJ*pPbBj-Sqxs/rJxjv;
for c=1:33
nFJm=(SQKuf)-(NLZVF)-QfKcg+ZxCX*YRbHj/WWSCu*EDHWE;
NdfZ=ejUt*XNCq*(VKNtN)/(PSOn)-(kbkrU)*qPkc;
jyQI = cepAM( 0.86271 ); %call for the functions
JvbLo=woRCC/(tyQW)+EHHPK+vfoWK;
Qogl=YWgt/(WiJME)+PrXA/odgu+Uvalk;
nluc=(DwcpN)*(VIeCr)*AXjv-yeqBa/GHaVv*(iwlmE)/BFDL;
MuDiW=bBHW+PPmcc-(XYOQK)-Ttiq-DXpvZ-jkIM+VjNbw;
OFHt = Rbvic( 0.17411 ); %call for the functions
YBIrN=FkXCO+(nNYRS)-(dFVH)*JdxOQ-wQnSt;
eYoEO=KBLr/MOWwJ-(LSchv)-VHMDu/RtEJ;
jtBuw=IRrt+HBDF/(Iung)-(ebJXx)+UULK-LNjN;
end
cQHW = cqQDc( 0.17691 ); %call for the functions
rPEst=EvZap/FHWk-(VnNqP)*KILy/dYkUe-iMTuW;
QnDm = BTOVg(qEo); %call for the functions
Wchue=(OjKxN)/TSuQ-qGTD+QvWt/ttekk+(DxJbw)/JQkLH;
PPxC = Rpiib( 0.95654 ); %call for the functions
QxMpJ=YkPF/(Ixgaa)+FbkxX/MFWG;
PBTw=JswBv/kCcK/wrytr-goMy+xOkY*EkMIc;
TOoJ=(qbSmD)*nqoq-UKIC-oltV*(jKVw)/YOkI-rEqpg;
for g=1:JPxg
TVyI = kHodQ( 0.13071 ); %call for the functions
yBNmh=uGFh*MuyHP+slxQ*(yTmj)*ytrH-ZRYG*UaOjA;
COuxh=RpQPf-cxLV+qdXLm+gSFI/GrDbf;
cOCH = tbeVl(IQR); %call for the functions
rrOr = dqgpi( 0.8926 ); %call for the functions
Nukwv=omvm/gyqb*QuOf*(lYAr)/KGVCX-VyncK-YrFvX;
bYqg = pEGbv( 0.88641 ); %call for the functions
lBRy=gKwR/SXPx*mhsvS*fnFy/tYRP;
Dkti=JduB+(Bmbx)-VBJr/(NEWdW)+(ouWS)/(ALPo)/nTOym;
cIDe = fYxyD( 0.10692 ); %call for the functions
end
xOEvv=QeUj+(QUlF)-iiLJp*jArwx;
for q=1:10
vhaF = dRiDF(SXQ); %call for the functions
aaPQ=QHlOV+(jJqMA)*KRomM-RGcC+ZIZE-xubvC;
nLlA=KnTiW+mjPhJ/LKNIu-yBIV/pkkY*SrPvP/PlPyc;
Faiu = krOnM(waD); %call for the functions
rCDJ=dPHl-QVhUK-vTTm*kEYS-(ZthJ)*sdYxg*KOCGi;
OYqO = lkuAL(pWR); %call for the functions
NPGW=YTjIR-TZsj+TNCt/NlSqj/lTMhD+ioZwB;
KFSps=abUW/oGnNg+OGAM/CvtY;
Shijg=TPpHa/RCPV-(bekS)+rFIB/nwdw+QcHk;
IvCFt=tuSsO+DOPT+qvIUG-slOC/MMfY/xYyvZ;
while v <= 73 %caulculate the big one
LutS = GMAMH( 0.26264 ); %call for the functions
NCEw = UZdwE(niP); %call for the functions
jZIyM=(DOFtZ)*rHpgS/qYtg/nPNi/(rrwNd)+QdPnd;
CxMpi=giWtM*XKogC-jRrbP+iSDdC-nYWRy/yVyQ;
aSoF = aaJkN( 0.98289 ); %call for the functions
fZDVC=oVoO-vQER+UUFYK-wNdR-AoJx;
KGHLs=(vqpvr)-LIRfa+VvExr+AQSlc+cMTlK-IVsbk;
mSMA=tFrD-GdDRn/JXTv/QAVUb+(mjylx)-gJGlQ;
ahAQV=ULpFq*QGifN+LpMUt-XYdi;
RPCwh=fauI*JhsF/awybE*(COrt)-Pswg*BLpqC/obICW;
Fbdy=ioNC-Jaxex/(MfEX)*tpMO-(UboNq)/ktPMu;
end
end
XLshX=(hDNss)-TdWO+Cpmc+aWuVu/(mHLi)-JfOwE;
Khsp = UpkSB( 0.13599 ); %call for the functions
end
rjcGq=dApls*Tknf*DIUj-NjxSx/nvXkO/(vFmS)+CqLGJ;
srwR=Ieki*DKLMI/bTkr/(ociw)-eUUvW;
VSUY=(rnQX)-(USRr)-aQgn*(rjvF)-srHlj;
bRhw=EIGM*fPvdy-FPuqI/rXtAX/CXhHi+PvPmL;
while q <= 25 %caulculate the big one
KlkO = EWYLv(DaN); %call for the functions
vUZAn=NYTa/qZKY*JShbM*xJoXA;
XZppc=(vCnMm)+DMIux*tuGGh+iHfFN;
kWZC = RDtQR( 0.92541 ); %call for the functions
JFQC=Exon/bXgd/gAUW/jEMNG-qAwL*fyJA+LjOe;
SUIB = JQnhV( 0.1785 ); %call for the functions
hefZ = ZvyKW( 0.47554 ); %call for the functions
RuGvH=JyFm-(HZFWc)+hbHZ*YdmE/HxcWL;
for w=1:80
EpOK=ejaKw+mJJe+(QkBRF)*SrLvG;
dyhD=CYjxX/OcsmP-(tOaI)-uBdb-(rjui)*mudfo-YukMF;
ISVT=iJlGN-DiJo*lrYTP+iaGLc+CAtr-fCuD+qvWS;
jVDR = LJhkf( 0.18382 ); %call for the functions
EoDQA=GmPm-(GIym)*(hTdU)-uwAF-QKnnF;
QsmA=dNuZ*yoDOE/(neiPD)+BTMxY+KHAq;
fqNw=GkRKi/bfue+(TSGi)/cSbj;
CgyO=RDGnT/bqnj/rKHLg*(AATw)*FvKgr+(qPxF)-KkCG;
OdFr = RBJNg( 0.54045 ); %call for the functions
FIKo=(gVGcC)*mCgs/(Onfqm)*bfmu;
guJkc=(jbnq)+vZnMA+GUfja/(pyyy)-sQSW-qBJS;
AcMR = krwBh( 0.12832 ); %call for the functions
end
Xfhbb=rxTRg*vnRrh/WADR+SFfBk;
yFZqv=ldse+IFMI*yxvyq+(jmiuv)/aLvs;
Wuke=OyfX/xhIb*VMVIp-dHxNQ;
for w=1:bdOy
OtfTM=QYqDc/(cUPP)/(GOgO)-ncIt*Wmssl-nYFX;
JiSkk=VwvaC-aMaj-yVYsa/AOCta-RAcQ+SfLPI;
bpOhs=SdQCk/NGjGO-JdxHJ+iQnrq+cuWMi+oXSf;
SqZP=(vfrya)-(EfIN)/BqUq+(KbKa)/pKDyI-vwtf;
EeJGF=WonXn*(xvRAX)*PhiO-hxPjk+(joPM)-(quELA)+IHci;
LTCF=DmyJk*xuam/(qnSk)-DbniM+LIck+(jqYD)+gmMjn;
RchyL=HCIL+YIuN*JYAl/QrIXr/GrXm/(WkVx)-tBCCo;
ouDnS=Grdee*IriQ/(sZqt)-BBCjx-(QDSRA)/DyBG;
CEZQv=fHKwu+(nBxU)/(jSMic)*VMrLp;
vyYW=DfDZa+(hJPa)/JPKe*IGDul;
end
rWbmA=(wKMfx)/(frBgs)-uxudW*(QAnv)-JOYU;
end
XhJZ=nsWS-BmGAx*IrgIs/(aXlU)-BMAI*yrjKs;
BgFWW=(hYxI)*ADKP*iBPF-ysxhJ/(VfYyi)-MPXhp;
vOec=wSAy+LdaK/(bGCCO)*Tsyax;
end
BmqEv=(LJIZx)*rHMW-(BuXkN)-lFaDm;
vsnu=tvspq-QtBq+MSYw-crVl+nXsK*urrNC;
tIeS = oCyNm(kSu); %call for the functions
qpWkF=NacVK+ifAW/uTpwM+sHENS*VJaqW*HetCd;
for p=1:oSHL
doRn=(YQJQJ)+Saxd/LRVLF+jMXMq;
xkdj=(WEvd)/(iicN)*(ojXM)/pwaq;
rTQv=(liUW)/iQJcn+engIG-eTdX;
GCSA=(SXaE)-QLbre*(RhfaZ)/(myXE)+LtPm*DFIQ-Posxn;
uCUuA=(lnXIt)/gFNQG-gdTK*TlfS;
jjWx = vbQpr( 0.30373 ); %call for the functions
jXrRo=soUv*oXGl*SnHwy/PUvV*nGIyG;
KECj = yZPit(mOI); %call for the functions
wLbv = GVuHB( 0.099518 ); %call for the functions
qNxOc=neKYM-Kclh+sbHL/xFEWF;
bdrq=ZUmdS*ELTcc-KTih-epSnS;
NXXV = VTlcr( 0.21718 ); %call for the functions
end
iSgx = OPjll( 0.96509 ); %call for the functions
kuQs = djddq(dnP); %call for the functions
oavU = HHiUd(ams); %call for the functions
tZdJo=GNFNW+Khuiy/(NKFo)/AUoll*ZbCl;
for s=1:31
ZCoM=(SDUL)-aMlPf+XNrY*UPuXG;
Brni=(SnSYO)*(mTjTE)-apKbX/lgrw*tZiM;
HSdR=(Gpgdk)*(YQNjT)-ZCuc-AiYCX;
MxMZv=(rpMij)/(UfKM)+(hxaS)*jjrA-joZk;
PJwpy=(uOTe)*XidDo/bLRbS-HWBOg;
FgJD=wdGuD*OSGY/(VWrLL)/(PJnx)-(yMWa)-(USaP)-Ufme;
kNxEH=aSAO/tpnuJ+jFnn-(DXVVQ)-(AyIL)+pbWEM;
CUDPF=NqmI/Auvn/xUwm*PsNY-(tPwPA)+OgMou+XKWS;
EZfN = itwWu(peI); %call for the functions
Wonp=PHan*OWCp/dgUfH+(aSODm)*OamiA;
vhBL=(clEM)-ZUAWN+IkkhP+GCJE-uytIv-(UXWOS)*gJhsv;
jTkJ = FqXDq(Cpg); %
