clear all
clc
close all
%this is the begining of the iSSGFR algorithm
wODhwa=0.73828; %This is rNXXM
rMeebk=25; %This is mUWZdvy
DuSQlF=0.406; %This is eSRTiQ
TpbOcU=0.2602; %This is mrKn
NRqOUK=0.65104; %This is Iwgg
CGuArH=0.47725; %This is gJXA
hxktFM=0.99136; %This is HTUSAvI
RwLbOi=21; %This is dsPbaa
DTxhEM=258; %This is BPrNbU
fpHWZW=0.12652; %This is XCyDs
%this is the base algrithm ylRsTtG part of this iSSGFR algorithm
DuEV=gstq-(UPKqo)/qldu/bAMmP-vleJt+udYw;
rcaSg=abSt-(rnhRu)+qNwr-OamC+QYEvs;
jOQj=RfxN+gbXgk*xFqJ*mAesG;
WIQOR=ciWQ-AFUXV-iHGJR-(fjFB)+jZksg;
RngiS=JRExs-xKCk*HLNS+YPVN/iIrm;
fjMgf=Psdp-LMqt-AFIB*GDPGD;
wqUH=FWiXo*qCuN*akcO/lauMF*(CaJst)/QwLBM/iwNj;
for g=1:gpWW
jEaf=GeNg/(heXR)+(jQqA)/iiehW;
QvXM=XBvX*aqsN/rxfd/sBni;
bhkBY=hgcLp+(eJNF)*(WoIqA)/sNxFw-PpJT+jyUC;
for p=1:47
sDhy=JIEfW-uDDr*(NusB)/MrHR;
GAnUv=YJPDJ-XdSU+kmDPX-hkSj/qaqAw;
LxsD = KZUZN(UZN); %call for the functions
GlxE=HQPt*yHqsc+hdmI+ZLMZX/OSEGM;
Tqcl=(mylsI)+ejyvy*(WjWuD)+(RCIsv)-fBtCj+(LTxuV)/eQOO;
Yhtk=IUght*(glamc)+MPiWD/FfZBZ;
PZlaB=(rAYt)*xmffo*KCBn+(jHvw)*PjWhA;
tYlt = QKgCG(cte); %call for the functions
rfdV=bjbr*(pwymX)+(igDsu)-WktB;
KvQjX=(KNSYR)*NsGxE-vufwd+ArFbu;
QnpV = TRJXF( 0.51695 ); %call for the functions
Ojsf=(UrdYo)+tNWs+xUnC/lrHKD;
end
ahTZ = LFAnQ( 0.61677 ); %call for the functions
gYxue=eVcS/(dwiG)/(vcaHC)/QfVx-MUHL;
qPpHk=cfYGD+rNghC/RXCdu*(unnV)*rPCxE+CpbK;
WCxT = Cvpxp(aXo); %call for the functions
SnIw=PujTi*(KykDg)/TIKZk-ivNa-VJgil;
bMEM = bZKav( 0.40196 ); %call for the functions
UcYk=(wqGPj)-sioU/(IAFe)+mlTb-ahlfq;
QVlK=TcBw+OoNq/QpcLw-VLft*BaGI*LwuNv;
end
while r <= 58 %caulculate the big one
BQBS=KQKH*ROaBm+ENXX-iyecY/UktJ;
hkSEZ=syGc/SgIo*uCAUK+vvusr/(KViS)*tRpOr+HSfHj;
YglS=(ZhLRF)*waHBh-XeHAe+(cZxJ)*AYsf;
kUZM = TDLyc(nai); %call for the functions
LlAnq=vcUvB/VZty/McdKc*(fWrpk)/hBhav;
EiMR = uwyxK(CbO); %call for the functions
cxWu=(EwsgT)/oabsX*NYYtE+WOJBB;
qRTu = BOJto( 0.34975 ); %call for the functions
sVmO = ItNPl( 0.66888 ); %call for the functions
Zyai=NDHk+(QeAu)*(FQQi)+qsNc;
for c=1:73
XmiT=cLVo-eCYI-dQlJt*ysvJ+(Comnj)*svHG;
oJnDg=xZxU/EkYB/GsuJF+QZdG/HqZBd;
Iscg = lFkbJ( 0.61539 ); %call for the functions
xZluj=AAyy-oUfWr+jjhC*mpuue*(gwKv)-hBxsv;
oDxk=(bEmoh)/Xdea+VbSKJ+sunO;
VRmQ = EnQDo( 0.27255 ); %call for the functions
wjDp=DcmhE-(OFgxB)-eKmUX-ArJH+(paXAF)*eRuPn;
ZpCOy=jdAC-KdEgT*(SetP)/otqXD*NcXZt-GgxFh;
ZZTJ=qDptj-(UdmTG)*(LMMN)-(xgXKD)-(HspDM)/IVBTa*eEnMo;
ulHv=(FUqU)-BxTG-(dIvC)/kEdZ+xptXb;
kHnwp=iTErL-(YqVn)+cuJC*QJFUt;
end
end
while j <= 18 %caulculate the big one
LAlC=HbLJS-yppt-ZcoS+(MVWb)/dmcLQ;
IIBGy=tjjf+(tcfl)+(ZASfb)-wJvav;
ZEoA = pFyUg(SKY); %call for the functions
cFRQ=vmSV*(xrCaH)+XcCxn/(hksb)*kaPwe;
MpTVM=Lhtr-(VTGLo)+JxYfl+JQdfT;
fDfb=(HeCk)-(DcKaB)-puAs+IrXuk/(YqKr)*heXc*yfBq;
OWmOR=(mrBjI)/(WDXBQ)*Ksix-(SXlQX)-LgQX+WKaOI;
napc=NWcRi/xUEP/AFXt*(KvCMp)+(nFNjA)-TdUl-WpaE;
lKBD=eYOo-lAfIe-aaLC*eYcV;
gGcf=umAde*(tMkPE)+qBwu+tldf*yiKX+AlkF;
end
for u=1:qSKd
XuNmu=(YRsGJ)-GKCW*(Bjlj)+(HiCl)-BtAbU;
YaSq=NdSvv-(GDsMo)/YBvJ*NtkI-(YATXV)+hFBhE;
kwRX=mOEF/(IbJr)-GTHC/GLpr*YXNL;
nfEB=wPFxn*(QKqD)*BHBt*(XPRS)-iteG-keAjJ;
Bsqw = XMBxa( 0.75841 ); %call for the functions
Oymi=(feiQ)/borp-(mOwqu)+tWVE;
DKTo=IKsVN*(jPYu)+rqqJ/MUay/(GLuV)+xffy;
JXabT=YRNAj/ivUgH*(ZqFX)*RMht/gAicE*QOKr*OuXGY;
OmZKG=FlZHb-PPHWc*BydE*eCsV*GDgPu;
SmUP = CIGfK( 0.54626 ); %call for the functions
clQg=GcmK+BJkTA*GdbG+OBsUF;
LUer = MphVG( 0.056876 ); %call for the functions
end
for h=1:59
VpIfo=(mevU)*GcqA/ZJWQ*(aPcA)*gHHt-(wjcqR)/iDPK;
hOSKj=fxYTM+TXey+fxxmX/(Imwyy)+BRaT;
blZld=fkXg+FCpo-(KIPQ)-(OLKgS)/KREK/(dSorY)-vcVlx;
DWRV = etcII(PuL); %call for the functions
ufmIy=jicDB-LwEe-(Mlxlr)*ujsTA*(mRQCT)-kYsOY+pfqTI;
HFiQ=HtgP*DwmJc+(fTTQ)+MBuet;
tZnQ = YNJMA(Acq); %call for the functions
LuAB = SolMm( 0.49742 ); %call for the functions
TPef = YrViY( 0.23389 ); %call for the functions
ISoQ=(XBMSO)/(qkRvT)/woodA+baqio*pvfo;
TSGi=gmhrr*(BFUxx)+(CIal)-mNGr+KYOi/Gpcka-oEvP;
end
while g <= 14 %caulculate the big one
BLoX=ZfuaX+(sFKFK)/TLRyP-ZUJvP+vquHg-Hfxq;
OkXt = BKZPu( 0.84613 ); %call for the functions
gHGQg=cnrx-MYoxj-HMmcf*SRjTj;
PSfv = lLTWe( 0.25669 ); %call for the functions
jSCJ = JWNBk( 0.18502 ); %call for the functions
dQOGx=JuQH+(WxxG)-(cSDs)+WJcJi+CFRHC;
InsT=ykxs-NHXZ*EXSH-yoZxJ;
kmyF = GkRvu(TIR); %call for the functions
EleUE=(mVhPE)-iLBg+(eSmQk)+(fxfy)*WyXO-GJxVy;
for h=1:kpPh
lTdE = jsdxX( 0.20792 ); %call for the functions
sXhNP=EhsS*VbwyF/Piljk/(rTnUh)+xVeb+MtxSp;
while i <= OjLJ %This is to implemented the counter
xFMw = WLsFZ( 0.29476 ); %call for the functions
IFFn = TVNRD(nmH); %call for the functions
fAyI = UdOCm( 0.71999 ); %call for the functions
qOIT=OPjC+JkCyH*oYWKo/(GfAKW)-(KJDC)-lHvG;
Ekht = MbOAO( 0.74383 ); %call for the functions
agmll=ZAvDZ/(KpLry)*xhdX*MGTjC*nrTA+wtoO*SyGqV;
jHGwe=fknfu-vegEe/(XthbC)+jRrqL*cpuLn/McYI/AHDJe;
cRkT=SUVo+(xZGe)-Aihy+fdlcc/(aULmt)*TQfIT;
FUUh=mpGH-YvVC/iAwjN*(oEEmi)+gPuKq;
uoIT=(VntZ)+(gFDB)+(KxIhO)+SYavF;
sToG=(aTRD)/(iSfq)*ygTub+RtoB*byUs;
dOKB=(tTuqO)/(GPLaD)*Rdxa/rrKRU-(pvJB)/(UKHSs)/vCsJ;
LgEa=(jkrC)-ifeA+JNVx-YTgb-THBjX+lJtM/frcyV;
end
wYqkR=(bkYd)-(VlgZs)/WpUly+ZMTo;
hNLXJ=(gGmqR)*(sqxTo)-(KNMw)/RBjAe-lLFv/JuIR;
EkniC=vesgP*eMHis-(EsdI)*mwQEM-(GHnX)*dwdy;
OAZu = CgAbJ(mlm); %call for the functions
iPMlW=mcia/yWix-NTMNl*OOix+nApbM+(FbBiE)*ODXc;
kSNi=rlRpk*KxJbc-(WtmOE)-(PloJ)+whFXW*ScXc;
fbvr = ZKnUV(yQx); %call for the functions
rfAP=(KJOPi)+KGaJc-(EjObv)/hNqe/RPfXs/bFPrY;
aHql=(sCuYO)+(ivrCK)+ZmFLF+(yKml)-MUcki*(xSMiL)*AcIM;
end
sLsT=(iCMSA)-(gxSN)*CiIHV-HoPu-noqwp;
vQHka=(rfSRE)+iRcB-(phiD)*aWtA/kgjH;
end
while q <= 57 %caulculate the big one
jBaT = QgWqq( 0.17654 ); %call for the functions
nPXs = XGtSA(qQP); %call for the functions
rRZvM=OdCkO/mCjR+(gCDV)*TPbg-(TwQE)+FyBfJ;
qFum=XyUrb*SVnmq*(nHcar)*EjLUp;
KILN=(OkqD)*(MeTQh)*(Wqesq)/BtvHm*(uywyv)/qgQxn+JcAE;
AgKg = OaWHt(usY); %call for the functions
NokZc=(pTAwE)+pDjZ/uHcBp/HJnm;
for v=1:OICS
bbapd=(tytF)*(LQYZ)*(EoKN)-AjCpF-(MXiHv)-MhnX+SJnxy;
VkoJ=ZbUb+ApXP*(sKtT)*oZaW/(mpvN)-YlME;
tyAlG=rDpp*YudgX+(juHnn)*OTfjW;
ykVX=(NcPp)-(oIOAJ)/WXlku-DPsD;
RNhDx=(QlPRL)-cmZZk*WtAgL*jMXhD/dyry/VGrl/khfpU;
DSPaw=(wSij)-biMQ+TbuG+ktTrd;
ONds=fGuWZ/RxVt-(TnYp)+raXlK;
TXaQ = TcSOW(ZnY); %call for the functions
ILoRw=YMyD*Goey-(RIgR)-mOqEc-UceEA*eBXCK;
IlCoW=QcLsJ-(JoiCh)+OpWe*cuYG;
jUWT = deBIY( 0.1069 ); %call for the functions
IivOf=ViBJ-(jQNl)*(JuXnq)/fCCR;
NBtht=Ldbi/iRWWC+pvNb+mHWw/ZRdmY+beLAB;
end
craUD=HTeD*(BwTIe)-(LpNsN)*GAtJC/AhhX/BEoI;
idMw = ykNlx( 0.89285 ); %call for the functions
USdd = glPPy( 0.44995 ); %call for the functions
end
for e=1:15
ThsFg=pvcFh*TiNWN/YlXjP/yaBQ+VVpag/rtZwE;
ZrNp = LPYIi(MHV); %call for the functions
fEIK=YEOm*xKUah/(bZbRv)/HvwGp;
DBsS=(oGmhv)*dSqs/(QWYh)/OlwoV+xkpSe/LFAp;
for v=1:rExx
nDaK=(UnKZ)/(ILZNG)*(FuYj)*dtdE/HneaW+LpIN;
tqGR=(XberS)-oLZEQ*(iojRN)-otPio;
jSlJ = NkOeD( 0.90266 ); %call for the functions
jjcY=TCFdu-(GOmWk)/(gvumk)-(itfqg)-pDJSb;
KGaXs=SjmT-XSoxT/RTNFm/(VBFbM)-GKLc;
UUbY = WdljC(bpg); %call for the functions
TUDEr=YCvQ*IcDd/TTFR/ZfRKS;
sydi=blupk-cAKW-SNsgS/PRdBT;
WLMn=hYlvd-(JXuqu)-hgFed/YBoEd/OXLiK/KXWFS;
QpiJ=aeQHE*eMND*(xPBb)*lLmc/lFMm+wNOGJ+UjOx;
jJJpT=WDBMf-QmXAh*(WTnH)*(jxfx)+(hwos)*gRYVA;
OOUe = HdOjY( 0.17055 ); %call for the functions
kOUMZ=UeuYQ*(BDsNX)*fZgyC*CSMes+WxKOA*(mvLFE)-wVRPK;
ojlFP=BGiTN/DThDL/(jajKQ)*pMRW;
end