Maple 2015 Questions and Posts

These are Posts and Questions associated with the product, Maple 2015

I want to solve numerically the PDE:

u_xx + u_yy= = u^{1/2}+(u_x)^2/(u)^{3/2}


My assumptions are that  |sqrt(2)u_x/u|<<1 (but I cannot neglect the first term since its in my first order approximation of another PDE.


So I tried solving by using pdsolve in maple, but to no cigar.


Here's the maple file:

PDE := diff(diff(u(x, y), x), x)+diff(diff(u(x, y), y), y) = u^(1/2)+(diff(u(x, y), x))^2/u^(3/2); IBC := {D[1](u)*(1, t) = 0, D[2](u)*(x, 1) = 0, u(0, t) = 1, u(x, 0) = 1}; pds := pdsolve(PDE, IBC, type = numeric); pds:-plot3d(t = 0 .. 1, x = 0 .. 1, axes = boxed, orientation = [-120, 40], color = [0, 0, u])

diff(diff(u(x, y), x), x)+diff(diff(u(x, y), y), y) = u^(1/2)+(diff(u(x, y), x))^2/u^(3/2)


{D[1](u)*(1, t) = 0, D[2](u)*(x, 1) = 0, u(0, t) = 1, u(x, 0) = 1}


Error, (in pdsolve/numeric/process_PDEs) all dependent variables in PDE must have dependencies explicitly declared, got {u}


Error, `pds` does not evaluate to a module





Hello All,

I looked through the Maple help on PDE systems and pdsolve and the physics problems that appear there. THere are a number of single-PDE cases with initial / boundary conditions; but I couldn't find PDE systems with ics/bcs.


Would you have a (simple) example of a PDE system with its initial / boundary conditions? I am attempting to build understanding of the syntax and different options of "pdsolve". Examples seem to be a great way to learn how to solve PDE systems. One can then pdsolve the PDE system without ics; then add them, try different options etc.


Thank you!


I have an arrays of data. One for x values, and one for y values. How can I obtain a numerical integration of y for a range of x values?

I have tried defining a function of X using ArrayInterpolation(x,y,X) and then calling evalf(Int(f,xmin..xmax)) but that gives an error message. (I don't seem to be able to paste into this window) The error message says

"Error, invalid input: evalf expects its 2nd argument, n, to be of type posint, but received numeric."

I thought I was using a form of the equation right from the help system.

I also tried the 2D version of integration, but it returns the difference of my limits times my function name.

I aslo tried AdaptiveQuadrature, but I can't get that to work either.

-Mike McDermott

Newbie Maple user



I'm new at maple and have a problem/question with the rkf45 numerical ODE Solver.

At first, my computer need a lot of time to calculate an analytic solution.
Therefor, I use the numerical way.

I have the following second order ODE:
ODE:=m*((D@@2)(x))(t)+d*(D(x))(t)+k*x(t) = d*(eval(diff(y(x), x), x = t))+k*y(t)
where y(t) is a realy big piecewise function, defined by me.

My initial conditions are:
x(0) = 0, (D(x))(0) = 0

With dsolve, I get the solution x(t) and the first derivative x'(t). I'm able to plot them with odeplot.


Problem 1:
I need also the second derivative x''(t).
On this page: is an example (eq 13 and 14) where the second derivative is useable, but this doesn't work with my differential equation.
I have add 
(D(D(x)))(0) = 0
to my initial conditions but then, I got the error that only 2 initial conditions are required.
What could I do, so that rkf45 returns also the second derivative?

Problem 2:
And in addition to this, I want to calculate with x(t), x'(t), x''(t) but I found no way to use them.
Only plots are possible.
If I reduce y(t) to a minimum, I can do everything with the analytic solution: plot, d/dt, d2/dt2, +, -, ...
I tried also to convert the procedure to a function but in this case, there is no way to derivate it.

Many thanks...



  My question is related to this one

but not exactly the same.


  Suppose I have lists




  The aim of sorting is first, arrange by ascending order of the first variable, then the second, then the third.


I could program manually by explicitly comparision, is there any robust way in maple can do it (possibly a few build-in function)?


Thank you very much!





Just trying Maple again after 10+ years, so please bear with me. I am basically trying to do a symbolic integration, where the output should be the antiderivative of the function I am integrating. However, I have not been able to succesfully get the output I expect.

I defined a function Cp(T), with constants C1, C2, etc. like so:

C__p := T -> C__1+C__2(C__3/(T*sinh(C__3/T)))^2+C__4(C__5/(T*cosh(C__5/T)))^2 

But when I integrate the function using int(Cp(T), T=Tref..Tsys) it does not output the antiderivative as I would expect. This is what I am looking for, but it just gives me the equation within the integral sign without symbolically solving integral. Can this be done?

It should be,

int(Cp(T), T=Tref..Tsys) = C1*(T2-T1) + C2*C3[coth(C3/T2) - coth(C3/T1)] - C4*C5[tanh(C5/T2)-tanh(C5/T1)]

Trivial integrals such as int(x2,x) = 2x solve okay, so I am assuming I setup the problem incorrectly. I just cannot figure out what I did wrong, and it is driving me crazy. I already wasted more time than is healthy on this. Any help would be greatly appreciated. Thanks.


I am trying to use implicit plot. The plot is OK but I want to put labels such as what Latex produces :

$\frac{\Omega}{\omega_n}$ for x axis and $a_0 \mathrm{(m)}$ for y axis

How can I apply this in my maple code as below:

plot1:=implicitplot(a3, Omega_r=1.5..2.5, a=0.00000001..0.1, labeldirections=[horizontal, horizontal], axes=boxed, labels=["W/w_n",typeset("a_0 (m)")], labelfont=[SYMBOL]):

What I have put as bold does not work for me, it is making everything in Greek :) . I want combination of Greek and math.