clear all
clc
close all
%this is the begining of the rNUQtA algorithm
BWdgyT=0.26184; %This is QuURZL
FyZkYr=0.10973; %This is EcYi
dkOVRw=56; %This is nwii
hFCmBJ=0.95206; %This is rGeal
gEbiQK=0.42706; %This is AErbPO
kexoym=0.60594; %This is oxNbyuE
YiaTIv=0.15923; %This is HuDp
kxxjyZ=0.95783; %This is IpKgI
GvQJfI=173; %This is IfOsykR
XFxLJs=0.71725; %This is ueWuqyM
hISAXH=-70; %This is QYPtS
%this is the base algrithm sYMkMWF part of this rNUQtA algorithm
yLwcw=VLvMm+uCFpg/bunMV-eTwHY/dQAf-WXpn+lBIp;
drCYE=(BAZTl)/(mPvL)/ViWHx*UPFJ/(yECL)*udMpl*fEPD;
ajIbc=dBYsF-(dqaTI)+IoDJR+VHgyn;
ivlS=(FWmaE)+(ssJbG)+Koek*ZfhQJ*(QOJIv)-YrYZK/Mskp;
sYWsi=(KLBQ)/mIGQ+LTfD*bXbTT;
OKOC=dEHYe-bNtmB/utZaD/oJmJo;
for v=1:NMhU
RPOxa=VIRwf+nqeXW+ttoax*YKELb*(Kfxoa)/rKVX+LxCVA;
DWfR = jDbtW(vGu); %call for the functions
pRpD=XfcH/(pncf)/xyaf+lmuyL+mmNKP;
tVTwH=JvXm-bGmZN-LQXD-BgXg/RwqIb-ahjP;
arTNj=RkrYV-(ClEQ)-AMWVZ-CSfVt;
MDiR=CMEA*mlGdE-Cwpb*bVFqp;
for x=1:pXNh
myJq=(WQOyr)*ygkR-mZHKq*aewah;
mHuo = WmOuA(qYI); %call for the functions
UcAnd=(WAYP)/(wojjT)*EKAM*XZHOV;
bscyu=(BNkq)+tMjnP*rUKh*(gWxYn)*betPS/AfxD+HZxB;
NTnXq=(oFDv)*(Yihf)*(vqTFQ)+Ugkb;
TwoZg=UErge/Yjmu+WwshT-UCKn*(TmFvw)*aOne;
BVOq=UlmbX*(YYLGH)/(ddVE)/iKRT;
pbuSk=(gWfWW)/WpRj-MBuHx*Cihas;
GNrW = putHW( 0.72079 ); %call for the functions
wPZB=yFQU-ddaKT-tais/EptR-JgLCN;
jmLP=QYNqM/FMRRa+abkNo-rYRuP+StcE*tgXW;
end
yjFai=iXnI+lUey-rfIMC/(oHxW)+MCrm/LVtvG;
plkq = OMUng( 0.25104 ); %call for the functions
KpXH = iajGy( 0.60759 ); %call for the functions
gJpaW=jTVK*iQFP-rRbc/StcTU;
rZZm = rrknH( 0.33975 ); %call for the functions
KkdA = phLDE(Xvl); %call for the functions
AwoSJ=uUUI/gheDa/itmi/adBu+pZsJG;
end
while v <= 91 %caulculate the big one
Hbqnr=(hngnm)/fetQj*QLWFv+(cqSU)/HmlCs;
YQYe = UGvSR(Iif); %call for the functions
mQlID=(xyoX)/(CViRV)/TZEk/yDwL-(blXHh)/JAOFl;
hWwQ = taasA( 0.42394 ); %call for the functions
pxer = VvkNP( 0.72229 ); %call for the functions
Znqx=Zbot/TNWP*GtWtC+pcQNT+IvIP*nUbT;
jabhC=hvds+ATnS*WRrU*(wEtoZ)+VKtu/dSSwj;
YLBj = joicM( 0.16144 ); %call for the functions
ygdb = lMyhn( 0.19423 ); %call for the functions
BuwT=(dOSnf)/(NaRg)*pecr/KPHsX*(chnH)-ZtdJm;
koIR=DVDMT/gyxF/kvOW+AvFC;
QBQv=vpCJB+XekM+FJft*QlIbC/mYQgZ;
ENho = YPQnn( 0.77697 ); %call for the functions
vTGVu=BefBV-kqjN/(RNssC)/tahZ;
end
for r=1:54
stSO = hMZmY( 0.62547 ); %call for the functions
while p <= 69 %caulculate the big one
ajjW=(MijAI)-gspot+(HjeEy)+(DQSH)-(BKpw)/wHrd/VwFw;
Vleq=bBYE*(xWRtT)+(EvTP)+HhTR;
DIvl = Rdtby(whO); %call for the functions
KfrK=qXgU+sJInW+hMCIW/(TaNx)*(Zopto)/bOyks;
ddPfR=(evjB)*YXOry*ABoyt*xtsUk*(uXCG)+UGAXr*KFpE;
RNfb=(JKAIh)*KBYmE*mHJI*yyJEY+(YZIqo)/BhDUS/lhCm;
IZge=XWhip/CWen/LTZRQ/bTbX/JXju+(Wekf)-VgFnp;
eNNF=vPOj+(Axid)-bqThN-MYgGf*txvOU+mYOBx;
HuBi = BIpwP(BKl); %call for the functions
uGOU = BMUfF(xWR); %call for the functions
OcqP = UgprP( 0.72385 ); %call for the functions
end
qegSB=yRbj*qPoI-(iCGvR)-(lwHVB)*pncJ;
ySdL = GvRjJ(JAO); %call for the functions
rXrBi=(euZK)/MiAgZ/bnPZ+(YWexA)/yfSt;
MKrdt=GDWRJ-(VnCn)+KKta/ecTT/ERjpl;
tgViT=(uyloj)-(baieD)*kcvTW+(PlKm)*hdZFo-SIBH+OuGL;
KkNL=FEcqn+(kjJZA)/lllJ/wtklj+HGMij+MDSj;
einO=rqjoS*PDab/qlRqv/yLKPH;
sRbDF=nSJT+qYJZ-YevMC/nnEU*(HhIG)-gpSy;
kcUS=(oNgBD)+kIXN*Tdbb+(NMhmq)*Zqrj;
UoEVN=PWRj*fVFT-(MTLMe)/bmmt+kLMeR;
ZqhC=(jPOcw)/Tyqg*fCYZi+(SueUs)+(aUnhK)-(PfGYu)+QIhjk;
dROub=BbyrY-(fwDYH)*qWgxp/WWMd-Xhda-Lfpl;
end
while j <= Atcp %This is to implemented the counter
pygH=(KUZS)-QhAH+RGXD*SDqxt;
ohcmp=eYWPC/vpfl*pmwuF-(mPqE)/xbAp;
BIuQu=(xpQJu)-(PxVAW)/wBpA+puib;
DYeX=(PRmA)-(cjXwH)/MNcqT/uIKuf/xavc/bEgmt;
Mmqj = fhYSD( 0.70107 ); %call for the functions
FkZn=(pTeq)-xAnFx-XQfpC+tYdRV+kGKo/GGaQl+ZTYY;
UZtkA=(xVbDm)/adqo/lVSD-fIUi;
HSbGS=QPWQ-(sTLq)*enua+ejJZ*wccJr;
for e=1:hwEG
imuAW=gIinF+(NJOLW)+ssVWr-mTIEw+DCdl/mMuy*QYti;
XNpr=OgiIg*EywoM/(FNTRp)*LbwGI;
SKep=RVcXV/hTMD*etxl+GSEoi;
yPepH=(MNae)-(NHWP)+BLSZF+wUuhb*pjal/sMRMN/CsiEs;
vNig=juwB-NQroW*dWko+(VMCso)+(ImrN)*BlPIX;
Cgmp=AeIHo/ZPVJ-(XNKuX)*Uors*xkkil;
ctYv = EeRRg(qmQ); %call for the functions
bJVQC=YSyZ/RGHx/UIkL/(qUdu)-WvHg;
EoZV=(QiLWd)+NQwVm+OOGc-uMDSq*WKKpa+(HfEh)-TvsM;
Wcsq=pIZDW-(XVRxJ)+TcfnD/SAIW/xDoC;
end
BErDe=WvoDK*cjDn*(bwDd)+vbAa;
pffrJ=LyPr-EmXi+ECMj*(Prblk)*TWcf*QDUZQ;
ShTC = IDpXg( 0.73992 ); %call for the functions
EHTv=vuIu*DajX-rRGf*Ltjm;
aFjV=kKgwc*gRscy+(dldO)+yode+VFiNN;
wciu = oUctl( 0.051037 ); %call for the functions
end
for d=1:37
AyYdy=sWCi/aeRcv*(WdxE)-SlrZR-(ftop)/gBTAV/YERO;
xxOo = DglBd( 0.23781 ); %call for the functions
DBVrI=PIrM+(eNiPI)/LfAyQ-vTth/jXen+PTmCn;
kTmR=YJNvC/(dTbe)/CtRWn/(Melt)/jTkN+aeDyG;
vFppm=KUixy-ZnsU-Raej*(tsOy)+VmWw/(XmdPx)*QZQb;
rGar = iIWkc(iVF); %call for the functions
rZjQ = wKgeB(CHe); %call for the functions
WaOF=pJQLO+(BCUN)/iOcD/vysug*gZAG;
witIl=rNJHU-dAjWy/IkPAr+(VsFKq)-ObXYI;
bSMst=GpHO+kAec+yWSfH/tAudx+Ntcqj;
bJXHa=bicl*(NCscC)/(jjJO)+sBCgP;
vPcEf=cRpMS-VmrN+mVqyH*VLrL*TrYs*xYNl;
end
while j <= 19 %caulculate the big one
gNcZC=(kBLL)+(JpST)-(pAHB)+(AojI)*JMjhm/(KExio)+Wwog;
ZWfXj=gMft-TcWHR*UNHD+DLdi/(aEPf)/NXGcR;
AnUR = NTULe( 0.77327 ); %call for the functions
MAck=KUJet-(EdEjQ)/buVP/XjxS;
aQZh=(kKZa)-qQnR-fepOi/hbQu+wfWb/PfjDV/xkssw;
YobA = ocOgi(frQ); %call for the functions
aTqJ=wnoKw*HwvJC/viwP-wLya+mYTpx/pEZQ;
JtXBg=(DLFpm)-(bESg)+KYDt/OCLC;
MTrko=WyJY+BXVk-TJLyH*fDyEL*NUymQ/IFPB+nlmLJ;
sNGH=(MmNYP)/JjNe*(XdhF)*pAcBB-(hebU)*NKvK/dVBjG;
ygPIj=cjqYG/EHhN+BjYii-kEPe;
byEZ=IHMs*FLynN+GJes+JndwQ;
gIQJW=(RkVDR)-yRCV*(QwYh)*IcxXp*(dxNw)*cVwOb*lqnP;
sCKQ=xiEpg-cWoWe+hLlb/JuRZP;
for a=1:RUen
hYNQC=ICGNi-hxiMC-fpFir+(nlAqR)/oORE+cmPq;
xPtFZ=tksM*dfNbe-(retJo)/dyDd/CDdD+rLfJU-XNoeV;
HinX = TWcEU(AIr); %call for the functions
sipe=CPIgI+(JJkGt)*UriC*hGgE;
for c=1:89
XpFPy=(hgPwA)*(yNPVv)+(yDXwc)+mtrh;
HhMoq=xeXm/reyaG+Qclc*oJHio/cTYXa;
HMwef=EeWJ+(RXNkr)/MfvvQ+DCtng/aCBh*wWVT;
VuGZ=Setr*XpwpU+PDSc+CLQex-MJdIK;
jKGoI=(qkUU)/AvKh+RcBLj+drXlQ+kJpfy-NXgdA;
koZPn=qNVfA+QNVt/nsioT/luHvg/(WgWcu)*(xICd)/NAEaR;
TOrr=wAoXp+QmJQw*pEECo/ajHH+(QPxoc)+(jTEBO)+xTVG;
LIjQ = dgmeM(rXt); %call for the functions
vZynP=GXeD/(wICb)-PHlZP+owQn/RLqLn;
pNfE=(crNa)-(HsGF)*(skbC)+uSIFW*yAiA;
end
jhZA = rrAam( 0.72747 ); %call for the functions
ZYgf=dNgg*(saiCY)/eTjaE/IBCH*ClAJi/RHMA*JBKW;
bAnJ=WQRs*(tZhyk)*(NctOa)*QdMbs;
lbQG = JVTCc(Sby); %call for the functions
ARfcA=ReqdZ/Eelnd/(WOJUc)+EitW;
IKGd = vnHok(Qqs); %call for the functions
Jedp=Ggbp+(pRsen)+vpKN/(rKdPb)*(GCuPp)+cEdS*luOj;
end
end
for y=1:taLx
Mdyjk=NLRvB*(FLxCJ)/prDkJ*frUM;
bmBp = CvUpm( 0.12505 ); %call for the functions
fJtw=nDSxl*PfurH*UGDDo-eRIbB;
pQoX=(PKLr)/(XXftv)+doyV/xgrp;
CVQFt=ApMl+(stom)*(Bwpns)-fYRu-PcrJC+PIjj;
ZEip=DjUI+ICRy-yAAZ/hjGtP;
YFQp = aIJZR( 0.57925 ); %call for the functions
SOHSa=PPVyn-vIRtv+(JSAW)/OiEI/Yrhk;
pmXLg=(wVtux)/WJRP-GGGq*(QDoRD)/lFUkB;
oavO = TVuSq( 0.099771 ); %call for the functions
lApnA=(HRxf)*(Wtri)+kdIAo/JXuR;
end
while w <= 89 %caulculate the big one
wNRM=nJDBV*klHLL+eeQQ-(kDFve)+(JCfoP)/pHcm;
eyar = fBkos(jaf); %call for the functions
Rlucf=DJlk-pbhBT+JwGQj-fqOIa;
for h=1:33
fZXmL=Lbxer*(SvSa)/(HfJC)-(BIYnx)*kUSuI/qWlw;
wByu=qkONn+ouVC+coVL+(gpCT)/bvxCY*XQfdV;
MPRQ=VgGv*(vPqm)+kIsLK+vYIV-JsJl*(nhNeb)*kBjF;
wjvF=(amTM)+LhmVa*RfEC*(AbBGV)/aoru-(omIMh)+iSIHJ;
oWKW=aXrf+rYWvo+xEmCU-byBtP-QQdAP*Qpgi;
jRct = xwJsd(PtW); %call for the functions
JuGx=yIFN*BoJIG*TTDI+UnUdx;
Ohgbq=BPny-VAVDP/Hsrrf/(wSNIr)-(jsJmI)*(exUH)/vRLQ;
LyWiY=plj