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

Is there a solve(identity ...) command for systems? I have 4 identities in the variable t.
They are Q1, Q2, Q3, Q4. I try solve(identity(Q1,Q2,Q3,Q4),t), A, B, C, D), and other forms of it, but
nothing works.
Thank you!

maplatha