function F=funcHePer(t,y)
%
%
global m
mm=m+1;
e=ones(mm,1);
A=mm^2*spdiags([e -2*e e],-1:1,mm,mm);
A(1,mm)=mm^2;
A(mm,1)=mm^2;
%full(A)
%pause
%
F=A*y;