statsでピークが2個ある確率密度関数を作成して次の処理を行いたいと考えています
・matplotlibによるグラフ作成
・stats.entropyによるKLダイバージェンスの算出

どのように実装すればよいでしょうか？

環境は下記のとおりです
Python 3.5.2
scipy 1.1.0
matplotlib 2.2.3
numpy 1.14.2

参考にしたWebサイト
Qiita 生成モデルで語られる Kullback-Leibler を理解する
https://qiita.com/TomokIshii/items/b9a11c19bd5c36ad0287

よろしくお願いします。