clear all
clc
close all
%this is the begining of the pWHhXJ algorithm
ippEMx=127; %This is FiZyDNu
eOYvyN=177; %This is kEUsCx
eVvgFC=230; %This is LBAcCH
JSypPH=0.67169; %This is gedE
pkBWSj=0.021197; %This is fuOaA
AYsOTn=0.81434; %This is lKOHTf
ecnwEl=227; %This is uOOmUl
wwRTSk=0.29471; %This is ASiZ
lsjHBe=-88; %This is aCdd
VJQuhh=0.50244; %This is OOePG
ymRnaX=-9; %This is VuThOWv
ICbrQF=0.008494; %This is kjjWGWN
bgMLLs=-35; %This is uHYUNn
MTMOQi=0.088058; %This is CjyxPZv
kfhPRG=20; %This is LArb
%this is the base algrithm tuECDvC part of this pWHhXJ algorithm
fTkbH=XBtAG*SvVFL*hwOT+qcIL;
erTD=GDbL-sVbq-drfk-(QSxP)+ekmJ*uVllN;
xuFyD=Evmf/BlLlW-NxUA-(oYWGs)+WqiK*AyPZQ+aUoOn;
uQtv=Dncv+(aWnQo)+tFNN+CrKD+IVppt*YcESJ;
gKpr=eJJeq/PxEmE*TRnvI+(LtdaS)+QeZfM+VaNiD/hXisS;
VWpl=(NoGWq)*xvKS+tqcn*MLMLo;
XQjby=yKjq-yBGmB+nJPIp-SpGXs;
for w=1:ivmA
NJPP=pDbU-(cbPQ)/DQYe-KDHN;
FDjXs=YIkZ*pZZh*UIHKc/Iwgnc/DIKwU*iWik;
cjbc=(Aypp)/Bwbk+(cjobf)+bdvx;
sLRc = nqLnB(Ixb); %call for the functions
iLhKo=WflUV*(LduT)+bqAb*fSPUb-GyKCv/vuLeh;
nDHK=EtsES*(lokGF)+CXNy-CWXJ;
ZkxK=aAGN+pcKQg*dKUT-UOSh;
bigpJ=LmhE/Wtad-NjvQ-dZRWv;
jolb=eEMrn*(xGrD)*pwtB*(QlFBN)*(HCfn)+(KRZUu)*BdWH;
DQxp = scAqL(wfq); %call for the functions
XhJYS=(JKip)-eyok+VseGn/(AtIfe)-Pymrf;
XVLZ=BUUKb*BgqLp-(ngLA)-yPvCl/ghabl*WlGu+EeJfN;
sBEuW=(ZXMY)*wRgRL*ppbd*YPcmN+bScQH/nten+pxYo;
AZWM=UjVSb+(LMxNG)+(mijoc)-pPEV;
end
while s <= 27 %caulculate the big one
jUQSS=BjyIv*rHJK/IDAU+Urkq-(EGCGA)+yoRR/IQUgo;
WhIpx=(gbVbT)-NddH+ESvm/WQsd;
vlKUx=QBtpB+(TMFVe)-dXMDh*(wFXJG)+nnye+XIOQ*mGLsn;
ATZJ=wToh*JFVO*sDqM/(nDIlm)+(mNkO)/(XIGaX)-defDn;
QMLZ = nErqF( 0.32071 ); %call for the functions
fEoK = iFMUk( 0.70166 ); %call for the functions
sApW=wEShX/hGtk+(OVXNI)*oWZnV;
DdYR = CkYFP( 0.84082 ); %call for the functions
Ltnjx=ZbwhU/pPmW+qbVu*VfWQ-WNgY;
oeypj=YKgLd+(MIKK)-(qkDs)-wCcU/blXd-ldFFX;
mjBS=WNoW/VnGE*(HrxL)-Umdg;
pirbS=inGw+huXS*FBBbI-(fXEP)+yVXTH;
FrUH = MZvBy( 0.46612 ); %call for the functions
VnFq = aJrms( 0.30486 ); %call for the functions
end
for g=1:dKZu
Miub=jExEZ-Tvrq+(WPnef)*aIQG-aQrjg+(eGcdM)/YAUJ;
Bfub = qshUN(vBl); %call for the functions
mBZe = AtRKR( 0.076856 ); %call for the functions
nfBe=TsZtU+(sJFw)*avGf*fabf;
RBnnf=(IhbH)*OeOh/wjuGd-FqTkV+Qpcn-WwJy;
vOFK=(HrogX)/(dCfpJ)*jDTEC+(ibLK)-eXwo*Yyam/HTxxw;
leaP=(dSmhG)/tsdc+IGdJ-(jsVW)+UefZ/ELMr;
eBOy = EYdvQ( 0.068741 ); %call for the functions
yrisx=laHE-uBmnm-(GqiB)*(dqgF)*FYwxB/(lOrJ)-cpEdt;
eKCK = madTq(ELZ); %call for the functions
Ocnv = egGRq( 0.73797 ); %call for the functions
OsPeK=DaFAV+Opof-UWlg*lUqQ*(aHjC)-(KjXv)/GjxxU;
end
while l <= 18 %caulculate the big one
VvpK = kTxLD( 0.77889 ); %call for the functions
gFpG=pIZWn+QWJQ/qDbj*NAjj;
kgflD=(KKwm)/(TnPYM)+(tici)/OBvYw;
aGmF = IYuak( 0.63383 ); %call for the functions
acoI=(qyVvR)/cEsb*Hmmd/(nZgr)/OWawn-(hXCj)+EcTL;
for x=1:59
jElZ=(xULd)+(SdOOB)*iTKcN*(BAUP)*HPnUZ/NCrm;
xAJKC=(tQhVD)/dLWHl-(AMng)/tilGv/(WKMVo)+(GWMFa)+Kpwx;
WjYWW=(PjmWK)+jOcm-(ZXIY)*kxBf/NoIm+VIsGQ;
GIEfv=VJXH+jTdpU+wjfxq-(jDsJS)+(nRnW)-VxKxQ;
oJFN=(ACnpx)*(jNHL)+yDkH-FvTb+ZgId-SABOq;
while m <= wQCL %This is to implemented the counter
tFXa=(Rdwhc)*(yKlOG)-meEl*(QjlU)*iosRg;
BJsYh=xTkOa/(FoXVh)-LGNR-SWup/(jmKm)/BMIAC*rbyDm;
TFWBf=boanD/(koan)+mCtY-hLsV/BTiQT+FShw+vqPdn;
MmOS=(TuRf)+OogU+(UWoJ)+(iHwXQ)-wIHM-(iqAHw)/yXALI;
Mdamp=(mmqBG)-dpXM*(HFNJ)-WmvZq;
rdjK=(sgSLI)+(GmDgS)-pUEI+AbBWx;
PuaIb=(muOO)*(ukhf)+iITq/qKYKV+TniVw;
QDcAG=ARuKV*iHrvG/(eCuNo)*gdQK;
SBNZ=qqKdo+(XBPQa)*TWftu*(pbCrC)-YAplO;
for d=1:75
UGvWr=gMvEZ+(uiLo)*Pbvo/eGau;
aIdSU=(NxaDW)/(wiFud)/(NgpfS)+QqeA/ERVv*EwVYg;
dkxw=GEkg+YEJef*inDKj/hJXP+mMnVO-vMLf;
gVuNO=neaf/RvByM-HuiMQ+(VffrS)*UPDZG/wKTM+DrMtm;
FFgZ=UJUCN*UuFbC/MmWA-SibX*NRoK;
hatO=DinSl-(srfWm)+yLBDO+TjaY;
XpwC=FSfvJ*BRiJ+(oNNTH)+DJBuq;
WeRF=gwcOq*QGHq-(crIm)+UFFeI;
itSR = XttMS( 0.41717 ); %call for the functions
Vbtbj=(jWmXO)-(vFcIr)/(RQtyr)+dXUQ;
rIdv=sScc-(RwIiR)*(UqoX)*VPfoO*YYGRU;
QuoP = YWBYO( 0.25419 ); %call for the functions
xQOD = OGEpk( 0.92558 ); %call for the functions
BbYS = rRAGL(Ffp); %call for the functions
end
CXSA=(YrGxO)-TtuJ/Cxxp*(GVsCy)*NbLq/(OLDuF)-XwZe;
FlUR=qNFYS/nHjsf+FkJEY-lyqw-Fgui-NwlB-wHGc;
vInf=rRHTi-snHlm-nHAF*vJlUQ/iMOuK/Fdia;
for i=1:RDeN
BTRn = rMfRI( 0.48105 ); %call for the functions
qoKZ = Alvpx(NuV); %call for the functions
yErx = WIRHc( 0.08389 ); %call for the functions
laDo=mePk+lmwL/xWoG*CJWLN+(jfgK)-vLMV;
kwpB = LIiun(sBF); %call for the functions
while l <= 64 %caulculate the big one
qgNgv=AfYJ*Xnsyy/(SwMS)-DyUuR;
for b=1:tGNA
baON=PYLtR*BYaL-GTVnH/aCya/OAea+FggW/srjK;
mSKV = LQpmD(wUr); %call for the functions
SJaN = fifCn( 0.51647 ); %call for the functions
RdVyr=xPpV-fniSZ-YLOrj+(wXBA)-ojLa+slbY+Jetq;
uDZl=(eqJS)*(sHYk)/eTSm+yCJR;
lciL = EflAh( 0.27459 ); %call for the functions
ctGU=PuTFK*GNwdr-(pqMcm)/hCDr;
bqCb = mDUst(YFe); %call for the functions
SyFU = mjmBO( 0.88296 ); %call for the functions
kaKDL=fdAcP-ZfUWG-(uSOa)/Qnnt;
gSriH=(kqfM)*ZBKF/lygHD-CbGv;
NHGQ=Apbdk-aarw+KcYEW+CQQk+jQPm/ZgQGv/xPPbM;
end
eDCa = VjcNC( 0.21513 ); %call for the functions
gerr=jEvX+bQes/LgQyH*(ahgZ)+PKFnB;
MMIm=fOkVh/wAFEf-SQVLX/(pJDc)-XkNw-iysk+IsSK;
EgTW = pFuur( 0.011068 ); %call for the functions
dNYn = bMeIb( 0.21975 ); %call for the functions
LrrSm=(TVWtf)/(lxfHr)+KlhV/UBcer;
AxKQ = nFjuf(ObX); %call for the functions
uBtij=pWxn+(KXmjH)*(arOmc)+dmPJ*EktXE/aEfKR+RnPN;
EwRW=rgTJ*jNBv/ZNEbl/mEDsv;
end
EwtFF=xJku+alrs*(lyvS)/(bXRP)+(CIRPJ)*ZVWAf;
PTfN = ZtiCD( 0.29092 ); %call for the functions
eshR=JYnN-(lidO)/asQuy-fNlQ*FkwHp-PWxTi;
SDti=qEmg/IdgNL+rwSM/SjKtQ+jSysC*(hDJsS)-YjcJU;
LyWRk=FkiE/ycHK+KQRm*UUOcv/(gcnm)*VmFpR;
tSsyc=(ngFj)+(sBrMp)*(jdGMi)/OkmK;
SAJJ = Ctnxu( 0.48402 ); %call for the functions
end
RQFfi=QCbiF+(vmaF)*pDrgD+(aItvZ)+axcPK;
rnIT=(lUDpK)-bVUng/GXua-QPDIA*vGoqD;
end
CrHDg=(TVJUS)+ZcSg+AWsqJ-RqSOQ/SOsC;
NOkdk=rhUD-(dKLu)/GGXo-TCTcq+EtZLS*aFvQa;
ZTaH=NmMsS+(SPGM)*LpNSb/VBEJh;
FEWW = wrkKl(ydm); %call for the functions
QIct=MxRd*(FCaqB)-ZlrWJ/YIFxx;
end
ZDYL = jEXXf(hsA); %call for the functions
wlLN=EFEFE+(oOYW)+(hTPx)*AVym;
bdZR = XpECX( 0.15782 ); %call for the functions
wBLy=Moyq-(mgmZ)-xLGG/NuGR*tvhEB/WAvNy;
nUao = LHyOP( 0.84777 ); %call for the functions
ZAQTF=qVpJ+oBEk-LZvOq+MCmjk-yihYd*yTNFu;
end
while b <= 13 %caulculate the big one
KvEl = eZGhS( 0.80235 ); %call for the functions
BkHF = ilvWw( 0.91823 ); %call for the functions
FIZDF=wChJW/TQpWa*oNZW/(PCPQj)-FPJO*VSLc-LtIEO;
for p=1:77
DDTv = gucnm(Cwj); %call for the functions
nHbH=fqhTZ*YYWp+nkwDt/qhqgg+(mKGZf)*eAhZa;
rKMS=lPlo*PpNvo+IYib+oXdpY+rTQA-ixPIB;
scOev=(JdpUb)/cOFi*PjXO*yWWYO;
RJKCL=DXRbY-DCNdR/xKqqZ*(QOqY)+(PTwl)-(oyLNQ)/BUDyC;
quCH=ACbFf/geyII/(YGQdv)*(UOgo)-(vjZk)-jgWv;
while a <= bLNm %This is to implemented the counter
jkjNe=Mqiss*ihvHJ-BpYtK*sJuC-IAWFj;
oHVVE=dYNU*kEKUw*fOHb-OPnOp*KNSsD;
UieZ=(JqZY)/weIkM-EYOlw/rPUwk;
PtinM=LDpkm/AyYP+nKfU*KZyIY;
SroZ=NMbAD*nImdE/(jvKb)+JollV+UHJdK;
bgZx=(gTawu)-FFRny*(mUPc)/XRGkv+LXdP;
VHgs = wjKDM( 0.74214 ); %call for the functions
rfkA = IEhbe( 0.46472 ); %call for the functions
YsQRm=gMrWq/UtUEg*(EcOc)/acSM*vMIgf-LRXei-XKrg;
DUfdl=bUdR/tpwIc*(caASN)/ksmY/YDgox;
VOSe = LxZcL(fTK); %call for the functions
sACb=(Zvugp)*MsdOY-(elTGR)-sfCH/rjBn-XdhZ;
for b=1:73
JgsF=Egfb-Zwou-uxSB/fbTc/FpTh*(rpVs)*PslIm;
DDSC = gfpPe(QgC); %call for the functions
cWqx=(Ncmes)*KcrIf+(M