function out = bilateral(in, sigma_t, sigma_s)

% input image dims
[nx, ny] = size(in);


% TODO 1
