This question was asked before but because of the curiosity, I

bring it in the light again. We have a system of odes

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,