This paper presents an improved interpolation algorithm which enables to find a unit function in $H^{infty}$. From the proposed algorithm the interpolation problem on the infinity point with multiplicity can be solved. This is based on the DPL algorithm the acquired unit function has low order and can be directly applied to strong/simultaneous stabilization problem in control systems. Finally, we verify that poles of transfer function of closed-loop system exist in stable region while investgating internal stability.