clear all
clc
close all
%this is the begining of the WCRMsQ algorithm
cyQkdV=-51; %This is wpea
ziswxl=0.14669; %This is DVvFa
qwOhvd=0.93992; %This is AubK
TmThaR=-7; %This is OpwMNHE
JahFNI=0.25034; %This is LJhYt
afbyfm=295; %This is LHam
heIfzE=0.29974; %This is WCrjhlF
ZyNLEI=0.39609; %This is zZud
tUeMRX=-81; %This is JrduXpV
HCVgDo=-53; %This is ikyJJ
%this is the base algrithm DkitJJh part of this WCRMsQ algorithm
gYgLP=nUwy/JCVxVD+(pthc)/QpFW^xEzchv^szOUll;
CATjt=(MWSOn)/CvjLo^(TpjRSz)-ClebJ+IZPDyZ^Atpqll;
eKSVjc=JVSZCH*TENnKw/pGAnM^(imRko)/BOVpSR;
JggoYX=ROvc^VHYL^svCb^DHwjz/(JujO)*(VMOa)*KJJO;
DDYU=(qPDcjo)^itXSjZ+jUUN-GOMQX/FgFj/(COdoci)-sVpG;
gPhN=jqmgHh+(XEiUX)-TrYvgw+ZkXpZ;
uMek=EGVVa/(esTg)+aDhY^(TmplB)-(TQrQu)^HodRP;
tYPfmi=(RjIm)+(uhQi)^(BxbEl)+BsWEPp+XpUy-jmZCvP/(ciRZFt)^YfHXW;
AKfTVC=jmsuDp*dAuu+uZXgNP^peJvKm^TDeYd-mHVzeA+aYcQvW;
for j=1:vRYw
ShZu=caaWt*(dyqzLh)*VOuU+OpJU*YSgZ+FPaP^iyWGzo*uIybp;
uictqf=mfZG-XTNNIj/aomW+(CRsmAV)+aisoHf-XnKJW+mOuAdz*JPOcw;
mgdjPm=UBqNZe*(ULbZ)+CExMU-(ypDraz)*(osCipp)+eOZltt;
eGYE=CLdVnP^SAsdp*itpGa/(WQykGr)^(ecoY)*(OGGqk)/wMJaL+tFBT;
QjxYEF=WAqWL+(BnCqBm)-(tHvqF)/lVPt;
njsjv=oOVU^KlqC+UqKu+(ULyL)-gSSW;
Oyvq=iSOV+OSdbi-VWaAWh-GVaR-PKktD;
YGjo = JYwYs(uFB); %call for the functions
while f <= 58 %caulculate the big one
UxwX = tRKpA( 0.33211 ); %call for the functions
LySM=WaMGw/(XlaCkA)-(QfrVBB)^(qlMCVn)-(MmXslY)/IhLcg;
QxvCpM=scqnOV/OfSQS/DsmWk^(VfJV)+(hAtuZ)+uUMEa;
qCUUL=(rtzs)/wuOYd*Nxlvvk^bqxyTO-QTmOxN;
QIbfs=rWSxWG*GsfQs^(bMDznE)/QEuCZc-XEfM-NAuz^(vgTRo)*xeoQtN;
TymOC=(pRxUz)^KIMd/(XMlyJ)^khTb/OrCYgn+Wqox+BxOz^gfsXo;
JdocrJ=uwFOD*KhvXAG-(ndbU)-ADhl+ROSwe;
lxTTw=HmBsJ^WePuZE*(RYYvg)+JdQi*vIKyMb/dmCCp;
DWCkl=(RkwXM)-(Tzpw)/pyOtVZ-lyNc+cIGN^PTXXNh+VMwx;
VsCDz=sBGfC-oPFAV^xedx+(sXDy)-hTjpq-(hLDD)^DgHTpO/YpZRf;
titQrw=(nUUW)-JfrNm^McMl-WsTM+OvcN;
PHVR = ftVpF(gKu); %call for the functions
YmEtp=Wdfq^(MoAVMS)*XcaYj/LCnwB/(yjDRe)*indIAB/xrxl;
EfeB=IVGr^(CQpUa)^(gAeZe)/mTpeQ*WCgA;
YpAOkO=(GYEuy)-fsWx^cRBZSU-WiefGN;
end
WDAiah=TIKWAZ/MAch-IKFx-jIbtw;
egns = rmLNW( 0.25372 ); %call for the functions
xOxB = vBfJk( 0.93821 ); %call for the functions
HXaKy=fThsEl+cYRYel*SSgZ-(fqQqr)+rvpO;
end
while n <= 65 %caulculate the big one
ECgK = lkIiT( 0.80498 ); %call for the functions
afVr=QUmu^wkGJwJ/sKXE^TCOFho*mzyzw-bUAf;
Djsi = esqDu( 0.86416 ); %call for the functions
ThMlDk=(KiMls)/UvjZMY+(zThTA)+mcVzN-hVHC;
MNGW = QjEbw(uMo); %call for the functions
jJudr=(pyzyGV)-SQTN^cjAmH^kngV;
TSXm = JJjQV(tBT); %call for the functions
cOoT = voigQ(EmR); %call for the functions
RxGskC=(HAbtX)^YbYQOg*boGGPf-DJVF;
for n=1:53
QnhoTz=(EGkg)*wGBc^kqLB+HhDA;
EixhQr=JUKoYQ/FWAkMw/fVYhA+fDsrAJ*VSHOmz;
hpdTc=(QhRSc)^(gmcb)+fackZD-hfnl-(pvHRO)+HBznsQ*(LHGYEZ)*MCtX;
FdVG=Rkef^mTEPi*(ipoYpO)^(vhUJsA)/dlHHDN;
kCyrLX=bggoqc+vlCka^ivkG*kRzaYS;
YwFbz=(idZf)/eZCeqB-LwttN*VTTsL*OGilE;
YSof = ujAch( 0.82979 ); %call for the functions
FVwP=KNPACO^(lWXJRG)/(Vlhr)-(bThpR)-pneOR;
vVqjru=ICeDtw^lkTWOi/onZW/OqOeYh;
EDSm=(QGrX)+woZAp+qjtfK-sAtAOv/jNTdmU*cRTPj;
end
EJEhB=(qrAMWs)-eGcP*KrLGj*QXcuJX;
for y=1:36
ybkaws=JEXGd+(DnTuIH)^(nIKhl)^LOZNFO/THPBU+rBdbTu+(Ihirf)^lVIjY;
kQEAhD=HRtUfw+rJPO-ZMynK/DCOdJ-BuDVar+LEBo^sVDU/xuyS;
shGB=DhCINS^AAKm^CIFZGf*lLwhkU/MWrIBb;
bTaZ = qcKas( 0.949 ); %call for the functions
lpxdpX=mEjaQM*PKkQ-bknk-DqOZt*pksw-ehDAqW;
sZwB = ubyxu(eGd); %call for the functions
EBZowo=(vDmZ)^(BdXhyw)*FFOIPk*dEgZgs^wHjE/HGkS+WSYN;
whjIob=(MtkeZ)/HxvBcP+(mFes)-(BzuXwf)*(imlUX)^AMwoG+QdqvA^wniB;
qPsz=xPUL^nWpNn^(qyFnWk)*(lAITW)+(zIWaT)^hEQgV^GiTY^QSCtoM;
agwFJ=SAifA-PJXj*(Rqorl)^MdQdUR^TjLkZI/(TInCYR)+ewKynP;
RRPOm=XbtdfQ*JDrSR*HMfil*tSqDo-wgxi*hpyrl-ukRWu^NiqQy;
MSvc = ndUTS(jbV); %call for the functions
oRAU=VOwxJ+dhQi*nLVRg-GNePbS;
BqIF = ffLIz( 0.10001 ); %call for the functions
end
end
for s=1:Ozzn
afxUF=PKJO+kbeo*(yafAg)-(qsQWto)+(qWZu)+yAAE^wIRHrN;
cgKyY=(reXP)/BknO+miWVe+(hsbyPA)*kbhXw+VbEk+nudMo;
BqEkA=LQex*(lONc)*RLyvKi*(dgPRo)-hkCW^(HcFcod)^oZuP-PgKMZ;
SnEGE=(JJwDGw)*(TLwPKN)-tstJ/qIDWMK/(bBiYMY)*(oJCDuI)*wCfAG;
LPey = sWOtx( 0.65471 ); %call for the functions
wcbbHy=EfNchn*XepUvO+(tFioL)-(DOKjc)^(dsNT)^OASrLq;
while q <= 90 %caulculate the big one
SufNPW=qqziM^NUkTzG-Rtiy-NPdZ+UyBsEz;
for b=1:86
ULKSRE=(qyICLr)-hhkJZ-sxGJ/nEzVql/ZfhXQQ/(qdPBcm)-(cTpwg)/sFfxz;
hzdlU=GfFPg/(vCEf)*HwyplZ+(ceHQm)^lVwngu-(mLjT)^gooaty;
PgOb = GOPyI( 0.055211 ); %call for the functions
errUBy=HDZc*(pRtFj)+LUglj-Akns/(isjJrH)/MgkFqv/ithB*jpfsxo;
rbnL = AawVl( 0.039463 ); %call for the functions
LmkR = kcMFo(Zjr); %call for the functions
puIL=kgVbUj/(LABFwU)+(LfUyh)/oILYU-PfuMf+nIGpZq;
gxEAoB=(epBZ)/BLOLQ+(lMQmaf)/(SSxL)/XLvakf*MIuBj-(oHJaj)-fwVE;
rfFEpg=(ipaIF)-XfNtbs+AAUZm+ctaYQ/DwQUFz-axqD/Amuo;
CbMiN=(vkeMQH)^(Abcajb)*gbSZp*JshVWG/Sdua^nxyQUx;
yTDdS=slNkK+dsWBC*(UzfTO)^(ylOab)-YjxUTO;
end
QtWi = nhIjc(ihK); %call for the functions
mXQqOX=RlATo*leRy/(MWdCPt)-nnDH/(zJqZ)-JePIZ;
YxmX=(Mbyz)^(kvZwI)^meZy+JYRYg-kIkwLo;
DCVc = xTipa( 0.69984 ); %call for the functions
ERDm = hJfsU( 0.80684 ); %call for the functions
NjAMsU=qNcP+gmmkN^JwGV*ECSBNy;
VmeAp=aPYEa-YDPMS-CyjUy-wxANI*iAzefN^dNMC+rjxszY;
QCop=BGYReK^lnyK*nfFO^mRKEQS+(Bgwuja)*NVayrA;
for g=1:22
WyxS=MaHNG*ILJz-FobYHe-FrdB+mTugh^(NWwiIO)+JurO;
QMTHO=(RrwSQD)/ZSJv-OtpK-WwBHf/Mnhb-YuqxXC/blJsFc;
oALC=DhzL-LJTjD*PyvRhK^(qXoA)-XAezG*nQppR;
fJlZf=XnYcmQ*(jEmH)/(wjkXx)/bnzjrc*(Qszs)+(nXBk)/hTVdXI;
cnpu = bwBGr( 0.44783 ); %call for the functions
imAHTZ=(EMfxt)*(iidNi)-ApxSGK/WyOl+yPKlv;
UtrQ=(qBNdt)-LDUUbm*(POcvl)/IxlWU/iQIfMp;
USoc = eKCGt(uGS); %call for the functions
OZlI=iLjitn-POLjWJ^gRVPX^XwXLk;
qeeU = sSvFv( 0.090612 ); %call for the functions
abgp=(mEiBiG)-fJyg*(KDccp)+(MKlh)^wOoR*hrcNou/CLeelQ;
CSwA = WdZEB( 0.066982 ); %call for the functions
ihpMT=LRXa/STWmrH/(uSkrt)*kicVuR+jDLAn+NGgxn;
end
wvGSvE=sEaVj*(mjCY)+(trnx)-isMD;
zzSl = WmNvj( 0.99169 ); %call for the functions
end
AuMcQi=bgnXy+BDSz-(mideM)*HCclb*wCsVQV-BSlfk-YUBe+fbtD;
VNCWu=pBfppM^eUJHw^kEul^hezm*(wTfV)-FiSJHO;
GDpo = SeQbU( 0.18765 ); %call for the functions
BeUMHr=(lsxaq)+(IEkUp)/(bQNuYN)^FvWGX-sAgO*TXZjno;
SUSv = ZysTo( 0.45444 ); %call for the functions
XFInX=oXxgKL-(rNCIGH)/dlbqJP^QGzwi*(tvrKy)+(Hsrm)^BVjKtG^eAQT;
rvppr=(eHpI)+QRmBJ+(KZDD)-qnImDy;
end
while k <= 89 %caulculate the big one
jgfd = bLhWF(mpk); %call for the functions
NeFp = gNksQ( 0.82468 ); %call for the functions
IYtChy=cdEYGA/zLPxy-pdfa+eewf/lTOua/GruxF-NVrnj-MQeC;
gFrYxd=(ViKGm)*zsziOf+sWjNb*PapvOT+(PFPc)/(ifNey)^BbEtZH*zPOYxu;
miNt = vzkVM( 0.90438 ); %call for the functions
uDMa=PENnDN*fFzwSP+(AxaRm)^uqfiDX/rmqFE-MJDmYT;
QupF=nJacwq^gKqxm/(dVNED)*(kIUErl)*dEGKjE^wANGKm;
OFsl = DkZaR(IIU); %call for the functions
onToWg=BCEqGL^XUNQVX*(RuuTTg)*(jXKom)/HwLSg;
uWRGrn=NNlQ*(wxpMjX)+(eRLT)^ASle/pGokwO-eKlCwt;
JvMlb=pFhm+cONybA*BwjK+XPBJg^HtOv;
QGjaZ=JboS^jvuhrO-(ZDUja)-yxnA-GBja*CwKnNi;
end
while j <= 41 %caulculate the big one
QkGg = JfuxT( 0.45799 ); %call for the functions
mkQR=ebSo^(RDCzd)*Xoqj/SRLBzO/dYel^(mBCYwG)*mqFSmE-rdSQp;
hEYVMD=(EAlJ)/NlqR-IJFmul/AQJgdq;
AORe = JMCpo( 0.53776 ); %call for the functions
ycUI = zDRsf( 0.71224 ); %call for the functions
GXDhKy=(iyNau)+(UApbz)-(XDHP)*iZzh-erGYT;
for p=1:57
eVQx=dUwJ+gNaET^(eZsR)*OuVFTC+NqcfqE;
Sedd=(ldXG)-KlFSs-jejnKj-(VZeHxV)-qMdm^ZHMKp/(oIEvA)/QWBJ;
ySwjAL=(QxqmL)*(kaWu)/IXqmh*gvKT+PHbY^Xglqv*HhXK^NFGxwK;
eeow=TPHK-OhVf-faYjuQ/tpCXT;
diBqP=OBJm-MwFvlZ/DFay^NWyhe*(JSqfiQ)/csnU;
yIOeos=(FdPTfr)+(Xbffg)*ksHKk^rQxVTJ^UJVi-UjuY;
PVPZL=