function m=pdeprbm(t)
%PDEPRBM Mass matrix for PARABOLIC.

%       A. Nordmark 1-20-95.
%       Copyright (c) 1994-98 by The MathWorks, Inc.
%       $Revision: 1.7 $  $Date: 1997/11/25 17:06:13 $

global pdeprbht pdeprbMM pdeprbdt pdeprbp pdeprbt pdeprbd pdeprbN ...
       pdeprbb pdeprbe pdeprbB pdeprbpm

if pdeprbdt
  [unused,pdeprbMM]=assema(pdeprbp,pdeprbt,0,pdeprbd,zeros(pdeprbN,1),t);
end
if pdeprbht
  [unused,unused,H,unused]=assemb(pdeprbb,pdeprbp,pdeprbe,t);
  pdeprbB=pdenullorth(H);
end
m=pdeprbB'*pdeprbMM*pdeprbB;
m=m(:,pdeprbpm);

