function value = u_init ( xy, t )

%the initial temprature profile
  n = size ( xy, 1 );


 value = zeros ( n, 1 );
  index1 = find ( xy(1:n,1) > 0 );
  index = find ( xy(1:n,1) == 0 );

  value(index1) = 10.0;
  value(index) = 10.0;