% 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 uncnowns

function [A,AR,AI]=Pade_parabolic_matrix(n,h)

L  = OneD_Laplacian(n);
%L  = TwoD_Laplacian(n);
II = speye(n);
A  = II + (1+ sqrt(3)*i)*h/4*L;
AR = real(A);
AI = imag(A);
return
%
