Add below steps in you Maple code

res:= dsolve(subs(m=0.5,{bcs1,bcs2, eqs1, eqs2}), numeric,output=listprocedure,maxmesh=1024): 
r := eval(f(eta), res):    
p := (eta,x) -> x*r(eta):  

pp1:=plot3d(p,0..6,-15..15,contours=40,style=patchnogrid, grid=[60,60], orientation=[-360,-180], lightmodel=light4, shading=zhue,  transparency=0.3,style=contour,axes=boxed):
display(pp1) ;