function [J] = J_fun3 (h,b,w,a,alpha,epsilon)

% Computes value of TV function J.
% For descent.m

[temp,N] = size(b);
x = [1:N] ./ N - 0.5;
g = a.*exp(-(x/w).^2);

bg = conv(b,g);
bg = bg(N/2+1:3*N/2);
J = norm(h-bg)^2;
J = J + alpha * sum(sqrt(N^2.*diff(b).^2 + epsilon)) / N;