function op=simulate1(ip);

K=1; W=100; l=100; 

[ig,N]=size(ip);

for i=2:N-1
b=ip(i-1)-ip(i+1);
roo=(b+sqrt(b^2 +4*ip(i+1)^2))/2;
op(1,i)=roo;
end;
op(1,1)=ip(1);
op(1,N)=ip(N);

return;

