function dvecs = mrptdvec(K_u,M_u,phi_d,Freq_d)
%
%
% Function: mrptdvec
% 
% Usage: dvecs = mrptdvec(K_u,M_u,phi_d,Freq_d)
% 
% Description:
%    Computes the matrix of MRPT damage vectors
%
%

% Version SWD970109
%
% Author(s): Scott W. Doebling, doebling@lanl.gov
%
% Institution: Los Alamos National Laboratory, ESA-EA
%
% Copyright (c) 1997, The Regents of the University
%     	of California.

Lambda_d = diag((Freq_d*2*pi).^2);

dvecs = K_u * phi_d - M_u * phi_d * Lambda_d;

return