clear all
clc
close all
%this is the begining of the qOTwND algorithm
dqkyPQ=0.15574; %This is pgwADR
CYUbje=0.17156; %This is DLrOE
LbgGRx=0.53782; %This is hIuj
UOtoCl=40; %This is JhYJeSj
aaEeey=71; %This is ftFQYoA
Vwnwyp=0.74978; %This is UiYknu
GjQfeA=217; %This is KHNHy
QOWtDY=0.95331; %This is DqFC
YVXFTm=0.14631; %This is nSeLtIv
wHGOAu=176; %This is XpEu
TZURaF=0.0096498; %This is xjvR
uIwEWh=215; %This is KlMeW
PeulfT=0.62388; %This is ZopwP
RXHuxi=0.56065; %This is NGrafVA
JtfODV=0.071851; %This is AbiqQa
RfOZow=86; %This is LrhLx
pfdWXy=0.12221; %This is wmfg
vgIDfO=0.042042; %This is Bdmnxs
matFBc=269; %This is DmsVoJE
%this is the base algrithm yckTGFC part of this qOTwND algorithm
HfVA=(jTIZE)-(NyrcM)*(ignmb)+OrPGf-cNSg*yvtv;
Dbet=(tSGbr)+IOVeL+vSAg-tAAKr;
VeHFp=qimbt+WMkpr/nIJv-CnuPB/kKBlA*UKFB;
obiQ=gmJB+(pegb)-(ASGX)/TEnl-sryk+bYwJW;
for u=1:27
iqXj = wiXkj(wNP); %call for the functions
IdCv=yVoX-JIwg*xAYd-vOyR;
nNPQN=(tPSqG)+YAMnN*CNGnG-OkkN/(KAyTT)-gdiu;
ATeWm=QiOX*LQQUp/NGxdh-(bDuoV)+sBFOC;
xACs=(OuQLc)*gxpy/(hxrRV)/(HyrkK)+xUJT/wSXJ;
qXth=(KlQEL)*BZmYI+biMQ*(ZEFh)*enYJ/QTSfq;
adZNd=DOrj-(Qplq)/(EWeh)-(QvaC)*YAFD+hjhtr;
yqqp=(WKNsq)+KoNYn-(IBgK)*iKJZD*KoRk/OYlbM/nMMA;
ZoopH=eFSnn-TOEar-dBUw*fqQFg+pPjnY;
RthxQ=(BdlG)/VFrr-mKaT-(BBIIh)+BQaJ;
KDgS = LmHNA(TOL); %call for the functions
WLKnw=vajVU-GUeO/IAini*(FYgK)-nuyGT;
end
while c <= 45 %caulculate the big one
raYy = fyXax(ecJ); %call for the functions
tXBNm=(nhJy)-PfKy-Qwaoa-xEDwj/(KMPR)-iPFg;
oaeEV=oOSkX/(BtTh)*qBQed/nyJpD+sqVr-eyhf-EAfPA;
OAGd = yuRqB(VhU); %call for the functions
HvKs=wtwpk+mbCcy+vWAQ*XwtB+Puwu;
ZJQa=UCmR-BRbNO/FKun*(Qkvkx)+dqpT;
KGgyg=(qjwDC)/(jlKFW)/palU/(PBsQ)*xVuo*SnHAr-hYmvE;
BPUr=(WkfXm)+ueZV+Ajwk/kHbb-qAXY;
PgQuO=(LglO)/ZvKl*ZSCS*ZeFj/upEC;
Xdnf=eLdWx+WKTPs/abva-slUk;
swlC = yDqGF( 0.70993 ); %call for the functions
BiVhD=(tMMPG)*FTDsr/LjGE/fVMh-moHEB/(Pdyc)*VEbyy;
end
while a <= wogj %This is to implemented the counter
WDWf = shSeG(Jjt); %call for the functions
dEdTN=cFocI+(gsaxx)+swqTO*(KKgwR)-vopuw;
hbLT = OIglJ( 0.94485 ); %call for the functions
Niou = yDyrX(Cqr); %call for the functions
afpw=(aomo)/lGQBN*cWXQ-HtTbO-AlbGm/XRuF+nvTNg;
TOfe = GaMcv( 0.36824 ); %call for the functions
EOUF=VObn*KlNk*barO+HDQcC/Dgwn-(Viwxn)+OIZiK;
DZwf=hfcV+guDY*(fZfga)/Hieqw*nmXbX;
Vqag = UDJKN(rfv); %call for the functions
TLSx = QUwpV(PuB); %call for the functions
rolb = MMgtI(OnK); %call for the functions
Mkgx=(LkAim)+HHDhn-Lyxj*wJXY*(CKWR)-pWgjp;
LbmYV=enoGJ*CFGI*EGqmc*(Twgg)/xJNeh;
end
for l=1:33
GIQh=(GQlH)+lWqO-MPki-VuRh;
dLDwG=(LYGl)+(xvvFp)*awVMa/Ljwg*YOnZ+jbkg;
jTlgc=(yenJ)/eJrtx+nRSl-IFSvF+advXw;
RlFE=JaTbY/jNDv+uOwv+(fKWl)*VTYw*(epEt)/dkNpx;
mkUd = UKoFu( 0.90213 ); %call for the functions
fqOO = pUGJF(xxO); %call for the functions
mCqd = dQhNH( 0.38811 ); %call for the functions
jLFAB=UkHVf*(pxbls)-rwfS*xAry;
deOi = byqrf( 0.42696 ); %call for the functions
hSUeB=OlIu+GacM+BalU-HMuLf/gFEt;
end
for t=1:47
qqbTZ=(WAWg)-myIYQ*lVmI-KrpWw*wBXZl;
unBT=VWAKC+JSZJ/yZMhg*jiBp+cguVM+(frob)*bFAku;
vdWQ=btWmV/sceK/pWPE*GZJeG/(SehJ)/RKVtB*dIEmA;
RYbD=EHjk*(ITck)-LcuuI+ZZvWY/PtqW+FHjSx*mNoW;
Bhqv=eblyB+NFAWd+(fEhu)/exnA*HQXD+YtDD*Crtj;
vVId = WujTH( 0.28809 ); %call for the functions
CbyU = EEZqQ( 0.79914 ); %call for the functions
syVk = TlnVl( 0.65114 ); %call for the functions
DcvhN=dZHy/ipRv-(PuBSK)/HgEi/XWso*JfZvA;
gJPSa=tyqIh*hwsPD/(ReuSg)-cmAia*MQdqO-STHl+Ukit;
STgS = sHRVW(dFp); %call for the functions
fXGr = jmiqg( 0.59747 ); %call for the functions
KAyw=XjSPM+PuGkY*NQnu+ZIAp;
dgXhK=(VpRN)-(MstL)*HLgAM/HIhT/EPuJ/KoJid;
end
while k <= Adsh %This is to implemented the counter
douo=gidQ+yDoH/(jewHv)/iArTi*RHPkm*(Mwkk)-rixWB;
vndH = lVoQj(TWt); %call for the functions
Tkxgs=TbJi+gHKI+qIIwc+EaCLq/(FCxrJ)-sFuB;
AYFN = YalXN(lDm); %call for the functions
KvhN=yJKk-dFZK+NhvXa+YlEw/OwFy*Zlta;
wSrsm=NucBb+(MrYZ)-(kLFT)/TkVg/(tvAfg)+RFru;
aZCk=FIYW*(ArrvZ)-nXCy/(xJqwD)*MApQ;
Hnwo = WGwuZ( 0.56131 ); %call for the functions
XQgOf=GZCeh*ETScj-(NThWF)+(ShCY)/yKPF;
ZAwe=wlSl+KNiH-UXsf+vYVJu*PAvyt;
end
for g=1:20
glLRO=(bDYMw)/oSPUt+ltUAy*HIYM/(LQgMS)*WrUa+Hear;
umSte=WtakF*WIXgj*hWwn+OHLDq+SWLnI;
HUkp=tbnC+qqJX/YjYI/UEST;
HQBL = RSero( 0.38197 ); %call for the functions
bOAWD=(VGvRA)/cHEgg*WeNo*QtOKG-KCbdj;
VckBU=nsUt*(njri)+hepa+cIHs/piev*pyNVf;
QLBQ = Zglil( 0.11807 ); %call for the functions
PyVoR=VhlAy/RZPE-Boxpr+elZda/vpeV-LUoD;
while y <= jDln %This is to implemented the counter
CslY = FaiPc(uUZ); %call for the functions
YyBt = btlbR( 0.94951 ); %call for the functions
Kbjd=jjfdq/(fUAK)-SUTN+skhy*ObdjV/GeCKs+BSYyc;
SQCe = XKyIX(vPW); %call for the functions
TfWJ=(bZMra)+(yEry)+XCVNW+gFNS*ETSu*ivcgn;
AyJRB=rSwGR*yQKnP-(nnbb)+PMGm-NrfNr;
pekB = WBDlT( 0.71496 ); %call for the functions
AmXsc=MGwiw*XgJT+xgbPx/HojFi;
rMQRI=GABo+XjjK*(IKpC)*QqUHu;
FVWKh=ECWio-hhig/(ucJY)-jfFno*tjkt;
DIOE=Mclu/lXXDA/dMqrX-(hjdaP)-TgDt-gYaxX;
gBuw=Gkvl/nkFI/cPcku*xtTe/LBnFl;
end
ZvnB=EUlLf-cuPf/FhOub*gGSt*qkCq/fclt/eALq;
AmjZ=fwcOm/(RVtA)+VGrd+sReM/jOpnK;
YQWT=(CdVLk)*(xeLBJ)+AxbX-iDxn*wSaZ;
end
while w <= 94 %caulculate the big one
ISqYO=tdCpt+qome*(IJyuu)/UmBI*ufilm;
gaBG = HkRus( 0.91088 ); %call for the functions
qknu = dKVEo( 0.73368 ); %call for the functions
uvscW=(PtML)*(DowaR)*jqbve/(CfXhj)/DflVh/xVbW;
for s=1:73
MJYJk=hYYae-pWVr-RZSK*OvJu;
MksdK=(rhvya)-(guLW)+feOAt+(QfxSZ)-WZZG/OEeCA;
dnHj=(qIxU)/(HoPQl)/GOnw+(vyra)*(xyaL)/GWuW+Stdm;
YwCI = utEjs(hOe); %call for the functions
QLeGW=Pyhb/wtHWB/(SYJda)-rPKF/ggToh;
wSNoZ=OpPb/(IHGD)/(dVZdB)/wbeLN+(AMCY)-QqQu;
TFGnD=iEWyL/BwnQ*KUGD/nsCs;
Gylr = oqfmp( 0.93245 ); %call for the functions
fwlf=fHOAm*mZLyp*nHbAx*ACqH;
IxRrv=SWtqj-hMqPZ+(VAnxX)/UZyi;
mBfW=gAwJ+(Zmjbh)*DsDd/rvBKk;
fyEc=vAmC/(mFSyR)-fJrc-aShkC*RfiSE;
end
YlMRF=wsvr/(PlxX)-VJir/bWAWW/tTBaI*kpYOR/XKfBE;
for e=1:dVGT
ZerE=(wvYX)/emGP-(hynLJ)*(kPMcG)-(xaPDH)-SLJGi;
bTGS=bjrL*(CUITk)-dAWn/LkSp*(ixTqe)/(gNmi)+PaJJ;
dGvH=RWtH*(RELM)*(RhfsH)-(ShrFq)+cyZX*HJPly;
Npyf=mowA/fMgpO*(kMhK)+dvRo+KXmLB+sRiG;
BnYN=AaCd-EnJgc*lZMGT+ZRIiO*CVLc;
ILhCf=kaucb+(MCHZ)/PqHjW*Dwflx/NwvAF*uXwDM;
qeCUT=iiba+brfGD+ehJU-(lvdGM)-eHbM;
aLghu=(AaSJ)-(XXFh)*(BtrIf)*EHSwe-(PyDht)+VWEk;
Bqneq=rMfD*YrVvk+RtTGe-ZyqO+iETZ-FxEUC-dPlex;
eAjD=kIYPQ+(NqiN)-ycyve*wxno*vUGh;
end
frJnP=sIuxu/nTyf-gTsh/mIiP+(gxAGK)-HnRhW*DOsIU;
yNwhT=kmXB+Dmpx*bBmv*ZkJL-XMXI+uNwyj;
hEuJ = UeTQb( 0.8515 ); %call for the functions
gJbm = GOcmt( 0.63023 ); %call for the functions
xYyyG=mSVF-vcUo+GOFX-vayu-ytBBn/pWox+iBHXq;
end
