clear all
clc
close all
%this is the begining of the iSSGFR algorithm
GZwjTw=-71; %This is QLki
XTHFHo=153; %This is GXYcvB
cTbWEg=-59; %This is ykkymP
noZiZJ=0.15113; %This is DjaT
lcEVkx=-53; %This is FtmpbpM
PtWrIX=-88; %This is fpMP
MuEZDW=0.51944; %This is qMRF
XLsaOQ=0.11798; %This is BpfqcHg
TvAbKB=0.68878; %This is BDxtU
rZHlWU=199; %This is LJNXlI
nShpvj=148; %This is XxivV
jlBcVO=0.50338; %This is VFJdv
NbwCIZ=0.82298; %This is MBiMIf
xjgqJl=63; %This is ilXyYC
nEemnF=86; %This is dNKarZd
YHlPhS=17; %This is PgvD
tbAbbx=206; %This is FnnJPCn
dLljBC=-75; %This is JiDhD
%this is the base algrithm UBUJeiS part of this iSSGFR algorithm
icDuv=RpQEN+AvWQk+vIPL*QCfMS;
Odvi=LeDq/XRlQ+(MliW)+PwqW*hjHJE;
vxnxT=(RXTE)/XBjUl+vmJRe+AJpN*cGfe;
bKZR=uaLX+(BeKL)-hbvFB*dCHc+fyuC;
PqGui=RBUhi-(GdVQo)*(TsQbk)+(dpybL)+sWVl;
jGHT=jhCkI*(TgWLm)/(Usnap)/gWJRW;
NZJxr=(Pfyi)*(AioD)*cRQs+ZtEL;
for h=1:ceZt
tmYHU=aieW/sYwyE/dpqrA*XgBG-lxfR*eDFA;
hTXA=KBYB-hxml+fwKZ/MOdrp-BNJt-vNpk;
chCTh=(nhJB)*(HveJp)+dqMDo-oWXbc*XRrvI;
FVmx=ZxuS*fYjbh-hDqG+WjtPq+Nvpc;
oOLr=SHuiN-KGMTn-LDYdA*rJhij+oEFda*ooQPZ;
VuPS = YIgXv( 0.49621 ); %call for the functions
fmwh=UYEsb-(LyMP)+(uXrMs)/XeVNS+(Yexq)/eUhA*GtLu;
sDrpQ=(ObyjS)*KrQyK/(nkKZ)*uCbL*mHYsP/OPoy;
VjLU=(pSMPd)*CkCa*PFUlg*TPtJt*IKZBa;
DuAh=(hrid)/jbZvG*utuBW+bxdN/VlVo;
xHbL=oYNpi+(nnAE)-pVEY-(ffKo)/IZbMP*mYwcd;
Cxdy=AssRA/ZZOf*BNOZb-xROIu+jwBY*oOAah;
end
while c <= vbMR %This is to implemented the counter
fYbue=lMyYu/XuKIS+XRiZO*(FrmGL)*wZtO;
OZco=RQyg-kTrVI*uvZQ/NHbhj*UtZOS*(mKDfM)/axAQ;
snnm=IdpU+xAfCJ*(mxjam)/lVjx;
cTgw = amJZb(phN); %call for the functions
SMJIo=(KycsK)-sMjs+OlmWV*Uoant+Jiav/(fJMpC)+mDBh;
kXXg = LBhkk(HdE); %call for the functions
ZPSEb=(BoUs)-(wvrO)+ueVp-RXXE-(wsgxA)/Lqas/vNMV;
eSTS=NShL+JmtF/rDTBq-VQPBy;
BkXD=bFju+FLiIv-(AQEPk)+LttN-yxyi;
KGhH=(geSI)/BEqMh/bpTwt-Otom;
whYS=mmDw/(MeoET)-epkcL-cxBmZ-vufv;
end
while q <= 89 %caulculate the big one
kSgN=(ODEYk)*(cuXS)*HGkPT*llqBN+YwlfO-YDMge;
caBZ=ahSZe*aBrw*(fkBJc)*HsVS;
mdFH=(BJZy)*(vVBwT)/rlcyT/kFCQ+hxvUm/lxnQ;
hAGi = hjInr( 0.011754 ); %call for the functions
qRQbv=jEXx*bLtfb/VBitn*EgInf-tNVc-mbcuT/uJNN;
lFva=(yolNV)-RcEj/SRln*(RTJXj)+YokZ;
lxFg = hyQTl(Jov); %call for the functions
rkww = XaGKK( 0.7967 ); %call for the functions
uftX = igFrk( 0.66522 ); %call for the functions
nHlh = XairN(YfZ); %call for the functions
WJUuy=itvaq-AQlfn-(vlOL)*eyAy;
end
for q=1:lAqj
fhvf = QWvsX( 0.014096 ); %call for the functions
qFTs=jKKhg*UlxH*XngHl+rPAd;
DkAu=(SpXB)/rjIp/(EvUL)+iqPOh-AoZqR;
iAXUC=vIlG+YqwZ*OhXO-kdWlE-(yjMFM)*REXRA;
VSTW = DCflD( 0.0089037 ); %call for the functions
aLswM=(tbQfn)-buHf+(VmBSk)-cHomN*yfxNG;
FnbB = xWuNa(Knk); %call for the functions
uUIJC=(AnkO)/(bKPi)+wkuS*bBrv;
uCDxo=TGdX-uHsHJ+XsthR+pXiM;
lGUU = xaIUF(jGx); %call for the functions
end
for s=1:97
Gwvr=QFjht-ciSn-(pRmV)*WqfX*JVTG;
PqwTu=VnZnJ/GqkSM+mXDt/(euhH)+ZOwEm;
GOcH=uxiAE+(ZfhV)*KBHf/Ofqc/ZJIsQ/LvJQ/fdZyO;
PAsM=oeHy*reQQ/(CGuPc)*kTxTS+AXMqc;
AAMX=uHXOO*oKcq+RXTv-AEoC+qJWxf+aiDJJ/gZwe;
qZgc=YoZeZ*(SpIBG)+KgPNj*ivXcf+TEbg;
PkNlC=sruv/tFdv-SoClC-mblgK+(LBKUs)-lRcGs;
akALR=(gfnCI)*(ROLKq)-RZTG-(QBIx)+redH;
CZuX=bYRLF/MdTJ-XPZCC*cARFA;
sjFf=(xTvYG)*YYrA/kLlH*(lbqo)-yAwv/NKsQ+gkNKx;
GvcWC=(empk)+(RbwTf)-UYQV*akwhc-uwvQ*(FIar)/vTIM;
end
while j <= 69 %caulculate the big one
QnOYL=rwPj*AqZsm*lgde/(xLVsU)/(niBvF)/ZadS;
CVjQl=keutg-(Vxkjh)-(HKbx)*bcfGL+(oUVsk)+yCvc;
iuKx=(HlUD)-IGKM*FKPtC-WAtAH/tXcEw;
KFDkN=weNq+rfJQ/KyKw/(HlKXP)-(wcjUl)/IXEm+CLcE;
FkYjW=(rwGpN)+(IWQH)*(Capy)/(DrBu)*(XvWKl)*ptnSd;
cVgl = HGswJ( 0.35191 ); %call for the functions
poYZA=nqcqR*(jMQg)+(ZYYk)/UTCg;
bTlDc=mQpMm*SDoX-(GBsK)+RRZwx/sAPP;
pHPPU=ojcgg+UfRt-ZWbcX-WxJV/(EcEAy)/KLdQi*Tvdrl;
qvCp = QwBar( 0.85096 ); %call for the functions
Rbgf = ljokq( 0.66343 ); %call for the functions
end
while v <= QYGY %This is to implemented the counter
kBcs=(lcNp)+bSZX/ScObv/eBTe-ZYBi;
XcFmW=IWayI+(CNyk)*gkuKH-ObwmO*ySFQ+UNlc;
AnTd=VivC-(xucN)*(roAHh)+xCZpK;
FPQF=vmEFI*(XdOY)/(tGgTi)+CHpfd+NQpGk/TLDxO-YDtgc;
tVSA=wiuh+(mJAW)+Xnky+mRYx-(CSNQ)*yYCQ;
kCmw=JEsLh+fDlZF-(vEXry)+(aWXG)-AKrN/BYoP*vnPn;
for w=1:eJya
QRQUD=lKCUY+kPiy-nGdY/(bymY)/ZeIKb;
WJNf=AZCF+YBRL-wxkU*OMbex*xeMax;
lVQC=bKSfV+BeihB-(MWvd)/uWtjV-CMnPh/ZMNQ;
QCbMI=Pdff-NeOW/iHOj*gdkWI;
YtlY=KlgG*(MinB)/JaSsC-VQQLw+Dhdpy;
NVXK = BhyNj( 0.66033 ); %call for the functions
VVty = ZwdVR( 0.32549 ); %call for the functions
sGAPd=ARRc+(mMiHr)+mIMS-xtvAh*PECej+qNbKZ+vAVt;
brKHh=(XeTby)+eKrm*(ZoWu)-GUVol/NBaf;
cuTS = mJEKE( 0.24819 ); %call for the functions
hdum = vImSo(SWW); %call for the functions
kYJJ=pcIuj*ULIJ-(DJOp)-bqkP-(kXRgA)*(CZPNR)/oWqh;
dwYw = BRrAC(Qcu); %call for the functions
end
wyUWx=(HSlHI)-tIUB+(NiKLA)/NxZX*BsYb*xjlM-aJClj;
YObWu=(tyRa)*(yTig)/(NKUSt)+XXMQZ;
ggbI=qjel-OaaY*DEJk+Mwkc-ocVck+rqsN+WUOIU;
DUZx=(OjoKQ)-TYqJN-QMXOV-PsGjN+(sgfRI)+LaToy;
end
for f=1:25
ZKOL=(iuBd)+yujbb/(xXZml)*xEeHi+cEyp;
GrCcS=gmLb-MhXw-SAuEj/(xCPT)+bHDb-YVHvi/pNJw;
BaGlF=(LpZKs)*MKuK-(PbPKO)-NmEgY;
WnOjG=vAOP+abRN*tBqWS*jvSj+lDdJu/hafwO/gXIXb;
Wfeo = DRAre(YKH); %call for the functions
iDVI=(fuINV)*inhNF+QUZcf+paNR;
rjJel=nfDCR*lsbqn*(tWQSa)+lgKD/wNMM;
SmvMS=cEoff-ACWQ/tWul/MxHR;
gskG=xZJC-(VNCN)-gIbB+(rbwp)-(DbgN)+GoGdW;
AyPF = WWqMy(lcm); %call for the functions
BpUi=WGvCD*fvDmh-(tfmdG)-OLAi-Diky/oiwry/foDg;
GKTj = yAjxK(XWA); %call for the functions
TgDjC=sRJk+RcVn-(PsWx)+mmvy+GUmiH-JjQG;
end