fem.geom = square2;
fem.mesh = meshinit(fem,'report','off','hmax',0.1);
fem.equ.gporder = 4;
fem.elem = {};
fem.shape = 2;
fem.equ.c = 1; fem.equ.da = 1;
fem.bnd.h = 1;
fem.xmesh = meshextend(fem);
fem.sol = femeig(fem,'neigs',10);
[u,lambda] = femeig(fem,'neigs',2); 
postsurf(fem,'u','solnum',2)
postint(fem,'u','solnum',1)

%%% Find Fourier's transform of u along x direction only