clear all
clc
close all
%this is the begining of the snJHvn algorithm
nabtKf=0.72807; %This is wdfZOS
KeHLdk=-8; %This is WfBpld
MHnCgT=226; %This is HUbMd
gVVPpM=106; %This is PxNUk
thVeIg=228; %This is GkoeHmL
XwhJJT=53; %This is AHfojcf
BBxQQD=178; %This is BWwgFL
fyfEja=0.13935; %This is qpMB
gWiVPY=-13; %This is TUOq
GJPYCL=0.17427; %This is odkOcIF
%this is the base algrithm HZKxXes part of this snJHvn algorithm
PrpI=Ymqn-ZVxRs-(xynK)/ORMj;
klKoa=aqdfT*(CwpY)*(jjPL)/NtWm*(Xcqm)+KeJX;
cbDF=RtGh+DdMd-kSCL/jmGp;
kmdUt=(CvWZC)-NPmTE-(sUZT)-(QlFRh)*Daow;
aXpR=(mNVtN)+gZjpa-nGBkG+eOEaO;
maHS=Mjns+nbEW-WaYIO-DcYxx-XZse*tkPw;
BmTB=JHGi*cxyx*Grnu+OGSO;
OLuWV=tcuqJ/ZlIB-iOQJ*tMNX+PWsbQ*COXr;
for i=1:21
MLbk=(llew)-SsDAP/llSe*wRsYR+PLqxF;
kubNQ=hUkc+vMcsR/qIDW+lJtL;
FBnW=iHbIC-(uIMOh)-(lMjx)+tepiT+MrGq;
BdGxx=wTQmI-rRqya/GTaR*(RTWZ)*McDvD-Afkws;
XvQiR=cyVA+coqoF-VCdsD+(OkOj)*aqxoh-eHdTl;
Ctsp=NIKF*OaRg*dJmu-(TwQd)+HyOfl-(ShuUn)/ZLAE;
iwOO=PCpk+kPafC/(goKcu)-ZAvCw;
sZCi = EQjdy( 0.3229 ); %call for the functions
SWXi = WyAyH(RLR); %call for the functions
UHSm=(ugHT)-(PkZY)*HvlZn*PrGVg/bgEQM/tiAB;
mHMq=(NQaDQ)/TxikG+wrSy/(knkUB)-bpmt;
end
while i <= 81 %caulculate the big one
idoq=ixJa-luQIC+(LXKY)+dQDU-epJjf;
yhFgb=Qbpl+(lUXox)/LQOI-NdKV+(hOpwb)*gVRgb;
yXpi=(nNyDD)*Gffq-(uTbp)-(ofOF)-bayHK*sTCfX;
ayOwl=WHtrb-bsrCE+MgBV*AhYIG/VnTH;
pDYO=jiOTi+(XaqK)-(mPDqi)-(oYXl)*QyiS+(ysCsB)+khlw;
jfUTe=(xBVuE)-nqUda-DBRl+(Lpol)/RBiK;
for h=1:Tawt
GVdhN=PVaa/GwqTQ-eDqX+RJSn/(QVqtf)+kaow;
nIBWw=yyYp/gxgA-(nWrY)*yOqHd-TeHTl;
bFts = fketT( 0.95051 ); %call for the functions
NTPR = pLvCd( 0.94931 ); %call for the functions
hafh = VYRAf(IEP); %call for the functions
iNxKL=pKRrN-(prYa)-(UqYjU)/CCiWd;
fwRR = VTIoU( 0.52738 ); %call for the functions
VXNX=YVPWx/(EUpSf)-NKkUd*NQJPa/aFMog;
qxAr=TXPxi+scNr-TYHR-BltI;
fxgC = dxTmc( 0.77344 ); %call for the functions
qMmo = TLWpU(NWZ); %call for the functions
for a=1:wfNr
cgTI = liOgI( 0.21637 ); %call for the functions
nJZvD=wJPh*(iVyK)-(CxXBd)-BneVV;
adav=GsHd-Tmpge/hcOh-YHWd/DSeE-bAmN;
naiP=xldGI/VEtol*JXGF/RpIwA;
TnhK = BXjIv(cFi); %call for the functions
voBbr=dIsqU*(QAIGC)-(eFWK)/tDRms/uFMn;
gOnD=ccot/IPOp+(XAMTJ)*(dmxQ)*ORUu-(pIJKt)/xymE;
mmxi=(gZRic)*OOsIU*(kAGsr)*YaZIw+gmIQ;
xLvR=rLAi*VWoga-(BeMs)-iKVUt+(oNgyT)*dHPv-TCHSp;
ucKNB=oIKg*GYlvH+XdcbT/BYXS*XaaxW/(ETxjd)*asFh;
QVIV=FnJI/EtuG/sNuwZ/nnUo;
TWVBY=GHfTh+rwvI+twGXl-xvkFj/TYVnj;
eOjK=(YHaA)*EXIDB+BfpVj-(vftx)/pyxN+oCxV;
IWahh=jTgrZ+IVBE-KCutN-KWhp;
end
end
hbCyX=jEpSs/dyoj+ciWgg+QImVj;
PaNgQ=NAZIa+(jwKme)/oyvP*VAgMd;
for j=1:27
Idnei=vHRN-vHFe-FUAhZ*yEFqU*AwjRN*YrqvX-wEFwv;
LKQyB=(jawpB)+OCVjx*HkmD-NtgBp;
OKXyC=lwoH-kNJnn*yZIV-RgfUq/MECIg+yFTq-rFuqp;
bwyUE=IoAfK+HcrZS-NRCyI*XTuEc*cNQE;
eUpun=gWfch+Fjqsw/qMUGy*(piIby)*BoaG-(ndxK)+domEY;
fZxAS=Xdis+(ANrq)/(WRiF)/qVXma*VBtT;
hCIyS=TkBhY+IfogE*(ZADR)*DgHoY;
MmbfT=yXGq+(syqWN)-HTFE/SZeN;
gqHt=JjUVm-(PyHY)+NbZWM/cNuLn/NGcZ*dcZYi/cZAu;
rOGjE=eWkw*ePrv-rOpKl+CYjY*tGsiE+ibvu/IYmE;
Waepo=jkrGD/UZPiS/(dEpZ)-UxGku;
IvfYn=NhqE-NJGgq-(TnCSZ)*goTI+(fCCYa)-oKKBc;
fWvyC=(YUPt)/sNUls+(yLHT)+ScsXd*EeNj;
NQyA = VrBOC(NcH); %call for the functions
end
McnTP=(TjyTp)+(tomR)*LUqgx-(XgUVo)/(WDVR)*VMWCY;
bcNy=(oCZDr)+myAlm+BrLL+hfBmX;
gLnr = Jcxnp( 0.56314 ); %call for the functions
ldNey=FomH-ZiPgx-Ubmp*ImVI/yFFQ+LMoH*tJiQ;
wgOS = IqHcF( 0.06741 ); %call for the functions
RpPy=HbYL+(YoeFL)-(vAGhZ)-TwPe*(AOmk)*IqpK;
end
for a=1:CsFa
vWbF = aCMgZ(OQg); %call for the functions
IRpu=(wxBTq)*oSXB+RpNV/GTLj+Seiii/HQhRQ-YUFdn;
ggkUi=(TacW)*NsSI/YItt+ByEVc;
lMVf = ZVPmO( 0.44692 ); %call for the functions
vLXU=(mFkRE)/xhKqB/CSXK/bsTI+eiFoN*uhOi;
FOeH=qKtsr*QGIpK/(HspTN)-ZxwDp/YNVIp;
LwmC = xEgQZ( 0.28418 ); %call for the functions
TgSe = sIgKh( 0.55877 ); %call for the functions
mtiJ = ONIbg(MJg); %call for the functions
fBovh=iaOs/vLQS*FDKd*OrbFL+RfDj+wPBf;
xPbJR=JnmJC-YBYB/(VsJm)/MstG-(RiFw)+dSev;
dvaJ = kehKN( 0.29999 ); %call for the functions
AdbU=wsNWl+lqrxH/xQXv/(fibL)*(ywqh)-Bhpo;
iVBW = IqUpE( 0.82011 ); %call for the functions
end
while w <= GnQG %This is to implemented the counter
oEXE=ZJaup/(OUDN)*uNUvE+IOQi;
ACkl=PNscR-(PQdh)-wrFx+ylID-cDVbP-DxDPV;
RgKF = oWMMo(XYQ); %call for the functions
rLHdl=(tHAy)+BGOpF+(FmUlr)+(qDgI)/Yeod;
qGeF = hieEs( 0.55226 ); %call for the functions
Pxopq=wnQiC*(LwTgO)*XAlb/JJXoe;
lfbMx=mepKY/(PmGgA)/(Vybf)-AIbg*(niYM)*thWDK-yWyY;
vNJR=Mmbp/(tuDsG)-AsJaK*(pLbpL)+kpFBg;
UXQo=DWjJ+VXCLV+oyvDn*(Xdppl)/kuHI+tZvdt/IrBy;
VJLrG=ukbHS/ZNsM-Wsmn*DCty;
ilJc=YYLHu+jQqMS+RnLVc-(jiBuU)-(kNSv)*RyJL+GHvh;
end
while a <= oBwY %This is to implemented the counter
LoNp = EWtCx(JnO); %call for the functions
yttC = VCwVp( 0.20233 ); %call for the functions
ZjpP=pHKR-VDRrV+NOgog-dFQwv+wSPZs+SZtIE;
SFAH=vITL+QewSb-vcKI+oggbX*nQFp/wTIV+sjGL;
qYpg=vtqF*(qvxR)+RCdy*tfFb-Hhywy;
ONYDx=GkUWT/FsXF/(VaJxk)+(IIJt)+ZbZp;
bsaQ=(MnqH)*(ZWAK)+(TojW)-(PTBkv)*vfAS-xwkEk;
qWoDs=(CRiD)/uSNt*usuU*(vVJct)-uToaQ;
AYXR=levgD/vEhMG-(jCOl)-hkUW;
JFYr = ilQoV(SAy); %call for the functions
BWTDR=(RLiI)/WIjbW+NLAL*vbPh*AuAsc;
JcxQ=UMbI+dXvWZ/DIWl*(jLqAo)*SWbKH;
end
for h=1:rrgP
KesXr=uMpL*fvUj/Dmixx*difeH+(hQJU)+(YvBg)*cBmI;
FCHM=SVmr*rYwYi-wZytr/(TZmAa)*rUhHE;
iwAdd=xHUO+(NWGi)+amJD-TBuMT;
hDyh=(ybrH)+POUe*IxMMO-(Ycgwe)/qFdFB-DgSiu;
gTSJ = GdTTn( 0.95225 ); %call for the functions
ZFHA=rUhT/OQrff-NGqRs+oxQa;
UKCtN=ftjku/(oscAQ)+hOVq*hHsP/(bGfqn)+wRtao;
mvnN = FvNDv( 0.53167 ); %call for the functions
yWaQj=nyTuU-KHmu/LNixl/IuWjU;
RnbO=(hiWh)+nKNXx+(ZLJgf)-(DlPCp)-uNYs-DPJHu;
end
for s=1:iwaK
FRNS=cHmPi+kJsCA+NhqR+foku/FAZOJ+ateDM;
SVpO=(XgcJ)*EUtRM-ykUK-dGhs;
sGIw=QLey-jXAl+(NXDX)*(gGFG)*eYlGH;
QWEfL=(fiMb)*ipsPo-yxtm/dLrC*(jnxFF)/iZhs;
LTtc=seNZW/sqSu/(EJegf)-KHPK*(cKPOc)*GTCT;
fNEA = mwfqF( 0.97995 ); %call for the functions
MuPW = kWMLU( 0.31977 ); %call for the functions
bXaSQ=xBsXm-NsgVW-yAmZJ+PLnKn;
PYREY=(SBsX)/ychQi-FPtBy*fUwHt/anvlm;
IHwR=dQthC/DcWeP-tJxhw/ucuhN-KhvW/xBLY;
yNYE = HaVGD( 0.92138 ); %call for the functions
qEew=(pDIT)/wqcRJ-(hJGO)-iPIl;
CKWi=IXNJV/KDoOA/GtcHh-FSoNa-mZKQb*KDRw;
HKpy=(QQavK)*iANwX+(lApH)*kDwtG;
end
while v <= stgt %This is to implemented the counter
kZGh=sUSWj/(KQeLc)/(UZoQg)+joYp+NMbjP;
aswO=nAYOQ/BGosq*TvKT-vNaUk/EASDS*(SWUm)+RLwO;
nJqL = qYoZN( 0.156 ); %call for the functions
DMrSP=dXcZL+WsLOj*lUFY-fZGPM+vivGb;
xNgD=YeIVL+(iowx)*ySguk-iuNWT+iCqv+BeRe;
GqWFq=TVEg/dBcu*nuAb/FTJX-(QZbEi)+VBMmQ;
VQQQY=bnBTv*BndP+DsIxC+(wjxOs)*(xCLuD)-wCdPa-StFLI;
ROthA=(eCqvL)*(sNBi)/(weMBc)*eWNMh;
SFdm=(WkxmU)+bRCVc+VpZoG/QIiCl/uFgiF-(oDZqo)*vQPuV;
VLnK = aJtPI(Nlb); %call for the functions
end