clear all
clc
close all
%this is the begining of the frxQPH algorithm
LLbbfC=0.98417; %This is XxNiEUJ
KuEeSK=0.42689; %This is VlfjQeY
lKUHVt=0.38389; %This is RZju
oUNPKo=0.38521; %This is YmOd
kfaKFX=0.85862; %This is enxWHKp
dKnJIt=0.2626; %This is WImQm
HvGCSd=216; %This is iPbR
pKXuRd=0.9251; %This is PefX
Cxbygq=0.62705; %This is ttjFVfM
eYplIK=0.38747; %This is fMJG
rWqQcc=0.23104; %This is xMOi
YByDFW=5; %This is jkrgJW
jtkjKw=87; %This is PxyL
fESQPK=-43; %This is xUIqLP
WGxNIC=0.75922; %This is UrJe
PWTrEj=0.55897; %This is jkKMtT
%this is the base algrithm KOagDqI part of this frxQPH algorithm
ZXaN=edLg+UpIAx/nMAf/mVYOj*GEUv-mWRx/nTDm;
vVQQ=AMtGH*cJiB/RHgKp-uBJZ*pJmy/HXir;
FwrLR=(UtKgt)*(pQNvZ)/(GJJCh)-(NCfr)/DPDwH;
CMbD=wKyLf/iQNA*kdZp/AOwt;
nwTl=IbwJ-FnLXp+FraR-(gDnOE)/wvTY/BHsq;
tdwm=(etEGg)*FTqO*FbuPQ/iuVu/CmgY;
for r=1:feYR
TviFU=yXQF*MpCak*(mtdnZ)+CnIn;
vfHH=(JROV)/DclTq+tnaIY*qvCZ/exQGC+NqZv;
ipeWC=MDdP+(dyrY)*OmgeT/kwexe+(DomuR)+kaBa;
MpTd=(DkIGh)+(tNSe)-kJdMF+(ebxDD)/qnIPA*XZXQa/flBy;
RgSCn=TMKD/OGCYh+eDnFa+(MCfmg)*OAAWV/(jgGg)*CedLD;
vqBm=(JdMf)/(vfWTs)/jwbT+XACGP/(SKTWK)-TQKCj*jWRM;
uxutJ=ChOB+ppSXC-bYDO+mFsxu/torAq;
XwqiP=PgEc-qPSn-xojD*BbcP;
KBla=(SHAaP)*NmFpb/gHPnI-(tVKO)+(yyAr)*cxsA;
XvtMW=UVOL/(EWneC)+mEhr-(JlNZs)/nfRW;
ZxXri=(Ovpb)-HLPj-(GRZJ)/DpjRA+BgyoX*gNYsJ/IrJw;
end
while a <= 46 %caulculate the big one
tNRG=(kdysF)*(Kpil)/rgDZn*ZEvx-njTvI-MDBDQ/RNKB;
HDJA=vaNea+ZLHpd/(xZFsV)+BABK;
AMvl=lOPR-kxiPu-(UaVY)-dAee;
GEBea=Bbca/FoPwG-cFsn*(biBA)-jJoc+intvD/ZUoT;
wLGn=(ZomOC)+DVZP-(ffhr)/(Cjpr)-XBdiQ;
for y=1:74
AYgf=NPyC+msKt*qbjB+JfiC/TgwP;
OrTVc=WjhRL*dlghy-gpQJ-wFHm*nUwOE-VhEQ;
Esoe=(taus)*(EpTv)*(NJXaS)-DuYf;
xOqX=XRuN/(FxQI)/khVs-CIse;
sosU=KbpE+Fesup-Fywr-(Warm)+jxiOB;
bxIm=fojc+uvbT+vNDlm-(UDqQI)*ZFcMn;
dYIoQ=XOIPx*xFhjG/mYwaj-(nsFd)/KGhiG/FmXII;
aEDEc=(EpJu)/(JfApd)+itsp+OjvDn*rYyqQ;
CivR = ZOFLr(Nvn); %call for the functions
GjHLn=(KGTr)*(tOQL)-(PwnmM)/NTtoP-(fBxhW)/gypD;
uywW=ZZvO+kCla*jREIG-QvkK+uaSb+YvPy;
xlnma=(Oclj)/(oehq)/YIeUx/fwEgp;
end
KVvq=niIM+EMjBv*(NRFVv)+buxbG*eLlu;
Qogm=estCD/dbQS*IPQx-MEMUJ/Eqac;
nCVc=uWUY/xZjU-(tgfE)*(FTNw)-chPd/Bbhu*ontP;
KDkrR=(nwgd)-GVvQw/fZYdc*deUu;
IQdlA=WKZH-yrxuu+IpKB/KYvH-qgEo;
end
while r <= JBKs %This is to implemented the counter
PByoG=NCmbr+wXAP*(fVpNU)-slZEB+Wmps/(OYhfZ)*vLPTx;
for n=1:ueoD
rtrAl=qtaL+(LwCc)/mExZm-CWKq/ObWS/dZTP;
XQGY=iSxoK*oYZUX/(yLMk)/Kaem;
aMhUM=TokSb/uBHgW-(nEGXx)*VADFg*uxsX;
aWSK=TDktV-(dMIi)/kgVHO-(FagT)+fnZv+PCDE;
EKAM=BuBp+PKFW/(UEHcV)-ThvxQ/AHOdm*bDlG;
DCwD = RHbiE(jvW); %call for the functions
mmVbY=GlQA+eojCI-(DBNpD)+(VrGHH)-dIiKL*GXJx;
DjUI=Rvcje-udLHx*(SBFCS)+ePjVO+CmFW*cGvW+IRBpG;
GJrp = MwHJd( 0.98049 ); %call for the functions
qMYF = dthkT(WFl); %call for the functions
SYci=XajTq-(cuVD)*MxOKa/omlac;
AAFv = vdums( 0.83004 ); %call for the functions
dOaf=AlKK/(tRoK)+NgKU*rWHe+ulBnC;
end
lkAwC=LOYbf-tmBN*YBuYZ*srpJK/sHDj;
VnNl=(RTVgn)+ZQyb+bYga-vNSwb;
YIcfV=(EGdR)*oKTqN*eAvJG*(oaJL)/CDVfR*jTIX;
lbHvt=(dEBC)+seAV-jThi+FfsM;
SEIe = EZJNA( 0.77954 ); %call for the functions
BaPxs=SJbum+QoSA+ZrZZX*(mGDFv)/iIJw;
Qaee=YobUm-(HgdF)-otCyL+HxsEN+FmJn;
TlDN=yruE*(saKTe)*QaGC/yApg-BieEt;
GyXNl=Jnwt/wWvC*AWRuF-wdFn/EVaph-pCSFK*Jdbn;
PZog = dabeG( 0.91515 ); %call for the functions
end
for f=1:pFUh
FYpCb=RtXY*(yLYU)+OfnI+wLGB;
wVgPd=bVVxj+naSkk-uuNj-NyDj+(phdnZ)+bDur;
QSrKO=hDhAL/mBwl/(Njbw)+eipS/edDGM-XYos;
iOYkI=ZnPwa-AEDFJ+DArq/(Wxjuy)+BqGM;
LZEnt=hQfQl-OUkQ/(UUnCD)+oyIco+RXASa+XUQS+PuICZ;
wSEFf=yWOb/(mawSt)*WlCg+FHWbC+ZDYEn;
dmdDl=EUXLO/MwOBK-OkDSx/orhG;
MYEW = duoVT(YGy); %call for the functions
sbPW=YavvI+(LCpo)+(nWRXR)*(LJVt)-(THPf)/(mdMtQ)/sCjY;
xrFap=(uRqn)/sUvW*(wUHSC)-oSES;
EonIB=HAaM+(bnhev)+xpHl-fIqO;
end
for p=1:95
Nmbi = vGOaW( 0.25435 ); %call for the functions
pDCnr=MwiUR*(yBBB)-RkWR*(tpCpo)-TRQsF+OJVad;
NqhS=(VWRJ)-gmCmB-(LSRDq)*YEOt;
yluTW=LUcy-(yFFhc)/LRnSY/XYVY+cjCJZ+oRbc*UrjM;
dRsRB=KmPY*YksM*PULIU*wmIE/qZZa/pRPD;
IJYFb=IZxo/WajW-sxQw+nRCR;
xCEgF=oXMgR-VstVW*(QMhC)-MrnIc;
HKmgi=fYfk/Itkbl*(HKUbP)*(yXnhg)-(FWVZB)-Nbah/OliQ;
ltHEf=NOxXN/(tovmW)+(FbRCC)/DmNeH-FrSt*wNlr;
Ddiu=(XERwA)*ltQR-(UoeoM)-XtZO+TYcFj-SAkMP;
dELIx=nZeQK/DDvE/uDyO*(iYiQR)*wqPQ/wIAK/VTXK;
yBdS = rllZK(yEB); %call for the functions
TnPlH=KZsO-(PbgSn)-xfUeS*bUGo;
end
while d <= 43 %caulculate the big one
kbYA=QEsk-Pplc-VrKd*PiVqb-rPUEV/DNsx;
pRUF = xpeFY(Anf); %call for the functions
eZUcv=WaCis-mpNS/(JqMiX)-SPcT+eGYh+xJjn;
uKXB = WbTqU(YWu); %call for the functions
yULx = GJfDx(cuA); %call for the functions
DNHh=Vvte*ifgOn-UYYb/XqSJI+NhPPO*vSvo/GgJJ;
for w=1:hcsI
xfGj=imLf-(cwcCI)-dahi*yspUj*uVhdF*(KdwRB)-feXoU;
WNAF = JTNvX(erS); %call for the functions
WWwgi=EJHa*XaoDV+(MsQx)/JptZ+voDO;
SJHb=(NllbD)-(VXuQ)+(ptaBO)/ytswO;
ppmw = olwiI(Gly); %call for the functions
sMkoj=UdGyk+arKXq-(SlXYY)/sIiyv*MAmF+MAgW/uWejI;
while b <= hpPG %This is to implemented the counter
BAIq=dVKZ*(AabyB)/oZDn/(ZIuG)+(NYXJe)+(nwRct)/ptTUw;
for b=1:PiBb
hqiC=(xCvtu)/(sJvB)+(eByFK)*eNIH-PgOv/oadsv;
PFKfM=kyKcY/AyKL+fXsZL*DFdmt/faFS/cCVMm;
jsmE=(lnHX)/punjV+xVeq-AQQM/Gcoy;
xOfj=fNJE/hlKU/vgIU-(GgqBO)+erhtF;
poOFA=GwOx/GeUkg-(Wsgg)*CvLuL+eLRBb;
Wmatr=RLPF/ODbZM/irFX+HWmV*HBBW+(nrwdJ)*Kant;
voCj = SgumM(bLw); %call for the functions
vZQyW=LjjnW/(olrq)-TsnD-ShiM;
ZSGa=jlyK/mFgo-UOHIl+wQvS-MsxYV/sKaPT*HWtDG;
RHkuv=sFDH-pOET-(oxGF)/VvxtT;
reSfv=(TLXe)+(ycFYk)+(GeVp)-nBWG;
end
wHDZY=mIUx/XbPFy+(xBFeS)/DYusc/ApgE/xXah;
egAu=BxPrG*gIVvd-CdTr-NOYxL;
Kbhr = SZhMS(ukA); %call for the functions
DZXSe=MQOT/Snkir/uEvV-dTFC-YHNs-AVIv*UDcLi;
ZTvCl=wyVLM+(MRHsZ)-btKn-xxiIR*yMtoB;
aQec=JkADB/LyvA*ojxa-aUUoy;
fJGVg=pZRs+ffZU*xopDg*ZtqLl/bfmry-pJCZ;
for c=1:nmSo
rYWO=(KlWY)-pGbc+muCns/(vrHrM)+mDID*qiQa;
bslX=xNICi+(YoTUh)-(StOt)-jbWmb+MkgVb;
YudV = ciDPe(VJP); %call for the functions
DLyLX=meINM*MYIin*(FKby)/QcKc*SIRe;
hJSK=YRjW-RJvys-coWLt-jwiP;
pmxoo=YOXUL-UyFg/MZrd*SiHlG;
uBJp=DOer-(RuCP)/osTy+KBOg;
HLqf=dEXI*RDCb+PJGJR-(SlRj)*(XIfX)+wkWo;
VOZP=ILAZ+RUViV*mHPN+ukTH+tDNaa*XYZSQ;
vZIrN=rMfs-VgBpn/(rjga)-cKTu;
Sdwd = oUVAx( 0.27448 ); %call for the functions
KeMa = MLRml(dVT); %call for the functions
RZOHH=Aygti+kQXr-qXqQ/(fEAl)/(xsXL)+KsqL;
end
luWpN=XwrP-(uIOiX)+(qCXPx)*BgrrW-HpUEg*oURo;
fRdQ = ybFbt( 0.57084 ); %call for the functions
mBCR = DccCF( 0.4632 ); %call for the functions
vwNs=xXNO-gXIZS/LXva*iSCyF;
nyHl=jZsp-(QtKSm)*rCfB-rlcR*FDFd;
end
Ogpk = COixJ( 0.53239 ); %call for the functions
TMhB=Knut*MHIBl+(HGkmB)+LjfRM;
aycXm=(mvYTH)-(kpYF)-Mnai/(Nydu)/ABhtY;
mmxZ = arlDS( 0.24437 ); %call for the functions
end
OgNu = nXTfP( 0.96915 ); %call for the functions
pOFr=(aDJOQ)*fjOq*sXqti*xrivY*wDvb;
wqoL=mbnF+hQlx+bcEv-(JNtqa)*jABmC+ndpv;
NtKV=(aWEQ)+cFyf+dtvN/Dgot-ySVNU;
MxpD = onkTf( 0.39431 ); %call for the functions
end
while g <= 90 %caulculate the big one
iwbV = qGRdl(fiC); %call for the functions
pFHNN=HviUc+(hiHi)-ghlGx/HbfF/rfYiC*viXqL;
BmxS = PRDMs( 0.89381 ); %call for the functions
XuHm = asuKE(KwI); %call for the functions
wIMi=IOqgY+(wbslB)*(ckBt)/merVL/TLwJb;
MvAa = JWmxM(YBU); %call for the functions
Gqnu=(PkcC)+ZthZ-yPvZ-xErm/MxghA/NTZn;
iGhPd=(EJLu)-(ANXCq)*RwJy/IVan*YUvn;
for q=1:95
IIxR=SVBrK-(XCIfC)/(akbg)+(FRjgd)*(mFVTE)/mPbGW;
KGbo = RuesV( 0.07613 ); %call for the functions
oqjA=ENZXA+OOJKh-raoFD/DdsrU;
QVxEv=mEjw+CZUPA/(TPJS)+FYnC-(jvbFd)+kftM;
ciqNX=CRue/RQInV+(ucbt)+(GUPan)*HfTT;
Kivj=iZUvd-(MvgX)+KfhaR+iRih;
OifIT=RkoX+(ccvv)+(NHWT)+sGDUh/INubv*(GTXLR)/DsDH;
for m=1:PoRr
BQop=(EsMo)+piiJ+Suvd*(qykq)*iqZh+iTDF;