function [phi1,phi2]=TwoPhasesLevelSet(Img,phi1,phi2);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Img=double(Img);

Img = (double(Img/max(max(Img))*255));

sigma=0.9;                              % scale parameter in Gaussian kernel for smoothing.

G=fspecial('gaussian',5,sigma);

Img_smooth=conv2(double(Img),G,'same');  % smooth image by Gaussiin convolution

[Ix,Iy]=gradient(Img_smooth);

f=Ix.^2+Iy.^2; 

g=1./(1+abs(f).^1);                      % edge indicator function.

dt= 10;

mu=0.19/dt;  

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Define where is the background and where is the volume

VolumeMask=fillingProc(2,2,Img);

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

tic 

[phi1,phi2]= EvolutionProcess(medfilt2(Img,[3,3]),phi1,phi2,VolumeMask,dt,16,.1,.1,6,mu,-3.5,1000,g);

toc