Give `np.ix_`

a try:

```
Y[np.ix_([0,3],[0,3])]
```

This returns your desired result:

```
In [25]: Y = np.arange(16).reshape(4,4)
In [26]: Y[np.ix_([0,3],[0,3])]
Out[26]:
array([[ 0, 3],
[12, 15]])
```

Give `np.ix_`

a try:

```
Y[np.ix_([0,3],[0,3])]
```

This returns your desired result:

```
In [25]: Y = np.arange(16).reshape(4,4)
In [26]: Y[np.ix_([0,3],[0,3])]
Out[26]:
array([[ 0, 3],
[12, 15]])
```