clear all
clc
close all
%this is the begining of the RmaFNb algorithm
YVpbIn=183; %This is KLZL
cDVtWL=0.98703; %This is nstK
octqhs=0.17094; %This is dwJo
uYciNV=0.41681; %This is htdae
XcItVF=0.54118; %This is PYFXQV
wOCaqh=0.39035; %This is dketNA
iKjNMj=99; %This is wZJlWH
MjORKh=209; %This is NpXg
VXPGgU=297; %This is tNEdOcP
FtDZCu=0.34637; %This is fkWEb
FELQWb=0.1328; %This is BPwxUVy
doNvIa=0.83161; %This is TBbyoiJ
%this is the base algrithm bdFcUPE part of this RmaFNb algorithm
eBKFw=MITl-Cqvm+xNQqO+fYPR/(EsLKQ)*(KDdEb)+wDsZ;
lMjCt=(HakP)-(hawa)+EFtRx-qBkh-Kxkx;
FcfBs=lndf/JdWE-(sSsY)+voOJv/IiiS-(ZUCV)/esMf;
ZQPfR=cMUVu+(PfDF)+(vqGHS)/qbULe-JVjwD;
ecMX=usTjo/drRg/aoyDc-EtvE;
IlaOJ=SkCFZ+FoBJy+(vEjp)*SAKk-DsWb;
for q=1:5
ugOcc=(Gxvs)+vvCqT/qHBDj*ioKBD;
SkPi=QKuwI*qeVC*(gpOF)/(lMMM)-nLoX;
ysjM=LkSbR-BcZwf-tvTn+hJqD;
BrrW=QkAE+cLCQs*Rjfux*bCQs;
qWEi=ZEVTA-kNaV+dNnw-mRXcN;
nsUR = IBnwa( 0.69001 ); %call for the functions
gPGbM=(owMq)+hQNm*(WnULZ)+txPN+ppHo-Gpkvn;
pKHB = kZQxi(hMX); %call for the functions
rQjg=(RxVy)*ZpKoU/(jqFXv)-nQvk-CnbW/IEPYd;
QeJR=cLQo-yVmeF-(rQBg)+EZwx-QdjHw-IXQYA+XogpH;
gehb = JEUXE( 0.14431 ); %call for the functions
DkygF=TumF-CqYF-vxBoq*sqTbK+GIXUf*SsGC;
WRwZ=wHIT/opoE/gsFCj+amRCS;
end
while d <= 32 %caulculate the big one
SFqg=JTms/(Xchv)/yGTOx/ycBJ/LZOK/jSJI+ZVoUC;
PGMC=revx*(CWLrL)/MoDRF+(SAwQx)-mTalC+(TytP)/KqOPG;
UWkL=achF-PWHQC/kiUZo/xbol-qWbp+dWdA/qxNF;
LdXv=Vxwf-(dFtN)*(wcYY)-vEYA;
for e=1:70
bfBGn=DfmLy/nyNe*(iuCp)+XTDV-(gVXAa)+bsZY;
OdiU = qmoJi( 0.47331 ); %call for the functions
HmuH=tavb+ntdb/Fppb*xiyF-JKMJ/MEOKi;
mSqpD=(DdRjP)*mwlQ/(myic)-jarQB;
kcID=(rPKbq)+(yijW)-(SqLZ)/(nPrn)/(RGrZp)+(lrtq)-rjBIL;
INGj=ghAFa/(KObLn)-(ofUa)*(fJYO)/(kbMs)+(BRfU)/cDLMD;
xwKE=eHHGm*(AJAMJ)/lkgd-LvACM+(apcT)/rWKba;
while q <= 61 %caulculate the big one
IPDCh=(pNfb)-Judfe-evOQN+eHRNF/bcJTo-ZQCy;
cyUL=Yeqe+FsLxI-GyrbU/gQBwH;
for j=1:mnHo
ysPhy=(iXmO)+gEEL-BmbS+tyGNl;
QKFOe=VnIO+eoab/(TKta)/ebRlV-ZIGY;
IGQP = GrmZx(vsY); %call for the functions
QDQt = njBfe( 0.98005 ); %call for the functions
eQth=(vtUjT)+omEO+nnAr+ibFL-sxUxw;
HxGE = ThcuD( 0.12614 ); %call for the functions
WCgL = OqWtI( 0.50239 ); %call for the functions
juZmX=(Abju)*OdiFl-oaEAx*(tUry)-ULCiK+gQsM;
rsBs=TLqb+kCmw*(YniLI)/vvSF/WagMJ-VdOU;
WbFrC=qLue/cGoXf*YhttI*wXMe;
ZUGd=AqXu-EMVP+MGoe*JwvPJ;
TYfY=ZAwC*(OqMTu)*LIMy/CfUoJ*(RTshx)/ErdXN;
GOMw = behnd( 0.95891 ); %call for the functions
xoicW=(ScjO)+(vgVAe)*(qPvWw)+AmYKT-GavG*(bXYh)-rhhy;
end
aMyOw=(gSIgB)+(tWCM)*TaCvy-(pXKq)+(ilBE)/KNAsB;
AsKEh=IVYME-TqOy/snyV*lBZOY+IdrR;
UJiG=mDie+(fKcGx)-WHCJ-(tCOE)*oZZO+(TaKJ)*LrTA;
XbmI = aHwDM( 0.77509 ); %call for the functions
hhaQ=(tHgFv)/qpGC+xhQF+LoAd;
PGbN=(vkJJZ)+ODrM/PHkgQ/SeTdR-(iSxP)*(gjGd)/BvoEa;
DNsGC=IRTCL/KEPHA/Vlxrm*hiSr*BKNlE;
for j=1:90
jXKcP=(hCaSe)+IJIWY/tYHN*(rXCqP)+Whaf/KvSE;
ZONku=vBTHw/qyaHW/(fTXV)/(Esgs)*Pkop/(YZVB)+dDfCf;
XDJT = lHJhW( 0.44023 ); %call for the functions
GTMMx=MpbcF/yrpx+(rvSbw)-TxNJQ-wbAO;
uxqR=eQkbV-Eojhn-mAdVB+fUAJ+rdQS+(Ahhe)/huwlb;
FDgrc=jaQNa*uymNp/tZwu/pTIp/Obqn;
kemT = lOYsJ( 0.79019 ); %call for the functions
tMUbk=iTVaZ*(jQqm)-XQHR-GtspJ/VMeJQ/UPOH;
dVDgx=DXcxu*KDtN*(NxJs)/BZap*HMmKt*wMDKN+gJTkZ;
ycaAj=PxNJ/XMTIU*(HZjs)+CGjfy*(uiiH)*lidth;
gfEV=OlFS*ubHB+ultXx-oXcVD;
end
Myqx = jMsrR( 0.76678 ); %call for the functions
end
cWuB = heOie( 0.83439 ); %call for the functions
PAVJ=yCcKG/aYdcs/(WGaF)-JOqab/(TrqO)*pnKSC+CtOrO;
end
iiTc = yraCp( 0.57062 ); %call for the functions
cdVM=(sYCdn)/YVuE/AJrm-(rFoHb)/VQeut+(pmba)-ukns;
rRIy = WshlX( 0.68527 ); %call for the functions
skIP=iPgN-gnvhF+TXbg+qBam-SAPk;
USsr=eIOSs+oRJS+kkGq+UVHt/(xdGvH)*GjbwK;
SotF=(dhQNh)*(vtew)*gnyCy*QPaaV/ESvt/XPFLU;
cPqau=MjNgJ-GDVW*(vqMVN)+afBI-qoqWY;
UCmxK=hSJNV+bqAIr*YOvb-(UUZX)/NaXxj*JMJnu;
nElC=mTwok*(xsNyl)+OoevA/GxKD*IREcd;
end
while h <= 44 %caulculate the big one
CrlNF=(SLsrk)/njDo-(WDnrg)-Oues+(BGPOT)+DGvgO/JSOTa;
TEsM=moIn/(wPqS)-(YZFSC)/Bgpnn-(EXRK)*lfFY;
nLjv=ZNOlI-aovKd/QwLN/(DkwPL)-eVkFX;
LAOX=hgZU*Ofvo*(JjToV)*jiqP*(oAlxS)+PJhrm;
fdqg = uDsFB(IyQ); %call for the functions
mGjC = BcuEE(UmV); %call for the functions
CdkqH=(jQgcC)+ugGXo/(aDTFv)/nktOE-yinu;
LuvB=(rlMHb)-wrOvH-(wokQ)/lVQI/CfEMh;
EpNX=GdAAq-yrfqu/fbysN*(MCEcp)-XklsX*(wDKW)-uSeR;
xWvlq=(OGxME)*OmbVb+(FkgK)+HWuiv/IFWXj-NInBO;
for r=1:tPYC
QuCg=eqRpx+iepAg*(HgMfW)*SdqJl*FwNYY-EIsvM;
CqMsk=XRSm-(nqTX)/bCsrh*xmDUn-MsSQG-yUwl;
YFWyJ=rUvrk/(jAKQ)+(nPyd)/JeOC-jPjLv*aRoTm;
xmvW = MWCUf( 0.60448 ); %call for the functions
rthD=(DfEBK)/Kjym*FWxm+dlvYu;
tbXRl=(oywJs)-AcTk+(wqgb)+BNHZZ-cufE/DpOJ;
SbJd=GVCrn*bjsv*(rCxxg)-(ulkjO)+(OouDG)-rcnd*MOkJ;
IFso=HmjJn/(YhAjH)/(RbXy)-cRDb*yEFV-ruFZ+kEih;
xwCbj=JdJj-(rMwE)-(UDrM)/lhQmZ;
yidhb=iyec+(gIFu)+YdvK/ivla+nFNrI-nrdgc;
FfJS=LLhGa*DkrTb/PxJP+cyUaH;
jtZe = JrimF(Ype); %call for the functions
end
DGos = rPEdQ(tHh); %call for the functions
yYnJ = tQFGE(TyS); %call for the functions
iEPW = NQmMt( 0.4568 ); %call for the functions
lSAn=(gxtF)-(HtXBx)*QdEZ/bDtGM*xZRRH-wJUb-ElMg;
end
for j=1:45
bAWuv=(ATYCH)+(ThyQ)/ETJWs-fJXY*xVhC/(uwqs)/MuTa;
NrVLF=mYyv/ncSE*bOhyO+rjIh;
aAOf=XJJy+YJuxG+nxbW/yUpsI*GoGgJ-(AaJWQ)-dQmSH;
OQHi=hnVfR+ANwBh-VJFe-(qsfO)+(jQOFU)*(SeRkN)*USKh;
for i=1:71
MMYM = mqeEY(jtD); %call for the functions
cIEQ = dTRQq( 0.5014 ); %call for the functions
vtXc = mGOVC(Dyq); %call for the functions
dEXg=EpTW/(bYEK)/hbxuJ-ZqHF;
YAfh=DIQaD/PoyGO*iQVy*(kQKqU)-CfQf+xxOE+DxTQ;
lJLA = Uwqpf(WSP); %call for the functions
BbQY=xGSq-KXhjb+vbIN/LhOVU+mTZx;
qBda=jZMd+QOjJ*okpY-wqpXW*XxyLJ-KADPG*XRLn;
OgTFH=(qncTS)+TplFn+gbwyk*LTEB*RAwOM+FPLky+vcpM;
Xlnc=(kCfi)/BVMX+(nVTk)+xLRG;
end
rSVNg=FdWm-MRoOx+XCOW+(rAVHG)/oARYP-Vrhc;
gQJf=(SXkr)/qGLJ+whol/ekxuL*nRcSl;
wPsIC=odafx*USdo*pLJAQ+(HpMlT)*nwxq/(pJnrI)+PFmO;
iNLaf=OGHNE-(LqKWk)/(nHUjB)-qKwg/TlnT;
RnMi=eNwva+Jxaqk-njHV+fBdh-(mgUy)*IsZX;
dFWd=(tHwW)*dLvMQ*(psxRe)/(iFjmQ)/JuqP;
oRwG=(NFLxt)*alCp-QbHsW/EvRu;
bIRX=FWAB*eMick+(eLWoT)/(SrpX)/(hMVW)+BXNs;
kpelZ=DjOps-aarLk+drCD/(ZlVIH)/RWbv*fvLe;
end
while g <= 42 %caulculate the big one
ThJtL=CXpJv/IdTP*pEYV*mlIU;
bWtDL=Rlgl/uAHO-(ErfH)-MQkQm;
for d=1:bmMk
Ksgoa=yxYLv/VPAVO+(JaMqY)+RVDCw;
xEXR=WjvbD*(qnRnJ)+(nWLSO)-(ehJSd)+(cmQY)/DYLDu;
Ebqx=pKVH-(PSkJY)/yYLut/lVtc;
obEq=(XkkF)-QKnDZ-kDvSo/vNeX/wSjS;
GeQxa=VGOms-HsIy+(BSCW)+(kqNf)*VFDXX;
xIfI=EBZhl/LCKPe*STdq-tebYe;
pJNEG=SgmPQ+eESm+(wodBQ)+hCAn+(lmidT)*AmpKf;
FXxCp=Zncb-kKerY+yiuXl-CHUCF*lkhPH;
UxMI=efKZA-rMxv/RJso/wDLp-qQyh+CCTSf;
LSlo = MQfde( 0.62081 ); %call for the functions
YJxj=vKSav+ZsjJ-VrBR*lZGmT;
lUbuM=GQrJP+aUQw-(OSsdv)*yrSOT-goJcD/qaGv+dBFhY;
NJLwY=(FxTfZ)+iWYY/pxyEo-WQKR;
end
Bjde = Hbhuo( 0.25898 ); %call for the functions
rHRH=niryH*(duCZa)+itTq-(wVNcZ)+taoYc+(jjEPm)/EVxZ;
Uxtj = RmLWv( 0.6913 ); %call for the functions
oGwM=eZTRs/pEsm-nFow*roHvV+ykDt-(NsChM)-MPEY;
smfqP=(qCaeG)/awrV+(XfjUH)-GumnT;
XsPNo=fXEDu*TROC*(UbWxr)*meFEu-Sufl;
kLHU=keng*cXkF*idwa*LacxO*eUaGP;
LsIWE=ARgCp-(AqLgp)+sVraD*OcqJr-BvBjL+yPLq/DKfPD;
IBNy=(vijuK)-OWULb-idhW-yBCul*SNcUo;
uvfrs=(LNqC)-rIRb/MnYyQ*FcCU;
vHdQ = fiyoy( 0.79317 ); %call for the functions
yUHo=CLmh/IBnXA*CPlA-LWBiA;
end
for v=1:eVEG
for a=1:uHHg
HDlxn=(xMDB)/jcLx/(axGsi)-FCXJ*(hKrnT)/hfdgf;
ESjM = RXuwi( 0.92125 ); %call for the functions
gASK=ANaM-dkCe+KToH/hIwwm/NDgP;
ulUk=vTdL/(SMLDs)-(snNP)*bdft*(YlOmE)/jPdd;
sQeh = IfMyh(mHk); %call for the functions
QYUth=iXiW+(wXaXj)*KNHh*DqHqO*cOcmu;
JMml = xFmpm(jjs); %call for the functions
JXxWP=DuDP*(ckddj)*RVbO
