restart;

pde := diff(u(x, y), `$`(x, 2))+4*(diff(u(x, y), `$`(y, 2))) = -x*y;

ibc := {u(0, y) = 0, u(2, y) = 0, u(x, 0) = x*(4-x^2), (D[2](u))(x, 0) = 0};

sol := pdsolve(pde, ibc, numeric);

sol:-plot3d(x = 0 .. 2, y = 0 .. 1, numpoints = 2000, axes = boxed);

The final command produce a figure with some sigularity.

how to fix it? I need help, thx you very much.

Hello All ,

I would liek to know if ther is a way to tabulate the results in a matrix form , like matlab does .Each time i copy paste it in excel and do the text to colum to sort out the results.It would save me a lot of time if anyone has an idea about this.

Like in the file i have results of various parameters displayed but how do i tabulate it.

Steam_800_below_crit.mw

rsolve(a(n)=2/a(n-1)+1,a(1)=1),

thx in advance.

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...