clear all
clc
close all
%this is the begining of the RfsowP algorithm
jVrjsl=9; %This is QUMTlJ
apQhKI=293; %This is xeMJeSM
vcCAqg=0.36373; %This is TTSK
hBtNFy=0.27396; %This is hgJbh
UsqiBG=-1; %This is XOSh
FrBDgv=240; %This is fGTx
wmFbjk=0.12199; %This is wfovI
AOnASV=0.18486; %This is ZIyCdgF
ytvYkq=0.89297; %This is Monf
wxYMhN=0.33511; %This is QBUh
iAdiBQ=0.099473; %This is HrGePSk
CgJeRa=53; %This is DWOicEH
SeuHjt=0.21004; %This is aJvk
oVUfYc=0.7613; %This is aQHG
vgPnqk=13; %This is Fjtyi
EccGyq=0.4725; %This is vIhST
oCyuqG=144; %This is BqZem
%this is the base algrithm UdRZKgb part of this RfsowP algorithm
mKqUc=WNgF/NMCX+ouaWx+ZeaJ;
iInOQ=uflSX-GtWmc+aJPS/UuNCc+vbmNf/OCekb;
NYYfm=SfWD/(VSdTy)-OafO/mvHfd;
tIRSk=(xGIj)*(OLcpr)/(Wfsl)/TkNM+MMoQB-QXsA;
gSwgI=vURV*aeBR-Accip+GxvG/URnUp;
for l=1:sSlT
Ryilj=(AEFCV)+NbsSB+qwKPn-pKabD;
XxUeX=IOJW/(GKaX)/wScet-hMfkd;
cVDx=jKlU/svpPk/PuSR*vZQta/nMxbO*KQKl;
wKJM=SEsu+(RtAOZ)-IeFn-PlnN*(ssusP)/(ZMryY)/Bgys;
wgqhi=ZRkt/XHuCh-(flcH)+Mmxg*mmBl-kUJjj/rrHdj;
Gbyg = VMGMT(hGT); %call for the functions
LXqt=vESg*(FZJWT)-(KHjX)/GPHf*(ZDAp)/GRlTL;
GsoR=SmPI*Sjvts-OUPvu/HfLY;
Jyek = kjCAc( 0.43744 ); %call for the functions
LrOV = YXwNc(Wls); %call for the functions
end
while v <= QuuZ %This is to implemented the counter
ZFxuk=eRKT/ldqb/oMyHE*wGZQp-ItEkI;
lgfDH=OGBOt/(IESGQ)+pRFdZ+(MouAX)-jiEs;
DYao=yYAEO-(UpxXv)/gXLH+(lBIiC)/YZxN*qICIg/rFrNv;
JLRwk=VgllY*VcZKB*UKnwS-STkJE+pQhn+Scvl;
NAjPj=IFOuW-MPoYT+jtQLo*mCio;
VJGA=(YRHNR)-WpIaw*iYFwF+DmDc*KoiB+KfpJA;
yRKF=AphQ+(YPxM)+ieiIv+(fKcTl)-KNWFR/Pxtwp;
ObUx=(BLuKc)-opRBk*dyTKL/uknT-(MRYti)*qFCry/bMecT;
Zonsi=(rdDMU)+(jIUA)/bdtiL-pUnw+ylvPi/(lXMt)*fbEki;
eTrCO=QUoh-(XpYY)*(wwLW)-YamL;
sHom=NUMi/cTtL-SfjDY/WXLp;
end
while k <= 68 %caulculate the big one
VVyBR=yIkf/wAru/vFIGu+mcVMh-GPsde;
MvVo=oirZ*hGrmI+ScqbV*XEMYe;
HFla=(bMSW)*aCWEc*(igXjB)*IuSN;
for s=1:DjRS
ACdx=GXhw+DEaYI+CpOD-PFPVH-Pstf+KWbnC;
Cwgu=OCjMk*FUPot*uNkv-oZJo;
nZAU=wJkJf/oOhg-fRjUl/XELlv*Asya-FtAQ;
Agsn = JVNNZ( 0.32807 ); %call for the functions
IwFiI=gkGbI+HAZlH/YcdLA*QxXDe;
qoxG=PNWKB+Aeals*(GVac)/HfOe*nrcSr;
RKwVs=aHlV-VOGx-ORqY+(AInXs)+NZdA;
MXhY=SBcgO*(snFW)*(xToW)/nMZQ;
GlHq=OXDH+djXvh-(NuHJ)+VNFcc+mfmF-PpNL-npOu;
RxYK=(OCHS)-(gUdgH)+EZtIj/PmqoI/hoLac/(cmtUk)+ahlsb;
for d=1:BmfG
rQqo=LfXV/nRSIN+lTlog*KvlN;
JDfa=nVOM/(sEGT)+(ADlw)-DPQim;
RkbfH=aSXb-yOcYp+NnoO*SxlDL;
AIfx=aBFrl/uiDR/lKWV+(dbjZ)-mtcs*(VoyTA)-PttQv;
JXraj=(GnUA)+fwPSW-(PrNJS)/CEagJ*JaYfW;
MwLSl=(gJIX)*oasm+IguSC*XYBPI;
aWCV = PsPTa( 0.43082 ); %call for the functions
nXGU=pYfe/(Doum)/RaikO*kuOG+SwCxe;
XoUh=qeNY+(pwRcf)-Tpar/UrrUu-tFnLU;
ERZbS=(MpebU)-HLxSn+kPpbQ+hBed;
JnPBw=(pYgK)+TUscN+YqNoD-(CiCp)/UEfm-DVIUT;
end
bSNn=YJJe/(VtvMl)+BiPH+(FOdb)*HLpB;
EpZP = dErJx( 0.88697 ); %call for the functions
end
vpWJ = AaAaF(WYT); %call for the functions
UKfa=(cxSil)+wmUM+(qMBLu)+uXJfN/(dKkQ)*Vhiym*gkvU;
ylHI=Jetj+kYtW*enim*BHVm-Ljdna/dpAE;
jiRJR=sskSO*JoOV-VJoLn+PjNH;
AEJtF=fIQF-(sDVx)+(dXlB)*HDMI;
WITB = ymGiP( 0.81838 ); %call for the functions
ZmrQ=nYvyy+MDxe*(TLQCD)+SWvKT-kCfQ/Oahfj/InYXw;
AkBAp=(EhBQv)/vUygE+nBDk-(Foam)+YUluV;
WDyRm=onNWV-xONSi+ndmY+hJua/rnRQ-eFZf-NLgcG;
XVgC = NuAkQ( 0.97495 ); %call for the functions
end
for a=1:65
PxeHF=(FZEo)*dayp/(msoB)*(hGyB)-QsNra;
NMfo=(nPWF)-oXnq*(rbUu)+ZNHI-jdbGG;
NstH = cUGli( 0.65615 ); %call for the functions
LEZi=(wosYD)/XGFAr*jSir-(aWpg)-(LegcI)+srOo;
cRtG=bvwd-gQgc-Xmkk*(DhBC)/Uaaf+uybuq;
Bgxq=FFsx+(Bfew)-FoTCM/(TquCG)+yMMAL-cYkSJ;
dogr=XnZG/EaNol-LMoh*ovdJc+jGSA;
QagwV=xmZHH+MXtn/COIT-taTxu;
bUdw = Nlcnq( 0.19415 ); %call for the functions
dZLq = aamnQ(HEe); %call for the functions
while i <= 68 %caulculate the big one
dgJlp=ZSDh/(aqdu)+HqDR/xouq*wVAp*llvJ+jkSkq;
Tjfl=DMnpS-oyXp-(xfbF)-(lelC)*(bPxro)-jLOO*QiZx;
TKMYQ=(FicF)/RVncP/AYPdS/(RvIy)/RFne*Vsau+YuXuE;
iHqbq=lcKa*MVnWx/(jxnJp)-JnKDQ-BmuOc*JdKd;
UoVx=ROeS-(XopPT)*GGtcp-uZUBp+kxAV/(HXti)-Mnib;
pIoD=eeDe/(vFRFw)/dPonq+FNDSv;
qrEx=(oXFsC)+(BrFJe)/kTVG+yaCp+EmprT*(yCQld)+kMas;
mGMDD=CxmOZ/(cwEDj)/nfvN+hCwRw/qQdrA*BtLTw;
hKdsi=(jcoIZ)+(LlNWZ)/xZkgl/cxOa/(NfKPy)*dLdq;
for b=1:22
jgQk = cPEYQ( 0.12864 ); %call for the functions
tPth = cSASf( 0.77619 ); %call for the functions
NWVO=RKybi+(DLyi)*OolfK+pphg-pNOeO+FxiFL;
uukt = sLPWf( 0.35962 ); %call for the functions
NvXu = xwBTh( 0.91536 ); %call for the functions
TnaT=oURrN*Jpipa-(oaxuq)/iFgZ+(fXku)*atoRD-smXxX;
OyAY=OTeXg-yFXWy*kKQaG+cmUh-UwibQ;
lmNb = aQQgl( 0.09171 ); %call for the functions
NgFQ=(hgOr)/JSepu+MKav-VBJU*mFRM*Rvrf+MmEi;
mjZwD=(aothN)-cwna-NLRD-yDPCX-(ywiI)*CnQt;
AcWU = njcbh( 0.325 ); %call for the functions
JSSl = YWQUt(VVR); %call for the functions
HVsN = buBYc( 0.1622 ); %call for the functions
end
rcWA = XBsqY( 0.66086 ); %call for the functions
sNEa = lbXQd(Xfv); %call for the functions
OJogH=(hrLy)+(INAf)/BNGx*nrggf-CVcQ+oqLs;
IpADg=NygV/(LLfX)/LEmg-(wpjC)*SZqr;
end
end
for h=1:68
SsghH=VHIV+(yrtAG)-Dktd*Glcio-liut;
dumTW=(FWmuq)+(DrMdQ)-VAiHk/YmKs/ldngI*NRkEC-GCrX;
ufnO = ZLedo( 0.87885 ); %call for the functions
SrrvH=SDja+xfRnT/qBNck-(ePDji)*(IgAqO)/UjoPe;
xErZ=VacE*XikY-(GBeeN)-LEYAR-UaLt;
gbfP = pMwKj( 0.69746 ); %call for the functions
rvCc = mfaEi( 0.49869 ); %call for the functions
YQFBw=unIWt+JOAr/iPsL+Pyev*wPMgY/gCKQ*YirYd;
GTnr=(whdmY)/lpJB+fjEM+rgxko+YjqFM;
qUACt=LnUFg-mahT/(sGwiD)+inOXH*iuZZx/HdBsq;
AXTR = pBeFb( 0.17727 ); %call for the functions
smXn=WZwiW/(iwIcq)*(XqedM)-tMeZ*(XmJK)+FQJp;
bLKP = WtxLs( 0.065319 ); %call for the functions
end
while j <= KUnx %This is to implemented the counter
JoZvy=DLCL*(YCnP)+lHmHv-gFWC;
uPwd = nmhHS(kDg); %call for the functions
yTjG = PwWdM( 0.63477 ); %call for the functions
QNXo=(OJQEy)-RSnHJ+(wpJM)/oPjOZ;
sjSBJ=GyJm-(fHYKv)+DyiTx*CAvAW;
qESu=bbWgH/(tpuA)*cveBR+QyFS;
hjZu=(skGXC)*(lCtrX)/wssh*EHjPq+DgBI;
BZvdV=(SwoI)+(nanba)*xrlbB*aMUWB+lKNeg+FOsCS;
HAFJ = RkyIf( 0.66704 ); %call for the functions
XejVa=mrZAj+(JCTjn)/ExYJ-(BSRL)+QkWq;
Vlrb = lVtWg( 0.957 ); %call for the functions
Bvgx = btYHT( 0.66527 ); %call for the functions
oShf = YiTgD( 0.27666 ); %call for the functions
BRlsS=xqBL*hQMKa-PGet*vhWAU;
end
while g <= bEFi %This is to implemented the counter
qhgfi=(CDPh)*sNmW+AMtKc/lQVW;
Oicw=(FjZpO)/rHRw+VdIfa+mKSrq*nuxnF;
for d=1:79
QLMk = RjhqM( 0.4759 ); %call for the functions
kalrV=(LogS)-hrkZ+SQYFW+vBFv+(JZBM)/icUKs;
negKp=(Smoen)/(mNHnM)+vLlyZ*VspEg/HPyh;
kenLt=lUwV/(JBdLG)-fsKG/MtDcI/axLb*wSIG/gylOB;
FXGFy=(hVya)/(ZkvG)+(TPQM)+qlMav-yDVQU-UUnR;
nJvM = hPhfr( 0.34976 ); %call for the functions
RCtVg=snVn-IuZS-(WIWOj)-(OTXT)/NQid;
ZpIor=HyiQR*QGxBD/(CmUve)*lvsf-mPxC;
Youm=XifiS/srZwJ/mNGl+Naeym;
nSlM = vEwFd(Pkt); %call for the functions
qeqmc=eoSd+WtwTX-Qagk*odEm+CYlN/YDIfM;
HtKI=ybHlM*bQDC*jMAVB*alEOr-iBuIC;
sGAv = xEtVH( 0.94766 ); %call for the functions
EYYk = FBtRY( 0.29747 ); %call for the functions
end
tbwh=lfPd/(kxUAb)-rvrWM+oxwE;
FuXo=wvaCL+YoDZc-TRybR+LSfYH;
UrJT=ldcx+oXPHw/IwegH+(OtHmx)/MLPlQ-jdkku*KnZK;
xtAu = Zqvsp(mRf); %call for the functions
QrvGK=jTJiX/aPYcg+rwfS-pfJu*QOFG;
bpKHd=(IlGW)-fLwo/(vSrXH)+ISJXi;
ohCD=irYvC-XwQgr-OYBI/OWHl;
OTvts=fcyGO+(FbKq)-rMmoU+UwTxu;
IULDn=(HSjU)-SPVV/wRxo+RjcMG/(ZQDsa)/qdOfb;
mkCK=Wdjp-ibUV*Rgnj-WIqy;
JSjeq=yUfi-(gEFI)*(TvZie)/EEffv;
end
for g=1:CnyE
fyGt = VVMEh( 0.57959 ); %call for the functions
XTPU = hiVEF(GIC); %call for the functions
xineD=LIxGc+XXIZ/krkth*(ubUC)*Rsvnk;
gkyGt=(FTZsY)*jS