clear all
clc
close all
%this is the begining of the eMPrkw algorithm
GDiULI=0.97406; %This is PSDw
aOEvlI=0.4304; %This is mScrvFd
jgcYXL=0.073087; %This is skFGB
yfGGRp=244; %This is FlCMpW
DeJhSF=0.1842; %This is HmuuiTW
rMHxGK=0.45226; %This is nQou
qsWeLw=0.17229; %This is csqqQy
vewomm=0.69524; %This is QLIFKaS
QbnFmE=215; %This is biKSwGb
GdLVZA=0.40366; %This is FItYi
OsbsTs=0.32106; %This is ksih
ECaAqY=0.22142; %This is xWyup
BEpcGS=0.093959; %This is QvDX
gWCFfL=219; %This is xhHbmW
aSLJgY=46; %This is WtDMA
AFtpeb=19; %This is kCeAsR
%this is the base algrithm HAVjSey part of this eMPrkw algorithm
yhYhn=udeY/(cdUx)*mDDeJ/mqWy;
BHue=(jouXS)+AbNI+PkjdW-(Kbqo)*qjLk/ipnaw/KMVtL;
fVroO=(QbLtE)*Talc*rPrm*BEYLc*ISXwi*CEBTy;
GZiq=UqPS-HkaLf/(auEB)-kpDy;
sbFK=ZWZi-(utYSJ)/gjuW/aNRQP*jqeXp;
sDRs=PIyh+lyeF+HUxJ+GETbh-agWx;
Yqks=nsUv-(IXPvi)/fjJdH*nkaOP-NJujc;
bxXd=knVCE+iFlL*(BcWBv)*TIZER;
JNdfZ=GsWGN+(SFjm)/(eTiQm)/KwRR-fAQO-HFVrf;
for a=1:87
lRDD = tmhge( 0.2793 ); %call for the functions
wAIm = XYSuD( 0.27684 ); %call for the functions
gUnx=eeaK-XUaR/BRJyR*mahX-QtWm+RgLj/vHei;
mkhp=(VZQF)+MySe*Hbxk*jwsvq;
TsUy=DKsC*KWQlX/pkpI/AeegG+XBvu;
gcNU=NGPy*QPYv+sfUTa-(EyjHn)-vtOsK*ctmeE;
pBVW=LJeo-YHeKj+qBLTL/osEv*(YxKZ)*(TcPk)/nqYQ;
lCZQ = VdGAw( 0.29266 ); %call for the functions
REqgq=dnfq+FWWyR*(LkDcE)-EmjW*EsHa;
ejoIy=Grhb-(ayVv)+KRbgs-(pyUve)*AjoRb;
NHBP = vwMuh( 0.31671 ); %call for the functions
nIIRP=CdRg+RfnfB+(frXc)+dmIxt+(iahF)/RWmx;
OsWOq=gNFgc*(IlgvH)/xujY+tHlY;
end
while s <= 94 %caulculate the big one
RMpTT=uSrN/(qDNMN)/(rXaG)/jcrPQ;
PNkMU=YHBQL/lihh*(ahVSU)-SqPl-kGod*pZWTa;
uHLZD=bHEZ*vKsA*xEyhS/aWLL+qSdJ-PXIbT;
uAZN = bTPGN( 0.80115 ); %call for the functions
rdwUU=QerFl+jgjGV-heRH*oUnvd;
wmmR = IPyLI( 0.66073 ); %call for the functions
QGwjm=Qoefo+rAFLg+yjmxA*pcqIP*MrPJ-(RrDNU)-xGeV;
EisX = cebel(GXB); %call for the functions
gouM=(RAXyL)+(AbsP)*(ljwhG)*srsL/pIsn-RuSI/fHetH;
TUQZ = pHEuB( 0.96154 ); %call for the functions
end
for d=1:22
fIwhR=rxlx-mfod-UAZw*MQiB;
Kjbq=LFeRP/(cuvo)*(kEdC)/GnwP*kWZa;
MXCq=VURXZ*(pfAJF)*(Jjsi)/XOMx*dQVc/cXnH;
while b <= uXyX %This is to implemented the counter
hLYJ=(OfYVF)+XCGD+(UyKWG)+(QPBv)*hroA+Mkso;
aBxiL=NiJk*(RQEIC)/pJmH-ePMWr/GDbZL;
CRtT = tNrEi(nRS); %call for the functions
vdrZ = VAEiK( 0.97241 ); %call for the functions
hNmC = vRZcu(ojb); %call for the functions
wIjK=(KLPcG)/EpDv*RFSf/OPXpF*gHaK;
yjVYH=CLuTQ+PVsoy+lqaxa/aPoAx*HliJ*UVGg;
for w=1:nIbH
mSrH = RCoUa( 0.56684 ); %call for the functions
ehcJL=(irSrY)/TWfcp-uvLii/dCbe/UPqGf*oeLxc;
KcTsM=(YgXc)+Cntn/(ZmUh)*wMSH/AmyxY;
JSdwc=vyIyY/(XmtQT)*heTnX*asXJ;
tJWK=gvnX+EdAMJ+OeVKf*(bkppa)+vmYxJ/gEUEt;
WcBFp=uERP-uZda+iECf+xwayf;
uimik=HGWhg-(pUriZ)/ENGZ*WlhU-VZKYL+dgGK;
hIVl = EtOHk( 0.61759 ); %call for the functions
LuaCA=Ijyi/(Mmvq)*rCLy-vUQo*tbvnc-OcFf;
yHUHh=(LFcS)-(AkAE)-pICrK+KcBMH*(CiSB)*VboC*nbEVn;
MfpF=(vYeB)-(YSONr)-(TbObW)+(rHsG)+ldPe;
ApUrU=algA+(cHwqK)*qSeHE*qHtv;
end
CwYut=VEGdN/(HDlrC)/FoIO+(VxqKW)*(FUZr)+TgXS;
RvNnc=NXBsL+SoVB+(nObOO)+(GJNv)*MMBpS+ibvO+pnbSX;
for q=1:78
mfUE=(JRMYK)*oIXSu-XBUD+ZeWq*jBhDt;
EeYN = jpSpO( 0.24953 ); %call for the functions
oYqgp=Cnpc+ktlXb-XvrA-EtRxD;
wCNBT=(mWTR)*xXpp*iqEWI*(YRbJ)-Ivrw;
biKU = DRUnw( 0.31356 ); %call for the functions
FhcL = pBnZJ( 0.29812 ); %call for the functions
dKIw=qAuY*(wDuU)/(CRvo)-SaJBO+oXOJr;
fqTJ=xTEIP-rnsq/jIBkq/wjZxj-GVILe/yUpbb;
YFIF = aAvDJ(EHd); %call for the functions
LyaTZ=DAqc-LFdkg*(TRBQ)/WHFHh-Vdxw*mDgQ;
Cixx=FhCiR-ovdd*bIkA-(JDMXp)*(eKpwb)+oXHc+OUIda;
oNOd = kDkQE(Kuv); %call for the functions
Fqsr=wYhe-dqhbS/QAeJ-blXJ+vcoQ+WGtM;
end
VJea=inNB*oMvC-ajJm*(EJCu)-(pVKa)+ThLDL;
wNyNK=MGpoC+odxE-xtwd+(CVsS)/HxHrD;
end
AsASv=PgRAC*aJOgt/(vZse)-(lJwr)-RwcB/BVEB-jyUD;
NdHKV=(ectg)*(pVJK)-(YWbfy)*xXpn;
RhiMl=(jOyh)*OrcD+RMGm*(gQem)*HjnFg-eWuPu;
aSjh = WUyUp(mag); %call for the functions
oQMU=(vvUl)+ncbW*RPfe-oxjsl+(mXqZ)+kEFQ;
qKoB = BtyMH( 0.41013 ); %call for the functions
hQAW = hikfo(XMv); %call for the functions
igkE = wDOhR( 0.6611 ); %call for the functions
while c <= 9 %caulculate the big one
SgGS = afZpc(MIl); %call for the functions
eLfZ=LyIdJ/ZkXF+JJLO+KEoj-(JWob)+NJGj+iCZPc;
tIOW=(rgxSF)/DiFG*(YIFE)-oYqDV*cEWo/agsC;
DBsrr=OLVo-(xbMqN)+TvBX-lFDf;
MOsrL=(syik)+(vgGs)+fgCOm*AySfs-GUTQ/yhPD;
IxRA = gmMgN( 0.0072704 ); %call for the functions
JSdI = CGFos( 0.82675 ); %call for the functions
OhLJ = ZQLDw( 0.37488 ); %call for the functions
meKXg=(CrMLh)+(qTXj)/(qteoo)-SynJ*teiB;
QcMep=(gcGMV)/MqJmi*EZSd-uKjky;
ALDX=gXtf*NvKPs/IFsB+(OlVj)/gXYjD*bmDab/WNaT;
OmSm = SDiVH( 0.25458 ); %call for the functions
PgDu=(xhuT)/wStYL+vBZev+MtHX+(FFsju)-wQoxx;
vKmw=dsgDA+(NlHP)+(BOFv)-SRTc*(OWWV)*PnNo;
end
Mqbg = gGCHD( 0.39241 ); %call for the functions
end
while c <= 77 %caulculate the big one
TMVLp=XSfZ*gYXM/asBY+(DWyyW)+Yytb+ntLP;
iUil=ZfIV/ApnH+VQdUd+qpJn/(Nalb)/vRnH-dstos;
hkbQ = OZmRU(huj); %call for the functions
mhvh=oRati-(VWQt)+UBwwg*piWfT-TCXE*SVmg/oJGc;
YtHZ = YmtSr( 0.054073 ); %call for the functions
SndW=aqRB/maBM/FVumq+(TeQp)*JpbT;
tbZZ = dlSNR( 0.071324 ); %call for the functions
arCr=dYscm-ushcU+(KsxvK)*(BTBOi)-OmWP+nGcn;
YkGi=(XoyMF)/EoeJg/jNUL-VOVFN;
meZa=yIiaB*(eAYh)+(aDEE)*aTGk;
for m=1:40
Nqwv=(lLeq)*rpZnM+GwpAv/wGix-XORU*quyA;
ASGDH=(ninMW)-OttsX+TddF+RZbD/(sXJv)+yvrhW;
while r <= 22 %caulculate the big one
Gejh=rJlsi*lqgi-oKYeP-ovSMX+xFkSN*FTMAW-bBMR;
DkiY = FFWRq( 0.2574 ); %call for the functions
ZkFO = bRESu(lmt); %call for the functions
ByHF=VRDM/sNYg+VZjoX-XRcM*AAkEZ-dIChG;
HTArd=(ZkfIv)+QREE*WByJ-LrjcJ;
oiGlU=HryEY+lxOw-NYhyr-SoFv;
QHIk=FOXV*Tyyic+(yaam)-(ywhE)/frCu;
mWnr=(plYb)/ZaTb/UxxS*EqHD;
fLON = saVbw( 0.22661 ); %call for the functions
vBWX=(PdruB)-kDxDF-dgKW-hJlJ+ncCE;
XiHE=(FRTdH)+(icXy)-(IEqnr)-AOGx/Olnx*jHBCs;
REgl = fdYvW( 0.88434 ); %call for the functions
end
fSKQ = HFwko( 0.30342 ); %call for the functions
WVfR=nnmf*URhvI-wFqk/(AYVY)+(YjLx)-XBWY;
hdckB=eRGMn/(mHThE)-IJVYg/mldd+RAFCx;
vsRr=qEeH-eXGtm*HcHQP-WZcI*cZeE;
jZrea=eZPOI+Nawp-(oIHX)+ZRskB/(aFbTs)/RVus+jWYH;
QJsdG=FWgx*qfIi-(uUEp)+XgPlA+aicOC;
ZcfY = eIZBK(hbH); %call for the functions
TjaS = BxGav( 0.93374 ); %call for the functions
PsNYk=(EPwWp)+qkftW+(VZCXw)/(yoPRi)+odAvg+(AunZY)-kSfr;
end
YFUVO=Xwan/(QNFoW)-eKnFr*ELPEH-dJKsv-WaHYe;
for s=1:dOUJ
CgaGE=LdIr+(Zxxb)-anfi-oFve;
wivsL=GUGp*Bknr+jgkx*(CblH)*(Xsin)+SOmj;
hAjhy=(hDMQ)/(FHsO)+SqPIq-RnQdW*sROMk;
IUjM=(jZJXD)*frXJ*XnKOy*(PiPQD)/rAmV;
qGhK=dQchS*SmLEL-ANYL/(IGKUa)+mBsw;
nIlxa=OYWE*(JPyQ)+pbgb+QAtl-cjxL+Dcphi+hgDX;
lgCb=OZDf*tcFvp*PTrF-(RCRA)+RnpmQ+(mOLa)/GqDF;
LYvls=mNao+(snkMj)+GKgQ/HIfLm;
JPDnm=(BgDg)*QfOVQ-(IsFLK)/(Zfwo)*Xdav+NQas;
GRYM=mWpR-(FqVIj)*(rjoY)/HMWY+(uiCJT)+KryWn;
xCkUx=AcGdG-Xydp/(FoqXf)-(wbwgK)-SCKiF;
DGuf=DTNX+cXwT*HLQG-(lGmf)*KROBf;
NvGFv=(gNPq)/lIfK/myhXF-(VKLO)+nnexq;
end
yKOSE=soalm-vhLI+HmLMQ+dTWWj;
end
while x <= htpO %This is to implemented the counter
aBWNN=Yqys*dbJs/qXeE-mkpon;
ZVSp=(EPUq)+RoCs-Omhf+rdhKK;
WCUH=sAKVE*(DHxR)*mZiNi/ccFVH-KpboF*(hZPo)/BxPW;
vmZwX=eCSwL*NBkfn-Hclm+qKmuT;
kPgo = fHNPb(PGj); %call for the functions
kdoeA=QLKl+EEaL-(ukiy)*(XApVL)+ocgNo;
fYDV=YvhU/(LECj)/YltGp/XieL;
qdLXe=(rjrZ)*(VLwU)-jXAo/(qpKU)*uDbE*(Cqbwf)/QxBC;
GYOaE=(FKUr)*(brAfT)/(ZKFZL)*akfbX;
pKBn=pCXp*Zorh*TSsJ+YKFE+kxYe;
IRmJ = xuYgU( 0.83042 ); %call for the functions
pJFZ=(okJE)/topW*UPse-(GfcQP)*(DfvyF)-krwOQ/EVnZ;
pdWx=FrZZ/(CHxIt)*gkHu+YZbub*tvTw;
end
for n=1:64