Please solve it using phyton, i don't need analytical solution

Problem 1 Compute the steady-state concentration (umol/um) distribution for the channel shown in the figure below. The PDE governing this system is: (aca2C) D + - KC = 0 2 ? where D = 0.5 um/s and k = 0.1 s 1. The boundary conditions are given in the figure. = The solid lines correspond to impermeable walls and the dashed lines to open boundaries with known concentrations. Assume x(in um) varies in the horizontal direction while y (in um) varies in the vertical direction (the lower left corner is the origin). 30 C = 40 10 20 10 C = 100 Justify your approach. Plot the concentration distribution throughout the channel using a 2D colormap with proper labels. Include contour lines