function f = lineaxy(x,a,b,A)
  [n,m]=size(x);
  f = zeros(n,m);
  intervalo = find(x > a & x <b);
  f(intervalo) = (A/b)*x(intervalo);
endfunction