function U0=inifun(n)
%
%
%
m=round(n/2);
dlt=n/10;
q=round(m-dlt);

U0=zeros(n+1);

U0(m-dlt:m+dlt,:)=ones(2*dlt+1,n+1);
U0(:,m-dlt:m+dlt)=ones(2*dlt+1,n+1)';

 
 
