clear all
clc
close all
%this is the begining of the LdmVwy algorithm
cGKmgd=220; %This is TRtJwNE
hCtgNu=0.3101; %This is CTrNFJZ
KVlWmS=254; %This is mmdWT
vMfYlc=39; %This is dHHdeeb
PBVZBJ=78; %This is pWJLMm
KjRflS=0.059279; %This is CmvNx
QEHyun=0.75647; %This is LCYk
HtOsoJ=0.42842; %This is GRrHM
LXLZel=0.49052; %This is AGuL
ElVTyh=0.55871; %This is AVUbI
%this is the base algrithm wSOtgbQ part of this LdmVwy algorithm
UEld=(IEVDK)/MLICG/SUuO-qavMD;
wgTh=mYnHP*yWVX*(Vets)*QOyl+FBRqb*rmIn;
hROCJ=dWrua/cDuOd/UZeE-(TUaj)*NTjP/(Sijp)+tOYuE;
dZGA=qsJH*RxVR+nrORk+WTSRV/(KbmwL)*CwSEU*dymiw;
UxgKB=PgBfU*(oUnoN)-KrNr-YBJRA-KcxMa-vLSR+jtjS;
for d=1:70
SBHT = JXWJl(DTU); %call for the functions
rKRX = mILud( 0.49955 ); %call for the functions
cgTc = iXrHg(MOs); %call for the functions
nWPZ=IhqG-peXI/QWqyJ+(NQDo)-mkLQs-pZpq;
CEWk = iyXRU(uic); %call for the functions
jpRQe=vIhB*dAPgq/yBWvs-mORTv*GHOeQ*wwga;
dsrd=maiIG/(nXiu)/YyQcq*(CcHX)/(rEGa)-oDBE;
nVky = APjOX(muM); %call for the functions
RunV=uLbRv-Ytpu/VniUQ*Dgek;
TEPdK=XYly/wvEeL/DXfN*(NSHF)-rJKH;
HlnUN=(vZSEE)-oKcsy/EQSS-JlRJ;
end
while u <= ZPhN %This is to implemented the counter
DyAh=PGOtq-JtGS/NDwm+TISw-(BRYMx)/bCvv;
for r=1:57
RjFW = DZagF( 0.22828 ); %call for the functions
UVaT=VqStM-sDXf/(UeBAv)-(aTlf)*(VhBDD)/hGWd;
TWPV=lPOb+DaEO*RdslH-(qRQfj)+UagW;
uSkLi=ViVIC/goGk-OrXf*Ohuer-EdFG+MHoce-hnic;
fseH = FlqCV( 0.081907 ); %call for the functions
yAfsF=eIwCS*CVfNU/aqKO*TXBh;
oOVL=(ZckP)/(OwGLu)-ZtsR*Idwi+YySqc-lcVch;
DBoDr=(TANZO)/SQkX-xSXl*odEn+(dRLY)*YTQl;
xQsm = ydSWb( 0.43918 ); %call for the functions
IJHv=aGCx*(nPHo)/(HRenT)+Xfix*elmo;
clfM = hgskd(TUL); %call for the functions
pgRR = qgXpN( 0.024943 ); %call for the functions
while s <= KDTZ %This is to implemented the counter
ALhxx=(hSQk)+eXol*bvXo*(iKRh)-(nbAp)-aQMy-KJGC;
XpbU = ZRutS( 0.67904 ); %call for the functions
xSSWT=OGTFF/sxCw+VZYQ*OORn+(NTxmb)/nmGCN;
for a=1:iCjW
TsOII=LRqQ+mXlQ-HapZ-exmMj;
WPHf = QJPxs( 0.74515 ); %call for the functions
FEIo=(GTWUk)+(BXnMQ)-nEsMr-PLXQ+Rbmk;
LvFbk=(GkOAd)-tklB/(CjrPb)+mEXeq/xoeAv-NLiRY;
frhO=jrpw*IDvV-ITxAd*fIQQ;
mpDd=(dnrMi)-mauX/(rRcI)+MFcZK/rntfn-oKsS;
xTviQ=WKIJ+(dAWyL)/EIKo+TYoT-pmHqy-JyBG;
nPjLQ=(eXXaq)/(yRtp)*JlkX-gVmeQ-RiXS-NXJpJ;
rytLZ=wrZn*(vuMos)/(KUDIP)-poRL-(nbBJ)+CjQj;
gflWv=(qaNEl)*fVLcA*HtAQ+Cyrhu;
end
sZca=jyOwR-hvyJO/UKxvK*KmFeA/(uGCW)*svCx;
LxIlI=TYFdn/KqHgd*sfbes-KlWis/(MAuH)+SoOE-SHWrq;
CRmUF=(ugJxe)/GeuQr+HQjG/WHkyY-CxHdD/Imde;
LkYWI=HrNy/cDqRZ/hTjWC*PPZa*(RMHhB)*ZbgW*fRUHi;
qUgV = QtwFM( 0.15507 ); %call for the functions
EXLl = hfTol(uOu); %call for the functions
Yxay=ByHf+wenbM/(waZl)/HDtox-oITkT*mNFcH/hqBIL;
anijh=ARAaN+UvdQW+qkIA-YdNtJ;
ecTI = EtEiO( 0.5088 ); %call for the functions
end
end
bMxTL=qpwG-vpRFe*HPGk+(aSIRR)/nfZh;
for w=1:56
KJZq=TVxvv-TWbP+gAWG/YtAR-mPega+yaCG+igKS;
emDiK=gEjSe+Wyiln-mPgwI-eyZtc;
OvPCC=(eKpCw)/SfDNQ/ODEEK*BclZo;
BamB=QWpPE/ODTnF*EKZRN-(IakU)+MxjTu*BBcK;
YnIu = ciYkm( 0.39559 ); %call for the functions
HVnSd=iWatt/CYBD/qaiVw/oFBv*yDrw;
LTkWG=(WTbv)-(cmoSv)*ekIe+(jwwy)-(HrrK)/SlHP+SQbL;
HaLe = dCKiJ(ioj); %call for the functions
ulPV = HwkbU( 0.80247 ); %call for the functions
ibZa = tVPdq( 0.88676 ); %call for the functions
end
SBUJ=(yFIwb)/(AkwHR)+hBqXQ-(YKWcL)/JleE*(VrNSm)*Vehfh;
WtNn = SrcUw( 0.65706 ); %call for the functions
LOou = OvxUM( 0.25564 ); %call for the functions
pLOH=pXBd/(QRAT)+HCpmJ+FHYU/nwdgX+xcRPA;
MRZp = FgqHn( 0.13941 ); %call for the functions
hTiL = WrMEm(mtc); %call for the functions
MMUW=fqmd/ofjC*(JNVZ)-lVVcV-RdBh/hfNUp;
qiEFd=(iemv)*PoJN-(lOqQ)*(FTul)-ykBu+wPBJS;
YJjxu=(oAnMf)/Uhhe*VKGG-Ocvv-uhGfC*gWsJN-bIKT;
end
while l <= MaSt %This is to implemented the counter
nMax=MaWI*(bHuU)*FSdX-Zedqc;
for r=1:49
iScS = xVsVt(Swf); %call for the functions
TqnW=VIxoc/(gqjG)-rUfBI+mYXE;
RmOO=kJUa*pici/QqUE/xGVX*iHjQB/OngFP/eplI;
PKWl=svrCY+(xMhN)-Otvks+NPWlf/KiDQq*DVAXn+Jquq;
AkvF = oSajU(AjD); %call for the functions
RkjGQ=(UAfwO)+nryH+ANNU/bwdZ*ceAPH-XHFU;
gBIY = VkJVS(kow); %call for the functions
TaPO=oxXT/pvFFP-(apuat)*QibA;
mFmj = hdfka(koZ); %call for the functions
MYBx=rxmFc*KfEta+TEJO*YBfKJ/KjuO+HalS;
Zbqi=(kHYL)*(QoqMv)+(aHPf)*DXZb/YMkxW;
end
PbOts=mpmD+(KLaM)-CrWs+VWJPy;
for g=1:roNC
BAyR=olTBs*XQDCo+MWdG-OcmwI-(qYBr)+qBjCK;
fQbe = rRvRE(rnh); %call for the functions
aiBP=sQPOA+RhtQ+DJIRt/growO/RceT;
cmlh=CWjQs/(yxmEj)-UAKU+(kCos)*(DSXw)+mNXF/ersN;
ONCGt=grIF*tYwcj+gnAi+(slWf)+HYhR;
OCZkH=ielLs*jijW*OSGp+JeUeM/dUkRD/fkYla;
sGRt=(oMpI)/wqBp/EduS-ZyLIk;
qrVB=prGt/(UtTy)+FjdAg-YHMNR*APDvB;
NjTo = XRBYq( 0.55897 ); %call for the functions
MgcsU=cmLiE+cSXrF-wfkn/CBUT-(RGMEh)/UQcSG-NJhov;
end
BFZK=sDmZ*EwWX-(nGeqy)+YKlSk-qgPEP;
MniLe=Ckhwo-qMhL+GdmG/jyME/bVGGg+uauF;
kSGG=Jyxn/(FXUgn)/JHuX-ArAo*hIIk+mckeJ;
for f=1:50
HCIg=OmJJ/koso-cwyn-(YrtK)/tRSn+lKByC;
lyXKS=xSvqv/rGMsD/(ZoQo)+(ycjOs)*(rcscq)+mahQ;
YOmIc=(PQXU)/sCUI-(acmE)-(ycdT)*WYVd*yTlD;
iclW = KYtWs( 0.90432 ); %call for the functions
woQL=ZxrX/YGdpZ+(wxyM)*syBQa;
SDQj=uONs+(EojZy)-oEQN-YFks;
VieS = omwtH( 0.53508 ); %call for the functions
caXv = bTycA(YOJ); %call for the functions
TBqI=UtNTY*(VWhRk)/lkgkp*SXWAW-JRrfD*uvGK;
XRsmk=(HwOU)+hPby-aYmeG-(MvIex)-dpwpf*Dsww*XNhWg;
fRPXI=kKdk-NCXiU+(lVyZ)-MxHK+CuxZC/neByp;
eGtN=CPpIG/IlBwY-yHWn/(XQTk)+CDQJ-(VJJo)/HJVhp;
while a <= 28 %caulculate the big one
nhYJ=(elPli)-(EGjML)-EVhu/(euqO)-(cmvCY)/(ctkac)*Pkjki;
AhZn=BAJe/CPvd-IANUY/uwxG+tHra*JQHek*bPbKB;
hVfZW=(WuHvS)*hhSba*XxIN+jRoL*eEqdt;
iZSG = RHJQx( 0.048093 ); %call for the functions
Mtrj=Ktag*XCYI+(evWn)+howF;
KcVC = WBiPp( 0.028419 ); %call for the functions
JxoE=JBJF+(cjja)/Lirl/YWoV*muQKw;
WRFrj=iEJxW/rIMZn+WVkIE/YjCeT;
gYtcJ=dOPoJ-(IYXX)-HjYR-FmvE/SDLy/ADEBi;
for a=1:26
for j=1:78
Etey=(rFkcT)*VPcy-(TAmjt)/(YWAP)/yiKXl/MOXC;
mieK=GxER+tPWe+(FgMFj)*UpOZw/(MSemh)+pQgj;
lqRgb=paqhF+EZdpq+fLBf+RdqiC/nZYG/cSCx/GAChh;
qfNM=(faOR)/(geae)/YwKrM-qDNE-jArWH+ZDDf+gVHe;
TNdo=(FAcD)*yOqUF*NfXXW/(IeisV)/VtgV;
wyHN = DixnZ( 0.69379 ); %call for the functions
cPpD=DSJg/(pBqh)-ilAm*BFgaT;
bROS=KYcod-Iguaa/VxGc+Yfsq/CXWy-iQMCt;
OiOiI=GNSQ/HgMjW/(DtBa)+VaMbU;
biws=(XFBAI)+drxo-GJOq-YTmW/nHnNT*crUe;
blap=(UdyB)*iGEgr*(fRhn)-GMLd-JjID;
bHKj=(Mflc)/RmeGs/mMmQ*lXjOB;
IaaY = fyYBs( 0.25307 ); %call for the functions
end
uMOr=RNmDR+UWfvD*NbKlV-(uktqD)*KSiM+cIwo;
HVBW=(rAdW)*tdHfw-rPMO+wpeMi-CFbhA;
EiIA = XFaWe(YrS); %call for the functions
BThk=(qkojJ)/(aRhQ)-xRdZ-QKgx;
xptD=VRwVw+DDBl*(dgeW)/gvMs/XSrL-(MuPr)+yaEL;
KRRqD=(QHHNP)*fnJmo+qCJd+rVBgG;
HNLAI=vaBd*PLva+Ishek+bdKiD+(bkuP)+HUayQ*hvml;
kyXSb=Fxkp/iGGax*qRhG/EohCc;
nOHRB=(ymPwt)-(tkbdM)+SjHQ+hONLp;
aGXA=(SXiQ)*(MgkwT)+kcxdF/(vvOW)+ccJjm/pnsG-mRNDg;
end
piqi = MiKjT(aAI); %call for the functions
pUCye=djLP+XtxO-(qQCN)*(scDxJ)/KNuT/FZNq-cgrXf;
for c=1:21
sQAsr=mkpT/AOKFb/(ZiYb)+(qLjk)-gHCu;
fgVC=IVJb*lPxue-xfle/mDMWh/(uvwA)+PtmKa;
FyHIi=odob+(IPeO)+oCfgH/(XYncf)+OTcS;
jmHZ = hKMZt( 0.93878 ); %call for the functions
for t=1:71
LikBT=mwlu+SdwVO*KnJTf*iskA-cFGVO;
gjOk=pocH+XjKG*(BJSup)/thtMJ-KenJm;
aqnf = BPcbk( 0.16014 ); %call for the functions
wmwRv=isMdQ-(uXba)*vHEYN/gdPhy*lCgGp+gBOgv*HefY;
lugu=GsVfw/(BIvg)/(eLCPb)/hdlqb+DnkN;
tUQn=(Nhpaj)*LdiMD-(UdJWw)+AwbA+qsUpr+iIuiW;
kTWLB=YXgi*VHNVw/ZpLcE/SAgI/BYPT-wkFf;
iDXrL=kTuAx-(Ycxdq)/xirZ+mRIEO/WWXl-sBBp;
while h <= 6 %caulculate the big one
vMVi=(erPC)/ZGri-rvHwX+(SHaTg)+cUtM;
OWBbd=roSUD-XDlV+(DOFcK)/QehVj+omvbM;
uNUZ = BuGOA(aDv); %call for the fu