A numerical code for steady state temperature variation in a 2D channel subjected to parabolic velocity profile. generate solution for different combination of velocity and thermal diffusivity.
This isn't the kind of question someone can give you an answer withing a few lines. Not clear whether your problem contains transients or is steady-state. No clear whether there is thermal coupling of fluid parameters. Please make sure you fully understand the physics before attempting a numerical solution. Once the physics is clear to you, please read up on numerical techniques for parabolic PDEs (these are the easiest ones to solve.)
This appears an homework. You have the analytical solution for the velocity, then you have to solve the parabolic equation for the temperature that is linear. You need to supply the BC. I suggest to start using a simple FTCS by considering the numerical stability constraints.