restart:with(plots): a := 1; b := .5; d := 1; omega := .4; h1 := 1+a*cos(x); h2 := -d-b*cos(x+omega); F := Q-1-d; de:={alpha*(diff(f(y), y, y, y, y))+G*(diff(theta(y), y, y))+B*(diff(phi(y), y, y)) +6*beta*(diff(f(y), y, y))*(diff(f(y), y, y, y))^2+3*beta*(diff(f(y), y, y, y, y))* (diff(f(y), y, y))^2 = 0, diff(theta(y), y, y)+Nb*(diff(theta(y), y))*(diff(phi(y), y))+Nt*(diff(theta(y), y))^2 = 0' diff(phi(y), y, y)+Nb*(diff(theta(y), y, y))/Nt = 0, f(1+cos(x)) = (1/2)*Q-1, f(-1-.5*cos(x+.4...

Is it possible the contour plot of these coupled equations

restart:with(plots): Nb := .2; Nt := .2; h1 := 1+a*sin(x);h2:= -b-d*sin(x); a := 1;b:=0.1;d:=0.1;F:=1-Q-d;Q:=-1;beta:=0.1;alpha:=0.1;G:=0.5;B:=0.5; eq1 := diff(alpha*f(y),y,y,y,y)+G*(diff(theta(y),y,y))+ B*(diff(phi(y),y,y))+6* beta*(diff(f(y),y,y))*(diff(f(y),y,y,y))^(2)+3* beta*(diff(f(y),y,y,y,y))*(diff(f(y),y,y))^(2)=0; eq2 := diff(theta(y), y, y)+Nb*(diff(theta(y), y))*(diff(phi(y...