function out=Einc(x,y) out=zeros(3,1); out(2,1)=sin(2*pi*x); end