Please code in Python

(d) From mass transport calculations, it can be shown that the concentration 'flux' (amount of product flowing per unit area per unit time) in the two directions at each point are: qx(x,y)=DxC(x,y)+U(x,y)C(x,y)andqy(x,y)=DyC(x,y) Here, D=0.1cm2/s is the mass diffusivity constant and U(x,y)=4y(1y) is the given horizontal flow velocity in cm/s. Determine qx and qy using finite differences and show the 'flow' of reaction products using quivers or streamlines. (e) The net horizontal mass flow rate of the product at every x slice is Q(x)=y=0y=Lyqx(x,y)dy Determine Q(x) between x=0.5cm and x=4cm in steps of 0.1cm and plot it. Comment on how you could use this plot to 'design' the length of a channel necessary for a required product flow rate (no calculations necessary; describe in words)