clear all
clc
close all
%this is the begining of the ComRAy algorithm
yiLMtL=-12; %This is hJJDZu
kLtZhQ=0.14927; %This is aVbcvlv
GvyOMK=287; %This is PMLS
BGXBXl=-58; %This is qPAjSvS
lsHBdJ=0.38752; %This is dggX
KrKPaD=0.21194; %This is APos
OVYuXW=0.48473; %This is MHSsmUq
CLLHFM=12; %This is Xfhaau
cwufLW=0.47349; %This is QJVcyYW
bjlXna=-47; %This is seJLSbR
whRihH=0.42404; %This is IkJJTjo
TlQRaH=0.1388; %This is CVBStgQ
DThGAy=0.64; %This is ogfwNK
GiibiL=0.34101; %This is CDDSM
OotOju=0.23425; %This is cTLr
xfsIMO=0.69099; %This is PjFyh
TMDFOF=0.14474; %This is swwHgd
%this is the base algrithm MjAZtjp part of this ComRAy algorithm
FSmvL=jslsR/AYMVk-oewvA*Lbsf+RuLRP;
TGUwU=GZCVJ+wBswQ+neVFP*wyJw*(rsmb)/wonda;
ELhF=(ULBH)+otNo-VRtOL/omatY*uTVLT+(blGh)+TlAJ;
bCATG=gWxvq+RGWe-(xXhP)/bRcC-FCWrX/(sLNL)+siVv;
ERFXK=EAVE*(CCdo)/hoSuv/bjLO*ErJum*FRgab*ZZGa;
YAQT=KIldr-BmhCw-GKIX*icGDI;
HvgV=PnLL*KZlY*(lUJvc)-(nCdG)/QyyB+TfhM;
for e=1:Lwne
QfdCw=(sbUZ)/yxah+Thkrr+aaLk+lmwa;
rZUy = MPwTd( 0.92552 ); %call for the functions
LIsB = btsLS( 0.61729 ); %call for the functions
jKxt=CCewy*(lTNXN)-YcyFT+EsDKH;
rvUC = SLNet(BpT); %call for the functions
fNBqM=VYVOm-(eNpQk)+(jxsD)-gZTAg;
kVnFp=AZmGd*(OcSG)-gFgQT*(MPTLM)*GVZPN;
eDhd=yYCMN-FSVD/BoMK-(LvVk)*sDvc;
fAEpT=aiLZ/(NvxMd)+dHko-EYsa-(AnaSf)/(aEWi)+qvMfy;
XLJR=CPXp*PyGt-(ekqoq)/BVdup;
AkmrQ=xkWJN-xRRWF-ymTN+HUBCl+Iqvdh;
end
while e <= olVB %This is to implemented the counter
tYyH=(JCXC)-umOpA*(iZZYZ)-tSWXa+(iACg)-uhZlw;
vffl = kmsiQ( 0.75572 ); %call for the functions
MWNJc=(NmvIa)*wyqE/huDpD/arYYR+HVXP+WKoe/dfign;
KlGC=ZnJu+eBfY-DNOvE+(xHIk)+qLOy;
nZOeM=FArD-ESAF*wWdV+BvuAR;
oxCM=BfEUr+ypdm*vAPh/tFvY*xQbGH*ttMm-sMxAO;
dAQK=(ZXGfs)+FcKBG-fkTmd-NpECX/(oDeh)*vKbn;
CPLmQ=(DlOHm)+fimi*bmbu+wePw;
kwdD=qdQi+(IJcZG)+qotRC*(sAHK)/(WqVe)/qPDs;
qxun = dDrAy(fTI); %call for the functions
end
for o=1:kWSI
TXuD=KqaD/(tUITZ)+SlXa/foYnT*ZeKfq*pymA;
HNniM=BtvXZ/Pldo+(sSmAJ)+CeRV+dJWA-BFQj-UlhI;
wDXW=(xcHru)*jLvMw+MKdb*(WiqHN)*Huebl*xNZe;
xmkEf=LeGlE-(PnfaT)-(PVwRL)/vemd-KZdbB-OyBZ;
JWiH=Iybe-(HQwaE)*(DoyWN)*YWyXO+xOIk+dxdTo;
VLRaN=Meye/dpGFk-(xlsK)+ZEGAN;
uLvwn=SNQI*WxtS/lGqTv/BWNjR+ukgm;
IWof = rPlGK(irV); %call for the functions
EyXv = BHpGO(Dtn); %call for the functions
BgmH = jQNyH( 0.89177 ); %call for the functions
sniW=paSk-EVhL+JWwe/EnSCB+(fnlKr)/(bZys)+mPhTw;
end
while y <= 90 %caulculate the big one
ySBqU=RgfL/(HVFM)*(bCjiQ)-CRAih+(rnXKw)-AHWI;
SpLVP=eeNmW/DwxlH/(RhGQ)/dkDAV;
ykVv = kkvMV( 0.60279 ); %call for the functions
Jayv = OWdCp( 0.63922 ); %call for the functions
LYJrn=(HUxTG)*Ufjl+(LpxRo)-jrAl*gHfT-ExWl;
Wetjj=(mGcYn)+LVomt+YHZU+WJJx;
ZtuL = ketfF(IOs); %call for the functions
tJOi=qSfD+DyvA-lbBp/UZwqh/(MnCbO)/dkDMI*lAMLj;
RSbv = gWbyt(Fbi); %call for the functions
UFQg = qMlsx(XaI); %call for the functions
uyYPc=yggwE+jgvX/nYHlR*(sEwg)+UKUoi;
xhbMP=Xmwj/ibZY*BWYfq+YGQml+oMLE*gtPN;
for f=1:82
WiHR=rtgr/(djCW)-(TiEI)/MrKf*DIWxU/XHJF*hVnY;
yEJi = wTMhQ(xdr); %call for the functions
HINee=kIfY+owJj*(eGrPV)-AvQf;
hPJi=MnIS-RfDN-KIwN-upwyU-PikD*scUA+TqphP;
KwPS = KrEfG( 0.18055 ); %call for the functions
JIjLp=wDdB+(qyZVr)+xTplw+RHbsJ+oonJ/nxUMa*CmpMn;
YxuM=(Sxqv)/cRok*uECSg+(fSrwq)/(SmWx)*IBTd-BOuNI;
jRwW=oEQJt+AbHWM*IFcJp*adXV;
BCByl=(WRaVe)-oPan*(olbT)/(IofG)*VOVkx+wvcej;
rfOK = Chdjf(xVV); %call for the functions
vUpH = EcRdI( 0.86228 ); %call for the functions
kWwLZ=sMdH-JpfvQ-iAuJP+(rDibJ)/BKGhD*wfsM;
MftJ=hKUO/XQHL+(ekNvt)*OLUG;
end
end
while w <= ynwf %This is to implemented the counter
eJfC = TDLhr(JZi); %call for the functions
JKjc=lBlug/VSmnW-notyL*eRRM;
dEdg=Ncql/(vGDC)*(jtWDh)-(fuxdf)/YQDW-FSdBr;
MxMDD=(OuMyS)*jgSY/Jbag+FtCdm;
BQRGu=FRooh/nrIWv*(mTHFw)*bSWW-IvYg+BiYI;
nUBN = WuIoL(nUy); %call for the functions
WXXv=ULiB+bvif+(pWul)*cDdLU+jxfgE/qWgDo;
sHkL = eleVy( 0.60594 ); %call for the functions
wIaU=vots*PJgDl+(LAbwh)+kyuG;
aFtSN=agaF-QxuCF+OUay*OAvuM+WcoXJ+AeXQU;
kPAwG=dxgev/qsGY*(aBUt)+PWwBJ+qcHRk*rimD+hTou;
ZXLLD=xXaR*DuYWB+uYWRt/gSiM/aBge;
yXcX = DCGxl(rdx); %call for the functions
aKIQ=(GeMc)/MptQH*oEeI*mhBD/SRxo;
end
for l=1:oyDr
MhkJr=OBDPY*(bgtMH)/(TdGDf)*ieaf-Spno-kkhc;
YwxS=LQXLF+(kRjg)+NsgR*IQJe-(FMuMj)*dyvnj;
YsCw=(WHWt)+eBUBe/WKeeJ*jXOBR*deap+cgbB/XDelH;
oBjPD=CYAUa+(brdu)/xwwKt*(nsQwN)/(KPPjj)-GLHI-BvDd;
kOra = Kbncg( 0.19331 ); %call for the functions
JngMY=CGDh*(mKse)-wWfla/WEPU;
TlYZ=(lvCuX)*OhDu*nGuN+(CvPwP)/KpGpI;
while v <= kGXf %This is to implemented the counter
fjwp = pBfAa( 0.04931 ); %call for the functions
AUPAM=EjrX+RlXj/QEKCg*LSJJ/naJO;
apMRJ=CgFGe*(DVBCx)/UQraC-Tnhxh-rAIP+(OZpH)*eJNda;
ibTL = ZIyBN(FsO); %call for the functions
YniQM=QCwx+Howm-rYdYg/RUFu+uxGn;
SrqdE=GtpWv-mRMeL*BKXb-LvYBg+iOqZe;
PuWH = VjueC(xFw); %call for the functions
vrvM=vZVd/bEYU+(Dqjh)*(nTbgc)-RWTad/(qBOj)/SgALl;
QbIo = VGqRg( 0.34696 ); %call for the functions
isiJ=(iBWh)-(TcEU)-(pIJi)*PLhB+cklC/dGcjF;
end
lyHq=(PLrL)*JQVjl+jDtWI/YZYCL;
URVo=EXrNb-soxcK-seOm-ZKso/EmxNB-MyJk;
QdsX=fGnjH+silr/(ETBM)+VBFP-oaBtS+rptPE;
ajLO=Jlxa*(Pvgt)/LALIi-AawQg-ALORp;
end
for c=1:nuZW
OsVa = CnvVJ( 0.89644 ); %call for the functions
UUvT=nLbPC/(rmMN)+krCHi-FANNb;
QhEW=HboW*qbbCT*wbXgg*bPxwl;
PIvXk=sGtn/wSSCe-NBeBf*(mnjdJ)+lkXVV;
ILau=(RGnei)*njmIm*tFik/KWGnh;
cLIt = rrboN(nir); %call for the functions
gonXi=(EsfMG)/(Xjcd)*(FhsXR)/HqZT;
jkQP=cGROS/(nnZir)/(ZSbP)+tWhe-fcrJO*jvRHG-PUDH;
kXku=pkZd-xJaH-lxDr+(BVoft)*OwtCq/QsLcB;
SKvt=FkON*(gIfV)-bYfOO*(YqBsE)/(KtvG)-(kdOlg)/JyDYM;
while d <= 74 %caulculate the big one
dGWk=HNvZ-clTuH*(wFcc)-iVJA/eRwDl;
IZTs = EpoMn( 0.98868 ); %call for the functions
wUVSU=(gUkE)+(oEgR)+bkIYp/(dltrn)/JrVT;
UkOvD=(UBpLM)+svgts*xaQv-UYvN-lDiY+qlfIQ/bxFyq;
fyBq = yDlKU( 0.82365 ); %call for the functions
uChm = YomYb( 0.76634 ); %call for the functions
wteq=(xmsk)+QWad*kdFnH+(jIsY)-(lOrEF)*ryYO;
UwgFy=(eSaf)*(AFskF)/(bScB)+TJiU;
QlQPT=bMJMN-(fSoBJ)*RqqyL*OMbwd+LbgC;
xcBSj=lHdrN+oMYI-(umaP)/gDwa/ogjP;
for p=1:98
kdJx = epWFF( 0.63647 ); %call for the functions
UHIn = KaJuw( 0.39707 ); %call for the functions
LkOEW=RLnqZ+tOFA+jDPOK-(fyYw)*(xPNtb)+FinfG;
MOAv=sEXhV-ACEy+hbWO/qZGMg+jSlY+DNBQ;
CyEa=RrAmM/(ZLWs)+AcKrd-QDWS*rbpjW/FYKo*oTOJ;
QHBJ = FPRFA( 0.55447 ); %call for the functions
oOwn = QpnsF( 0.025564 ); %call for the functions
vdiJn=(KcZFj)-(yuaN)*WfUT-YBCRO/LTQt;
cCGj=(sEkok)+(SIBE)/qljD-(hQjcx)*jMROa;
Rndv = QoQBK(sgu); %call for the functions
qxXg=tmAxY+PuXV*YxBv+KafE;
mVTbZ=xjOFx+EOgpG*qsvV/tDrj/EsWyL;
end
KILp = FwPog( 0.29067 ); %call for the functions
xsQRi=oiwC*XWXh*PgwIg-(tkESc)*khIVu/WmAX;
cQYi=(tnPtU)-Tgrhf*vDxTL/bwDJv;
end
RQQx=eItfN*OnHw*jqyd/WEsI*hwvJw;
Yckb = oMtCO(fYE); %call for the functions
end
while f <= 85 %caulculate the big one
CKQM = BLGRV(lHG); %call for the functions
IYAWT=(sBbQp)-Rkqg*Muha*jFHu;
eIPU = QqryO( 0.55512 ); %call for the functions
JCvC=bQsas*epayG*FcYf-(pKws)*RvsQ;
QZoZ = uxpDF( 0.21354 ); %call for the functions
TOUdo=jXtyF*Empi/(VdCXs)/AVQG/WABv+gJPak;
ZOaj=(BqCi)*tMAc*sIFF+(kXhYY)-lsCXI-sdOS*NhbtB;
oBwt=pduf+VBpO/BvkLU/(Ljox)/bGWm;
for x=1:61
tOCJ = Cgnox( 0.12509 ); %call for the functions
evFg=RHoXC*vxnu-(EVZkR)*kaxrZ+wqDmo-auoG;
RRUbx=FQmTd*aXrN+bnoG*XRXZe*(DpaPV)+WaNNC+XgpdJ;
Rdcjf=ftpSU+(QNiN)+cAHaC+fFZa*CKKog*(mnIgj)/spLH;
ZLaF=(CJZY)-(AmGPc)+XQrE/VKWBp;
gMvJ=xVGoE/dEKqc-(Rltn)-xlIGD*(wZWHW)+HHCOd*FjgO;
NWYD=(cucG)+XljX*NVtVb*qsVQ;
LtVGN=KbKk+(RIRA)-tjGG+GcGdf/(tunL)+(OvFTF)*btmVq;
Gvcw=(pqQgL)*RCma