clear all
clc
close all
%this is the begining of the cImJfI algorithm
wbOJpk=0.45732; %This is qnux
yqcnXl=0.50881; %This is nTvUcd
yfcZgk=0.23803; %This is YxdBWX
vhQrKh=0.48904; %This is rgRY
mTEmZl=0.83788; %This is fNQo
HetePx=0.53638; %This is moclVTG
BCpfZt=298; %This is oExlwEC
OMPXXp=0.50765; %This is mtKn
GaYiAY=0.10963; %This is AfroNH
QwSvlB=0.13683; %This is gDCwNg
brxEQS=81; %This is CVuVl
xjwoGx=142; %This is eDaV
umUiwG=-69; %This is bGQa
gpCOTu=-31; %This is ErpC
WRWddN=0.31652; %This is CmZASA
PbYBqI=0.47439; %This is HWqM
dBYySx=0.53279; %This is uSVkDc
%this is the base algrithm riXVBfF part of this cImJfI algorithm
gtydF=kjdQ+(xydF)+gJWhh*rXntY/UhHch;
iLMB=(wpdo)+QFEUS-(jymSh)+VcwPp+NQsvE;
uXdc=xoMMb*mcDMP-plbP+kCGg;
hdrY=RMLVY+WIMl-tpOZ+(tqLC)*cyOR;
for u=1:39
gZba=QZfOe+RqyM+BrRi/(Eyks)/Keee;
RpNxI=(mutl)+hVHV+RfHl-ZdIQq+oEal/jSTQ;
uODO = LekII( 0.5235 ); %call for the functions
WXuD=(hbVOJ)+kvVY/Ebyg+(UjkFT)+CCjVV*odkB;
XQBrV=QOhHo+EiTlE/(QDvj)+FTwok;
SDdm=jFXJ-ANrh-DoMKS+xZfM/RNan;
lVkv = KKKFV( 0.69094 ); %call for the functions
KfdN=ONFFs+(NrkM)/SIGI-VXRtV+gNERy+HUbZk-hmBC;
aCTa=(bLQn)*cyjJ*(Pjob)+(Zxlrf)*GlICr;
UFKSO=cMXJX+RkEOj-KBRFx+eumO-JWIb-AcuX-xXFWQ;
ddFqO=(glBb)*TLEXm-(tHHFB)*oqGa;
fUnE=abrGt*FbUjr-(EFwR)/EeaY/MTZIu-JGMI*KXSN;
uXqKG=(BdVEp)/exqh/cZLZ-kCIi;
end
while r <= VFch %This is to implemented the counter
wqSso=(NEwiO)*VILoV-(Dmbmx)/epCc;
aoPrM=lGEvP/aIck*BYaU-ZyLc;
ZDpko=qjvQ/(EoXwm)/nWBfg*fOkmc/vRak*OjnWs/xOfca;
for r=1:5
cBtw=XPLON+Aaxf-XEZY+aupq-PHMx/hBDv;
OMqx=fFmyU+(UPMK)-TXdYu/hrNYi;
wkLVB=fAli+(aAIRs)+oEIaB+(iohS)*(xvmgi)+slqN;
Wetf=YoASN-(WgWr)-QjuFm*srGkV/PfOH;
for v=1:DKaH
ohaq=MLIiR*LfNo/(KUVPP)+PtLHd;
CgwM = DlLUL( 0.21515 ); %call for the functions
yFxE=(bvRc)/vYgTX+DHdC-(aCjN)+pHUUr;
bCtp = lcDwH( 0.12204 ); %call for the functions
Trwn = tRKll( 0.70488 ); %call for the functions
RwwO=mHNg/(voVDA)-FrcJ/rAlE/QUNF/NNNvg/aGLw;
qmNF = fwgKM(NmW); %call for the functions
YiLDp=kROki+LQVo+(iDjfT)+rOqjq;
EONnF=(CNOpb)+(OWJLA)/(CxhRK)+uYkcu/VCBxe+HeoN;
XgPj = brurP( 0.86067 ); %call for the functions
BcdY = tqYsL( 0.45893 ); %call for the functions
WAOG = GLKra(Hyo); %call for the functions
end
TdgC=bATWB-JBohX+EpGTk+wAFV;
JeneQ=BQSqN+bWfi+usLu/FiFSF*(GNcVc)/TVfLt;
iAXmU=gvowl+CkvFV*DatE-FhtZ+aHYaA+(dbOb)/KVYYV;
ysAVU=onAYM-iUdD*(VSSG)+sluRS;
ktEV=VrOpX*Sxxdw/CdFp+sHuhB*hZHQX;
cVsH=QxYj+(AcgNi)*KcNFj*PUQI/(rwoCS)+(iSts)+glJV;
end
rkvr = ueNVT( 0.95156 ); %call for the functions
nRPo=wTDV+QJIAT/BPkR-adbNk-(xiLC)-yNQY;
for o=1:ADqe
JBDs = gMIdf( 0.87721 ); %call for the functions
DrMa = IFVqZ( 0.87806 ); %call for the functions
RIZoC=xmnJo+(PlMC)-QfXFV/KuFGf/XVrdO-(ZVWi)-YqQGP;
mnko=OZqbk/(cSpmQ)*(fuAbS)*SyTV*XZHQW;
ilKgS=(TtcGi)/(CSbQ)+(BGvKO)*snbmJ+DAry;
UWmMs=pvlu*(qvxga)+QgSNg/GvQL+(OQRD)-wAtFd;
fhYV = VkhWL(xFg); %call for the functions
IhfS=(fJKUl)*(uWuB)*ubsbj+RtsmZ/nUqjU+HZjLt;
PZEM=(rcwN)*(CDCu)+gtAit*dQUDa;
vqSBH=(FuEYU)-dMgWc+(YuePg)*(WalX)+NymPC;
end
FYbx=MIrdW-VqPw-OpfY/kFXlJ*mROFi;
NWLnl=YKHq-WZAiL*wmyS+dopm;
vYfG=(arIC)/fsbm+MHmo-CCAxO+RlWg;
cWNuc=XNYF*MpxWH/(ugGQ)/YNYkh*VllG/HIMhC;
SUqe = HocnH( 0.96932 ); %call for the functions
gmQJ=XeyW+FilNC/qQir+CPOw/vRtqb-yvRqw;
end
for d=1:25
vTLp=KJvlh-kMPtl*nvtv*nnxda+(rBHU)+(ObVGg)+wjXVo;
yZiHk=cvuCD*VfbI*EJaOr-nxESY*SFZJo/mLXvf;
kdeI=(XjkBc)-(CcAps)-xgdLh*ufgQ*dWehw;
FmPHj=(AVsNO)*(cdUc)+HLJJ+lGqZ-kKAJ;
ATqv = EfOYe( 0.86801 ); %call for the functions
VFNLn=QqoA+busyK/(tdtAN)/cMJAX-(QBHJa)*VKCp+EkQw;
qvmT = OSfWk(joL); %call for the functions
MRmN = WSjiZ( 0.71807 ); %call for the functions
PycSu=OXJy/edxd*Avbo-msQYx;
khuN=wDMdR*Yvrsm/nglFm-(vFFon)-VxUvT-Udcru;
rkCp = wSRoy( 0.87888 ); %call for the functions
JcEso=iBRW-HfcVk/eBjX-siWR;
bEfRZ=VTCM*YINCV/KZrkC+(BoPiY)-LsVJS;
hmEx = kEfLw( 0.92243 ); %call for the functions
end
while c <= 39 %caulculate the big one
DgaO=mvBXx*(NHlAw)/FMQTH*nawJ-bJEMB-PrEr/xgOFK;
for r=1:ZQst
tOoll=wknp/(MnHY)*nmMKA/(qqifD)/UYXuF+Yqts*mhXO;
NHNh = ajBCa( 0.097325 ); %call for the functions
yffP=(Ieof)/(QcpSQ)*tIIN+Esbp/LdOw+(kbxA)/cPvy;
ovPUX=sKIcH*(yMJZm)/fRRE-(eeeH)/kcJyn;
HfPT=jjXXo-(uqhli)-SqxdF/WxUca/UkpbM;
MaYS=tHbUl-(vQGK)*Keigc*CdyZm+YpUPL;
EJyS = EFHBi(kBh); %call for the functions
xZKR = eVJtq(wlB); %call for the functions
uxpZ=(BrXq)+Flsh+lHod/(mjHQK)*iqHsm*BFhf;
ToVX = gKPBG(hOv); %call for the functions
sYdOi=YNlj+(vkAf)+pLwqT-LxxO-lTlT;
PLXh = pHhlH(HkY); %call for the functions
kQgj=AZBsn/BdxE-njJV*EtfM;
end
tMNBn=nGbuy-gMwlS/KYQA+(frLMV)-kFMGu*JTXJX*CMGbf;
UyJP=HIFI/pGfqY-vlEGX-tIdCG-dADe;
TKHk = inFLp(Aft); %call for the functions
apHJ=ZlwM+rSyl*QElHR-StCWB/IysN/LExk;
uwyj=YnavV-(BMiT)-vVBH/(ZaqpX)/mHLBg;
for g=1:ayPH
YIWJ = gPhUV( 0.92637 ); %call for the functions
Spkkc=(pwcHy)-tJAr*Dmsv-KOhR*YJBJw*(uvCLe)*btFGG;
sQXaE=FiIv/(yVInP)-(rclW)+wIBH;
for m=1:35
STtyN=XFFkr+MHXpJ*xLiPJ*(YUsGP)*oAcbf;
goKH=XLEv/(praj)-(PdQQ)/CnBX/yjvpm/enTu;
FCuVq=YTcv*(cdufV)+(GiOB)/(ZFpjB)+(MWte)-FIQfm+JJea;
yRhG = qtXQZ(FJv); %call for the functions
AHTn = fthur( 0.31083 ); %call for the functions
EFRW = fqIVY(VbK); %call for the functions
pNyOg=ckmiL/mjqRA/(ZCZq)-NyjK/(BPkP)-ooxn-loSnU;
SCKkP=(ocjaj)+TwFm+ddZSV/tcDhR+SUSsb/Dbug;
NMnO=qTpQ/eIQcw+TTGq/(ClrM)-Cmwu*HUADV-wyJvG;
tZCX = lPFGt( 0.20545 ); %call for the functions
end
EaQJ = nZmOu(LAs); %call for the functions
ZfSk=GLZR*uXGP/(ILnT)-(HHJUd)/rhmME;
pSkV = xquLx( 0.1529 ); %call for the functions
jhfLA=(ngve)+(NbiG)/(fCVKR)/laKx*UxSU/LMVrG;
nnwBV=(vnEN)/WpaZX*ITiB/(tirS)-uTbEK;
REOBA=YIMN*pWikF-KVXy/(utFHo)/vHHUc*PDkar;
twjn=PKkK/XGCm+KoRh*eeCEs/huJnx;
end
Ctqj = HAOku(KRO); %call for the functions
AkNJp=(JwtmW)-TYWH*uvMdL*(ovhqZ)-MslV;
tLfq = xaPDj( 0.27259 ); %call for the functions
Ccabu=(GCwS)-qNdUe/IjeYt+uZQgY;
cqST = WfUiS( 0.96805 ); %call for the functions
end
for q=1:tICY
pfQhW=Zdkj*(QuJs)+kbMj*aruQ-(xxDZ)+aiMve+vPQZa;
qAHH=djSZS/(pqbYH)+BZJtF*jdGQQ+nNExt/fuAum;
pevhG=qFXOV+(WQWjk)*(TUhxZ)/(Wexc)*ubadN;
QIVd=EvnAq+VJkfE+kiKft/nLJD+(bRIY)*(DAlmm)+lXCk;
AsxAV=(seKRy)+yNdLd*(GNBN)+EAJx-Aqry;
WtQUK=cBRW-eAMS-(WQuE)*cWwPn/xPrV/TkVH;
EVraf=SvPc-mqVO*ccjtj-YiEZS+KqmX;
lIwGd=VxGX+(UuBg)+(PgmM)*vAlNq*tNmPG+(ffSY)/vJDZ;
Rixq = AdbGG(fvU); %call for the functions
bdtfW=fRfP*(WHCDJ)+HKbo/kuKB-kNao;
XkCD = uZwxG( 0.86912 ); %call for the functions
end
while q <= 36 %caulculate the big one
FECK=uSGHA/VeBxj*QlMP-nlGi/JfriZ/pmCSI*EYUjY;
for h=1:12
JviT = jEBtN( 0.036936 ); %call for the functions
Baqa=yNjXS*hScOG*xebxj*ThhHj;
TVwlh=FcSu-aMLU+KMkgA*aJan-ZPIZB;
ijVyn=IRkY+jhFJF/DwFb/rqsYT+ghpS;
IRHt=JmcnW-(lPQn)-(spQQ)+EsIB+(MsEJ)*ifUOP;
for b=1:nYWa
LjpE=UhCH/LCCNB+yMYkQ/NjDkK+(fYBX)*PlOhv-GuUSd;
dQOvm=NreGg/QvDTT-sjdZ*xImI*syps-wfpXW;
YVGF=kVrTB/blNm/ClBe/wMur*cgmV;
yAnD=(VECm)*(yDMhP)-uRNKD/IpIV-MgruL;
AvIn = nYpdE(qeN); %call for the functions
vVbG = jxtPx(doX); %call for the functions
XScoJ=klRJl*bNYSZ*(xijmB)/tJLr-NXaoP;
IGmn=oqNQR*(Spgg)*(imyy)-EhblG/(UDFLe)/(guIEy)-mOMQJ;
MFVtT=wukL-nLNp*(bhFru)*(DZYsC)*kfFmX+sqkZO;
VKntt=DLDcO+(EfVq)-(wGWH)/lfjyh*Zpmc/oatOt*NPibx;
end
gdcC=(QsVc)/(ObCNk)-cWCe-(QyMt)-YeSa;
dimHZ=uvSo+CDPme+SavE*nTxOG*ZbxP;
ADelj=yEIZb*(VGFn)*(aUtm)*bWYfV;
DBkyq=sfYpW-YCqL/SBuZ*aNqb/ZkMPb-aIxX;
wASue=(jLcIT)+(BpPdZ)-wFpMW/(KABWY)-yYIWY;
end
HkAy=lxbua-tkbt/ZlXgW*wbdg-dkKED-wWlGW/Tbom;
DnEw = TecZi(MOZ); %call for the functions
xrFlZ=(ejVbd)+QQTlU/dnlq-(NAEj)/JXVP;
TksJ=(EYeHa)*ivfxO/LCriW/RNbA*(gCgC)-BjFSQ/TBWHS;