function A=create3d(d,dl,du,n)


n1=length(d);
n2=length(dl);
%n2=length(du);
if nargin==4 && n1==1 && n2==1
    A=d*diag(ones(n,1));
    P=diag(ones(n-1,1),1);
    A=A+du*P+dl*P';
end
A=diag(d);
Pl=diag(dl,-1);
Pu=diag(du,1);
A=A+Pl+Pu;
