fprime0 := diff(BfT[0]*z^ls*(ls*p[0]*z^lk+1-ls*p[0])/(.99*(1-(1-.989*f[0])^5)), z); 'fprime0(1)' = subs(z = 1, fprime0)

 je veux savoir comment je calcule les valeurs de fprime0 s'il vous plait je suis coincé depuis une semaine et j'ai pas pu terminer mon travail le resultat me donne 

"Large sortie de plus de 1000000 noeuds" 

I want to solve two equations with all roots.

eq1:=-7.506556000*10^(-23)*omega^12+2.243366854*10^(-78)*omega^34+3.651941282*10^(-50)*omega^25+1.145341958*10^(-105)*omega^41-3.837716260*10^(-102)*omega^40+9.952750000*10^(-10)*omega^5+2.520327960*10^(-113)*omega^43-5.961864252*10^(-98)*omega^39+2.934254874*10^(-60)*omega^29+3.877088344*10^(-109)*omega^42-3.482296310*10^(-82)*omega^35+2.846747460*10^(-123)*omega^45-1.089582984*10^(-70)*omega^32-5.927854900*10^(-48...

Hi there!

I want to  plot phase portraits of this fractional system in maple please help me.

                          (alpha_1)
                        ...

eq1:=diff(f(eta),eta,eta,eta)+f(eta)*diff(f(eta),eta,eta)+beta*(1-diff(f(eta),eta)^2)=0;

eq2:=diff(theta(eta), eta, eta)+pr*f(eta)*(diff(theta(eta), eta))+pr*Ec*(diff(f(eta), eta, eta))^2

+(2-beta)*pr*Q*theta(eta) = 0;

bc:=theta(9)=0,D(theta)(0)=-Nc*(1-theta(0)),f(0)=1,D(f)(0)=lambda,D(f)(9)=1};

Para:={beta=1,Ec=0.5,Q=-0.5,Nc=5,pr=3}:

How to plot -D(theta)(0) vs lambda?

Thanks

Hi every one,

I am facing a situation here. Whenever, I sub n = 0 or 1, the following ode has a solution. 

But when I adopt other values except the above ones then the maple gives me the error

"division by zero". 

eq11:=(diff(((abs(diff(f(eta),eta,eta)))^(n-1)*diff(f(eta),eta,eta)),eta)

+n/(n+1)*f(eta)*diff(f(eta),eta,eta))=0;

resx:=dsolve({eq11,f(0)=1,D(f)(0)=0,D(f)(12)=0},numeric,method = bvp[midrich],

Hello,

I'm solving a system of 3 second order ODE's for 3 unknown functions, using dsolve. But in the output, I receive various sets of solutions which contain only 3 constants of integration. A friend solved the same system but in Mathematica, and received more constants of integration. 

I guess this has to do with the issue that Maple omits additive constants of integration almost everywhere. 

Is there some solution? i.e. some way of using...

I am looking for a numerical scheme which can easly handle such type of problem.

I already had the analytical results of this problem, now I want to know how to treat this problem numerically.

Any idea, then please let me know.

Thanks

Hi every one,

I am trying to find a way to plot a dependent variable against a parameter. 

The data I want to use for this ploting should be from the numerical results

obtained by using shooting method.

restart:ShootLib := "C:/MP/Shoot9/":
> libname := ShootLib, libname:
> with( Shoot ):
> with(plots):with(LinearAlgebra):
> A:=0:
> FNS:={ f(z), G(z), H(z)}:
> ODE:={diff(f(z),z)  = G(z),

Maple 13 does not calculate the definite integral

int(x^4/(4*x^5+2), x=0..1);

although the corresponding indefinite integral is calculated correctly. The reason?

Thank you.

restart; b := 6; de1 := diff(f(eta), eta, eta, eta)+f(eta)*(diff(f(eta), eta, eta))-M^2*(diff(f(eta), eta))-(diff(f(eta), eta))^2 = 0, (1+(4/3)*Nr)*(diff(theta(eta), eta, eta))+epsilon*(theta(eta)*(diff(theta(eta), eta, eta))+epsilon*(diff(theta(eta), eta))^2)+Pr*(s*theta(eta)+f(eta)*(diff(theta(eta), eta))+Nb*(diff(theta(eta), eta))*(diff(phi(eta), eta))+Nt*(diff(theta(eta), eta))^2) = 0, diff(phi(eta), eta, eta)+Le*f(eta)*(diff(phi(eta), eta))+Nb*(diff(theta(eta), eta, eta...

restart; with(Student[NumericalAnalysis]); with(PDEtools, casesplit, declare); x := .5; a := 1; b := .5; d := 1; Q := 2; omega := .4; h1 := 1+a*cos(x); h2 := -d-b*cos(x+omega); F := Q-1-d; de := 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, diff(theta(y), y, y)+Nb*(diff(theta(y), y))*(diff(phi(y), y))+Nt*(diff(theta(y), y))^2 = 0, diff(phi(y...

Hi every one,

How can I get bouble Curl of this equation

In the above equation, the first term has the Laplacian, the second one has the gradient and in the term "k" is the unit vector and also u=(u1,u2,u3). sqrt(R) is a constant. 

hi guys,

say i have an expression like this:

Hi,

I want to solve this  fractional differential system in maple please help me.

                          (alpha)
                         D        y[0](x) = 0