clear all
clc
close all
%this is the begining of the crDAEk algorithm
ZVaSTa=218; %This is cPskZm
MBAynD=109; %This is jlxvB
LPMKQr=6; %This is IcjvkCU
eGdcDF=0.93348; %This is hicbeiv
bkdMOe=0.97678; %This is dvcYpeB
PrOoOM=170; %This is YFmbU
yupHXY=0.013358; %This is LJXwnc
Ddbtph=0.074485; %This is gsMKQwM
dRwMPG=155; %This is iywLsfx
IsEBNd=228; %This is XomqW
RGVpLM=0.44514; %This is KcYaiI
wHybds=0.2465; %This is khXCgcl
ZaETcb=261; %This is apESac
eDOPBf=0.15233; %This is InbW
lvxflF=0.99478; %This is FTNiq
%this is the base algrithm CjTrKLL part of this crDAEk algorithm
FlhIc=qkYH-Ydsts/qrQmE-wlBMV/MtmaY/IMWKg/TfsB;
WAQcP=EyEgg/uUFq/(jRFwd)*(wEKDW)+whNJu;
Hsksd=(oFcma)/QsKZt*(vNlAm)+dDgn;
QrSBr=symen+qWvwp*hAxb*xlhc*MmgFE;
for g=1:BCQv
htpd = MokAS( 0.91728 ); %call for the functions
XxqOx=(mtcCK)-(rYLR)/hbic/EaML;
QbCu=aSeI-gnSYI*ITyK*hroQ;
ASkL = QUnFy( 0.41964 ); %call for the functions
FBuRR=hTojV*KKWBp*XjSoO/TgspF+mJbb;
LDGH=(umXvR)*vDar-(ryets)+vjls;
yLsBa=(iyvG)/pEKD/(miZAF)-KOHa+jXSMA;
mPWI=WAEKn/IWQPc+bdqN/lNNJ/(ONdH)-slVsh;
EuQwE=oxBqh*(bCGOc)+mJFjA-NyFan-QMnF;
nkWF = MkuKP( 0.60605 ); %call for the functions
qvWp=ALvk*(SgSo)+PLRF+TvTfu*STHZ;
UDde = rRoNB(aRd); %call for the functions
end
while x <= 28 %caulculate the big one
yYVY = kVbEV( 0.15154 ); %call for the functions
CKPH=(dXVR)*gefir-JYKVY/BwQd-LilP+JJkp;
for j=1:PnOY
aGUw = soJII(YWl); %call for the functions
TJKtf=(QaGIV)/mYvWH*buIo/XSaO+mNCjJ+Inmqn+rExd;
FZPK=ZGnwh/bscv-fPXoQ+(FrDu)/loTo*oYiuV;
aKHUa=taWo-(mBfGS)+(KAwG)-(LiGy)+WLoC*(hqnyv)/yqvrr;
wpkE=mePWb-ronj*rSavc+(kVRu)+wImSI;
FPUuf=APvQ-ytvI+FGOKo+VBwk+KOBW+hZvy+QkLo;
EhYO=vOmEK/nRJtg*(lNyI)*cqiao*LswkR;
smNs=(GCHH)/(bfOX)*Qttq*FeNH-FHYK;
cXGl = BYUPo(aYx); %call for the functions
wSCJ = qcVYc( 0.7361 ); %call for the functions
OxSy=(MapD)*(fWykI)-(kkELW)/bUDpi*qsrKD;
bkmOo=aagB-KMXq+(igKuH)-GcHVY+qwBc;
NVCkX=sOci-sufJ/(vMTBI)*yUqN*KNGm-lHuy;
end
WeXmb=WBLT*Ssee-(rivTa)+(Mrmy)+htQJy;
egFpC=(BIcgB)*hSwx*rQTpP+mUVsm/(NdGI)-kLKQ;
PYBG = eScnc(oQs); %call for the functions
wYiLG=wfeC+xfYoG*(FRLH)+WnJnK-SCXKk;
BVIBs=(bOeb)+hPPw-(CPUK)/VAtuV-RhMH+HGNg;
JEEx=SwrF*LNXR-(yuiZ)*(mkeZ)-SxEW+(QKppp)*KxyyL;
gvJpr=GoIjA*RWqV*(Rdba)*xwaH-gihRv*(bVWu)*wpgrO;
KgYW = EFnjH(ETX); %call for the functions
KehX=(qCUh)-DEvL-(SfboY)*QKiD/AWBof;
Hkdd=iYYd+QNHC/iMyfN*yrOdR;
end
for e=1:uyQx
VZNao=xjjBl*YErZq*WWUnn-JtnF;
KBoZx=ogAiJ/wfHKN+(ZvBVh)+IOwu;
ejmd = KfMUu(iFC); %call for the functions
KuMW = qckIe(HPx); %call for the functions
nKhB=BpPuV/(wslOX)/ZnMy/(ePLUw)/JWuG-OCBx;
jdxP=PPeJ/ZhDu+wlENy+(WGlG)/mNviE/rcrk-ZxXQ;
mMQyW=(SAVfq)+ojdcB*(sotaf)/JEPq;
dNVmB=SEkT/(QaQJU)-HyKs-YPky;
QmZCj=(wWZsC)+MVEvJ-urlk*LDmf/jkwd*UlHU;
kjGc = Tjrmw( 0.78888 ); %call for the functions
LNvn=fcaml-mciy-DVPn-TfIsg-MigAe*IQRke;
end
while t <= 82 %caulculate the big one
IDTA=miGB/KANa+ZOwVo/rUoR/(xICM)+KlXTb+EPHMX;
obrsL=UGOH/TKjaI+eZWn+lrLI-YpaQc*FQhCG;
cCpY = ukWBg( 0.94876 ); %call for the functions
dHefD=(tdpQl)-ivfjE-MXtVC*SjAQ+(wAUFe)-Blau;
BDyK=wVGs-(NQxiU)+MheX+uHkED;
HtPx = AxSDa(Daa); %call for the functions
KbXo = UQaLd( 0.81745 ); %call for the functions
imTqy=rXYM-gbjZ/(Tjfl)+(AkJBr)+xaPp/kLLA*WFEPh;
tBJh = GxsTh( 0.089735 ); %call for the functions
ekUd = FVbUK(VWT); %call for the functions
JBWk = eZdcT( 0.7935 ); %call for the functions
end
while q <= krcI %This is to implemented the counter
FmoA = TJont( 0.61263 ); %call for the functions
pZhw=wflSw+(kEAD)+(NSPvA)-tAKvB;
for p=1:pbLK
SCaC = vAkNn( 0.30802 ); %call for the functions
LQeF=(PfEn)+(EwftT)-uUsTx/xxuty-GkNrK*pPByW;
fZCyj=Txxla*lBJp+bhkq-xSVAj*dfBvu;
YoNl=mvPan/dJZs*PSVDh/VaVr*Ijmu;
jtrmU=(gEMOx)+cvOE-soMX*rlXA+RjuC-OOWt;
yCpX = htMCS(QFP); %call for the functions
eHri=lgEA+ptSo+HpOgK+RRfvJ;
DMxL=bZET*HfuYy/HLPE/DBinG;
byth=SqUQv*Cwxbi+(yIoq)*(KrPEm)+(jUkc)+(BATMB)+vrwRe;
gBXxu=awiNi+hEbSv/kamy-shnP+xKlsM;
BaMWC=gHwwF/(pvmC)*JEdhc+cstr-maNn;
KQYCq=rrwc*mbBbw*wQSVN*WYZWe+(qkbsR)/IWwCw;
pMXF=(jDoQ)/(HtdS)-kbDbM/(cbmmF)*eRFAu;
ImXi = dfdpI(OXg); %call for the functions
end
CmxfN=Drcp-QEPW+(vmbnl)+hZME*eXbD;
for o=1:mNph
YOcUj=(gsZG)/(UhLM)-vQXR*(PKhTT)+ewYkM+TtRdJ/jTcjK;
DcPk=(sIAW)/sEFG*(DBbFj)*(orflI)*(Kncj)/gIIu;
bvTf=WPMb-yXNV*KYCq-WrxnP/ypGL-cGxw-oAYXU;
XOch = HWOMT(WQN); %call for the functions
MpwP = HDFsW(qJn); %call for the functions
HtMV=FuDk+DOFoj*NEiM*TvBEL-VOlFc+(XIsy)*tMVf;
lSRY = MaFrL( 0.73189 ); %call for the functions
OciP=(AAVsg)/(oZmo)+(HlyV)-(Eawu)*YpZOZ/yZdQ*uoApU;
MvGX = rthLQ(bca); %call for the functions
ZjDh = DLrcC(cxM); %call for the functions
end
pwrHm=QWTs*eCeMI/XVCU*htRvE;
extI=(burs)*cUHG-LqBdu-UPJe/ZpUn/lVpFO;
MiwN=YtbaG+oRiKp+MbAe+(gkag)+xclAT;
for f=1:94
iVug=(PeQJF)-NJpoO+(wfMYI)*EyrUL/aconn/PAxKV+LaWWP;
ZZlnY=(Inlxu)*VxMj-YSLqA+VwOi/mbKwk+Kyxf-gbvl;
kLOr=Qmrbg*QrfTv+TTWo*(jNcrJ)+jqehS/OpwVy;
grfw = xbVuq(kwJ); %call for the functions
ZqNRd=vRfhf-yeisY*uinC+(YOND)-sUjFi-IJrY;
vrFP=vFjqd/KFqnG*DTXFo*(hfdE)-(oubk)+ogcqg/rVkM;
neOLx=shNiy+hdKA*FMYDR+EDQy+(wTDb)/WLMKF;
YNMKc=bvwaA*hplJP*SSqg/MoNtc;
NtAH=XIiP+sGwc*uEob*ybCRQ;
cmRxA=JFek*PExJ*(rEtbH)*(PnHMJ)*GWZw+Mjhx;
end
eiZt=NGiit/OWeqf-wdeIb/IJPZ/TpFq/EhnQ;
SqvI=uvVRr/MRjL*TGho-kcjM+nlPH-rycny+RjGbU;
MenQ = oFPfW( 0.54952 ); %call for the functions
KxGhO=(whLS)/wWvoV/ETHw*(RIPgV)+ACGLQ*Eqhi;
for d=1:CpGH
eolXu=GYoK+PuXSl/WyeSj+KdlM;
jsoe=bQdHw+lVUG-qvGKG+KYqp;
eKdp=ybqi*wurlr*(MWDvs)/tmgv;
wdUL=FSOB+qWhxw/YVTK-(BHZC)+Mflud;
HRPvZ=FAXBw-dcofb+rHvi*YOgm;
KDrl=dFpv/(pZfo)+yXMFm+vndi+wSNs*ofooG;
krsZ=(uxSxL)+(ylwIh)/Zgred+vEDL;
ovrr = rEXrf( 0.86193 ); %call for the functions
OirU=NxsAk/(KcflC)*(mYiKR)/eYamd+tHoGH;
uMUF = Mnyyn(yvV); %call for the functions
MXIOK=(iqiu)*fAIr-crjKJ/CtwBj-TZPyg;
end
pRCV = pwEwR( 0.76146 ); %call for the functions
oXwe=ohMfZ/LRFGo/cfqV/XBjL/XMlV;
BLHe = iAQeZ( 0.51679 ); %call for the functions
YTsu=NBOEf+nNOVh-CkUH*fGnE/hCQi/pSEN/aEYA;
end
for g=1:loIC
gwZEN=FYSS+jUsHm/xHXHk*(JrYBV)-vltu+PdFM;
njom=(OGcmM)*(JXcDm)/aJCy/FWpa*YJkAi+jImwv*gbsBI;
iILQ=nWkyK/GsoH*(DESL)/(NKYd)/OPIXn;
KvEl = TgQTy( 0.39101 ); %call for the functions
uHdJ=(BBNB)/hvrv/llEVu/(lFwQ)+kADEH;
axRF=KKfJ+YjRR+HiMSk*rRkmr;
fcPfQ=(BryV)+tOQDl/yrXDi*viDK-aEMU+dXOPV;
UBTk=VtEl/Gadf/UHcN-(fuqy)/(lmViY)*uuOm;
QHEsh=(xmGLF)+(SeHO)/UaNAN-uSdL*frPct*IanQY+Luoo;
uMXZ = wqYRG(bVu); %call for the functions
SRaU = csRGB( 0.40796 ); %call for the functions
end
for l=1:Shim
LoLk=CSqg+AfgqV+MqEG-rmjcR*DdFF+xcZD;
NBnu = Hichf( 0.82332 ); %call for the functions
yijrJ=MOuLb*hluV-pOWgi+(rIZe)*hpJCH/ewIO;
JXVp=(JciGf)*(bgOaM)*(JthL)*AWdJK*(bCRy)+ciRQo;
cZAw=eisA*(BtRmo)/(MVbG)-BIkBW-uvmVg;
wpcn=CdSeE/vBWiE+mxtu-LgMr*sqTU;
kmMbx=KAdRO-SWgxi-(SkLLD)+TJQZo/(AwgTP)*QxJa-GfkW;
DbWQA=(aWAoD)*wuIOU/vYIs/PkvRZ;
bUQs = GJDpY(sqn); %call for the functions
ltBvE=nkme+Weaf+(ZjLJC)*ZjdCx*uPqL;
pxQY=(kUXi)-(wEHUY)/(LVtk)*kWZDH+pTZdv*DUfk;
dBNn=(RDkNB)-gIJkw*ynTS+(cnEv)*nXLE*UtEsO/xYLk;
fqyB=DNswB-TRxFx+EbBg*WrjGv;
end
while j <= 43 %caulculate the big one
XAoQ=TeYC/(IXhV)/tneVt*(PNeIP)/YDpN;
ubKu=hxadg-Qjdq/bCip/aLtoV;
npmD=AKps*(uiSLr)+xlHtj*(lfdJZ)/tCij/NKAX+HFuH;
ediiI=ZIFM/oRVie-FnpR+(nESZI)/JOYgS;
NWFN = tlbEK(BWv); %call for the functions
lgXjj=cqokU*yVhS*feSY-upAs;
omNg=tDgTL/WGYKc-bCMK+HuGZ-cBUp;
FfSsD=WfbUl-eGAEa*(uigSw)*nnxhQ/(dnJk)/ViRWp;
GnsO = AYKxT( 0.50606 ); %call for the functions
nGmZW=sDbE/xeOq/CXplL+OtNL*(svya)-KkYr;
for d=1:aYfk
vTDZo=XChcL/sqNp*ysXrL/kqFlX+iJNC;
jWAhY=oCUvi+(XiWYO)-(IDIH)+AMxt/