@tayyab You will find both explained in the help:

?dsolve,numeric,bvp

@Axel Vogt You are right. R should rather be:

R:=map([op],indets(u,algebraic^fraction ));

The procedure was made in Maple 14. If you change identical('all') to identical(all) it works in Maple 12.

The unevaluation quotes in identical('all') are superfluous in Maple 14, but gives an error in Maple 12.

@Axel Vogt You are right. R should rather be:

R:=map([op],indets(u,algebraic^fraction ));

The procedure was made in Maple 14. If you change identical('all') to identical(all) it works in Maple 12.

The unevaluation quotes in identical('all') are superfluous in Maple 14, but gives an error in Maple 12.

Did you want w to appear in the expression equ1 or did you want to solve the equation equ1 = w for different values of w or what?

@jean-jacques One parameter - one window ? What should be going on in the M windows? It is not clear to me at this moment.

Incidentally, if you are not aware of Explore you may want to try that:  ?Explore

@jean-jacques One parameter - one window ? What should be going on in the M windows? It is not clear to me at this moment.

Incidentally, if you are not aware of Explore you may want to try that:  ?Explore

I consider it inappropriate to start your question with "In the name of God". Who's god? The Christian god or what? It doesn't matter. What has your question got to do with God?

You are sure to antagonize Christians if they think you might be a Muslim, Muslims if they think you are Christian, and atheists in any case. Thus you are less likely to get responses.

You are more likely to get good answers if you provide more details, in this case the integral in question.

Here are two animates inside each other. I didn't try three.

animate(animate,[implicitplot3d,[eval(f,{a=1/2})=50, x=0..40 ,  y=0..40, z=0..40,axes=boxed,caption=typeset("c = ",c,", b= ",b)],b=0..1,frames=7,paraminfo=false],c=0..1,frames=7,paraminfo=false);

Here are two animates inside each other. I didn't try three.

animate(animate,[implicitplot3d,[eval(f,{a=1/2})=50, x=0..40 ,  y=0..40, z=0..40,axes=boxed,caption=typeset("c = ",c,", b= ",b)],b=0..1,frames=7,paraminfo=false],c=0..1,frames=7,paraminfo=false);

I'm guessing that this is a physics problem of some sort, KE being kinetic energy and PE being potential energy.

Is PI meant to be Pi (i.e. our familiar constant Pi)? Not that it helps though. The integral seems difficult.

In the situation could the integrals be definite integrals, then at least numerical computation could be done for specific values of the constants?

I'm guessing that this is a physics problem of some sort, KE being kinetic energy and PE being potential energy.

Is PI meant to be Pi (i.e. our familiar constant Pi)? Not that it helps though. The integral seems difficult.

In the situation could the integrals be definite integrals, then at least numerical computation could be done for specific values of the constants?

You could use ghostscript.

http://pages.cs.wisc.edu/~ghost/

You could use ghostscript.

http://pages.cs.wisc.edu/~ghost/

OK, something like the result of this?

int(int(exp(x^3)*sin(x)*y,y),x);