f@2: function [db] = compute_features(experiment, src)
f@2: 
f@2:   if nargin < 2
f@2:     global src;
f@2:   end
f@2: 
f@2:   N = 5*2^17;
f@2:  
f@2:   % set the appropriate options for the particular feature set
f@2: 
f@2:   options = set_options(experiment);
f@2: 
f@2:   % extract features using the computed options
f@2: 
f@2:   features = extract_features(experiment, N, options);
f@2: 
f@2:   db = prepare_database(src, features);
f@2:   db.features = single(db.features);
f@2: 
f@2: end