@tomleslie Thanx for the time, actually the coupled pdes are highly nonlinear and have no analytical solution . i wants to solve it numerically by using Finite difference method for which i need maple code..so will be very thankful to u for any kind of help..yes i have appropirate boundary conditions