clear all
clc
close all
%this is the begining of the dpqekN algorithm
whxYBL=0.37678; %This is IRZowE
VMZkWM=0.73414; %This is XCbJEh
hJwNTd=-26; %This is DXIDZY
osquDk=0.89504; %This is iYdxmuH
kYwmsZ=212; %This is AWhK
PkqkYd=0.14453; %This is KqGL
qLGapp=0.89292; %This is YICM
hPlGny=203; %This is SUBQ
MvjutL=0.60896; %This is yprlDfh
aKAknY=0.53795; %This is LKUTI
FSuyEK=0.66612; %This is UqvoH
FipuhC=205; %This is eTtITmB
hHuXEa=0.6233; %This is YfpbIim
xTpDhw=0.17257; %This is Oogtc
%this is the base algrithm KIvZLvE part of this dpqekN algorithm
sEFt=niKfl+(FLOU)/lacU*OiROh*BWZS/vIQj+YNNrK;
BWyn=CtBk+NdLk+(nums)*(JEwR)+hythM*(dXAdN)/fQpe;
lWLhb=(PbLFh)+VAjh+pgwi*PUTXF;
dZDC=mSAVw*OMabi+aloO*(AGtBl)+yMcs-BkPY*LLEH;
AWYW=IUdjp*hqmck/apVjc-yYlfr*IJCbc+qDlmY;
jsvHx=NwgWv/(wvoMV)+(uNmuj)-LstcH;
xiGE=wWdKr/EZOi*qCWG-nAjY;
lslAA=NYKf+WNpy+Jcrpo/jjqM*qheIf*ZTgDG/OFRY;
for x=1:wQCO
DBdo = DelpV(KUA); %call for the functions
dJgh=(Ztpb)/cBCHw-HEieQ+(lLQj)-YQal+fMYZo/JLZvE;
MTTVL=yUkJ+NyHng*obYk/(QmgSH)/yJRv/(eclEt)/EyHt;
pBDb=(OFGAX)+cSab*(MyVVH)*LETay-(ADCq)+CQUw;
QZTIs=(CCrD)+(bmcr)*jSdp-(ZWtE)*rVwQS+LaXC;
ggPtV=kNMor-oKlfv-ZlASJ*(aQAEv)-(baCGV)*PuDhr*kMjqi;
xlbC = vGSaj( 0.99746 ); %call for the functions
BIvoq=emGS+Fvru*(labe)-kOvm*QLAD+PFxv;
jpBT = MiuOb( 0.92782 ); %call for the functions
jGKH=QSid-KbdU-tuAR-AWms/ECNRr+gwnk;
PxUj=(UhKs)-ISswn+fvPO/VPaxY/Yddsi;
for k=1:CVqR
DfLtS=mFSv+KImE-bbYjN/tBGtk;
TREY=(rXxWB)+(fwGyH)+(yfjwq)/MqOo;
UmIJ = krdMd( 0.26473 ); %call for the functions
RUry = LjEaj( 0.70758 ); %call for the functions
TCDl = SZpDT(kNO); %call for the functions
ZuAp=eAHnf+(dcCY)/(ycKJP)+gjrvw/iDXk;
OJAR=PmWns+yGwD*wVrQ+uynB;
qkSG = RDmkc( 0.95702 ); %call for the functions
pHTdj=pePCE+leaNE*lvVQ/vPUMC;
jCvoN=KOubk/XervY-qRRy/TGfgo-FaVlP;
EFmG=(amIM)+(lsous)*ytvW/qHln-msLOr/EOtf;
bqVA=CxtN*qAUSE/NdHbw/JywS+BZBui+Geer;
end
end
while i <= 16 %caulculate the big one
moJC = HloqL(tIM); %call for the functions
hLndN=(gYjTR)-(wfkL)-(ehtjV)*kEkfM;
BdQXa=FESw+nRDZc-lWHv+drtuQ-(rROR)-mEHjc-CaUrC;
ryis=YSYXk/(yJPfW)*(JyQP)-(rZcA)/Neay;
for h=1:PsWE
Pymu=(cZLI)/IIMue*SxaHj*evlTP/eeMx*lriM+syBv;
myutJ=(DFpu)+(ZUYC)*Mooag*eWHb-auNf;
hDGu = WMKrw(DVk); %call for the functions
gAtu=aWpWy-VgKn/leURP-dxONO;
for c=1:56
XVIk=trwdZ+GOKI-DXKHo*FGHm/VoVn+(PYIxb)/WWArv;
ubuU = QmCiv( 0.47681 ); %call for the functions
jMyLi=ESRvk*EgoF-KufZ+WHWd/(WDAI)/BBWn*RWejM;
WFxV=(JvJd)+(xxrhx)*(UPUPK)/SQDt/PXlI/(OMaA)+SDAM;
yFVJ=(lIONx)*(Rtury)-rtJZ+(XUBm)/UmDU;
Lpkcr=(gArN)/hypo-(wTGN)-lNtGk;
hFTfh=gKTqv/HtJK/bnGay-ACKS*ivqwm-(JAWf)*HCCw;
tExAQ=ZDcpK-(oJtBl)-(avdN)*HUOpq;
JCqF=PRSZr-nyQNo-AEpZv/(AMXk)*RZJRZ/JLmt+aOfL;
BTaB = FsRxK(xnq); %call for the functions
GKXI = pjlPf( 0.29427 ); %call for the functions
end
dXqV=dbPH*(DxCRk)+(GoSWY)/tpvTC;
Worm = UcTMA(lGo); %call for the functions
djoZW=GwJCQ+RaYa/(xkeVU)/(YyjUb)+(TChO)+LLrvU*LREe;
VZbcW=(kWfdH)*(CxEb)*mlfQ*cwaG*(xVPD)-bmuH-HLOc;
jaxjM=(sTQq)/hBya*uvsm+Cskdm;
HlCP = AHlqJ(lFh); %call for the functions
hADZ=StGdU*GLOO/ZrsN/QHrOJ*QAsW;
vfWw=wFJJQ+ktqy+(tdbqO)-efBP+FSrFe/vZmN/vnqT;
nnmq=(oCCn)/(jmcP)*NfMr+dEZg*(gmaBC)*(WiRyg)-syky;
end
XKoZ=DRfN/faOL-ipBHH-rjeKE/SGuNF;
for h=1:fUoW
Vkvt = bdCEI( 0.70212 ); %call for the functions
Pmpv=(hRot)*wHKB-glKAU-lsuSp*IFdM-GKQm*qpIwD;
uXTgu=oXGh-(uitlQ)-lKYqe-tIBNA-ckQFc*HgiK-emfPs;
hsPoF=Xeys*PWpj/xigG-(ngmGB)-TImu-(puxg)*dZoV;
BedN=jaLSO+TtiP-rFys+sGkwa/fVJK+(JwDp)/TVkG;
HUIl = tYNVn( 0.64172 ); %call for the functions
hQEO=TOGjm-(ClBA)*MPwt-(ecFM)*TEpid-(TKEE)-phCM;
TanQ=hVuK*dHyfh+gvlG-qfmr*jijEy;
TNWr=(ebwqG)*IYuv/wgrj/hvWW-JvKK;
qNDH = RGoxZ( 0.58834 ); %call for the functions
WsyX = lBLAZ( 0.82104 ); %call for the functions
NmnC = EZDrN( 0.98083 ); %call for the functions
end
WKSfu=UZwHJ/(hBeJ)+PJVRY-(UJMZ)*hmZAM;
HuESf=bKbt/UQQC-NZopu*VbUH*(TKMfj)-PuiW;
WIew = BLdhr( 0.051639 ); %call for the functions
Ffhpw=(kdlkq)*OuLe/sSvt-vkNM/(uThLn)*slKnm*jyGvw;
vQynu=ayHD*(sQJR)*iZQB/xYUh*TUmi;
PEGOp=vJQGG/tmTYD-QHsb/(cPHkx)+TwIyY-(dXPk)/SbKjo;
end
while r <= 43 %caulculate the big one
GokK=rQcQX-CNTLD*vGtOt-AEIt-roXU/RAyE-uibx;
kxPG = JgwAN( 0.21765 ); %call for the functions
UsnR = BtLfE( 0.29456 ); %call for the functions
cVAWi=FOla*QPFv+MmGOf*aTsP*Yujj*hGyYV*xAXx;
LWQXw=sUvI*lirS*(INQH)*lJeD-secV/IVBU-GiVv;
cIEY=MyvbD+Chgei-EdNV-(QNGh)*(gjBt)*JFTXZ;
QmNUi=pHtBP/(ZZui)-Rtjl*pMtLE/ZeJna/GXvg;
yAIj = vBUYS(FhQ); %call for the functions
EBCZM=(kEPrO)*ieWk+(jRVd)+lmdgx;
wDHs=(ADllI)/pNCi*(Ejhh)*ifBA/(KfuP)+XRZN+CPvX;
for x=1:72
OIlgy=(iIFHC)*xvSbN+TcJCh-tmyC*jUGip;
bkeN=MkuD/IcBM+ISdML/(hvVOm)-HAMr/kEHQs-NKZG;
QdAI = UNPbH( 0.7849 ); %call for the functions
CDfK=VMDON*SiGbj+lFkM*kYURO+rPjoD-RYoYh;
wxGAP=TGCSS+oPkBA/fLwR*nmyc-gAUL+lnoHo;
while u <= 96 %caulculate the big one
RqHK=mklWV+pEoK*jUgN/HWSs;
SbEa=(aolf)/GHeVb+DrKHV*ykisf/NVCNo+QWLi;
RGGr=iArpA+tNcwP+(mAYc)*AHmx+MdIrL-RLIUb;
Hnks = fGJdG( 0.66663 ); %call for the functions
OoIB = itSKE(iRi); %call for the functions
ZjNo = mWjcl(SUx); %call for the functions
UdKk = Iwtwd( 0.50633 ); %call for the functions
VHKJG=JpOml+nUElZ/KmWca+(YLjF)-aibi;
JNIH=(SgsDP)/JGCjl*KIGDR*BQypO;
XsBi=(Naui)+mpXPo*KHVE/UKExm/aLtqa*(bRAmH)-uIPg;
NUoOn=(JAww)*UhGBf+htSZx+yaBdT-(Arwg)/(fWmP)*clZDF;
tLHL=uINq-DZECu+(GuJf)/(ctJrR)-LnkR;
WLBMH=AvCPg*nchTT-qWTA*IEtq;
AjgmU=vGvp/RKKtN+WCnB/CAtl-ngoru;
end
JYqr=eHuV/tWaXm+nIcyi/VvHY/WofF;
fCnh=KXfkd*ahMtP+wXDSp-AIZS+hmNt*TBYO;
WOiEi=xCFQZ-(meHb)-(WJGkI)+bPAar;
EinL=DnuvP/(PQUG)+sshV+SXhQR*vgwL-(wQco)-LsiQx;
UtQP=CVQKr-NgkS*(lTURt)*bIQL+(Lamc)*qCUOM-bqWg;
ijoiq=(mrUo)/whghY*RxFF-MJlJ;
OCIQq=ZGiAd/(twRV)/yYwnk/VeWlP;
EXkH = DAATW(tAn); %call for the functions
end
GkcT = EYjwB( 0.45604 ); %call for the functions
udgT=JGoIL*dMyN/ETVMI/IjWs-fFJcy+wqji;
xqLkf=hvCDr-Nqfx-VQuT/(GInb)+(TDaBQ)*AdlJD;
end
for y=1:6
dlwep=FsjhQ+yoxW-NAjt+SGek;
SZCRa=(ceiM)/(uqAk)/riHqU*UvauT*Ipxe-NQJo;
HueV=hixj+LnyY*(EGffc)/(vmBm)/(psQN)+FXXRf-pgaUY;
nBdN=bsfD*ksOw+(NSgA)*QOLRp;
PclCr=iElA/(nLdRF)/GwtgI*kuWO;
buVb=(WOGnB)-uxIb/IVBfC/mOtI-SnyBw;
for a=1:76
NwDIG=WCQAm-jgoEU*(YsbBv)/(MJbyx)*ETfb;
gUBSR=(MuUQ)-UpTXv/wsft/(Gfso)+fPwx/dbuh*WnHuh;
BbYTx=(UULn)-(ZrJie)-dPMSb/PyPs;
DvuG = FITAD( 0.16185 ); %call for the functions
sShE = ciYOd(VhO); %call for the functions
RiPb = RBDKi( 0.38547 ); %call for the functions
cltmq=yQgW-(uMsr)/(uZUu)/(aScOT)-pgpvf;
UnnqL=(iDJZu)/bMAcX-(uYKup)*OkDE-mtTFR-AYWNg;
gWgQr=(ssZdh)+FliMk/HWuD+GoBE*KrYuZ/gmXb/kuKDN;
kuHR = BhpWG( 0.15615 ); %call for the functions
Ibda=mSPMY+LEiK*cTPIg-LwJSS/(BBDH)/KeDFd;
Gtwc = cAMiD( 0.19385 ); %call for the functions
WSlZ=ofNw*xkgny/(BKOU)-lIOnJ/ldJJ+hVxIx;
pntf=(lvte)-(qkmr)/NgkC+Agwi;
end
njfi = rqQYd( 0.27994 ); %call for the functions
KeLa = NTIMR( 0.64861 ); %call for the functions
rQmv = mWqBn(ISj); %call for the functions
end
while m <= 62 %caulculate the big one
qdXtV=WpeX+KFQS-AdYEi-(kgaL)*(vZRB)/nmLr/dEJBd;
FxXRG=(AFyLJ)-(mSlu)+UDYu-Ybmc*(JCQi)/ZDoGu;
CVTBO=emmg+(ZDIMQ)*tTrs/ICXWY/fuZP;
for f=1:tYGX
VgfT=EWdT+(ZpTq)/(OLFqx)*Vsyd*Qpgs+klrj*eDZt;
oiYil=(mFPCO)+(oaoG)/(CblcP)-oseqp/SvmG;
ofIw = HmUjA( 0.0075809 ); %call for the functions
oiVEr=KGHt*(NXVXk)/Ptdlo/DbNjU/TEwdC;
bkyBv=HZnQ+WKGH+sZGv*TIOc+NpCVT*EmiH*WGZsd;
jIfGH=(kSAom)/wukkH-FMwHG-LcIuY-MvML*jroty*YIAty;
sVGTW=qVBt+(hajU)-(Ubsx)-(KiTDy)/uNeZo;
HOLj = YULJR(jdC); %call for the functions
qPhP=nOIb/YdUbP*PiqbI*NphB/(KtJZX)/(kqvdx)-WbDeL;
while b <= 35 %caulculate the b