Question: Nonlinear ode equations

Hi,

I have a system of nonlinear equations in conjuction with boundary conditions, as below:

(all coefficients are constant.)

dsys3 := {a2*(diff(f1(x), x, x))+a3*(diff(f2(x), x, x, x))+a4*(diff(f2(x), x))+a6*(diff(f3(x), x))+a7*f1(x)+a1N*(diff(f3(x), x, x))*(diff(f3(x), x))+a2N*(diff(f3(x), x))*f3(x)+a3N*(diff(f3(x), x))*f3(x) = 0, a8*(diff(f2(x), x, x, x, x))+a9*(diff(f2(x), x, x))+a10*f2(x)+a11*(diff(f1(x), x, x, x))+a12*(diff(f1(x), x))+a13*(diff(f3(x), x, x))+a14*f3(x)+a4N*(f3(x)*f3(x))+a5N*((diff(f3(x), x))*(diff(f3(x), x)))+a6N*(diff(f3(x), x, x))*f3(x) = 0, (diff(f3(x), x, x, x, x))*a16+(diff(f3(x), x, x))*a18+(diff(f3(x), x, x))*a19+(diff(f1(x), x))*a22+(diff(f1(x), x))*a23+(diff(f2(x), x, x))*a24+f2(x)*a25+f2(x)*a26+f3(x)*a27+f3(x)*a29 = 0, f1(0) = 0, f1(1) = 0, f2(0) = 0, f2(1) = 0, f3(0) = 0, f3(1) = 0, ((D@@1)(f2))(0) = 0, ((D@@1)(f2))(1) = 0, ((D@@1)(f3))(0) = 0, ((D@@1)(f3))(1) = 0}:
dsol5 := dsolve(dsys3, 'maxmesh' = 500, abserr = .1, numeric, range = 0 .. 1, output = listprocedure)

but i encounter an error as:

Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system

what is your suggestion to solve this nonlinear equations?

Thank you.