In this paper, the hybrid differential evolution and symbiotic organism search (HDS), is the first-time developed for general solutions of a piezoelectric stack in ultrasonic transducers. The convergence and reliability of the new algorithm are verified through comparison with corresponding data from similar previous publications and differential evolution (DE) algorithm. This study also presents and discusses the calculation results using HDS for commercial piezoelectric stacks. The Matlab HDS programs for a segmented piezoelectric (PZT) model have advanced features including its applicability to any configurations, thickness and arbitrary layer numbers of PZT. Using the novel proposed technique, there is no requirement for initial data guess, no limitations for piezoelectric stacks and the convergence rate is much faster than DE. Therefore, the HDS is promising for direct evaluation of specific aging or degradation mechanisms of ultrasonic transducers.