function y=f(x)
y=zeros(size(x));
[m n]=size(x);
for a=1:m
for b=1:n
if x(a,b)