Question: How to integrate of a numerically solved differential equation in Maple?

Hello Everyone

I have solved the equation below numerically and now I want to integral p(x)

ode := diff((b+5)*((1/6)*(diff(p(x), x))*h^3+(1/2)*c*h^2)/(b+2)+De^2*((1/10)*(diff(p(x), x))^3*h^5+c^2*(diff(p(x), x))*h^3+(1/2)*c*(diff(p(x), x))^2*h^4+c^3*h^2)/(b+2), x)+diff(h, x) = 0

ics := p(0) = 0, p(1) = 0;
                   ics := p(0) = 0, p(1) = 0
solnum := dsolve({ics, ode}, numeric, method = bvp, maxmesh = 600);
                solnum := proc(x_bvp)  ...  end;
odeplot(solnum, 0 .. 1)

but when I write evalf(int(p,0..1)) I get this result:int(p,0..1). What should I do to have the value?

Please Wait...