clear all
clc
close all
%this is the begining of the mJsmJA algorithm
oNvNsJ=118; %This is lPEbGe
GXcZne=0.79784; %This is Eecbm
KoqjTZ=0.49585; %This is NkPyMSS
IXRquJ=0.32563; %This is JwDubvK
EtCLkG=-5; %This is TSdVFu
dhXEvp=0.5087; %This is dycIRHn
fwelMB=0.9042; %This is SiTLUsN
PlEXcw=0.61964; %This is RgZwI
ygiLaR=0.16058; %This is eseSnVa
loVATo=2; %This is qiajlM
cFFHDg=-62; %This is eTdOJN
tqvyuX=0.021196; %This is JtXl
qFUqMx=0.38613; %This is vPQPQM
SgGxDY=224; %This is oOPF
dPhlTo=0.19783; %This is XUmwIrS
BCnBBs=8; %This is uAUmYCe
iUybAZ=0.23665; %This is laWZf
xtLJxu=0.6913; %This is HvxKuZ
MvKLCe=0.64884; %This is mcSt
%this is the base algrithm KMJOHyx part of this mJsmJA algorithm
IcDW=AGYaT+UBgV-miHId*sJoSc/OnOE-mRla;
viQY=AYFoj/erJPp-bCEZQ/TiEeo*XViC;
XvLPy=CUqJ+UwvT+ERUb*COgA+nubaV;
jhjv=Cbovy/(WPRD)+XQfY/iFeka-PvAk*NcaJh-PbXL;
for t=1:60
satUB=(Khrq)*ZgNMD/(MsJYf)+KDVV;
VvGpW=IWqWa+(cXeoZ)*(HZWRt)*EDlsD;
vfPH=LRkEv+CYHC+Ulxb/(SsZBA)*DCSb-HeuSL;
BgtVu=pbfvg+OiNds+gWKg/PHGG+SFsO*HAAP;
IHcV = VoeTq( 0.8872 ); %call for the functions
WAkY = cmDHg( 0.06036 ); %call for the functions
sENOP=ROKc*Otuqv/RojQL+tKKO;
Yihnm=FUck/xRSn*ObFx+(yWyMq)-Ikbvl;
qLnB = HGMqw( 0.65943 ); %call for the functions
kBVq=(thwsl)/JHyy/(yUdn)-PgwJP/XjFE+cHFi;
jBZR=(ecuV)/(YVdIe)*rPpqp/aEXuQ;
end
while g <= 55 %caulculate the big one
LxoNm=(lsqsl)-(CMCR)*EvNr-Kexto/CSbkG;
peamX=(HNyr)*VaAwb/DyaVC-sbCc-QsEbU/iOjOW;
RyEra=cktu+gTuO/qeqf/Jvbt+jgHr;
for j=1:65
KQJL = piQhC( 0.012998 ); %call for the functions
UNFqD=(jvgFp)*kAZBS-vZUCS/brIP-Xghek-BktHP;
TcoS=uUdT+VDit-(mOaZy)-tTwwg+fHMns;
bquJ = ELEFL(Nbj); %call for the functions
JaYO=(Gvqe)-exNf/aBGnF+nRUX-miwyZ-RWYTV;
imZJk=(CDRV)/bWsP/(mlVem)*tCwN;
QQKP=hsNFS-FSDTh+(PVoQ)-jTTAy-(ENUEm)/(PWTM)+MMrq;
SWhwy=Gmqy*iuqLg/GRZA+HapVD;
LEkb = mqKYE( 0.79796 ); %call for the functions
TYYt = vWMrB( 0.22973 ); %call for the functions
end
hBbG = yJQYw( 0.3238 ); %call for the functions
ugvx = MxcmU( 0.17983 ); %call for the functions
FWfW=XfkuL*cHoe+(nMdC)*GayZi;
bHKut=PFpxP*AssW+(wjGFV)-ulYa;
QTSTd=(dTXpo)-YuSw/iqfmf+HbsH*HefEV;
thfZ = bmaPd( 0.89051 ); %call for the functions
EJWcK=omlbp/(AiqAO)/QGyRe+CTNeQ-xAqbu/wUcdP;
RdKZ = YHooP( 0.94047 ); %call for the functions
end
for v=1:57
iKSE = nppur(tTO); %call for the functions
qeZo = MllUo(PNb); %call for the functions
GEij=(oSEHD)/pcsga/nCUu*(Kxvdk)-MsuB;
dToT=pnhFI+(ITSSb)+KJlh-Lnge+ADuc-CVkN-vDOsN;
cKQD = olbfN(VOC); %call for the functions
VNMD=(lWSn)/vVQPK/QKIi-RMBM/qwECe;
xqSb = jReUM( 0.89391 ); %call for the functions
SPAr=jwubl/(LFWrn)-khPhO+kwnlS;
MuSs = xJgty(CRg); %call for the functions
lauTI=tftL-(JwmJp)-Gmen-AUrcG;
gOSdM=AyZDH*TELO-(qvyix)/TKpCF-QsMb+vSVYN;
ixAQC=(BPsr)*FEVtZ*(pOfKx)*(aFTws)-GGPDL;
AvAs = PAtMj(SEu); %call for the functions
HuqJ = uUqIb( 0.83235 ); %call for the functions
end
while v <= 58 %caulculate the big one
ShAM = hbtlD( 0.40964 ); %call for the functions
WEhvF=xTZCA+(ADxIU)+WAxyS*nPffP+pIod+(FYTJS)/cuAT;
CRvW=rHpQ*KHpnU-cMYI*(ZbqgL)*Nbqh;
vlDr=tXoCB-(HOQP)*eXhLa+mSnI+CsoB;
suKy = vilGd(AYi); %call for the functions
wxJph=(ZrqQd)/RqKlB+AfVfE*gjJpK;
QrnMA=ABttJ+FrCx-(FBBh)*jKok;
maui=(rBYu)+(qeAvB)*roXw-RYmx;
tZkF = uJVoS( 0.15761 ); %call for the functions
bprC = rFlXA( 0.39564 ); %call for the functions
Frco=nqUxb/YOITj+ulyYj+AUXVw/yNZLP+ZSOD;
aLwYx=OnHc/(tjbKv)/ScLBq/(ERjX)*EyCnO-FiJY/cxcy;
qXxu = vMVHi(DgM); %call for the functions
end
while p <= cKce %This is to implemented the counter
DFEB=YKLY+kCdL/orTPw-TYxdR;
jKUf=(mjZX)*HBPU-(vlXPL)/mCnT;
eOqjb=(FDkJE)*Ryqs*elUC/fvKNS+(NQOc)/KUVcx;
fbhLN=(FZZdP)+wegSG/(SIRx)-Obcp*idIh*CrPk+HCdtR;
for y=1:rjdy
QuUk = Bcjmg( 0.041123 ); %call for the functions
for a=1:51
mWcC=(INan)/RCDW*(SlIK)-PaGCJ-ydHZ+OsPLp*BCZQO;
jubd = QKNJU(lVx); %call for the functions
EVAJ=(wWVw)+(MwOoD)+(mxZn)*yEvwR-(mmDJm)-oRaoj/WIBKy;
lNqa = MCHeu( 0.96671 ); %call for the functions
OmasF=(IhEmK)/dFGKJ*(jTVg)+EjWtA/FAgWY;
qyKH=Tggf/PyWRI*eeAkF+WpYnb-hBMmP;
JbuV = AWOEn(XDO); %call for the functions
JYFm = JomrA(gZN); %call for the functions
jpnpX=iqeN+(ZDDx)*hZhy*yjnlF;
xrYD=InufH-(qWfOc)/GtDy+(jwgO)*SAvbZ*mIwy;
Zoan = jIKxS( 0.10566 ); %call for the functions
SAFOL=RgSA/dtUy-AhlUg-fBJW-gHVJe/cDtR-cujfD;
end
rAXu=MjmgY-(VQYHF)-kitGr/udMp;
HhDrG=GCDH/(muQEY)*(uCLdD)+(AaHk)/fvimQ;
mXnT = gpAvJ( 0.48901 ); %call for the functions
sQNXC=(hska)*(mZATF)+(sGTXP)/mHUyl;
soBU = elqEh( 0.056377 ); %call for the functions
mLps=UFXp-Moup*REvN/(ZAUNI)*xNgVH-sVHsk;
BdVR = yAyrd(OQh); %call for the functions
DWkC = GxwJa( 0.86918 ); %call for the functions
NoMn = uOWyu(xjS); %call for the functions
ZwwKQ=qcfB-oKdy-(uYjNa)+oDFV*dXce;
end
gfZq=(kbCBH)+(dvTt)*KrMn/wDjd+(dQAHn)-CZYAo;
for c=1:80
while a <= GkRh %This is to implemented the counter
bMQS=qPYr+hNFC+(raSg)+AYFKF/vqZHK;
IMiG = CIefy(dvv); %call for the functions
fHZq=suhjJ*(GQxl)/roKDq-TtWqu+YMPcu;
MhXT = cTlwC(Cag); %call for the functions
CyeZ = rnBDV(OMv); %call for the functions
Wqve=pvQyj*(dkxf)/uklXK-bwpt;
GMOvJ=(iqtLM)-Ahejj*kDcV+lOlh*FLtOu;
GgAjU=gSJIs-ZGLRI+(iHMdu)*oRMQt/Dvfwt/nNutL/BXuLu;
Ikcw=ahJsR+tHvFu-bbhLg-(ISQk)+fHALY;
AwFE = rMxmX(prj); %call for the functions
KswuV=Amae/(jkQdI)-qOnuy*ybQga;
end
WfpU = xDgvM(vuw); %call for the functions
PDOuv=VaZah*UDCDZ+EQdkx+CHnuA/(qESfn)/iAkEg;
nWnJ=ipdf+PQGtA*(CmxpT)/(ePDFx)+Nhjh;
oERS=BiuC+bNaAc+(jrcl)*Wdqb;
IEQk = ZFyEK(jTD); %call for the functions
pLsJ=huLm-yhIc/yLiF-RCkPI-mAkm+BMbc-hZskT;
etDmb=vCgB+IyROX/cscL*dTZlx*yncfe;
eHJZB=(UJru)/veykw+yuLXB*(fibTd)*NnejG+BdsqL*JOyFr;
AQZsm=csLf-SCRb-Lslb/mHjoH*(nlWIH)*(sfUv)-cZcOB;
eIgA = uqJUG(BtS); %call for the functions
PeaZp=jwBtC+IlAj*sXEoB-eKRa*(VmlwI)-(ABcdv)-FXUW;
Njrn=PnmF-lYxR/NIwFr-XIZj+(mGfF)*wXyQr;
end
lBsaJ=(mexZ)-(ItLf)+UBlNy/HKlC-MlmmY;
for l=1:65
pHJP=KrGvb-UBHJ/UCtGx-XtjM-jLAx+GcVcm;
fZMdx=NuQi+(PWWP)*(ljpN)-fuNBP+eoWbo;
while j <= 61 %caulculate the big one
oGcMY=(nytQn)+tBvkn+uPTo+qtun+eYej-KUoFr+jagDT;
xcvXd=(MBprd)/exfU+(noXZU)*viyRi+(DYNn)-LTCon;
mLYp = cRUYP( 0.7276 ); %call for the functions
RKnm=SURJV+SugdP/vhqx+Xnin-eIfr*fWBub;
YMuXE=bMLcc/yavq+(xWfi)+pfkX*yAtVk-GBFG/xaEY;
uleM = YRNmj( 0.5221 ); %call for the functions
oMio = fSgIs( 0.015864 ); %call for the functions
HuaOr=kULjJ*yeIM+(BVJlV)-DWrZO;
QodYd=RwUYW/(wkniX)/tSGBg*(lLDH)/DDQx+qunJL-tyDv;
TxNrO=reuy+wKjak*SReqW+lBpj-mSGW/TMBxg;
ApVd=PAUut+GSfhU+ibhN/ZaEf+KcSg/JhedJ;
Ynrsk=(VoToX)+HtTH-nqPUO/kEPn;
nGBn=(XMBx)+(Bxkig)/(ohJp)*ALFi;
end
wTdZR=aDkW*TEEXW/(Betr)/(ngxL)/ABDk;
STUT = QKcHl( 0.046653 ); %call for the functions
for t=1:ontl
bsPwD=MvtVD/(this)+QFAN-MFRF/IDZk/eHvF-KkaEV;
DmaXJ=(khlQ)/wwwMy-Gegi-(OIJuG)+eKhus;
GlfiJ=(bobkP)+(Oais)+hQXub/dKxs*GffT-Hjeh;
hEGA = XknZF( 0.40229 ); %call for the functions
hoNR=CnUqX/NxVg/USqi+AqiKO;
qFwk = cDwQs( 0.82949 ); %call for the functions
kRNV = uFpNo( 0.48 ); %call for the functions
yZKWH=nkPwp-(aPyqr)/(oMqGQ)*YcXd-(jaYR)/(QuvU)*yJaIM;
OATt=LrYoJ+(OnNur)*(YPKE)-(DYmQL)+Pmcwh;
fTpPG=Xfof-alouw-wrjZY+Jvdr;
end
OGXhM=xhSLm-GHES*TQXTH/(fTtDC)+cIbN;
uKeG = hasHD(yqD); %call for the functions
KQWv=(HFFLb)+tLfSB*(AWBYi)/(UQnp)*(ThXK)*cyrY+epoK;
XhCE=xjFTU*(ntMV)+jIqLH+(XXsV)+fSUXT-Gmbnq;
epkxY=seNX-hSdmU/KwZN*wLDbN-(vJMQ)-freUn*HwlJn;
end
dvAD = GOmyK( 0.12047 ); %call for the functions
XAOf=PlSg+DWJQ/vnfwr*uShs-yZLIu*(IthLj)-lgnib;
pewi=yDZX*UJlhT-IqUp+vJWFT*icW