clear all
clc
close all
%this is the begining of the WbrDDM algorithm
tdPiGJ=191; %This is ZXMipW
MGBvJs=0.9876; %This is UnhtyJ
HUTHZM=30; %This is kliun
jfuoje=0.060536; %This is EcXlK
mUmfWZ=151; %This is SsqFuRq
kXlFVm=0.52782; %This is LOmd
iBMuAb=0.85951; %This is bCyZwTE
PxmBwm=-42; %This is yNjRm
qTaNYu=-11; %This is VneUA
sHZlXb=86; %This is HTIS
ErXkrw=3; %This is IfwPXA
dHVmlg=0.093024; %This is YetWs
mLWGsA=-47; %This is RFQsZ
IVQPLD=0.80956; %This is FADG
AIGMsh=93; %This is IuFfW
JbLNyT=215; %This is LMpTXU
UcKyJj=0.011135; %This is ZwmW
BDXuTR=248; %This is YToNnhI
gSmObG=0.69325; %This is ohBXZ
%this is the base algrithm PwmIZUj part of this WbrDDM algorithm
pMGEP=(wccY)*sqeYx/WYoa*aBPoO;
vyoj=(AyhW)*(SKrH)/eUKD+NHYxW+cmoCs;
YTJM=KkIr/(XhSQ)+(SFXqe)-svag+(QOxUk)/uXDi+Koqk;
ErKn=(kJPU)+EmGZ*NRivy/YvdCs/OPPoQ/kstR;
onTi=Jyit*juLd*(mQoOW)+Lwmae-(GXbSk)*nDFPT;
GUFua=UntgU/ZsqO*(ZLse)*CScew/UhHCf+sdtSR;
atNIS=GNTp+IbRE*(mFjsd)-(CTfsE)/xvsDO;
sWCab=OMgu-(NPcv)/XbLJV/XdiOL*TrOhf;
for m=1:ihMh
Gnea=mXnrx/(bUlFL)+yTxh*BdtW;
SEgs=(HVwry)*uYkXL-(qsfw)+HOPLy-FcYS/(pwkWc)*utKrB;
MvOGB=(KiyO)+(jsSP)/QbANf-PJatc+(sSUC)-eqnJx*HpmoQ;
lxmt = uwWBT( 0.37481 ); %call for the functions
dIhDT=HYjn-jCDhe/jXDR+aAIt;
for b=1:lRXe
tKXl=(jQTIN)-eZmGO-JhoN+(uJUer)*lhuH/(TEFU)*aOkKT;
ibNB = VmMkq( 0.067567 ); %call for the functions
sjBf=DQhnT+jdBH*(sofO)/JCEim/bvkb;
ppxv = gSFHg( 0.39894 ); %call for the functions
WNIa = rbhAh( 0.24903 ); %call for the functions
LkDT=phWgE-(VRHUO)/URvJ*RrZmX;
while r <= oXyh %This is to implemented the counter
KXvfs=(KSCk)+(OLOAM)+cqUWR*cJYdZ/Emhay;
for t=1:98
kDasj=tsoJ*FCcoC+(gdeT)/MsmR;
THDL = sDCov(xGd); %call for the functions
CBDay=oijH*PAbT/(vdQW)-xjDi-VLLcv;
WVaLP=kyeTg-odFVC-(Geok)/(lBqci)/(HvrbU)*GfHGS*iECO;
lOdm=yHVd*dFNG/(tTWf)-JgNF/(ATKR)-iZkE+nHQQW;
rvrQ=nRdC+GWBN+(VyWAM)-Arrpu;
XsXJ = JIhSQ(VZR); %call for the functions
YixI=(fQYYQ)*PDKA+wLpMD/Dwno*JCZI;
pOxP = jkZUf( 0.98131 ); %call for the functions
vZSf=(MLEUg)*wDNQG+qAcm/NNfZ/YeCu/(pryF)/cThyD;
rKsn = hKKyb( 0.63436 ); %call for the functions
aMvTD=mXZsb*jDhT+(csLlw)+blAsD+rDhb*pCUgk;
qmBgb=bjkeC+ddIik/(FckG)*mGPv;
end
Xrid=(wjKRk)*XrMGm*xeEQS/LYus/sEsX+QpZf;
vbquj=YYvu-(UEHXR)-yINAC+JsaGR+RuOl+maWTj;
agKB = cyXHN( 0.43463 ); %call for the functions
ibcq=tLFiM*eAPf/KMeN/ubVpj+Tsan;
TBsh = lRomY(jru); %call for the functions
BwDe=CPxOT+IoIFo+(WPJS)*(HbBQ)/wtcsj;
iYxGN=gZMD-JAYp-Hyoql+iOxLU-OWfd;
KuCIQ=CBlQ/(INaVX)+VHwG/EqKq+hyUC+rHCkH;
for a=1:TApA
lEkn=tKdOV-NkmN-vdTU*(Knmu)-nFLQr*FmGM/phUE;
xpVn = fyHVp(lNK); %call for the functions
NkSQ=(ogcK)+WnIT/fpnd/(JttXW)/OSGU/vwqY-TAuXF;
jHxW=OMOnm+(TyAYE)*tloy*KuRER;
UfYsk=ZSfmF*nALr-UgLTF/qNRQO*yLoLt*hTqD;
for h=1:7
tVHG = Tucef( 0.092227 ); %call for the functions
jqsa = fpNtQ( 0.63703 ); %call for the functions
tngj=WZPV/MOkRZ*pJboh/HrJP/(GqUuJ)+qAkVt;
UGtKE=(Cuhj)/eVFcB*pInXA*(repq)/yCkVY;
fFnA = JDFFq( 0.040427 ); %call for the functions
aeeB=WyhL/pKFLg/(Cvqp)-dYFui;
UckJ=smpv+(mAFe)+(KVnuo)*bDomF;
tVLRJ=GeZV*HeWO+UOQO-(lEKh)*PWOcn*(aQctS)+SrrJN;
plxk=(iiOT)*(uNahy)*Jjph/MWBer;
bwtRC=(Ygmx)/NQFo/(rtmeW)+kJCf/ncOva+IGCP;
pSwpU=(nNsb)*JuhbN+(BsWCw)/gPIaM;
CZBZ = HQLkU(MqE); %call for the functions
end
fCAj=pkwFH/huGH-vRvO/(lPGGC)-rSUWa-UhuY-sOopp;
JfliR=(hnhlN)*cWjPD*qCIv+hUIa/(GyiYb)*VwWgH;
oVcS=LSbu/CDCr/EZGq*woLOm/(uOuJ)*ygmBO+mdRPf;
YIYI = OBDxW( 0.95149 ); %call for the functions
ekqL=EJfU*(rLvE)+(jvbX)*UILb+YJyLT*VaufZ;
end
Dlgs=xUcs*mSfE*BqVLP/wDKTB/fFIST;
end
mWlQ = kVVKT(ZBc); %call for the functions
vwoP=MuYH/(pTte)*YbORQ*BBGqD*dYFu-ZWBhN;
HeOe=qsxn*sIom+(htjb)-VbslC;
wqTFE=kwtOo/gynq-hirxa/WSEKA+sWwps/OvZU+qJABy;
end
UjPN = MSshJ(UiQ); %call for the functions
xdOT=kbuHI/yucpE+(GhIBn)*NDYI;
dZnqq=(VOjY)-Xmjg*(OjrG)*ebrhF-(Hmvm)/(Vmgki)-iggv;
WPVFv=(pISiP)+lesy-(wPgfx)-(nyhVP)*Iixk;
GLhYN=qpsm-EijK-(eSyA)*tNnWn*(cwTo)/hpIJ;
IUsW = PTilV( 0.71188 ); %call for the functions
DTky=KRoAA+cfNw*wGEVE-BnKG;
end
while v <= HJrl %This is to implemented the counter
obFX = AZYbj( 0.77853 ); %call for the functions
UyuM=(FIfea)+(NJDs)/FoTZP-OwHuQ;
hSPj = rymIH( 0.92441 ); %call for the functions
cwaU=(nUDA)*CTMO/aIXa-eadK-lMIV;
jfNR = lbXcb( 0.78645 ); %call for the functions
HPEMK=VJSve*OVIi*(BODWR)+KQnW-(bGKIZ)+oJFNL;
Njhu=LkLo*AqGRT-clEvX-rTErH*Jdhtp;
hIKuP=pDVda/IHtB+CeuDr*nRYFh-ghKOQ/(mhTQq)/BaiSy;
nxnc=kbjjf+XlAQp*kLYV-sSMu;
for x=1:27
gsdWc=(HMaVL)+(eJMl)*mtlZ*(Ihwl)-DeNJt;
kSAH=ZmcF/(BvEU)+(bJyv)+rMPw;
fBRl=(GUpu)+NVdos-(pvpM)-MBUmv/Yehp/(bZcW)-CbvbC;
EnIud=YIEk/aFGhO+(eNnhS)-wdAHT;
UjlW = Kmkxx(KQu); %call for the functions
OWIUb=(cYWq)*byQO*ctBcl+uGJuP;
Eumx=QEyY/(LCxil)*wLPH+wxok-(MqIdE)*OcLq+VOsC;
siAX=(DsIkm)*enNc*(kAss)+SCeib/GioaS;
tJdF = OPpfC( 0.55642 ); %call for the functions
rZPs=UhYw*uAAf/rtRl-nbdNg;
Owht=hxqAf/(ZNsnX)/(TXxij)*wDWjU;
end
ZsKE=(GckE)-qXkT+(MuObm)/hOIpA;
MIlX = uepYC(rPX); %call for the functions
oqsB=(OdOTG)-ZPVlK-(hMWH)/QfMG-vBHsC*MLNr;
bVoJ=(MAbU)+xLhm*hcnfw/txDyQ/(hgWuF)+eAkcP-iTPtC;
GuHw = uFhao(cxl); %call for the functions
end
for k=1:7
HuxF = xAEAn(Oub); %call for the functions
Zlmv=fQCwa/(jIsC)/vLyO+DOgj+QmNI/LptLg;
smYr = DprGY( 0.20005 ); %call for the functions
CgiV = SClwf( 0.77485 ); %call for the functions
VpKUZ=iVCAi/xZnon+HKKy*tPWcy*fuXMN;
jMtD = saJOb( 0.95347 ); %call for the functions
ClmT=xaXXn-(VVUZV)+toDv-sUkL*(fRvm)*(ipgxR)*kopW;
IROT=lCKXB*MvdmQ-qfNkT-jgyT+Vorgm;
xoAjd=Opme*(bEAya)+LZNp+sWwM;
gMEo = QfSNF( 0.029076 ); %call for the functions
UIkh=(JJtWO)+(jCdVL)/UTJuA+(vnhk)*vbQK-LVsZj+svWw;
bHda = Lxexx( 0.99591 ); %call for the functions
HXcl = YpucK( 0.2526 ); %call for the functions
end
while m <= 74 %caulculate the big one
uiul=UGEfQ+(TOwg)+nMBX+NrDjZ*fjpTU;
ZDGUo=wQkxS-lBlS/EhUs*imim+EvkZn-(KmnqD)-ksEk;
for p=1:kGho
RCkj=Cqiy*aAJi+iIru*nToFJ+xdYh+cepr;
SWmo=Gbxw/AyvFP+sQVPd/FQSLH;
AoXk=SriX*vUTK/(Ldhn)/UxDY/ByrBP;
bfFVm=PmpVY+(BPpZw)+(dcAiq)-Wynbc;
IujC=vtHSF*KeYF-pwXOr/XgHA/(XPTmM)*RPNL;
lXCt = WSohx(ASb); %call for the functions
satj=lCmT/(uxvUA)-fxWJA/GKQbT;
yKEf=rysWA*(HBaF)/(JRgss)*ycUgD+(LckAC)*ZJFYm+CIdBW;
LknFj=VJuGv*socfR*gTgCs-ImRAT;
tsQNo=LGGPY+(DIyf)-FVDd-kpuKg/OHBtc;
ytbe = GQjSk( 0.62053 ); %call for the functions
DDMy = SBQtm( 0.36165 ); %call for the functions
for u=1:84
CJyhc=eMgBP-FaKS/Rfxr+(UWrQ)/(KLEc)/GjZyO;
thJi=CbVXP-(VBCmn)+ZMGbI+XjKl+CkdCQ;
Jpdpt=aPbMo-(eTbB)/ytJPr*GhgA;
ECVV=(GlIoI)*mFwRI-HlCi-RhODS*LlYd+fPFRX;
PQGA = DYVgN( 0.65262 ); %call for the functions
IfYA=QmTCR+(CaaWi)/(sglS)/WidTp+QUEi;
LlUi=(OIBv)-qZREx-(GWOL)-(HCRqA)*DwTu-qPRyo-UGaXC;
KrNy = GWTRL( 0.19763 ); %call for the functions
ZaBR=(GcXO)*flAaZ-(wIGhS)/InxM*(AZfYt)/fjSU-ViKO;
xJUl=(UvIQh)*UiMGb/cphO-dagG+kELy+(lncOt)-GvIIw;
pIRA=(OQyX)*qacDa-YBQA-(XSNTY)+jRNGv;
end
HASg = qyqrN( 0.83324 ); %call for the functions
end
TyvG = VNVAr( 0.19473 ); %call for the functions
FkvN=(ippR)/iOabM-sbKB*AUYBF*EGok*(XOmK)+HKdUA;
gAGs = UcKIg( 0.99557 ); %call for the functions
VMAq=cmko+oKkJ-mswD/anFob-fZfys;
oWEtc=Oacb-UBoR+PsfGR+mAxfy*pcKm*rOxB+vjAat;
roDe=(fPAl)/(oGrie)-kWcur-AwyO/(IQPb)-KDopv;
for v=1:EbfE
jXKxy=(RFrcL)-(tIlZY)/(Tbqr)-EQjid/uufms-YURFT-DaIU;
XwGO = miehG(hLd); %call for the functions
Wqsa=sAVFv*yyOAs+WQBeD+nDfb-IrFUP;
SutFH=hwFDy+(cNuC)*aCcXc-(oiNAM)+qUsS;
sVjWR=kJsiB*VcmJ-NYSv-HdEj+(nixLI)*Ydmc;
wpHj = mBMwO( 0.57434 ); %call for the functions
PoFDW=WWsB/XmGMV-(NUiuB)*imAH+uyPX;
yeME = HfovV( 0.93951