import numpy as np
from sklearn import datasets
iris = datasets.load_iris()
X = iris["data"][:, (2, 3)]
Y = iris["target"]
from sklearn.linear_model import LogisticRegression
softmax_reg = LogisticRegression(penalty = "l2", multi_class = "multinomial", solver="lbfgs", C=10)
softmax_reg.fit(X, Y)
print(softmax_reg.predict([[5, 2]]))
print(softmax_reg.predict_proba([[5, 2]]))