Trying Maple on a textbook problem to verify my hand solution. 

But Maple pdsolve hangs with the mserver.exe tallomg almost 100% CPU and over 10 GB of RAM!

I waited for almost 20 minutes. Tried another time, same thing.

It is no problem if Maple can't solve this, but Maple seems to suffer from too many hangs when it is not able to solve a problem. I've had similar problems with dsolve also.

This is on windows 10, 64 bit with Maple 2018.2 With Physics updates version 218

restart;
PackageTools:-IsPackageInstalled("Physics Updates");

218

interface(showassumed=0);
infolevel[pdsolve]:=2;
pde :=  diff(u(x,t),t)=k*diff(u(x,t),x$2)+Q(x,t);
ic  :=  u(x,0)=f(x);
bc  :=  eval(diff(u(x,t),x),x=0)=A(t),eval(diff(u(x,t),x),x=1)=B(t);
sol:=pdsolve({pde,ic,bc},u(x,t)) assuming t>0,k>0;

then

* trying method "_Fn" for 2nd order PDEs
   -> trying "linear_in_xt"
   -> trying "BC_equal_0"
* trying method "_Cn_cn" for 2nd order PDEs
Trying travelling wave solutions as power series in tanh ...
Trying travelling wave solutions as power series in ln ...
Trying travelling wave solutions as power series in tanh ...
Trying travelling wave solutions as power series in ln ...
* trying method "Wave" for 2nd order PDEs
   -> trying "Cauchy"
   -> trying "SemiInfiniteDomain"
   -> trying "WithSourceTerm"
* trying method "Heat" for 2nd order PDEs
   -> trying "SemiInfiniteDomain"
   -> trying "WithSourceTerm"
* trying method "Series" for 2nd order PDEs
   -> trying "ThreeBCsincos"
   -> trying "FourBC"
   -> trying "ThreeBC"
   -> trying "ThreeBCPeriodic"
   -> trying "WithSourceTerm"
      * trying method "_Fn" for 2nd order PDEs
         -> trying "linear_in_xt"
         -> trying "BC_equal_0"
      * trying method "_Cn_cn" for 2nd order PDEs

And here is hangs. Notice that because both ends are Neumann, there is no unique solution to this problem.  So the solution will contain arbitrary constant. May be this is what made pdsolve hang? 

No it is not. Trying with only one end nonhomogeneous  Neumann, and the other end Dirichlet, it still hangs. The problem seems to be with one end is nonhomogeneous  Neumann, which is a function of time. So this hangs also (same place)

restart;
interface(showassumed=0);
infolevel[pdsolve]:=2;
pde :=  diff(u(x,t),t)=k*diff(u(x,t),x$2)+Q(x,t);
ic  :=  u(x,0)=f(x);
bc  :=  eval(diff(u(x,t),x),x=0)=A(t),eval(diff(u(x,t),x),x=1)=0;
sol:=pdsolve({pde,ic,bc},u(x,t)) assuming t>0,k>0;

And this also

restart;
interface(showassumed=0);
infolevel[pdsolve]:=2;
pde :=  diff(u(x,t),t)=k*diff(u(x,t),x$2)+Q(x,t);
ic  :=  u(x,0)=f(x);
bc  :=  eval(diff(u(x,t),x),x=0)=sin(t),eval(diff(u(x,t),x),x=1)=0;
sol:=pdsolve({pde,ic,bc},u(x,t)) assuming t>0,k>0;

But this does not hang

interface(showassumed=0);
infolevel[pdsolve]:=2;
pde :=  diff(u(x,t),t)=k*diff(u(x,t),x$2)+Q(x,t);
ic  :=  u(x,0)=f(x);
bc  :=  eval(diff(u(x,t),x),x=0)=1,eval(u(x,t),x=1)=0;
sol:=pdsolve({pde,ic,bc},u(x,t)) assuming t>0,k>0;

The issue seems to be when one end is nonhomogeneous  Neumann which is function of time. 

Is there a workaround so it does not hang? The complaint is that Maple hangs, and not that Maple unable to solve the PDE.

restart;
#PDE
PDE := (diff(u(x,y), x,x))-1/x*(diff(u(x, y), x)) -x^2*(diff(u(x,y), y,y))= 2;
BCs:=u(0,t) = 0, u(1,t) = 0;

