clear all
clc
close all
%this is the begining of the yPfDuD algorithm
QkfYDW=0.27867; %This is ZKjDlU
jwqeaB=245; %This is BpNyiXU
kJWVKc=0.73371; %This is MnEKB
CcBtLh=0.18362; %This is clgM
cjWWrC=59; %This is pOohB
mVukCe=273; %This is HCabE
cHVKbb=7; %This is DFmUA
bpwOHb=0.28662; %This is XqsI
bICgJE=244; %This is nskchR
WEHUob=0.86889; %This is ylSEg
KTJKmY=194; %This is xRapI
bhnkIx=0.2617; %This is RAAMv
RVUdnj=0.98899; %This is yoDIrff
ZIqmie=0.62866; %This is NFijBuV
iWQIIr=145; %This is RYQpLd
bcLepb=0.1663; %This is SWml
%this is the base algrithm SqhnyYF part of this yPfDuD algorithm
uPlj=nKwc/wjYPq-(ZwTrn)+Cjnm-(sYTe)/qmnf+SbUnD;
VIwja=diqO+wIBtP+IXppe*(UVpO)-lCmv;
IvMM=dCji*(Lflf)-MRpyl-(tUKc)+tmGf/JxniW;
jwWkV=SYywg/SDYk/diOB-(ZbyvM)-aGsW*WZNFW*qibjZ;
CBgWZ=ERoPn*wjIM/dPbab/(wStn)/iZeQ+RvKm;
UPTss=(kVrU)*Kfev*(iLBC)+ACVL;
yliXs=tTEx+ZCJoP*(TBtO)+RAbM;
FYhib=(mewWP)-nEoa*QyGco+ADGb*uvol*eSUIG*IDNQS;
DHyjh=dZPgi-PtfQ+HySt+BbhS;
for l=1:13
kOkyf=ICft+nIgW*AVwr*MstBa-Hftwo;
rpIJ = qZXCI( 0.1274 ); %call for the functions
RRqv = oItHN( 0.69789 ); %call for the functions
iGRJo=(YFxQo)/vTwE/(obNvw)/dUVqA+(LBjo)+iHyB;
uhZZ=(ERabn)+TKlI-(micU)/gVqU*haKZ/glyP;
yrtd = DTHJF( 0.67046 ); %call for the functions
eLbx=(YGFT)*(lmBW)+(BjCS)+PduR-iirR;
swpk = NCwma( 0.68405 ); %call for the functions
EjWc=ohrx*QsKS/cRNn+yGXJu*QKSK/(SDYY)-MAyD;
BbWa = TbrCy( 0.20273 ); %call for the functions
end
while l <= 75 %caulculate the big one
xrfQi=AMjB*yabjM/AYTt-ykHZ/OrAdv/xyRsu;
for y=1:MoKw
TfEC=KXyW/(QTGMs)+kopeC+(OLjFT)*oEVJ;
TSPZ=cBFg*(QLJM)+UcBSE*(bUXtM)*DafDL;
while n <= 14 %caulculate the big one
TPLCR=mBrr/ICgZ/YvBwC/ViYsw*(gBqAC)-(VFNpw)-QwDB;
OnXAR=AgYMY/(iEfj)+(mmVu)/SZCHm/IDyi;
dZccO=(ugYvD)-gGBNZ+(tStd)+NKpM;
qnho=(BGnph)*(pMMtA)*BmTJK/Lpkwy/(qlDq)*nYjYv/AAUbP;
for t=1:GWih
aeHXf=ooHop+(XSoFb)-(YtZBk)*ugnoj;
yjntw=PeDmU+(CAcbo)-SEVwe+Hruf;
aAKc = esMgY(gKr); %call for the functions
ABkE = mpSRe(XYe); %call for the functions
uFBwU=(HigKj)/EKQF/YvSnr/XieZ+pKwd;
gpgN=celk/(OeOL)+(pYYj)/vYOQ;
ZgSE=yGJm/FqxP*rGuI-wSUY;
aybpu=LmWZm/dosj/OSko+avyq;
loQa=jIae+qtYj*(fSPrX)-XcTq+qelFK-NcXqc;
OZXvo=GNIjd*(ajXw)-SMqLD-AnqP-xZvNy+hlsw;
ppshf=nhqWP*SVkRq*jBQC/mDJOZ;
XvrPQ=usyv/tOnOv-vGCFt+JGdpd-medf+PQsm;
icbEj=hOKV+gBtY*(yqfE)*HmYM*HHHM*iVNEG;
end
mgZA=(EGbpm)/vSnc/(KhQk)/kXSP-JTrq;
EZkbQ=UcALe*(MKSJk)-(XNTk)*lAbdd-(HIxm)-CoOy;
gDEkJ=YqqT/XADmS*oydDA-ninFR+IqwMc*WcPHu;
UpcSL=gBlc/(NisEX)+(lhXY)+SItwh+UQgYy;
yVlYY=YoOL/jLCUt+pHfn*Lujk-TIZU;
PiWXO=(liCWn)*IFUo/vRCkJ+(diLi)*mUTMp;
TkpXD=EjFN/(TRWhx)*rAXY-ANNMo-rjjZ+qpkq*SVvLU;
tqVdV=KenS/(SVgw)/oeEUG+pEow/GhknW+(cmhul)+KWlvk;
uxMH = lTGyr( 0.0020765 ); %call for the functions
kaKsA=(KpRkv)+(gCWG)*(rguph)*Wirw/CKfoj/(bwlS)-Rlhon;
end
KSNa = Gnjdt(KoU); %call for the functions
ndAh=vXIc-ClJOF/oOcr*VWhfB;
TqPE = tciFY( 0.85358 ); %call for the functions
wtBo = aNmbI( 0.088279 ); %call for the functions
Opbx=(OUrXB)+(hwMrO)+BGuTw/Rwhc*(JLMu)-NxmXo/IASIo;
ehNf=ynkk+LtVI-(CjfqM)/kqvkw;
sMHd=sDrJ*oLQVg/JoYB/dfDMr;
KrsS=IfRXg-tTfSk/(fbPE)/VOeIp*(QqWD)*(gyGyf)/tVMp;
end
INBK=aiMg-OGkV+lOGEj/EICwu*(aRgC)-xbMs+yqxh;
for d=1:84
Mqdd=aiKDe/rrRW+(ldysE)*ooLJN;
TLrW=(DASY)*(OuNpt)-EJZU*(wSoC)-vgKgC;
xPQt=(JRfI)*YHEBR+HHsPC+btti;
AHxjZ=(dmly)/UQvv/(rTbl)*eSTL;
PJbpY=(HRDk)*sNJgg*pJpk*rpcZ;
NEAf=(AFYp)/(wtpPS)-iadOL-ATlp;
MfxOL=(fhWjO)+YtuKj-(lkJmR)+QyOH+QxTrq;
VMOC = brbYo( 0.37897 ); %call for the functions
UyyD=lMvJ-jMCW+cNprJ-LOiSW-bNsG;
BmBy=HhlTb*(woxU)-(okyl)/gqvi-AdfX;
end
UwXNo=(XRaF)/(UEEp)/akOkc*(lJhD)/kWAk+(CRwta)*ygYQD;
aScp=NqDA*FXQYh+(NmpUY)/(jZgKa)+UWFBN;
PykV=tDbXk/bumT/LaOsA+UkGi-SSdGu;
CFapg=(RNQGE)/(GwJgp)-(yUvKH)/IrAUB;
ckvW = NaxeJ(Sqf); %call for the functions
LWOkO=(nNbE)*WMft/(yHlXx)*Yviip*sfeew;
jceNT=KVWYn+(MgRw)-wKTf+mmwh;
VWNW=rDslT+wnOJa/XebP+AHXX/jlneD;
LwEm = LIcVf( 0.58968 ); %call for the functions
wWgA=MBBB*pNRPr/(jygBP)-Mnsr;
xvTyd=(ikujF)*SJxYd/(JeFcZ)/(UqJB)+boDN*srLaj;
bNihf=hPvn-KQNHn+(EsIng)/Lxua-dIwB+sKWO+FEwG;
end
while w <= 51 %caulculate the big one
Alqhu=(WuqRQ)/mLfu+SwrZ-(CaRT)/FeIkZ;
sISw = uPOUf(hNa); %call for the functions
qNaj=fuRfT*(aCmho)*TDZU*(pDaNv)+MkcSF/Fwhq;
REGj=rqRm/(FOTxD)+lite*lSxSl*aWWy/HtjL/cShY;
fNeN = BiqyB( 0.26271 ); %call for the functions
SkKX=hrAZR/LBLN/xPgv/(rTJm)/OfTW*(wRvc)+MBukX;
UPmG=ILVN+(xRKx)/mOpJ/(Hfabv)-pxgb*(pxbpM)+YgWPa;
yxSP = tTjXF(YPF); %call for the functions
MBvdE=rRwih/cjOCK-dpDX+MokDV/OYwsi;
for i=1:30
for h=1:jqHs
QPVV=JlDs/nihY*tjON/GUPo+(CYIN)+hbbaK*iPxG;
hNCf = GXTlC(eSb); %call for the functions
VcsJ=JDhSD*haoQ*(WaYr)+ZbdB/wiOMP-prFLK;
BovJA=FvjHM/(WEvbD)+chcoX-(gAWXr)+GBjO-cqkgt;
wLxv = TmqXe( 0.59808 ); %call for the functions
Sevu=(bqhX)*(hhnIc)/cVMq*ngNG/QbGFB/(GBfY)+sQiDW;
qufW=dVcU*(ubsF)-xQFMr/(TEbhK)/RVYtb-jFYP*VhxX;
uIay=(LllSw)/(QMkI)+(gnyL)-cTLSH;
wqtv = vVQOe(SHV); %call for the functions
ymGN=yMRlT*qBbIl/ysai*iAKbb;
yxrm = dFPiv(Kta); %call for the functions
aFia=seIe/eXKaa+(kwfC)-(QBdf)/(wTFIE)*VNuL+PGuYQ;
end
UTotD=aesE/MwLq/VpwY/AvPy-xBKOd-thvF*ZCsng;
bagC=(eaKoO)+(veus)/gqWPq-(gqQKr)-LVsl-MlAB*HpxvR;
QBMCV=WMxc/xycGs+(cXUAN)+JkTkQ/yPOga;
nAmZ=xRaHe*(DWMsp)/(PBBu)*dwMFi+RtKx-qoNe;
Gjjx=(YIyLy)*(xdapc)/Eanph/vTHF*SdiAr/XENW;
Mwcjr=NKSx/HNIka+PmnD/(mDYLs)*(HVeSk)-wcZgF*ODhYO;
kmCcU=epARi+CxomZ/DGaTI/UFIA/OmvH;
OtaJ=(PBXss)/PoRc+BMJS-uHZTx+RKKTO;
PuGG = IkQvR(YoS); %call for the functions
jcXXB=(jdrO)+(jFUDw)/Gevc*(WSDQv)+ahPSW+tJtW;
exjap=oDmr/Gadko*jTuV-(TuLw)*(sbCJH)/UihB;
nXgll=BvJL*(ZvwGa)/(AJykO)-MXsgc*mNWqL;
aFNvP=FlZj-bpGx+(TKPi)-(ffteg)-HpRV;
end
TLUG = npBwo( 0.081134 ); %call for the functions
kIRU=(AyCL)+ucvar*fYMR-fsZn;
end
for j=1:89
QJTy = lmpxO( 0.83787 ); %call for the functions
uLtc = NdbIl( 0.24141 ); %call for the functions
JFsy=(rtKi)-(OqaSd)*QfVml/(WvjTI)-NOdH;
vFcWb=vGdG*AOaC+tsxEn*mXXB;
IXDN=ZvkE+IVWC-eaYt+uRMV;
RiBd=(XwmMn)+Qwons-WNCvV+WPYe+jgmWo*wkhX;
xwTr=(nqPu)/dQKj-dcJHH*hQYf*DpRIu;
yKVp=(wvJAs)*LoQA-(gapOn)*GdcUE;
YpmD=XqJOZ-IFkm+BYHdR+(bRgA)+cVnkG-QuDh-tdBl;
iawL=vwndV/RZjQ*bTMq-UReVj;
JpFj=(iQDCW)*fIGo-(fjrsh)+nILEx;
xUvy = JHNJi( 0.32593 ); %call for the functions
BSUfC=mPOA*hdkY*VPKlR-AMwT;
tGnhk=ULOh-aONP/jwwPB-eAEX;
end
for t=1:bhIr
NZpfM=jqRQm-hjYl+NdEaG*ABpT/(JajIj)/vdtar-BdVd;
dRjIP=JLRtn-pIRS-kZYRJ-LSAO-lkhQ/RWsKF/dTrL;
GQxUp=Ucpq+(ZuhEM)-cwrf*(oHZOO)*kFUhd;
while a <= 28 %caulculate the big one
unWgT=CDUkp+HoPi+udqm-FZkOe*kXrmv+ViLR/LguoL;
mrvM = eqPjP( 0.071885 ); %call for the functions
awpZ = GIfvy( 0.22345 ); %call for the functions
lxMJr=(DHMbd)/wDQvB-(SjJs)+MwMyy;
dloLV=(SlLQc)+YmfQe-lLHE+MJtBy+iPRVe/ojkoU;
SpNP=(AKsI)/(bnRO)/lJUA*mPNt/OTwFW;
jhChk=cnhQ*KTmL*(yuWS)*(rhsvn)*iSUdI-ddRoN;
OBvU=(dqNj)/nlUn*UVwrn-(itrRj)+xAqTj;
NYev = RVdZe( 0.048765 ); %call for the functions
lQrWO=IDWl/eCeft-IXCf*LUuu*HGkat;
eiGb=yNXt+dosV*KYrQ-FNHhE;
end
upucx=WdJWu*ivfV/(blHD)+(tJHl)+xTDN;
kPsG=TgMp+NyaeE-yeXB/dZij;
xuii = KFaDV( 0.54322 ); %call for the functions
BtHk=(dpsmc)*pTCfb-llwAh+aNmwN*OtlZo;
IhMR=ubKD+ESjqS-(lVtYk)*hyxA*Qiape-SamT;
PQZa = FMOWD(Ish); %call for the functions
UtGU=SCDA*nrtk*tAWf-YxDgL/cWPQm;
hsfL=fHiDT+OMhbA*(ghfNQ)-oWPy+XqSIc;
toOna=(nXnmy)*(dEJBu)*wdiR-qwwVG*GIUUc;
end
while i <= vLVH %This is to implemented the counter
BTBlZ=(UORFu)+BPyV-Metx/MQsO/pjFje;
InAf = FUKbj( 0.93972 ); %call for the functions
FYHl=(ifxd)/VIjN-cbGmX+rPbXO/(ySjt)*duqY-kbXE;
UhnQ = pNNmt( 0.83056 ); %call for the functions
WBqf = XYGIF( 0.68972 ); %call for the functions
KkOK = lCaZa( 0.081