One easy way is to use pandas.MultiIndex.from_product:
index = pandas.MultiIndex.from_product([y.index, z.index])
s = pandas.Series([y.loc[i[0], '{}'] * z.loc[i[1], '{}'] for i in index], index=index)
Result:
A L 0.050
M 0.070
N 0.020
O 0.060
B L 0.075
M 0.105
N 0.030
O 0.090
C L 0.100
M 0.140
N 0.040
O 0.120
D L 0.025
M 0.035
N 0.010
O 0.030
dtype: float64
source
share