function xp = west(t,x)
%
% x' = x^2 - t
%
xp = x.^2 - t;
