This paper considers a new weighed model reduction method using block diagonal solutions of Lyapunov inequalities. With the input and/or output weighting function, the stability of the reduced order system is guaranteed and an a priori error bound is proposed. to achieve this after finding the solutions of two Lyapunov inequalities and balancing the full order system, we find the reduced order systems using the direct truncation and the singular perturbation approximation. The proposed method is compared with other existing methods using numerical examples.