% Generates Pade parabolic matrix
%
% ( I + (1+1/sqrt(3)*i)*h/4*L)
% where
% L - 1D negative Laplacian (1 -2 1)
% h - mesh-size ( h=1/(n+1)
% n - number of unknowns

function [A,AR,AI]=Shifted_omega_matrix(n,omega)

L  = OneD_Laplacian(n);
%L  = TwoD_Laplacian(n);
II = speye(n);
A  = L + i*omega*II;
AR = L;
AI = omega*II;
return
%
