function [ u ] = I( x,y )
% forcing function
%   for tests with NFE , example Hutt 2, pag. 4
%sig11=1/32;
%sig22=1/32;
%I0=5;
%u=I0*exp(-(sig11*x^2+sig22*y^2))/(32*pi);
u=0;
end

