function [U,V]=sample_field4(X,Y)
% sample call
%phase_portrait('sample_field4',-4:0.2:4,-4:0.2:4,200,0.02,30)
  U=sin(X+Y)+cos(X+Y);
  V=sin(X-Y)-cos(X-Y);

