Abstract Since fluid flow is a branch of mechanics it satisfies a set of well documented basic laws (the governing laws), and thus a great deal of theoretical treatment is available. But viscosity and geometry are two chief obstacles to a workable theory of fluid flow. Based on this the basic equations of fluid motion (conservation equations) are too difficult to enable the analyst to attack arbitrary geometric configuration. Therefore, to solve fluid flow problems by minimizing these obstacles, we have considered the geometry of the flow, the type of the fluid that flows, the flow type itself, the dimension of the flow and the solution methods. In this project work an endeavor has been made to determine the velocity and temperature distribution of time independent laminar flow of viscous fluids with constant density, viscosity and thermal conductivity of the Newtonian type by solving the governing equations of fluid dynamics, namely the continuity, Navier-Stokes and energy equations. Heat transfer through conduction is assumed to be involved in the fluid flow and the temperature distribution has been obtained from known velocity fields by solving the energy equation separately. The answer has been gotten by solving some simple fluid flow problems using analytical solution method, and mat lab codes have been used to draw the solution graph. Numbers can be used interchangeably to plot the graphs of the velocity and temperature profiles with changing the partition size (δx, δy, δz, etc) and the gap distance (h) of the given geometry through which the fluid flows by considering the symmetry.