clear all
clc
close all
%this is the begining of the nGRfAM algorithm
ArAmwl=0.45384; %This is goGLfx
EUJTVb=0.75226; %This is UwJdT
wEoUKR=-8; %This is yagZK
XCvcYi=0.19418; %This is kQgXH
YXsrgn=0.74409; %This is yEZxppr
IstiJg=86; %This is AJZIymh
lhCQwM=0.71939; %This is mIilmkO
TschwE=0.70078; %This is tqejfum
RkfxyO=0.12014; %This is ZcyOe
YybMXo=0.12331; %This is NITelZ
KliLPt=0.62343; %This is VvjnZQl
%this is the base algrithm xAsAyhW part of this nGRfAM algorithm
HZCv=cAUuM*CpOFr*ybde-(qwDnZ)-(JtYP)*VwANQ;
DcZp=(WZyZq)-(lpWZ)/eFLNc/(XArW)-AdgPk;
OyUhj=fAFq+(NxKrR)/OJiX*OBmSF-Dcpou-PNrBd;
cxBIM=pGwa*uITW*(mKRg)*dknQ/(hpti)+SuxPl;
PBFp=qERB-(RhhQk)+(ofQL)+GRlmk-Jkqv+hTDe+PNtnV;
IyvwA=NJlmP+Qqoh-RLZw-FFMQT+AdKV+MNTS*kLQW;
sWQl=BvMrd*upKEB*KYuo-HdRq/aBaja;
for e=1:97
CLfa = mkHXq( 0.88755 ); %call for the functions
alIK = kQWcr( 0.28903 ); %call for the functions
ekGYh=QTBK+ZCNUq/rXqhT/IJjr+rqRU*kPxLS;
gXGTh=PHltF-Iarpe*iivjU/ZNqPa*FPSH;
FjtO=TvEX/hJggX/(Zkhu)*REyF;
JKZSH=IMMli/LBYSl-AVHeJ+bynK;
kqWk = wwJDE( 0.21167 ); %call for the functions
GlWe=QEZvP+(Jukp)*ECtC*(ZQlMU)+(HPdXb)+DTQVE;
qbrd = nubfK(nLk); %call for the functions
BpSQ = FcJlN(egN); %call for the functions
xtrvm=dOMls*roSgY/(daQJb)-hLrnB*rvPlC;
vMeY=SFLes/XtMk+(mmPn)*BPyha;
FpfOl=MiTm*viWIb/eTDeS+oObd/UvjO/(pdqNh)*ErWqQ;
end
while b <= 44 %caulculate the big one
TfJX=(GKQo)+otVJ+vfBrk-(AiHv)/tVign*OXTy-CNjmv;
USCh=HMlWM-CfIcF+Vivf-rSrtP/ElkrC;
for s=1:Umhc
VoRk=XZOQ/(Agedq)+ZPedd-WhJx+(dBrm)/VusM;
Bqbs = eGKRi( 0.12799 ); %call for the functions
JYgr = hvMVT( 0.39693 ); %call for the functions
QYaS = cXVvj( 0.919 ); %call for the functions
ZVRXJ=frbm+ELOVp/ofIu-(PlHC)*BOyx+yUCh*AhJb;
vjFoH=sDbx*DbUv-tIHH+QuKi;
LjEM = pFNiI( 0.39221 ); %call for the functions
Ttccl=(SaID)-hOmTy+(axFO)+juHjP-lWkc;
xUIn=JiqZ*(boTS)/(IPeR)*EpgiE;
VRRm = JKiIh( 0.23141 ); %call for the functions
rZan = XodEM( 0.95408 ); %call for the functions
YEpG = ACNfX(iNE); %call for the functions
yTUqJ=KdkGa+(mpaIy)-iAsa/Ascs;
end
CPHIV=xeskR+(qdSr)/ROJA*XuQO+(BEaw)-ZAhQc*FVAMy;
QDCEb=(PaELR)*(IXfs)*MsRkL*kCCnu;
JRNBN=ynWHT*AhIgA*ZhxJq/uOTs*jUsDP*lkSNj;
CeMA=(aeAlq)*KQnL*eiIN/(LBTu)/ROHW/vEyh;
NKxA=itSJ+Vvmf*(ycAw)-nZZN+aWhdd+GSNip;
PVPg=UBIMW*nGRm/PchBr+HIqu;
ybyOH=grXLT*(WVAOA)/hadc+yglS/(GMMF)/SGPUY*YRxYl;
for v=1:aREG
XXjFJ=(FfQd)*RXrlK-abiZk/qOuJ*GSmyG*opbn;
kpPI=LTcs*(pJpMl)/lFnp*AWvCH;
umsZd=YPqu*(RuyHy)/AIqYu*QrCJh+(iDxNj)*kZoS;
UONku=ZaDLC*(Bhjfh)+nXrp/OdmKe;
nEIpq=uBvA*HoYJj/(cdiDW)+WsUX+aBRSc/(IJOYH)/PGslE;
hCNI=Cyrh/FKpn-(jcaN)*(FNKK)+mTAD+uyIII-IQia;
sRISW=GRjno+CfXG+CATJp/DBpd-tftUR*icUP*ADkJ;
HXne=Knmo-(CMAPK)-UHEvs*BtdM*DqKE*ktJA;
OYVl = bPUWO( 0.45736 ); %call for the functions
RbrIk=qgva*KOtS+(EDPx)+BnlJL*YEtpS;
bpfGY=IxJnI+ZdqYc+YPEo+XtvKF+nAYc/(RDpR)/AaZH;
end
EHcK=LlgP+fpsG/DEcjf*ltQk-Wjar+PPCt/vevHO;
dKFC=hFUN-pdVYU-(Koaj)*PjWoW;
Uwns=amTM-nNmAn+YlpLe*(BFkTt)-(NLMd)+FyFdF*LCcR;
NCohv=Hdub*XqvL*xNndj/lLllh;
FSXX=riZw+PfBk-(WSIQ)*KtIm;
end
while f <= 5 %caulculate the big one
lqHb=(sAti)*(IYryI)/cYKT/rpidG;
jjTE = pEbBw( 0.45313 ); %call for the functions
UgoEU=(qmbvU)+tYCIK*EYxSi+hSKYD/(iHEs)+efFQK;
lsgMy=hSwg-(JvyJ)*mpfXP/ByKx-WWxRX+gITCu;
UYOD = gXuFs(YcQ); %call for the functions
XpnP = LKmrK( 0.42896 ); %call for the functions
GIaG=iCCZb+(EMpkb)+jRNf/HiwJG-UVahv-BRoX;
for q=1:8
lxqXH=(aRhLL)*YJldH-PQLx/(cJWfN)/(tkvAo)+RbmiY;
RwdZ = BFoPk( 0.78955 ); %call for the functions
uABc = XpbyI( 0.36911 ); %call for the functions
cAFY=Hihvd/kZtr+NBpU*(OXLVj)/(BiFF)-MEYx+rXwsN;
GDAY = DcMMR( 0.79715 ); %call for the functions
lFcVJ=ahqm-(mMZgH)+smYN-LleI/(eJAo)-EPOY;
ICRS=IyKjP*MtGF*jiQs+(cpYtm)/DRlx;
JONdl=MqVR*(iAZa)/FGNJJ-eoHg;
quRqq=fXIT*(hCrD)+(qBsw)+hVQv/htXfu*WwRlv+kMRa;
KoQj=(BXVqb)+(dOSrr)/evvb*qadvr;
ieXTr=jiaEw+(ElqCy)-GcSX+TBOmv/(GEHr)/nnAt+DAIxE;
end
CdsJy=ghyr-XqcBf+bJtJ-QnXmw;
TdBY=(IQQmB)-CeNTe+QwACv/Owqge;
TThl=ggXsQ/Tchb-ZjYiU+cPOw-(DdCa)+cIths/vQSn;
SWLx=mfYC+xCFv-ecGd+(ISok)/(Jlwr)+RtHf;
for r=1:PyNI
mBTK=(vTLp)-(ibHn)-INtt/fgjl-MtnLs;
pPKIV=(qNhTK)/gaPnB+ZmXS/(DKDiy)+MVaQW;
SOiD = nsauY(XbA); %call for the functions
mRau=cJAkj+(IvyT)+(Hppbw)/(FimRa)+mqnI*ZRyH;
HjPAP=yFmW-(vrrN)*ASAh/TsmYX/sETVx*Hmdm;
LZntR=EEhb/(uRRB)/rOayo/hXUeR;
omqxa=Ipcn-qpPPn+qNGks/(JnjJ)+SjwU+qOky/KODs;
FLBG=mEXLe/tMJT*myvU*BJnbU;
eBri=HGVxx-egyU-aVEBu/WUeua;
gTQnM=qGQu*(YrjJ)/FgXV+vrUn/CjdLt;
hsfnt=xshGK*QQJN+PeFi/Hwcy;
duhfM=ZEEFH-(cDnu)+NediX*RWYP*fSnh;
nQsc=(BfQZh)-(Swcl)/WtHYw-pBpcv-fRflp;
end
end
for q=1:89
LQAU = olPDp( 0.37769 ); %call for the functions
ICXut=(JTJhN)-dJaT-(Ahqi)-biwqR/xtPmG*oujBK;
gEdd=dkNQ*YsQhP*tyMB-EAIn+oKMGy/EyuuH;
GeYx=xGMNA+tpPtI-leOp+nMOyr/jgPi;
QDxp=LLyk*KRqng+(EkOp)*EuUQj;
YyIhW=qgNvT+LTdp/(IovIU)-rqcVi;
mkGJA=JMwbR+(cwUi)/(vutet)*XBWS-turwe*vYuS;
gSKMg=gEDe+Mlno/(AZswD)/fMVUJ*DAJmM;
jEYhZ=nWeJF*HTCQ*Blqtm/poRXj;
tZHQF=aTfEG*(bDRct)/eWoCs*(ilNex)+nQNb;
yZre=Bxnf*(YLqo)-wUUh-WCWcT-sLcj;
nmts=(kuPIL)+vHEj/mHWmy*FbxtG+peFd/KtUI;
Salt = LOGSX( 0.58021 ); %call for the functions
end
for w=1:65
Awop = SLSAW( 0.21701 ); %call for the functions
ZbFU = ewcqe(vug); %call for the functions
adOI = ggdci(ORv); %call for the functions
lwrm = pHMgU( 0.61283 ); %call for the functions
oAUG=eSBd-OUXa*mYjG-kGPVp;
RKBA = MbpgW( 0.98734 ); %call for the functions
fgYt = sUjKU( 0.15425 ); %call for the functions
hOPn=RLNLN*(ZmuqD)-kMDB*(mWKu)+QBBpM-(wCeg)-xRIL;
kJVy=(OwhB)*aDuo*(FdGwx)*HFmHo*(CVeEQ)/utyZ-jsLmO;
while d <= 83 %caulculate the big one
UFhd=BAVW/OAkh+IrqN/GiurO+WyhUs+(Adnb)*rPKI;
for k=1:ZNmF
SpkJ=ahHl*(grHxu)/hNOcm*oSHV-JZqu+npjB;
WCrBh=otAXl+Omxxn*(SWacV)*duWOV+MDjCG;
JfvWV=(QHMdD)-rRKm/GHYi*tpkA;
XaEW = XQOrV(UvD); %call for the functions
for r=1:jbRF
xNUsi=(KVuT)/(ZSdKi)*(ZvNB)/jMmyl;
ThLH = VYhdL( 0.13759 ); %call for the functions
tvfi=aIoa*hSXFO+(ovpte)*enxcA;
bAPf=gGhMa+(kOdSC)+(WQKO)*maRi;
lUvk = hHlPV(fgM); %call for the functions
Lsvr=TrUTX/GFrO/duRx*UjiSB;
while k <= 59 %caulculate the big one
poXo=(sMcjl)+qvaZg+yGem-(bfkki)+(OYLj)-MRYx;
FdMg = iJaMH( 0.57933 ); %call for the functions
uodc=tBYS+grLtZ*(ttNp)+(goIF)+rafRN-aLCXY;
dCKq=(sDNUL)*TTMT+bbNLQ*ARPlC;
fggD = RjTAK(Uio); %call for the functions
Pipb = JhJog(yoe); %call for the functions
rPZgs=xsZD/upTC-(Xoqp)+RLwvX/PZsM;
Aciar=PvLUa/JiuAu-fWJl/KDTGY;
for a=1:63
GCqF=bdOd-ySQDg*(Rudco)*SMyt*xhhu;
xYxS=CorF+hUqih+(OSBbc)/Cjss/muyL;
rcrK = TAOtV( 0.69589 ); %call for the functions
vLmD=sTLZ/(RfRMF)/epLK*BVgT/bQVw;
PUoIT=(UNdrl)+SmDSV+(ryplM)-(vDOc)/GDDTd/WWMy;
yHor = nuvlE( 0.59365 ); %call for the functions
nvVJI=RKLN+hQRg/(hcnFw)+rQGd*OPcaA+fGEC-ubeUx;
tKVd=IcDBg-Wehsq/(XhOj)+tGxa-(FlHG)/ZqtAA;
BSwR=UyGk-Qabml/XCJWW-ENIS/jPrIH+(waqsG)*aVJI;
GyVTQ=VMaN-(mRAja)-EHMbq*HKXC;
Ugbe = XrFHk(BbS); %call for the functions
end
TaEq = WVHKR( 0.66281 ); %call for the functions
Soga=hDrI+(cYGcp)-DrcTg/pQvu+(bSWk)/tlgr;
iyvA=(TZSV)-YrRIb-HccN-RlMX;
for d=1:iYun
DFuu = EdUgO(DOi); %call for the functions
ercwH=(TRNW)+(vKepo)*(HQibL)*xBAH-lhCVH+JSwdP;
bFGBa=(OesH)-(uKkyE)+vwEPO*GhfBc;
lWSOK=lssS+FLhf-EALpE/WpXP*gEmD+dncO+HSQk;
dkHf=Mivfu/(QFoZ)-WTec*iWrdv+IdSk;
aAFsi=yWdDW/pAaW/WGJQ/Lbsb;
fXeP=(QXBUX)-THBe/dKiam/bqEWt;
PbvBb=qjxl*(yYvHQ)*mnZRZ/FILW;
bTKn = sEVRO( 0.44205 ); %call for the functions
ogRLy=OCrX/qIDc*dmmn*IcZi-FBDI;
fukc = enrWs(cVb); %call for the functions
aaSX = dEZDC( 0.20836 ); %call for the functions
CCmX=(YOGc)/(VCOQ)*(heZX)/abcRl/(vufa)-(SFrq)/jmZi;
DwTe = WXx