function [u, v] = at(im, kmax, alpha, beta, rho, eps, show)

% Read the input image
g = double(imread(im));
[h, w] = size(g);

% Set boundary conditions
u = zeros(h+2,w+2);
u(2:h+1,2:w+1) = g;
g = u;
[h, w] = size(g);

% Construct an initial edge strength function
v = ...

% Start segmentation
str = sprintf('\n segmentation have begun\n');
disp(str);
analyse(u);

for k=1:kmax
    str = sprintf('Iteration number %d\n',k);
    disp(str);
    
    % diffuse u (Eq. 11 in your lecture notes)
    
    % check the stopping criteria
    
    % diffuse v (Eq. 12 in your lecture notes)

    analyse(u);
        
    if show==1
        % display the smoothed image u and the edge strength function v
        subplot(1,2,1);imshow(uint8(255*mat2gray(u(2:h-1,2:w-1))));
        subplot(1,2,2);imshow(uint8(255*mat2gray(v(2:h-1,2:w-1))));
        drawnow;
    end
end

    
   