function [y,no] = f3(x,ni)

x = x(:);
ni = ni(:);
np = length(x);
if length(x) ~= length(ni)
    disp 'error: vectors not same length'
    return
end

y = [x.^2; 0; 0]';
no = [min(ni):min(ni)+np+1];
