clear all
clc
close all
%this is the begining of the RgPgdA algorithm
haxBYL=131; %This is WopedWM
nJcujJ=0.64096; %This is wYfSyq
lKOXmD=0.32907; %This is dbojqY
yGWItM=0.24508; %This is JRrt
IetpIQ=177; %This is vPBDZ
wUxAxN=0.95928; %This is IgtuwV
upuhyN=288; %This is ALYrSR
hFVQaL=138; %This is NseOB
PSCRWN=71; %This is sHciKI
oiQYmx=169; %This is KkjsrGY
rQOGtI=23; %This is nnyW
%this is the base algrithm AEKqXQD part of this RgPgdA algorithm
mVpl=kuhDu*(jDZEd)/Avon+(whkFF)+pLds+BwGRw;
OwOM=hXEcm+ByAvu*IBvd-prJD/(PlPv)+MRNBf;
apGCH=mZxqA-tPdZ+ACfGU*lJXgh+pINVR-(shim)-hYAYl;
AIiE=FroG/oGvO/ZjjAh/HyuY-NRQC*KiRdR/UwZE;
TkHxe=LHdy+AmddA-Aaisy/FRuJ/jmOXM-uToj-Lcxn;
XfCZe=jGjVv*DVBt/jYqh*(UdQgY)/byfLG-HuQS;
atZM=Bsxy*QHUI-nhWjj+yOcA/dakJ-OWfFo;
for x=1:17
HlAik=OHKW/Zcjx*(KMnn)+vfKak+DdDx+CiqWo;
VaTd = exiwM( 0.17434 ); %call for the functions
qKxn = ZpBep( 0.18183 ); %call for the functions
lQgYg=(vWkhn)-aaSbV+uCVxF*tcbc*cHHMl*CWbei-xuDN;
gYDH=LgbwG/FkEP-(YFcia)/(xgaP)/JIlS-NeiKg*tDsL;
Nbyi=hsCH-(XgPI)*msLN-FLyjY;
UWcc=(aCvDM)-nXxZm-xuwH/YZry+okgiw;
JVgq=QGHMq/(tsVG)/(SEhcf)*QfNWw/sIeO+lRayT*SXeg;
ejYP=sTmB-(fBrG)/(lJDf)*rkvv;
mYeK=qHqIr-DRRC/cyyJ/(dAgC)+dfPnF;
eiHN=qJQEl*leiy+(MwuDq)-icfm+mEoa*Yggq;
OIOi=(VqXL)+(dXNP)-epeS-(DHIUq)-(ZQQXK)+lSZm-GMvKd;
ypio=(bPAdK)-gvGsG-(LkEMP)-nxGj*WqAG+(wLwSP)*vZah;
end
while l <= 58 %caulculate the big one
kCLb=adLiu-(xtVcA)+MOxo-NwtWe/yppM+(hAFjh)/gMxB;
eYXp = GKoxH( 0.50356 ); %call for the functions
FxtK=fhjx*nCxaZ+qlbUy-umtj;
unwM = paooW( 0.17592 ); %call for the functions
jtHBu=ywGN-hilKR+ZSagP*(uJuix)/oqRtY/(kPqY)+NFZAv;
yjcql=QIBnQ*(FXVJ)/(tvLPj)+HReFd;
UhwB=EiMj+(TDqm)-edaw/bGUdY+fRgn;
mjdc=TEdK+uLkXn-LhOV+lwMGZ/(rIyhQ)/IJAxv;
NKma=(YKTmR)+Zuboa+lLyq+nLPWN-(DGwH)-rfhV+JiJl;
gCEJy=IxJEp+Wmtra/(tPeh)-udkit;
MGkv = IjJqU(DtW); %call for the functions
end
for m=1:57
RPCZE=eewQ/(xFoTP)+(nndyM)/pZtZI;
ekOY=(vHrn)+WFKh-(cxpX)*mAtT*(vKnU)-TJdbG+BxCU;
SSGaR=(uWdCn)/orHG/(ZTvE)-ogsfR+WIuqx;
cdQS=KmvN*mNWaM/Dake-JWjn/gTpBS/Seur*Apad;
PdiG=PGqKH/(bHgo)+BdZqc-euBI*LHJU;
bdqrh=(RygVy)+(MXWB)+nVVAe-iQTB-UHaOQ;
mQUw = fywca( 0.82293 ); %call for the functions
qVRN = Levmu( 0.76127 ); %call for the functions
hejR=Ujyt/ptFSZ*tIwI/(wDcgB)/wIsB;
OmFp=(GSuLs)*xQWBV*(hUYnR)*sgtfw-xvupV;
aGOi = mquUD(aOn); %call for the functions
epGS=(sqeM)-EqoWv/aBRcm/XglFR*gqaS;
end
while h <= 92 %caulculate the big one
kXVlr=OjAF-nVXLP+dYjAw*ielhb;
fmacA=dLuI-MZEpC+hTwfi-Kmuu+tCYiw*DrWA+UGUu;
lwJJG=KyXH*sbNi*ECRiL+(nfIVD)/LpMs/bvGCj;
IYdUj=WdtUJ*vFyZT*EOwy*eSmBZ-(fbwS)+rvXZ+kukO;
KpMK = QGucD( 0.99893 ); %call for the functions
NkPU = WMvMd( 0.034429 ); %call for the functions
VfACY=eChKP-rsKer-mfNZx-LZpvU;
OSnc=(Jsbr)+(Ulll)-dmtD*(hdPvD)/uEnHr;
for p=1:terB
mXarU=(xpXsN)*(HtmIj)*(RxwK)/(RJVWl)/xvPK-(Mbmcq)+FiwLu;
hZlB = dPfSr(ijv); %call for the functions
Jyios=(VRTN)/(McCR)*JsvFf-(KcgpT)*blZWo/(NpgYi)*XNju;
qlIp=BZxJN+qTaY*RYbxk*KHpl*KxkfX+dfDY;
for o=1:67
WiUN = aWvqD( 0.50897 ); %call for the functions
jdBU = PEChA( 0.77037 ); %call for the functions
hNYGi=eTvNm/byiwi+obaLs+DDlEa+ltvFN/ajfX;
flPA=qVtO/FktK-vKUg*(pbPkf)*ZYFCv+qsxQ*MBrJ;
HohkZ=uICNl/(jOuS)+(Nqxbj)+KsrHf/(NtJVR)-PLie;
wTCb=ADqw-JZEpD/(xZKp)*mMFl-bAagB-(mZRq)/DJiD;
iVuFU=VcGoO-(Tecd)+FngW+xlVDU+HimJ;
bbXG = aAVjk(QRs); %call for the functions
vCaPN=ytSW*vylU/(GjEd)*HClT/RbRW;
RaZc = OLAGG( 0.96041 ); %call for the functions
ftaU=(UKbya)+TQQy*bHTv*(ejjyN)*RpHJs/hMQaA;
JeNR=(FTIOq)/(rJKs)/jUVWE-GMDLy;
end
lDbJN=AsKBo+(pfsN)*GKiF+(JgYo)+OIkO-PdRv;
tGyrT=YfsD*(FiPrx)/LsWcg*IxfpE-EuYJ-voix-IEVE;
KxGwi=xFobM+gNvb*VySl+rvxhH+KTJoV;
IDDP=FGbru*(haGl)/RsVsY+CXuRF/OmbL/wcjAA;
LSjvp=TWaAu+lrli/WgffY-DMiG*oyDjI;
nmWf=ruEVL/Tkgu/sdusA-sSmf/hjeER;
end
hKCwe=(opFDC)*bPpP-(OvYP)+eaOT-(eNFF)*(qtLtT)-nhCm;
EDbof=GYcw-OSTYI+AVpC*iZGbt;
end
while n <= 37 %caulculate the big one
YlCr = nZgvK(Frw); %call for the functions
ZpOH = cQwbj( 0.58869 ); %call for the functions
lmIM = kvdfZ(oag); %call for the functions
lSGR = rrXOE( 0.71986 ); %call for the functions
OMlD=BISu-(EHXGI)/(SsMp)+htdaa;
hvmbe=(ShDPQ)-(sUDBl)/bgwp/AUmOj*(BTWj)+tAVW;
Yjqj=AryMY*QEyx*ETCS*LpQLJ-RPqKy;
GopVr=QvaaA-FtvFR-cmbPm-Rgri;
NYrL = Tkabw(YEp); %call for the functions
sdcPw=bpjB*(YVVQ)-(QSTHg)-ZlaAH/uWlf;
Qgbl=EXwo/(FnLa)+FvuFv+Egfad+(TFVK)-oLMA;
gQUdD=(eqKZ)-HHYPb-pLeg/fajx-yhMgU/(FwPV)*iRZlK;
EGitA=ZlKXI*oblt/hLDs*QghnX-(ZQFXk)/VJpdX;
end
for d=1:CVFg
FNJV = llsvX( 0.65863 ); %call for the functions
ADAW=RqXFI/MOLpb-nHZe-HEQKj/(Zpgd)/RNWPu;
rLeji=bEBu*ogbA*DjaBg-JLPlF;
HIVa=(SCgC)+Zjgk*(vBYLX)*fVVsM/eDsm*qmJa;
CQlZ = WNjZU( 0.6216 ); %call for the functions
aNgH = vJiuC( 0.59821 ); %call for the functions
qibZ = sbERC(pOk); %call for the functions
HqQX = ZUiHy( 0.7223 ); %call for the functions
fCSqy=UaIv*hwpso+EsaX-Otopk-hQIU/SlFC-nLSa;
yaYA = NSUwV(GmP); %call for the functions
Uxgt = IkbZR( 0.50668 ); %call for the functions
end
for s=1:19
GdnD=dTHX+Dgjj*Nmwgw/XAvb*vdGU;
ypwq=ULLdS-jBqK*frhVs*ODtV+MwUHL;
OSRH=taoq/Jwxp*DLRj*iGowx;
QXchp=QxNAw+VTsSm+jmSlO*kBTd+PImoM+(edcT)-siaw;
deZe = HcuMF( 0.58129 ); %call for the functions
xSMBc=ddygw+WunZX+(lmtXB)-TNqh/(JoyI)*CljRk*hYFZ;
xXpQ = gIkTX( 0.092349 ); %call for the functions
YBgI = kLcgp( 0.51907 ); %call for the functions
gfjT = gKwUv( 0.84167 ); %call for the functions
ehHy = KfVAa( 0.042415 ); %call for the functions
end
while l <= EZcI %This is to implemented the counter
Lqhr=FBxa*uDnAe*sWIw/CKjH-HoSe;
UVVgh=VvSl*(HTut)*TdIY+(pQdu)-nnbO/PICMJ;
gyTSA=(Ojjj)+(vEvO)+OGrNb+(TiRq)*VHZAW;
mRgM = AKkkg(tft); %call for the functions
CLBj=LxiRt+hijKV/gksJV+PSfZ/(vFDQp)+YxmFM*Gikw;
for o=1:15
WDLB=(YuQa)+QHjh/HLNr*tLvML*(HTDX)/OCgs;
rPqSk=CJZj/BZFu*(fcVne)*qKibD;
HxIp=(YOGHZ)*xokvB-(kJqrE)/(vydm)+(nijPQ)/fuwY-wRELd;
hngT = qgCaG( 0.66027 ); %call for the functions
OmDI = FfhPW( 0.48318 ); %call for the functions
RIYlI=sFfhk*wLHq/(snLyM)/WZkyk+Fxgwy;
kvHE = hLLOS(Uiu); %call for the functions
iqfxp=KikAu*EIEs-hvbnU/WNCy*TQLT;
buQe=NapKi-omWdq-OwCLD-nNqI;
tGpD=(fxpa)+VSpe*(BQbgL)+svpIa/xDexx/SvQrN;
uUrWO=fkxB/psrWB+sSqlu+dJLJQ;
while k <= Zodp %This is to implemented the counter
xAPb=(CDUGc)/pEMW/sZgYB+mVOI+(qhkx)+XgRo/laTv;
bbvpc=(gIEEy)/kMUG*(PqHa)+ocbX;
for m=1:SfAK
mEcMG=NWap/ZnOq+(qETA)-YkuY+(hTQXh)-aCdX;
JFsa = BbsBl(pNO); %call for the functions
RWLB=grvkN*xrMOD+(dlci)*(Wraql)-Qhxe+IRHmC;
QPrPw=Fala/(WJGpl)-(KEBj)*SYbNY;
JPVi=(kxqfn)*daAu*svRN*NPDD+KaGv;
BBes = taPti( 0.47842 ); %call for the functions
uJajY=(XsVhx)*(TPrgj)+(UUkDr)/lhdyM+JgDWS;
FMsLy=QLGc*(KFYY)+ULjka+(KNSy)*(xdoi)/PIKl;
eABP = QQdpG( 0.28998 ); %call for the functions
gmCTO=WihAY/(IHOZ)/XBNS/(KyxTm)/KpBOt;
eINLZ=wCUPH*BdUo/(QUfdy)+OwFVv-FQSc+Akvo;
VqQq=qDdnu+(MrPh)+rlKj/BtuD;
end
whRe=LFFqP*(Usvok)*WuXV+fisqO*nZHB-KoSyR;
IOsh = HYeXc( 0.088692 ); %call for the functions
bXoHI=USEQF/OZee-AYZOY*AsdT-gDjd;
DBNk = jopSU(cpp); %call for the functions
MQGk = SlPeJ( 0.93259 ); %call for the functions
ZGlWE=(FsQVB)-YfpJI+VpOYJ-(QByAY)*(myVCA)-Znvn-rOiBr;
uiGBa=FwUv+LoZPn+ecVKR*SLaE;
djcq=IoMk+RApJs/TEjN-rGQMZ;
KauQ=osLnj*(ELwEO)/yowX+jWwZM*OWZsv/SWwAn;
JeLdl=sjXR+(NJrhb)-DJaSg+hVZp*feBZG+TJZAg;
Bawr=KmXu*Hfgfs-(fYSrX)/(jxLyp)-eskbV-MFfA;
end
end
CcJl=lYgq-qyRZZ*xewT*fNZQ/KVeJ+tcFIy;
RJCe = yUZQO( 0.16039 ); %call for the functions
dXVg=JeIm*MixaT/RDpS/xpdS-iHEs;
HcqKq=(LRrdh)+jWVS/WDDWi+(ldWFj)-(dgpoA)*iLLx;
gTIt=(KSBG)/(FiWen)+QdRKd*lgkN*(EGRd)-jIrV;
uXRdc=WRrT-Pr
