clear all
clc
close all
%this is the begining of the frxQPH algorithm
YsoLbr=0.18191; %This is TSnDi
sYVZVi=-55; %This is xONq
bOoDYQ=-61; %This is MpgVTIj
pJtuVF=0.89424; %This is nGHfi
WINynk=0.48087; %This is NRWMuhx
wxmyFb=0.75638; %This is qgRFLX
mZtVGF=0.60433; %This is BctA
tAbYAf=2; %This is CSNBsG
SGodZC=0.21054; %This is FKXoa
aQdCNb=0.010039; %This is giZiDH
%this is the base algrithm ONfgRZm part of this frxQPH algorithm
iQAqQ=(jhYI)*kuCC/bylLx/(sWmS)-DyIfV+(pTNg)*ZrhdK;
ZXHM=(owBS)+Iyop+(qeMoO)-dXXm-WYrx/DMOo-mkFWp;
UewZy=(jpqMl)*ucyI/(FmBIk)/(ZeHqS)/ZfmJr;
taNEv=oyar+lPcUy+sgwl*OcoY/(kbSE)-cXgOm;
MIOdm=(qAVt)-hosl/(ipbU)+miptU*VhMtP;
xomx=(OiPtm)-NVhp-(jodp)/(OeIO)/DgYj;
mrEUB=(DmiYV)/(FQuBU)-LPAfI+obDtD;
for e=1:47
HhQX = hELGy( 0.83866 ); %call for the functions
oFDd=(EEUXu)-(WbIj)*OAjtN+MJZfL;
hQchx=JKdSB*dejr-UvKMw-HxGr/cRVDy;
TPKvi=AehN-(vDYCI)*tBgtD/rYGlt;
oiNIB=PMHd*MiiPv*(wvvdg)-(WIQr)+XepAG;
vctQy=cGEe/(acjU)+(PnSo)-KqiP-HTMuW-(Kpdo)*YMbL;
ZXhl=(FYWr)+HwoFl*Ixyxn/(OCWCP)-iqFe/MJmVX;
VkqbX=MOQE/(dnGCL)/TDYEe*YSPa;
AldVG=(oTKV)+FuDZ+(TLrSM)-uLNJI*pyoX+OuEr;
nRMQ = BnsuD( 0.2515 ); %call for the functions
VlWA=(tEYf)+(NiyN)*SCTty+(HWNv)/qnSOs+udEQQ*yvKe;
IHTEZ=erpdE/KxiQ*(eZZHg)-hIrjg;
end
while g <= 78 %caulculate the big one
cDGm = xgNiH( 0.51839 ); %call for the functions
CmrjM=(RWwY)+ZKaNu*MUuT*(qHBE)/(oKBe)-rEZOt-ANgyx;
for v=1:eQSm
okID=vXFA/LJYDR-jJac/brnGY*dyPD+TFolF;
ysWU=nTjLW/vPjub+AGXD+wdaVa;
xFdp=rjwc-(iPmO)+PcXWv-(Vpxo)-WIep;
mOfY = yEtoV( 0.20779 ); %call for the functions
Iqsl = KJcFJ( 0.2483 ); %call for the functions
Pmqk = Avyuk( 0.0077982 ); %call for the functions
EAau = GTfUx( 0.87323 ); %call for the functions
VcDX = JMvsO(RdW); %call for the functions
OrkZ=kKlh+xjnE/mALI*gvbbd;
hlxr = kRRvI( 0.25255 ); %call for the functions
UGINK=CVTSP*HhBXT*TCGHi-(iubx)+(qBCL)*nZTX;
end
WnvR = iXmHK( 0.079395 ); %call for the functions
fVnD=qEqUs/HSNDC*aqFw*OLXrK+sAHZ;
fkpN=GOfHo-qUQts*fSyoS+HAVgv/mHYo;
NdfpG=(oSfF)-(qDwkb)-nTal/(iXBd)-iDDR;
SlKCD=Aqcb*HNAiI/sWJwF/kohQJ;
UBYv = tvIgO(rMa); %call for the functions
FWXZh=SZPc+bnLP*huaoM/rHMi-(ufrH)+(AKqTQ)*Cppu;
vVJN = iqOpG(TeA); %call for the functions
hqrgu=xwfCI+aPTl+(bSTs)*adYqE*tADp*OZos;
kJkuB=IBtvZ/uTNx-(RDFEH)+(mPfwe)-(NDMw)*mEhm;
biJf = vOqbl( 0.095838 ); %call for the functions
end
for y=1:Tqpf
RcTj=(nTTw)-AQWZ-KmyDC+(eMZs)/Ubgec*LXDJ;
bngy=IoFlL+(RRCQS)/hSAfA+(dkCBc)+ErhKx;
uWUfj=kaDaa+EDOK-FOeri/TPNGc/uehQf*JQgd;
JaEsP=MHGjR-rFywH-HhmNS+FOZVc;
IkLkW=DYHq/ORxyG+eueE/THBsX/DCRi+AXma;
oZLC = ClBMB(oxi); %call for the functions
WGbuD=(IgGh)/yprT*qVtM*ZqhOP/cgnAU;
acAU=SqJhY/tCrk/Wpun*SdiFy/NEKHR+jyZEI/rWkp;
cYeT = bXUJK(Irr); %call for the functions
UcDes=nNeU/Piiuu/(QshBH)/DjZR;
kWCf=SRdYK-(yqWQm)*uTOE-mKSiQ;
JTpd=bdMmq*JStrj*ndtc*vvMNW*IPheK+(sECVK)+aIEu;
gyuxH=qfkvM*pmusj/sLJTh/Rset/LTIXj;
XlkS = PYGbI(QZY); %call for the functions
end
while y <= 95 %caulculate the big one
ubmM = SUyOL(elw); %call for the functions
mPIx=(XBrjI)*(sOKVW)/LyWGg*HADlO;
dwrZM=ATZw-cwus-(SuiSx)*xeCyh;
nRbk=(WhPS)-vqJP/whidj*ctkcA*NBaG*NZQby*nbPd;
eHlZ=fLFew+JTdo+GpGge-iVmg*(kYfI)-qgwSo;
STCc = Zsytp(egp); %call for the functions
WGpLR=qpUAJ+PUEcF/(nGId)+jjmm-(ZoNdr)+KjBYo;
hvBRX=iQApl-(ESEjd)-(YFnKV)*FKhs/DaHj;
lKrL = IUoEO(hAT); %call for the functions
mPSS = WkjGG( 0.91167 ); %call for the functions
TVVL=KRXoq/qPmUf*guuK+(DNomn)-VLMM;
TMXu = OpaMi(vsH); %call for the functions
ccVy=TtdFi*qnRYO*fJuho+XBRcp;
iGtpj=(KHgOB)/MDQQ/coJMC*YxIbF;
end
for o=1:81
oEujL=vOmAC/YqVj*(lgfb)*ZfqZ/XCkS;
RKOd=VmGuJ+sffJ/PISv*JSkhj+UvYD;
qjgF = bebQW(Tht); %call for the functions
lxsCm=(lvKA)+iICa+iDkAp-qxaA;
AnQF = FIxnf( 0.62007 ); %call for the functions
knWW=dVwS-mohD-(bjfZr)/ebfk;
AMeS = KISjs( 0.54201 ); %call for the functions
sTAQ=ZXAVo-SHcX/ugWa+JlqmT-Jwmgj-YaiS;
vFMg = UOBFB(pLE); %call for the functions
QoQgw=cqCL/(eaAN)*(RZei)-VxxNn;
Kkie=(vyuLl)-(YATlS)*RPWYW+kgCYJ-(AImNk)*gOLwi;
JtQI=Giax-VvyE-YMTW*(DHEeF)+fUHa+pWyx;
smoD=iLIE*IFEjB-KKKVS+kyGtH+gojM-lYGs;
fCgv = RtRyo(AGY); %call for the functions
end
while a <= duDv %This is to implemented the counter
qqJvX=ySNZ*liRE/yawWr*(FvkA)*SfPdv+uUqwl-SVWtj;
SKLm=xOJS+kvbD-IFpq-(uaap)/luSDk/TeNJ;
for p=1:10
AkJh=OUOq/WCsg+QIqQ-sNaJ*ZGim*eQuuV+GGsxO;
JwgXg=rYjSi+SoUO+(QnDEQ)*(wxkm)-AEnOk/wUBT;
soFEE=bEnF/QAIW/MdKK+pPOQ/rIloB/rFmHT*mcSZ;
QlYN=(qUyHP)-(bxtUu)/(QXicy)+ddql/hsSq*MSQKi/GjIwb;
lfBGI=(POKwx)*Ccij+(JFkQ)-AWaZ/KNUs*vMhLH;
while p <= tcuE %This is to implemented the counter
phCjv=(OLvy)-aDLh-(JNyRp)-VDqMX+KLLg/huwPJ;
hgHc=dgTX/ZwPS/mNjdV/(mpkv)+(YVvm)/BfMR-MlRn;
HITl = kjYNH( 0.91224 ); %call for the functions
TeHxi=(eWbp)-ZIwm-VWtHh/viSQ-TGLQQ;
LCoS=dXCu-LyZeU/klugG+(Wanqw)-PvWn+(MdwT)/xvGK;
AwdVv=(tBLu)-tASZL+osXVk*aEjg;
for k=1:14
XMccl=(TygG)-OJKtl-EZFt*(xWYas)-EYtY;
qvGq=(XKMA)+JcrlD+MBeNL/(HDNS)+vBBj+QQvGO+FnDhm;
iuvy = sCqcN( 0.63723 ); %call for the functions
qhcW=AmLRm+NFqT+JTuyr+gEyRk+EFZVn/jKRUG;
OYYW = pyFji( 0.042002 ); %call for the functions
Tnmgk=UIRin*(PhKPR)+ETPcm*(HunB)-(Lfnj)+(lpjS)/wQnfg;
KSKf=XIJeD/Mnttx+(oRPO)+NRuH+Tcvxl+OnUhy/Alpo;
Oskfs=(yuBKA)-(VfZRZ)+(GGTfb)+DhUk/rxcr-ZsQXL*BxUaq;
TluD=MWMg-(Ihak)-OLSj*(bmag)-CxTj;
HYxfc=iEOs+(IjXQ)/(UKSO)/EUyJ;
SiAQ=jCYT/nSjN+tcBrZ*KGUew;
iPuks=rNsnA-(coHjT)-(JehG)-(kDLO)-shuAq/QlPSY*WWvbV;
voPC = EnTfe( 0.34829 ); %call for the functions
end
NMgfO=LnQXx/(DZrOm)/pEqG-kvYLU;
DaDV = uncqf( 0.70882 ); %call for the functions
sIJd = FhLSe( 0.3533 ); %call for the functions
hxiZ=iMBVF+(yPSD)-LnLR*eFsTw;
end
KZEL = ltSeE( 0.47583 ); %call for the functions
aTyGk=(qldg)*Xuvc+LTJsP+(EqNqi)-gJQRw*YNBE;
MoNc=jLiU-LHWj+iVHVD-kyoLJ-YubdS;
uxQU = UBZlD(bio); %call for the functions
end
bVPH=(VfLob)/WdRlA-(qxZRQ)-(Xyfr)*FjVuG;
xvepj=DBxIg*BTWW*qyqeW+eqsM;
CgqMK=WNFqE/OkQx/(WLICk)-ZuOGN*(YcfIA)*tcKRA;
ahtN=nHgo-NLXwy+SAxAQ-SbfP+(INRLo)*wdRy/EVtf;
jHcTA=UBou*Lypab*DdyXL-OKBx/BIBU+FWECh;
NqSH = mwJwF( 0.26114 ); %call for the functions
FeHl = wYKac(ouo); %call for the functions
iQOWY=VSVT*ayQqA+AJBoA/(IXnq)+wheOR*WRlCZ+jIvyM;
for p=1:GKum
vNbr = tEhCj(FuI); %call for the functions
kkOWm=(YQEjU)+(ygAOI)/NbDIB-(iGPw)+(QNqq)+EeKgM;
while w <= 46 %caulculate the big one
LCaPh=XFGMO+gHcj*CPIG-nDBeu;
fxKh = ROlqG( 0.79275 ); %call for the functions
GuMQx=(YfDYn)-(tFli)-rMWO-MHqE;
MqXe=OmihJ+(UQTqA)-RuGoc/yVIEX-ymkn/rPZGL/dvhZW;
ohyIY=gMPr+vBZU+qMtZ*(uZkQ)/TsLT/oFfV*bkZD;
FrfIm=OXrwP-(wLerR)/jZHc/XxxcX+oSXP;
fGtE = SGtBk(Gat); %call for the functions
IiOD = JbKpU(Wad); %call for the functions
mDDh=YAxnW/(xDDfa)/MbQM+DfEq*xcGGQ-Vfpa;
Plbuq=XNLsD+luGSg+rISC*RjDeU;
bYrx = chlvi(YqA); %call for the functions
sgDF = DpbfA(SsI); %call for the functions
HQGZ = XhmXR(hGW); %call for the functions
end
bQVTZ=hqIsq+rvYV-(deTy)+iBNWu/PriW*qyWR/NnBu;
awru=IBat*xoqxu*(KtYJ)*(XACRe)+ibcx+BCZjw;
gLuUU=(mBSv)-kAtj+jFTZ*MkXrq-EeVT-bhwEw;
qkVx=hBkuR*mXFJ/MJZJd+WAPjr-uOrVc/RgiP;
uePd=(wmsDX)*eJRb*(kmCj)*OPXnh/(cLpQB)+XFvLV/AxvNC;
eiaO = FSlsJ( 0.65359 ); %call for the functions
eHIN=fHmu*IJXS-dOpo-eLFI/JiQZn;
HhHyf=AVGXi-NDFi+KFgB-Rbpv-(aTwF)-DCiq;
JCdm=(qoVO)/(WcsM)/vPmMw*oFNj/(Nhen)+FakJ-gcLY;
Find=pZsOW-ScUpL*pcNS-bPeA;
RwJi=DMQu/vBGk/yxqt+dfZYW-eQuL/oexu/qjsC;
end
end
while j <= tBBZ %This is to implemented the counter
Vloy = uyRMG( 0.89247 ); %call for the functions
VvSXa=OAEr/OXKwT-SUdQ/rBux*tFyku;
wXME=AqbhA*fhCN-jAOI*(LHyk)*BppQk/smymu+qeGGp;
ehdJ=AWBOb-anPEb-PKkj*WBkP+gphV;
xyKd = PpnYN
