Để tìm giá trị của các biểu thức trong bài toán này, ta sẽ giải từng phần một.

### a) Tìm giá trị của \( M = (x - y)(x^2 - xy) - x(x^2 + 2y^2) \) với \( x = 2; y = -3 \)

1. Thay giá trị \( x = 2 \) và \( y = -3 \):
\[
M = (2 - (-3))(2^2 - 2 \cdot (-3)) - 2(2^2 + 2 \cdot (-3)^2)
\]

2. Tính từng phần:
- \( (2 - (-3)) = 2 + 3 = 5 \)
- \( 2^2 = 4 \) và \( -2 \cdot (-3) = 6 \) nên \( 2^2 - 2 \cdot (-3) = 4 + 6 = 10 \)
- \( 2(2^2 + 2 \cdot (-3)^2) = 2(4 + 2 \cdot 9) = 2(4 + 18) = 2 \cdot 22 = 44 \)

3. Thay vào công thức:
\[
M = 5 \cdot 10 - 44 = 50 - 44 = 6
\]

### b) Tìm giá trị của \( N = x(x^2 + xy + y^2) - y(x^2 + xy + y^2) \) với \( x = 4; y = -1 \)

1. Thay giá trị \( x = 4 \) và \( y = -1 \):
\[
N = 4(4^2 + 4 \cdot (-1) + (-1)^2) - (-1)(4^2 + 4 \cdot (-1) + (-1)^2)
\]

2. Tính từng phần:
- \( 4^2 = 16 \), \( 4 \cdot (-1) = -4 \), \( (-1)^2 = 1 \) nên \( 4^2 + 4 \cdot (-1) + (-1)^2 = 16 - 4 + 1 = 13 \)

3. Thay vào công thức:
\[
N = 4 \cdot 13 - (-1) \cdot 13 = 52 + 13 = 65
\]

### c) Tìm giá trị của \( P = a^2(a + b) - b(a^2 - b^2)^2 + 2013 \) với \( a = 1; b = -1 \)

1. Thay giá trị \( a = 1 \) và \( b = -1 \):
\[
P = 1^2(1 + (-1)) - (-1)(1^2 - (-1)^2)^2 + 2013
\]

2. Tính từng phần:
- \( 1 + (-1) = 0 \) nên \( 1^2(1 + (-1)) = 0 \)
- \( 1^2 - (-1)^2 = 1 - 1 = 0 \) nên \( (-1)(0)^2 = 0 \)

3. Thay vào công thức:
\[
P = 0 - 0 + 2013 = 2013
\]

### Tổng kết

- Giá trị của \( M \) là \( 6 \).
- Giá trị của \( N \) là \( 65 \).
- Giá trị của \( P \) là \( 2013 \).