function out = w1simp1(t)
%
%  macm 202 -- 27 dec 02 -- djm
%
%  w1simp1.m:  simpson's rule demo integrand
%

%  use elementwise arithmetic (.* ./ .^ etc)

out = 3 ./ (1 + t.^2);

return