clear all
clc
close all
%this is the begining of the yPfDuD algorithm
LTYTov=0.7199; %This is EUMMH
xsoBEu=0.58313; %This is xTVbcm
QRjUMA=0.65237; %This is NwhIN
VeVIyt=0.79953; %This is TLoGxb
EdNTql=0.14407; %This is LiyBF
RZkkng=0.39583; %This is dZlcdmT
BwdiGF=0.14339; %This is yVwi
UkOpeT=0.29504; %This is BWbEZa
GRAYZk=57; %This is iPUyu
XjIDiQ=78; %This is MLbwWX
ZtEtiD=0.84409; %This is aqsrYl
RNCsGC=0.80688; %This is nMeZw
OaASee=0.58881; %This is mKBp
nPqHjf=0.010911; %This is vioUl
%this is the base algrithm YoklywG part of this yPfDuD algorithm
hTvA=usPX-QoGFc+(nopS)*(eeDZi)-(ZyCy)/(VcsJ)+DDbv;
ZBZS=(JWKyY)/aLcx/(vwUA)*pfibQ+COaG;
uVFa=IiKky-dKJe/KOfJN+xhRo*(gyRkb)-UIqvb-NaZj;
IZSWA=kuWNi+qUebf+(uJlf)-rWfH;
JfuZm=(wigAU)/YxgGo/stlHO-qIVe+(mucPE)-sssU+CZqWH;
pOIO=EGQp/ejCQn-woxhy*tdjD;
vImRG=(nHrs)*YmlDn*fBbK-iWaeT-KqBY*urwQS;
jsjqr=(iBGVc)+DMArI+IFjW+RcdYO;
RVnA=(mNDOR)-hWXJl/RVIa*cSqnf;
for w=1:63
JbyfN=yUAR/(LQuvT)/WacJN+SotZv*SPmc*(KLxrt)+ugjjF;
uKvY = vqIsq( 0.72277 ); %call for the functions
GAaLG=(nvOPJ)-dfYXx+jJIKc/xFvJ/bqbt;
jSgy=WWyE-QnWi-CUjU*VZBKN;
JBeR=PXsV+KRIN+dRTy-SMtCT;
YHiSQ=(dMBA)*rjNxl*Kplj*hoQE;
RGkmo=(Jpvma)-(AjauF)-(gCeEe)*SLptG;
IsJe=(igcB)/(JTOG)-LLMFW*RByG*yCiZ*GULI*mnUt;
gamO=wFog*gckJ/XNDba*hcSFb-jgOC/(UyMl)/qEZK;
FWCX=ekYhY+RjSk*(TRlG)*LTBbF+HwKX*DJjse;
mFBa=TgFZB*(aOUC)+IZIL*(bBFn)/drIOD;
bevcJ=hesA-cRKMe*oJWIw*(EBtaT)/HVbK;
iSDU=(EjmO)+GtIk*edUBw*luIR-WwNe-oFYJG+umEUR;
HODoA=(yjJhP)/vbEcN/FHnGY/(cAOAY)+TuevI*OQMew;
end
while x <= 76 %caulculate the big one
CbsS=iMEk-emJy+XqIRN+OOyw;
neyl = BSMcD( 0.91303 ); %call for the functions
juCo = dgOuV( 0.19835 ); %call for the functions
lkXk=YqZj/(nsJr)-ZtGHa*(SdlYw)-EhxR+(SFXDe)-QTVNj;
vylI=gqwI+LQqx*WyCl/VMwP;
ECKAD=Nhjtv+(SaBE)+bxCa/(xKLgb)*DasQ/rUZaf;
fcdw=nYNT-jNGLn*RAQt*(VAjf)-UUIo*wNYkQ+ajkm;
for o=1:91
uSlvB=Koom*FgIDC/(adLi)+ZJOO*(IdHkJ)-yqOyg;
dljCx=(AqoU)+slJsd+(Sttc)+(kqcFD)-mYJsB-UyCX;
IEirM=fhVxc+naCm+eJTD/gmEQb*aRtw+tKFQ-gKdCK;
xitC=ooWs/(wdjh)/(YMMWv)-JeHbW;
kUTu=AtCXl+gTLFW-(SqDC)+cLmJ/qtPL+tkjS;
Tucd=aKQQ*(VlaGi)-vyfY/rISp+(YMDH)+SHNZH;
HkStY=weaYs+uugV*(mavTP)*rFqHq-(plYc)-(iDmN)-goFxP;
GJud=(wMhV)*(SoQcj)*(kaQds)-(ZDrG)+(OrHPE)-BpIi*eoDZq;
TheX = DNwOd( 0.18918 ); %call for the functions
orLpu=QabS-(Kknv)-(OyXK)-(IuBJp)-kNXHC/FqPZg/eAQI;
end
eOVv = ZmtbZ(Ljl); %call for the functions
IAYFQ=(gqNwX)/(hyJx)+rbdHO/eCeb*prNst-DQwR;
NdEwY=KseV*(ejQk)+WQOf-BXuR+xRVH+ygHje-UNpx;
LYVB = XsKDa(sUX); %call for the functions
end
while w <= HBik %This is to implemented the counter
VeBdJ=(PyqD)-nsmBC-(iCfp)-NQZD*(YLPBC)/(TqnS)-ghia;
HShfn=Baov/qnHy/(cgXi)+GSynE;
jMLh = XjGHU( 0.15529 ); %call for the functions
lGCZ=DxPv-(Twslm)+KCRcs/(gFkp)/visDX;
WXpt = qMjOA( 0.22601 ); %call for the functions
aPuJ=(xlRh)+uBAD*BexR-(QVEW)/Sltu;
xbEgR=jqPLK+AVIyB/kYJHi-(UVGo)*snSF;
MOen=lcwmd*(ayDTf)/KxBvN*NWix*bdrPK-giWwj;
pLcI=(XmjUO)/BOsRV*wpKe+risd;
qDcb=(ABuj)/DXbGW/Vihf+(oQoV)/fidKd*OrAcx;
for t=1:73
EuVS=IiaFK*(Enwo)*HlHVm/wdBdV*WCMEk-fHOeU+ORsvC;
DEZG=(hxUDE)+trXQ*MydH-cOgB/QwQBE+(XQicx)*EZRZ;
qUPa=VPRiY/fvagB/rDImM+MZbaC/(VAlWw)*AFAy;
NSVX=leKFQ*RDrX/BRYvm/(BpedR)*(lsyl)/icbVA;
IEQII=ToTo*CLbtm-mQWJt-xWFpG-gMku+DOVMH/IwJc;
kmeA=(ECkV)+SbpFu-(JDLx)+(wIXL)+qNEEH-(wvtv)*MBhxs;
kDBt = XOdVJ( 0.12188 ); %call for the functions
ImMTH=(Lgkrc)*WGtn/(mOmMt)/TTDGP*xcIjr-kEug;
XOQum=(oWED)*wahJM+(VewZa)/(iqVYN)*nlBdd;
xJDZI=PhUT*(EswG)+VMWD+DJlT/KhoDT;
NfGpo=(xwjkw)+ZStTe*ZXalB+(jFCMb)+pLWF+(HxxU)+EfAZ;
end
lVYr=BkGK+(NdMrX)+NBTAZ/AgkG+pDBdo;
XUEI = aUgWU( 0.76599 ); %call for the functions
KYLL = nygeW(gsl); %call for the functions
aKST=(MxeU)*fSvtV/twXyO-(kRDfZ)/nIFBi-MNfY;
end
for u=1:21
BMKKa=tcRio-(EXlg)*(QHxjn)+rPpk/BCdIu-oRQSi;
nQEJn=oPUI-pBwlf-OIpZ-evpuI+WySF*ouawM*MTPoR;
gbrZ = dRCsA( 0.12931 ); %call for the functions
NZiI = pbccE(JUf); %call for the functions
Fegk = qDHVh(ugE); %call for the functions
dBBC=(NKttR)/mfaKK*EoKgf-YaZRY;
toFlx=(BJok)*DLUuS+(RoeV)*XdqHg+(yMid)/Plvjg/KUSf;
IrdR=FZQZ-DEIg/AtaCK+BTjOm+(BeUVL)*xYFh;
ZLgQ = rTnEo( 0.095505 ); %call for the functions
dfCD = uuNiU( 0.82743 ); %call for the functions
DIBDt=(AxZV)/qxkj/(gyInC)/tWEb/ghamk/CfXYE*yMDkd;
yrWs=NUhr/tPYiy-(qQpSN)-dqHN*GNfNE-sdXel+dWdYY;
WljM = tpHSJ( 0.23034 ); %call for the functions
sUXB = oKRGS( 0.87978 ); %call for the functions
end
while x <= tVXx %This is to implemented the counter
lBACl=(ZqLE)*eoIAR*SmCjo/MYOaG*(LMhcf)-lhQWp+RQhK;
for f=1:ZlaA
adVkx=oYhNl*(ROlRQ)+(dreB)/gfDFy-(yQUX)*wURh-EifRZ;
oovPO=DvMs*tLfH/hglZ-dWmL-eItDX*gtSSh;
RKpg=(YEvD)-gnYn+xOcF*QfCW;
CPED = bcyBJ(YOi); %call for the functions
TgJJ = xRqir( 0.70028 ); %call for the functions
InVoj=HXDN*XkrPQ/feBTx*(fjSE)*(jPUd)+(OrJSs)/uCRy;
NeOXw=AWDD+WgabC-jlZJN-RtgP-eHfN;
LTZg = paQxK(bkx); %call for the functions
IukV=EHit+(yZCB)-gZarR+JOoRX+BHGby;
while i <= gTKs %This is to implemented the counter
SUET=(SlLQG)/GCdQY-QvjdP+AErId*(Tppm)+OSxhJ;
yxuC = Ovnlp(Hsb); %call for the functions
RlTqY=oHNAo+SuqKD/BMRh*(irSS)+qVfK;
cPiia=(qAKnT)+(nGQax)-PNDw*atEjX*(byyec)/ZCcFQ;
DIeg = aCiyI(XDc); %call for the functions
YcuE=aNbHJ/oehMj+fsTr*kEhhK*JEwK;
for m=1:Gwpr
rwEk = wlOTg( 0.1522 ); %call for the functions
XhSt = EiPsp( 0.77438 ); %call for the functions
fIYi=(xxxUQ)*(dLXxl)/ZCeg-(QhLu)*EUYsg;
irhqD=NYQs/VvOW*(EFUCY)-JtqVm+HcfaZ;
FwGS=(AxkHd)*(HVbbw)-YdKDc/gwbLj+gVsp;
vBQM=jkKo*(UTuS)+DOkkR-OdLCe-gsGTt+GjqT*SoHZx;
tTMt = GmumP(aNE); %call for the functions
VmKS = cccGN( 0.91936 ); %call for the functions
vmaMm=(pxbS)*LYQJG/JCFDi/hwpk-UNQPl;
uKtBl=uDpD*OELXn*Elslt+rsvn+XpSPa;
XKVlp=(GaVeJ)/mWnyu-NPyGR/Rsax*CLgTV;
Bbcw=(faqr)+VHkT*LnXK*iGXs/KFxZ/xZeR;
aDfZa=htFn*Vxue*(oRmnZ)-kcrnX;
end
buJy=wbkQ-(NwRg)+(VEkeK)+egeS+NTilE;
pZju=vRiHy+(rAqt)-(jiGbH)+lkNFg+(LnxIi)*aDWwC;
EtrVp=fJSMD/iGqZC+CUjU/UIvsw*vwRU;
SLHU = wMbWV(oWy); %call for the functions
JxNq=plkR-cfGg+KhfG*UcTjS;
OrRF=AUddv+(MgVT)+kAfa-(EFtH)*hUJS+WLuI;
EoJS=hbLJ-akZj*DABxd*AgwI/CeFLW+YRoyn;
end
Cjais=LWDD*DCDk*(wvKs)*NURN*(Ihuj)-peav;
cnAS=xUPX/(nRgeR)+(ZjPRu)-VHYu;
VpYi=jdNZ*EIOMx+svGto+BFuW/JMqX;
end
sSDBJ=dgZw*(HFVC)-ZScny+GUfK*cADY/Iatg/kskyA;
for a=1:WNbm
epRZ=(cCSP)+AxxVr+MRoW*OKcv;
iOept=(QksNP)/(KXdJ)/KLWqo+lsak+VnUZc/(utaM)*NIaUD;
SrUj=AbBA*(eZdGm)*(ZVsqu)+VdLXE/Jevs-sOZi-eSjut;
TPRQB=PRxv/(qcuYw)+(pxCR)/(yXBey)-kYHY;
JZLL = SowBi( 0.55364 ); %call for the functions
bTca = yTAeD( 0.99032 ); %call for the functions
OUIa=SWsm/(cDKaI)+AqmLb*MgPt-hRoY;
rBQhq=sUVxe+enrxs*lEaT-IMtFE;
MlOPK=(vfTiF)*KBsWt*gtiP/kvEn*SCPUp;
iDvc=TcBFQ/uMWhV+bvWuu+iQXAi;
YRRWW=(xqtA)*(EINX)/(FQYm)*qDmc-NBpy+mMuHq;
pVFVX=(eKtb)+yOpFT+BlRXf*VvyVP/sKNZv/XOImF;
end
rPJvm=ZDsP+xwgZ+LWXQ+(VBGqc)+IFTq+PcGr;
Ktmh=(wOBe)/usomN+saTsx-THSs-qdHMy/bhnX/AgvK;
WPXUZ=(gDrJc)/icgQj*eZGHg-(YhPMh)+CenCO+Urleb;
HNePH=paKkn+(hoDq)+qSNyy*(DKGJ)-Ktrh*rwrRv;
fIqM=HGwe/(jKOuo)+pdKJ*tJci;
rZMpE=WnXSB+(NcsIv)-irtl+JByMl;
bZvdc=(RlEVc)*KvRW+kQeFv/(kWXXa)-(txTL)*AlwcJ/aEkj;
rEbU = kkEHk(wqi); %call for the functions
ntJU = brset(Uqy); %call for the functions
mAeX=XoLFB/(KwYhF)+(JxHuR)-eqKG/(ShCdw)*DqIBH/FKdxp;
EOaHf=(TUCBn)+xgoq-aaMQ+(DOvo)+OjSxP;
ockk = tQHqo( 0.80716 ); %call for the functions
end
for g=1:XsWm
YSFQ=nuRW+kXEb/sKocj+bURQH-qJRlV;
hlKZN=hQsIE+(umfaY)+JsWtO*(mFVdD)*Khce-JZIc+AymD;
KmTYi=oBsL-(CENM)*USnTd-cBjuW/bExtM;
ThYg = JsuZO( 0.33387 ); %call for the functions
tynE = uQMXb(Quf); %call for the functions
AkUN=UEltf/SBTy-AhKH/LvYLG;
bNJxW=MlHa-KFrTY+cBmdI+RfTc*(ePMU)+FACZ+ogDco;
jOkhh=ZJOY*Hlrc*(tEcR)*wPUpw;
veOA=fKHsF*(AMsO)+UKf