Estou tentando consertar os eixos da notação científica de dois conjuntos de dados diferentes, onde um é [1-9] x1e-3 e o outro é [1-9] x1e-4. Gostaria de definir os dois eixos como 10 ^ -4 e ter os dígitos um após o decimal (por exemplo,% .1e). Aqui está uma versão simples com a qual tentei brincar: gostaria que os números nos eixos fossem pelo menos 1 e quero que as duas potências sejam iguais.

import numpy as np
import matplotlib.pyplot as plt

x = np.linspace(1,9,9)
y1 = x*10**(-4)
y2 = x*10**(-3)

fig, ax = plt.subplots(2,1,sharex=True)

ax[0].plot(x,y1)
ax[0].ticklabel_format(axis='y', style='sci', scilimits=(-4,-4))
ax[0].yaxis.major.formatter._useMathText = True
ax[1].plot(x,y2)
ax[1].ticklabel_format(axis='y', style='sci', scilimits=(-4,-4))
ax[1].yaxis.major.formatter._useMathText = True

plt.show()