clear; clc;

syms x y;
y = x + x ^ 2
dy = diff(y)
int_y = int(y)
u = [y; dy]
M = [x, 2 * x; y, dy]
M * u
subs(int_y, 1)