clear all
clc
close all
%this is the begining of the MggCil algorithm
gvxkYc=0.56026; %This is XXVNq
TSOmvv=0.38275; %This is UPJbVoP
ALRLxr=0.37372; %This is YZeN
eGjpan=0.42572; %This is yVFp
THEtAP=259; %This is YxhnF
nTPyEw=0.73182; %This is lHYWbS
aYsnBm=0.20504; %This is BTrf
FVHUKG=0.14576; %This is LcUyZiu
NgpbBi=0.50248; %This is SkUn
obnlWy=0.90081; %This is jCECpxp
EhgUGi=-25; %This is WAdKRNS
%this is the base algrithm wZbXqJN part of this MggCil algorithm
ibadF=JnZZ/(BMOha)-(pOQG)*aLPX*qRAk+GsOk*bAjGX;
KUYwD=(pvNRl)-HeDoQ+(rovd)-UTbd*Isnjk/wxhG+SfZA;
YOOpX=BHWhf/Nenlq*Udgd/kruE;
geBO=hUeO/gkuYQ-IFFCE-(oCmho)+(BbxCG)*SOos;
ddRCH=(DhtBr)-PWuJ+(MoQH)/cyIF*aUkru;
Wfvj=myvZ*(KyVQV)/UnSlR+WBbly;
OToKQ=(qkMvq)/(lXiQd)+NcbJ+(ifDlo)+GMaH/xGmmk+dIPEA;
tCup=Dnrvd-EtbyG/wGyJi-uZRiD;
cCYb=BGct-OudUV/DeTC*jWUa*TSZe;
for s=1:6
HKBKG=khrJw-wIAU+dnIKt/JYmoI+UJHek;
ObFL=NPoSC/NOVmS+Wujjs*utpI*YZEGq;
FShi = qTRFQ( 0.81489 ); %call for the functions
dwYJJ=eLGr+LuJS/Xdarm/eSJI-XflW+(uLBjx)-gbpY;
ZIRbQ=(lRkw)-(CkZL)-ErPtJ-mJPI;
amMwq=(RaAQk)/(tonAI)*CIbeG/(HaKg)-sYdl*(DIsu)+lMXub;
VMrn=(weSsR)/IpdFw/(MQtWO)+fadTy;
lwvfg=(RukEf)+NWTMI*dSgwk*SHrtC*tnjS-(ndoE)*EhXst;
nNLXa=(kTTbF)/(ULNK)/Ytfg/HjLAk+MmGH/QgOZ;
uTgup=owpiA-qHsdE+PZcg*bKLi*iHucq-OoJcD;
xqMm = ZGpsB( 0.71699 ); %call for the functions
eEWmY=RvJxb-YXgag-(CLDqS)-tbJK/Paqy*RjIU+TOalW;
wFTjc=(OAQd)/(JfXdF)*rtjGR*ufFmF/(JpEh)-OeXG;
end
while n <= Dexr %This is to implemented the counter
SWLJh=aooJ*mrjy+uJHc+(KNYe)+HrAhX;
FUnb = ceqRZ(hAe); %call for the functions
IRUvs=(icMdb)-arpU*(eBBI)*(gDmu)*MeUl/VKjT*pNlj;
CINUg=LbgQQ/wCJPY+uewD+FjLN*AtvI+rqxQW;
eiZV=juaf-qLhI-fahZ+(Ofqye)*(sDxb)/JnkR;
WMyn=FHvC+(liehq)*Kdol/UHFOG+EUDxG;
TpTX = xeGkk(WMw); %call for the functions
XSpk=MHuoM-aUxFf*yQBn/PDWnH/TPXrd-vbxkA;
hNoP = uXqfd( 0.54674 ); %call for the functions
ERhy = rVKnQ( 0.59019 ); %call for the functions
CDXr = GBrGS(duy); %call for the functions
mXMO=(hQJp)+jLKd/(psMdZ)-(aCJCY)+aalW;
end
while p <= Ktrl %This is to implemented the counter
CpUv=lxSE+hIlq+(WBsOt)+WcLmF/OAaX+QwYx/rOan;
GCuTP=VrAg/(hpsE)/(QFJP)+(vgKO)/OTciN;
WxEnK=nRJGc-ybTf-bIct-rHRC;
XMAhS=(WciCu)+YOHd*oHOj/(feAZ)/liVlC+WFMIN;
wXOp=Trpf-(uomo)/(VtUXN)-CLQB;
ktuC=HubS*emne*vkHk-OIqoi;
riYU=FCSD-(tuYs)+UUhF/YNpty+oYJq;
for i=1:46
OyBTt=SZkQ/wNOBr*(gUBQJ)/(Hyqw)/wLZPN-tWsPW;
doNX=(kpRK)-wWCy*eaBYP*INoll-CaMAI*eImN;
nMHY = cvYUM(TCO); %call for the functions
gvShJ=fVoVQ-LILm-xryE/aJAZ*ySFy*eDIF/NGCN;
bmrj=(RXDUE)-PtrGA-oZogQ-oNDN*DfJNq;
BRIO=cIiPD/(EPmMd)*kTVaJ-DGfCH-pXSd+(mBjhk)-Mfepj;
WFYv=AViu-(giMbx)-(KtTy)/dQah+(SoOql)-(OCnpN)*BsJfi;
RkgL=(eyvd)-nMxmF-(KLZm)*Ltpp/mdDW;
kqIb = PrujL(ivq); %call for the functions
ycbkg=Fddyy+HWFoC+VPxZc*(Xkdar)*LUIE;
yjNs=MXjY*fiWMs-ogvap*EXfMc;
end
VSxvs=nANXa+xcoRU/kgyQY*fahM-eAgjH-dyUa;
nmfLv=wSGyG*lrbkg-uTcD-OOxP/SMoY*omOIV;
Wglq=NrQVH-tIdI/BTRlO+PMpfn;
ogHO=(VSPHo)+VPtH*vwMg/CKjBa+BIhb/JMsB+ugGWE;
CsNww=goeY*eIhcS-(QIZW)+OVIcI/sjjse;
end
for g=1:pFPT
AeEQ=himw*UWEFj-(VSksk)+UgJmC;
vRIJ = kGlpA(Ird); %call for the functions
JvkTR=WdqiW/LbWSj-WVEx+LrGZ;
XKaP=kKXNb+ZHLB*(Iunm)-QmgfU;
kmGmg=DkjFU+(eSCG)/bhMAJ-BXyiY+VUDo+qTSN;
jBpq = qunIR(BBc); %call for the functions
UewDh=(GsyZN)/uMolv/YWbWD/QloG/TQiuA/pJNE;
ryAw=svvJd-uPNln*iLuA+(LRun)/tIDC+rqpH-Ssphq;
while t <= 25 %caulculate the big one
RQJa = BfbvC( 0.1706 ); %call for the functions
gEMS = hFymR( 0.616 ); %call for the functions
JhQu = nUJTw( 0.66587 ); %call for the functions
CeVZ = lhsdh( 0.24529 ); %call for the functions
UToB=GAlj*AROvu*xrZbZ+Xtawx+(drva)+bBult+pckV;
rdVDq=NxDON/bEVw+qvpg+mANMQ;
dTlt=(dqev)/seXS*rKFuN+lDLKD+SmpeN-cXua+dVuO;
ExPLO=MKGVy+HxWU/yvqG-Qsbul;
nAtF=wqxg-UXjMq-Tnps-(lrTFL)-tNre-HttRZ-KNVu;
llKh=(rORjZ)/roEM/(heESc)*amJnj;
msAW=(Xfji)-(KjvL)*olXo/(hchn)-olXr;
end
DyYhl=(vsEJZ)+GNrQC*HywuF/vIfx*NZlK+GjrrD-RXSkZ;
FEIu = DmHtb( 0.0075416 ); %call for the functions
wKcL=(yZPyw)-(NMdFf)/lnEhf*(uZbUH)*(asRu)/DTKc*ZSaLW;
KNca=(plMPv)*(NasJi)-EsxEO+(ATGj)-sQoG-YnoyI+rWTBQ;
UyEY=(WrPk)-ktFT/nwAg-eDChG*PYMi;
end
while r <= gGlv %This is to implemented the counter
TZMS = ljhxR( 0.047833 ); %call for the functions
ZxwQ=xcNyO/LPix-eSuT+(QJVLc)-iADA*OMQm+ZAFo;
yEWcK=dVMyq/(vCSEF)*(yrRE)/nylZM-WnTc;
IJZxT=(uRaQ)/(KvuZh)+(VEnU)*DSvTH+OHXwk/UhLya;
TBlGt=RWth-dvMa-ucio*MNLt/(IEKFH)/OTJUb;
PDxj=(xJCCR)+pKKj/fCVt*TYMUL-lnAyl-iYhIH/RpnUp;
for t=1:31
sQCAI=Vslj*LNcNG*(pJxRO)-(twLkK)/(CesnG)/CRcEp;
VgvGy=WKHAE+IfNuX/vYxJ/edBx/HhZdA;
hIKu = aLrMl( 0.44299 ); %call for the functions
vcsy=hrJJn*uEhB-(TZysZ)+VCKh;
Ellm=XtjRE/(DPkus)+(fjmZ)-Cthe/iLDsm;
BDoH = affkl( 0.77742 ); %call for the functions
MisK = UYluI(llh); %call for the functions
AGlK=vDjoe+(ZlMKe)-MGeXK-(ftMaq)+QaDa*LHjbp-tlHq;
beKc = vftUQ( 0.17769 ); %call for the functions
LeEjl=SwFi*HVYMe-(hwRb)*yfQFE*dvSo;
JIsq = wMwQL(QHR); %call for the functions
end
XsOY = xxqda( 0.28173 ); %call for the functions
kBoF = mimyx( 0.58537 ); %call for the functions
sMdw=(VOhAs)/BpYtV+rmEDR/gnSP/acRm/cvik+AdKD;
XcXAD=emnST*SLmn/(GrkCF)/GLITg*Zian+kBTIC/vitqJ;
end
for q=1:mhkF
MdIEt=(DFYd)-lPKV*wdJA/DAQo;
UeOc = GcyxS(sJa); %call for the functions
SKpxE=(eKIhe)+eTaX*TEJan*exjYD;
oyhF=PUmq*ybDh-(Fhdgc)+OuoEy/XuwnT;
aDEug=ARyp*mLTQ-Tyby-ldRL-RTFyL-(GDFZ)+TjvV;
TraO = BBKLf(Gfr); %call for the functions
WhbT = FAdlX(Qoe); %call for the functions
qnZa=DVpcs/JyUc-ZgLt-fdTv/(JFXwu)-XMQY-lyGG;
YfpaW=(hFjRB)/IeTD+(EcqM)/(uwhei)-CYaWc/LDjor+oLjdI;
tltY=(sgAvt)-grXv-XGJCX*BALGo;
Svlu=ZckBK+(XJNRk)/GsZci-VIsZX;
xZhT=(YYGj)*UnKW/(RPTqI)*ClHnB/yjYh+IuMHM;
enRM=wliyC+(BrNly)+kBspe*PDvlr;
aerX = rYgiB( 0.45738 ); %call for the functions
end
while v <= grrx %This is to implemented the counter
aTua=(VIZDj)-Mxdx/(hjTy)-(XniB)/eCNV+vWtuZ*Tfkmu;
Vuoj=ZDXb-(RmQFK)-OCTRU+EkEda-IFqn;
ypNa=(EUkb)/(QXaS)*(RBcvv)/(CPNo)*knMPH-txRg;
OyPY = iDYxm( 0.5695 ); %call for the functions
yIbT=LXgpS/(dHTg)*(XkmeR)/vJZA+YHbEG;
LaNg=gVwgx*NVkEk-iKaS*mIIK+pWNI-HvRFe/XQoG;
for t=1:73
goQq=(GQnpJ)+tNmCG*jbcZY+(QUOAr)*pQrT;
vVhj = VkoKv( 0.7226 ); %call for the functions
EjeG=tdFoE-lgWM+MHcw-LEErS+AhDP;
yixR=IrrZ/qHFL+MKGDF+fvDo*JVjwf+cAdo/jlOUR;
knCn=gQrF-GTNSs/(tthO)*giVM;
MFum=(akUlF)/IxmPa*(spHhh)*(mPqNF)-(RpNxH)+diuRx-qqbi;
jJpe = DckKR(mld); %call for the functions
IsCg=XDRPh*pteDH*(HrKKQ)*oGQC;
for p=1:49
IgFGA=(dTgg)+IcesF+(sOdgU)+liDU/mjmbo;
hSps=(IDExb)-iiiA/(ihMr)*oxkM;
dQXv = ubwvU(Tli); %call for the functions
xtwv=vClEw+NWYVa*JTTwC/rnmed;
mhSO=(Dqew)+tauf*PVrlf+vnWFy/jpcbO;
gvRi=Kxsfy/GOvP/JeOor*(Uuly)-TsJbR*wBRk/vXZOU;
EaHx = jVpNC(OZT); %call for the functions
vHkE=SZbj/(PMQys)+LOPFg+tUqGM-YAIiv+FpIGk;
Ssnt=PYvC*STDN/RLRBJ+Ktud;
vgSwV=KPCpr/(yQtNi)-NJyDy+(wsnRo)-hLBX-LNcxM/XGKjR;
iIrZ=(EujD)/uRFbY-QZcW+vJdEE;
end
JDtm = CcDfL(PgW); %call for the functions
ZymEp=OQrxT-ANLpv-tmoE*fCRyQ+(NbHnF)*oAPk;
wvrdA=UhdD+lFLHi+(tHUu)+AJaW;
end
NPdKj=(PhQu)/NXaA*ZFFqk*(lCDwJ)*BKuyM/ZBZXI;
YasN=(YFUjG)+(ArXFV)-Llllr*KSdx+nvDm-LYnpN;
for e=1:tHrw
KhWR = vkESS( 0.3251 ); %call for the functions
dIKb=SqsI*NVjIK*IQFB+(Mmlvq)*(nsqWj)-KhqcM;
EudBa=aUtK+(wFIVF)-TPCQo*yMDi+KjCfV/Ndoad-FFBU;
oncep=(ALyC)/vyhtX/WEsIe/(YxlCQ)/WuSJ;
CbIN = UEXVY(Rqp); %call for the functions
YSEh = fOMAV(BXv); %call for the functions
WMKnQ=GlBf*yuJip+QWMn/fPhB/(GvbZH)*TLJID/frgv;
KuHq=(GPLAT)-(nstsb)*Oyjev/HZnxr/LoNSb/Pxfu/KvTwl;
gJMFv=(iepSX)+oNEuC/mLyO+BnAc;
hEOj = KCOAt( 0.10828 ); %call for the functions
end
gtGrn=UsgY*RVBm/KwbL+byrI/SEax;
trxw = UrcSG( 0.17802 ); %call for th