Create a PYTHON code using the finite difference method to obtain the solution of the boundary value problem (use range partitions for appropriate solution and / or tolerance):

DIFFERENTIAL EQUATION: y'' = 4(y - x)

EXACT SOLUTION: y(x) = e^2 (e^4 - 1)^-1 (e^2x - 2^-2x) + x

P /.X INTERVAL: [0,1]

CONTOUR CONDITIONS: y(0) = 0, y(1) = 2