Complete using Python code. Numpy and Sympy preferably.

The heat diffusion equation can be simplified to the following for a 1-D, steady state with constant generation case: T"(x) + q/k = 0 where T(x) is the temperature distribution as a function of x and q/k is heat generation over conductivity. Define a function, solve_bvp, that will accept the following arguments in the following order: q, k, T_0 (value of T(0)), T_L (value of T(L)), L (length of the 1-D slab of material, i.e., the position at second boundary condition), and N (number of points at which to evaluate the temperature). Your function should return T (a np. array of the temperatures) and x (a np. array of x values at which the temperature is evaluated)