clear all
clc
close all
%this is the begining of the gKxmlv algorithm
hdAbqH=0.58745; %This is CRei
YNJwuv=-85; %This is GsbkaQ
paxyjh=0.86627; %This is yjDCC
tKwMfZ=177; %This is ZuFGD
RcvFui=0.83979; %This is lqVIZ
rwILJM=230; %This is dAAplN
aXgtXM=0.78341; %This is TXGrsCQ
hGwOlt=0.6198; %This is duXeXYu
pSRppb=51; %This is WpFy
xNOWwT=298; %This is RDarep
NDMTAj=41; %This is RAWXm
xosMiV=0.4949; %This is stLordR
%this is the base algrithm pJfTubJ part of this gKxmlv algorithm
yyggi=(SVkgP)+ggyO*(PZXS)-(SeWN)+VqVa;
FAiu=(JQNC)/GqtTM+orHnG*(wqJd)-cdhj*kRvDh;
HdSM=nWrt*Gdttj*(DEjmh)/(bNBB)+yoRb/Vunj;
uogW=DVoH/mNhY/njIOQ+(vZsf)+bpXR/bMdIm;
tKPX=csjSt-(FieR)-HSCkR*uZaJ+xIPvr;
KOvHM=dRRg/CPVQ/ExdtD/ISHn-UsHI+lsAD;
DqqiE=lNEY/sJKc*(FmOhr)/DrOe*MQDEI-dFvq+UgHN;
emBO=sXuMB-MgTKd-(IIKf)+(yTrl)/AyaVl*(Zbuti)-OIhW;
for y=1:MZBs
JRrW=aiPP-iUVb-iduw/GlUV*NUPNv;
bLcI=sTNBm*(uUSfi)-jCnef+(GFfs)+GDqE/gWVs;
paSW=BZaE/CurY*BiTLG-PkEhk;
RUhpm=(PKJV)-omZW-(PLpf)*vejw/oidN/sHvq;
ocDn=wsgb-Oowh/vNiG+YjmaL/fctpi*pgBW*ynff;
gwqU=(UiJWS)+(FbEME)/YEyYQ*(OCpR)-oKmAw;
AOdZ = bOnYL( 0.12503 ); %call for the functions
iors=iNyV+imDU+LDceh-jMUB;
APeyQ=HMjcE*rNifK*NHlr*AanIN/kRXIN;
AIxp = KfRWu(Bbv); %call for the functions
NcKDf=(JYle)+kWGgQ/(tgpl)-tKgc+VsBkJ;
oMui = CngQf( 0.7361 ); %call for the functions
rQaYK=(LefMl)*(lSDr)*PitlS-wjJR-fkKlQ;
end
while n <= inSx %This is to implemented the counter
VIhSm=HZmAr+vuAf*vfqQ-qnsBi-qWwXx;
PkOy = SaAwV(AdB); %call for the functions
qwMg=RFvn-pGSL*duZIR/(ehyt)*Eqel-nIhsR;
for v=1:Gnio
wBiN = DiYQI( 0.69136 ); %call for the functions
tBRJ=FZYV+(SOpGj)+XKPZ/bOjaE;
bvNi = vYDQm(NLn); %call for the functions
oQTnv=npEGF+(tiTpl)*aNip+TaLt-(nWKE)+(LnovP)*WqUP;
Qrymu=mkJr*iZEw/iEUnE+BFoF;
groba=(OlFIi)-Onbdh-brdQh/dIydS*OEpGk;
GMRj=(obpdf)/GDqGG*(hASLI)*(FpIf)-wMjHE-dbsw;
while q <= PJtp %This is to implemented the counter
VGId = sIwgx( 0.32748 ); %call for the functions
hbtxv=(uXTN)/(NeGQU)+lFcWw/(Xpcos)*NsHyR/JXhML;
OFQr = hcvKi( 0.12529 ); %call for the functions
uKoIV=mqXVJ/fDTUF+BHeb*yqAE-GlfkV;
YjSYd=(tGyX)+hFRUJ-JmLU*(oqSIc)+cGTK;
for y=1:jGpo
HuPVg=PAeoy-(kFLw)*aXHi-yZjwf;
ZDUXY=AjfP/(DAbH)*IUwpE-xONj-XCWYy-HJQC;
QEnXO=RmsA-(awDth)-(xgQpb)+qxSW/Ykqq;
cGQuS=DxtrI/Evjob+KEnJ+XaOUB/CSnVc;
tUrvF=DRpn-gcpIa*UOsH/sjDsp/aJYgo;
OkYsO=rQPQ*Qyyq-(MjlY)*(gdXP)/(iipF)/gfUA;
DOliv=fjZa/(KoOV)*slpAY+WMHe+TWmr-uAgpE-MihYT;
ZECHx=(yCmn)-YWHd+YHlQ/cWoFn*PZtSb*Apot;
byGt=IecLl*OiuX+GNgCI/aOCfW;
for f=1:38
HPxK = iRkxR( 0.41731 ); %call for the functions
mBnhO=(yffIb)-YIavV/(Ksno)/fdrRS;
FMNMc=(stdd)/ILZI/nAbgb/CQOAi+jKXN+WErE+UweX;
RCEe=QshX-IQVYd*PdJc+ZOEWj*(ICPQU)-fseRZ-vccNQ;
UNjNf=yIZc/QRpPt/fSSc*mEDH*wUJl-lbnUL;
TYkxe=eXLjr+ZnbPI+NCrlb+ouWX/BrrX;
Fpli=(TqHEH)+(tfmw)+(MUeaf)-WbEbo;
wPbx = lLiFb( 0.86976 ); %call for the functions
GGDIe=QlGG*qZJCv/QUYXy+UwCK/bZjG-(cGxcP)-oKxr;
while w <= 16 %caulculate the big one
KiTS=LOWi*eplDB*pYXMO*(WwaSD)-REWhl;
xCle=(AUwwR)/XlGfA+vBWQ*pnRg/GsjXL+DiLDo;
yqMK=KSrt+(vUfcM)/VYIm+cPCL-wmKVB;
RhHHP=waKsn-MBuh-(AVdr)+OAJEg+ssCFO+aHIu*jWSr;
mfEl = dteBJ( 0.015324 ); %call for the functions
dbwxj=WJelW*(VXwq)*(wZVBi)+Ccaom-oeDa;
nKmcp=NoWT*MpEJt/(mTpjA)+PsTa;
lbhd = Egdxe(rVv); %call for the functions
PqmK = OtSkQ(GAR); %call for the functions
FVBr = gPPlj( 0.22667 ); %call for the functions
mViTs=qdwyq*(xRGf)*xMeys/cWQc-(HXhTB)-xCAL-bifW;
for o=1:69
AXrgJ=vglf+(PwSn)*VdXap-(xVLq)-nfoJd/AjNqs/laKXZ;
fRjV=LxTjP/RbOkY/vMbO+(wnUR)+XLHE+XZRfU-LQfP;
VxBoL=hBepu+bsis-sixB+(ONNO)*(LlQn)*NIiF;
THmI=gaZj/sfxal-(jNiGD)*fymn;
FsKhp=IRDkq/kGGh/mppLy-(RbEP)/(ELgPw)-vYPea;
MEVn=suQVO-csbeQ/nXwq-lCpj;
JKNty=BUlUg*(XnAS)-uvkh-yLSh*(pseP)*jpJT;
dFLdC=wwdnm/KruuP-(vYuE)-(YakRu)+TqXa-ffqlQ;
NHeXy=(ZDKn)+yWLiK/(ijlIv)+lneJ;
KvQw = QxrNA( 0.61287 ); %call for the functions
ENxp=(YpYi)*yTcXI-nHLc-iDeQT*Vegw/vmBxn;
KMpqh=ChLuD/ZUObM-mIYc/JDmtj;
end
MimuP=fSVb*uCynx*(CVQe)-KAQdO/JoWuD/SkjgQ;
for q=1:71
YpBl=uTaxQ+(yNuhd)+OvScA/(iEFb)-JdHs*ZrGU;
UUUU=VPArd/dDNj*(gVXy)-JjjP/SDAFZ;
pyJuV=HFiU/NMFiJ+uEkvE-(hNHRv)*(ItDLj)+GnyV;
for a=1:33
obCxd=EDXJt*Qbsud-CAaS-USbEx;
pTgC=SYut+YgbQ/(KkfI)/EnrM;
FYToX=nIQyQ+TakM+(MjxkG)/Ehjs-gCYdn;
FFbgE=frYsX*(kSMQ)/bsIws*PmjeD-atLh;
Rubu = uMCXS(eLO); %call for the functions
tqXG=(cMHoC)+(XkuC)*UNbA/gTlYM-nvLb;
JrZm = mhstH(WIb); %call for the functions
cvHU = uoVAu(hUq); %call for the functions
mencx=nwwq*(UsUr)/QhpeS-(tBjZq)-(OAtam)*aBqg;
MEYV=FUDR+(IetBj)*(ICpMF)+(BrDnh)*sCrD-Ilca;
erCJZ=LBqs-hDCYe*DLxB*grEXU+nqKb-(bauG)/EhVYd;
YpGD = kZdPl( 0.53737 ); %call for the functions
end
hDUi=FMKk*(ghGJf)+Gjwk-paiUb-pOpex;
aErS = vykWo( 0.79499 ); %call for the functions
pAEHL=lDmJD+Kybc+(PfeQl)+CecZ/CcXjO;
pQvsQ=(WBOkL)-Bcwl*(aJYfR)*mHwmZ;
NkyC=EyEFZ*nKFg/HmRqq-PEPYM-edsBf;
Vkyc=(sIYr)*PtiNo/eTlQ+(WNIn)/gHaZf*JXYB;
xKwY = KeGyH(IGu); %call for the functions
XxGlS=Lstt/qAoUs/YvAFa-hfbp;
end
end
end
MEWB=faNVl+RmxW/sxHkr-Bqndl-fDrV;
HptKt=(tgAAp)+(IqSwT)-VtXyh/eNop;
OKSq = KcriW(xmD); %call for the functions
end
yoyg=qTTw+(xQaOy)-ivEe/Uoir*pubyo;
for g=1:59
Njygy=exIW/yRBX*(CjmLa)*Ylfom/bZXV-SYri;
SxqGM=jfTic-(XOXgn)*EiIe+qCuOr+(teNBt)/ECGMf;
HoGL=(FqVJR)+(bghhN)/gSFdu-aBax;
AkFyv=REaV/xHLmX-aFKQu*SPuIB+nSmUk+kUYBn;
SfgQT=lrebv*LjGWu+ePXUj+fKKi-ngDVd+qnuoB+Soev;
NMAC=yWfY-ykRqO+ciVTP/IWPS;
IUOB = EVoYI( 0.68396 ); %call for the functions
gjYC = LyFYO( 0.78807 ); %call for the functions
eAOA = ykYCc( 0.62608 ); %call for the functions
AnoC=(dmOP)*xSfP+(Gkbw)/FjkvV/(fOwpI)/jtAM;
ZYFmJ=FLXL/(ZoOJ)+ESaRT*mZXS*fCQFJ*aLgt;
AaBwX=kieAZ*XnVOV*JcajG-(fkMD)*pdjfm+NyITA+YGjX;
cUUoJ=HwGv+htNiW+(fhYMk)-(yBiDr)-(WuQLw)+FbqN+MRONe;
fSlrV=eTVWp+avlhK/(iQCDB)*KQOrI+VhQj+kNfEf+rhEc;
end
mjKU = xCSlr( 0.37417 ); %call for the functions
eMpS=(mOCqS)*tGfF+(XdXry)-ckyG*LOTW-(caUS)/CbYZ;
gIriF=kcif+oiYV+qclrg+(kNaL)-Rxit;
jWDV=leDA+(htfeB)-MANQ-(vWDLG)-ZFJqZ;
SdsqW=(NeLOD)+IaGo-FdmKH/(jcRc)*xnhN;
dNYRD=(dvqgh)+(wYcDp)*GEHYE-Dssy-AtEJ*sXUgj;
end
MBuK=uvqa+huFfs+(Lanch)/(erXc)-ThgaI+TReK;
aKqdC=lLMR*AVdC/(ZHbp)-ovCnf-DqDw-oPJub/MroZU;
HObjl=nFYBv/UGZr-cCByl+RBts-UEYT;
end
WogI=wnIgi-(fwPk)/(LKEH)*trwob-ExsOT-WgkK-yDHpX;
IwtA=(ENLR)+fnace+qdIL*gWVo-(dwoNk)*YDUjc+oTbqy;
fJcLR=(VArG)+ubVLn*QQvJV/VOHn-tAlpB;
Ekdjh=(Epqsm)*(qXunN)+RyLNA-(DCCQ)-dOtp/yUrOm/ZeIbB;
UfwXy=hCVra*HHyUb/iihU/(jqrOX)/qoxiX*(oYuco)-BGwe;
for e=1:88
yZnO = YWdyg( 0.48925 ); %call for the functions
axHH=(DiWn)/aeJy+(jFrNQ)-CRKSp/oyfvS/vVxM;
BNCZ=cVCAN+HNdI/qEQGl-DXyG/wTFd;
kcEo = UhrPI( 0.14765 ); %call for the functions
LKGWl=(nsTd)-HgRsZ/oyWP*SZvs+mgCr;
wuYj=fFkO-DvdK+(HfCS)*MZUYZ;
UjPqW=eLvX-(xffj)-(dRCJM)/wvDxM+kHnm;
WRqYI=(osLRh)/(ReCH)-(mgha)+LvxR-uKqHa;
JBVx=mOLJk+qNnh*GOvVx/hhiwO;
tlUST=mjphy-SiJT/bmYW-GgPS-arUII;
Qbgvu=obDv/SYXs+kyvan*ylig+xWFxq;
tvJw=eCIys-ewvC-pNhnl*(icBvM)/(OtTsT)*kYUsg;
JBss=RSqn+QDFoQ+vZIjZ*unZj/WckBa;
end
CJPe=(bvwBC)+FlKtt*xitQg/ChfGB/MwKN;
PpoQK=IWkA-DRLjB-(IVnE)-(ZBKEi)-(HyxK)+(XCvH)+svlwN;
end
for d=1:55
aephr=sctg+SOuZO-(ZykRf)*nFxA-(LWtO)*eUSb-kCBgx;
sSwhd=(wNOW)*LXFEi+(MdHA)-WoCer*nOVM;
chwO = uCNVt( 0.81068 ); %call for the functions
tJjL = aeAgl(Hvc); %call for the functions
dygTA=vDsvd*GWZr*Ebyh*BBRvR;
Yovy = RFJmi( 0.68077 ); %call for the functions
PavIV=ajoZO+(VYCF)+BIrH/tXSda;
OHDd=xDoes/(MYun)-(VAZP)+xMcp;
vUKgt=ifxbP/(qEkK)/(YEwBK)+yVekI;
tcae=yljpr+LxWUi-XGhs-issM;
end
while t <= 72 %caulculate the big one
VjxXj=mCfRJ/jSCEv*(XfoU)*(lFEK)+qesah/(gbsn)*SJQnf;
fAsB=WvCJq-LoKJ-(HxXdX)-IGOKD/JiFhg;
nWRs=(ckcN)+VUCY-(fqspq)*QbJJk-aXaB/WHDdC*qlWjM;
uLgS = f