ICs:=u(x,0) =0, D[2](u)(x,0)=2 ;
num_sol := pdsolve(PDE, {BCs,ICs}, numeric); 
BCs := u(0, t) = 0, u(1, t) = 0;
ICs := u(x, 0) = 0, D[2](u)(x, 0) = 2;
num_sol := pdsolve(PDE, {BCs,ICs}, numeric);

 #ODE
ODE:=diff(y(x),x,x)+16*diff(y(x),x)=0;
 BCs:=y(0) =0, y(2)=3 ;
ode_num_sol:=dsolve({ODE,BCs},numeric);

Suppose that we solved numerically an ODE or PDE like above. 
How to calculate values of function y or u at special values x like analytical solutions?

For example; u(1,1)=? y(1)=? etc.

the_code.mw

Why  

eval(diff(u(x,t),x),x=0)=A(t)

gives

But

eval(diff(u(x,t),x),x=L)=A(t)

gives

I was expecting the same syntax in both cases. It seems for numbers Maple uses the first syntax and for symbols it uses the second syntax.

Does one need to worry about this difference?

I would like to animate the motion of a bicycle racer on a classic velodrome track i.e. one with varied vertical and horizontal curvatures along its length.

Is there a source which explains the math expressions which model the shape of such a track?

Let us consider S := [x^2 = A1, x*y = A2, z^2 = A3] as a list. Is there any command to obtain  SS = [A1=x^2 , A2=x*y , A3=z^2] from S?

It seems that the forum activity of the user  Markiyan Hirnyk   is no longer accessible.
Is it a site problem?

My goal is to find a positive linear combination of a vector vec over a set of vectors M1,...,M6. In this case the Mi are the columns of the matrix M.

My reasoning is: solve the linear system M.x=vec, which gives me a parametrized solution, then check if there is at least one positive solution.

M:=Matrix([[3, 0, 2, 2, 1, 1], [-1, -1, 0, -1, 0, -1], [0, 3, 0, 1, 1, 2], [3, 0, 1, 2, 0, 1], [-1, -1, -1, -1, -1, -1]]);

vec:= Vector[column](5, [3, 1, 0, 0, -2]);

x:=LinearSolve(M,vec, free = s);

But when I use the LinearMultivariateSystem function I get an error

LinearMultivariateSystem({x[1]>=0, x[2]>=0, x[3]>=0, x[4]>=0, x[5]>=0, x[6]>=0},[s[1],s[2],s[3], s[4], s[5],s[6]]);
Error, (in Utilities:-SimpleAnd) invalid input: a string/name list is expected for sort method `lexorder`

Can anyone tell me what is wrong in the code?

If I manually type in the expressions, and use variables such as x,y,z... then I get the solution, but I want to put the function into a loop.

I am very open to suggestions on how to decide if a vector is a positive linear combination of other vectors.

Maria

Hi

I have the values of of a function A(x,y,z) in 3d cartesian coordinates as [x[1],y[1],z[1],A(x[1],y[1],z[1])], [x[2],y[2],z[2],A(x[2],y[2],z[2])], [x[i],y[i],z[i],A(x[i],y[i],z[i])],etc...where i vary from 1 to 300 (or higher).

How to plot A(x[i],y[i],z[i]) in 3d.

Thanks

I am trying write a code for this question at here 

https://math.stackexchange.com/questions/170319/how-many-triangles-with-integral-side-lengths-are-possible-provided-their-perim/170325#170325

but I got the wrong answer. This is my code.
 

restart; 
L := []; 
for a to 18 do
 for b to 18 do 
for c to 18 do 
if a < b+c and b < a+c and c < b+a and a > abs(b-c) and b > abs(a-c) and c > abs(a-b) and a+b+c = 36 then
 L := [op(L), [a, b, c]] end if end do end do end do;
 nops(L); 
L

I got 136 triangles. This is a wrong answer. How can I repair my code?

how I can pdsolve these equations?

moadelat.mw

Download moadelat.mw

I have multiple accounts. I would like to delete them. Please tell me how to delete it. Thank you.

How can I plot functions Vc=(0.5+t)^n

with t=z/h, 

Hai any one help me to plot the curves and remove the errors  i am attaching the codes and  sample graph. thanks in advanced.
 

test.mw

Hi all,

I'm struggling with solving a delayed differential equation with time-dependent delay.

I have already use dde solver with matlab and find it cannot reaches my numerical accuracy. 

Does maple can solve this problem? I explore a little bit and don't find there is a document to guide how to tackle with this problem.

Constant delay cannot solve my problem.

Thank you very much,

Peter