?
function?[sttf]?=?st（timeseries）

TRUE?=?1;
FALSE?=?0;
%%%?DEFAULT?PARAMETERS??[change?these?for?your?particular?application]
verbose?=?TRUE;??????????
removeedge=?FALSE;
analytic_signal?=??FALSE;
factor?=?1;

%?use?defaults?for?input?variables
%?calculate?the?sampled?time?and?frequency?values?from?the?two?sampling?rates
t?=?（0:500-1）*1;
spe_nelements?=ceil（（250?-?0+1）/1）???;
f?=?（0?+?[0:spe_nelements-1]*1）/（1*500）;

%?The?actual?S?Transform?function?is?here:
st?=?strans（timeseries025011verboseremoveedgeanalytic_signalfactor）;?
%?this?function?is?below?thus?nicely?encapsulated

%WRITE?switch?statement?on?nargout
%?if?0?then?plot?amplitude?spectrum
return


%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^


function?st?=?strans（timeseriesminfreqmaxfreqsamplingratefreqsamplingrateverboseremoveedgeanalytic_signalfactor）;?
%?Returns?the?Stockwell?Transform?STOutput?of?the?time-series
%?Code?by?R.G.?Stockwell.
%?Reference?is?“Localization?of?the?Complex?Spectrum:?The?S?Transform“
%?from?IEEE?Transactions?on?Signal?Processing?vol.?44.?number?4
%?April?1996?pages?998-1001.
%
%-------Inputs?Returned------------------------------------------------
%?????????-?are?all?taken?care