How can I plot ROC curves for a multi-class problem

Question

I am trying to reproduce ROC curves for a multi-class prediction (Iris dataset - 3 classes) problem. This script from sklearn package works fine if the test dataset