Hi everybody,

I am using maple to solve differential equations.  I have been using maple for some time and came across the infolevel[dsolve] command or option.  I tried finding the answer to my question online but cant seem to come across a straight answer.

My question is what does the integer do that you have to equal the infolevel[dsolve] command?  

For example I am using a book that states

"An even more important diagnostic tool is the infolevel[dsolve] command, which will give information on what methods are used in attempting to solve the ODE when dsolve is applied, even if unsuccessful.  An integer between 1 and 5 must be specified, with generally more detailed information being provided as the number is increased. On applying the dsolve command to ode2, the method of attack is summarized in the following output and, in this case, the general solution y(x) given with two arbitrary coefficients C1 and C2.

What does the integer between 1 and 5 mean or do? In the example above the book uses 5, why 5?.

Thanks,

Matt 

Dear all;

Please in this code Maple, I plot using pointplot3d, some points. I would like to add a name of each point.

The name of each point  P:=[x(i),y(j),t(k)] is u[i,j]^k.  Please how can I add the name at each point. here is the maple code. 

Thanks.

pointplot3d.mw

v1:=[seq(tau[lambda,i]^2,i=1..5)];

latex(v1);

gives:

[{\tau_{{\lambda,1}}}^{2},{\tau_{{\lambda,2}}}^{2},{\tau_{{\lambda,3}}
}^{2},{\tau_{{\lambda,4}}}^{2},{\tau_{{\lambda,5}}}^{2}]

which produces a non natural output in latex, see below.

I wonder it could be improved in the future release?

'latex' command is really useful when i have large algebric output. But having to change this tiny thing at many places is just very time consuming. I'd really appreciate if this can be fixed in updates or future versions of Maple.

Casper

Hi all,

I will use the following dummy example.

with function,f

f:=(xid,yid)->sum(x[i],i=1..xid)*sum(y[i],i=1..yid);

and a complicated term, myterm

myterm:=(f(3,4)+f(2,2))*f(1,1):
myterm:=expand(myterm);

'if' i have some previous knowledge, or know a bit of the term, i can find the structure by doing this

repar:=[f(1,1),f(2,2),f(3,4)];  # Or with more f(xid,yid) terms

tmp:=seq(repar[i]=ff[i],i=1..3);
simplify(myterm,{tmp});          # This is fine, gives me what i want

But, can we go further, and more 'obvious'

Given the fucntion f, same as before, and the same 'myterm'

can I have this
restart:
iwant:=(f(3,4)+f(2,2))*f(1,1);  # as a result, straightforward

so I dont have to go back to 'repar' and find that the terms exactly are.

Thanks,

Download suppress_warnins.mw

Hi All,

I need to suppress the warning messages when i try to pring my Table. Any help will be much apreciated.

Regards,

Xcyborg

I’m trying to figure out how to find a basis for a subspace, V, of R³ defined by V = {(x, y, z)l(2x-3y+6z=0)}

I’m using the student linear algebra module for maple 17

I’ve tried defining the subspace and asking for the basis of V but I always get an error code.

I’ve tried consulting the maple website and looking through their help menu, but can’t find anything that answers how to find a basis... At least a basis from the subspace defined in my problem.

I know how to find a basis for the subspace by hand but not with maple.

Any help will be greatly appreciated. 

(1/2)*(-x-y+sqrt(-3*x^2-6*x*y-3*y^2))/(x+y)

the complex value is always a constant whatever x and y 

i suspected that this is a wrong function

Hello ,

2014_PDE_project.mw

I was using maple as a documentation software and was preparing some notes for my exam. 

The file now opens but most of the content is gone. Even the one which was properly saved.

I see Joe Riel 5676 has solved an earlier case like this. could somebody look at this and correct the file and upload it again please.

Hello: I'm looking over the Help section but I can not find a Maple package that has a command to compute the symmetric sum or the cyclic sum. I just started working with inequalities. Please could anyone recommend a package that allows me to compute expressions related to Muirhead's Inequality (see part 2 of  the answer:  http://math.stackexchange.com/questions/488967/majorization-muirhead)

I want to introduce a matrix of order M by M as for any m, M, pl help as show in file

Help.mw

Dear all;

Please I have maybe a simple question. I need to point these nodes in R^3 using Maple. Thanks for your help.

Tab := [[0, 0, 0], [1/5, 0, 0], [2/5, 0, 0], [0, 1/5, 0], [1/5, 1/5, 0], [2/5, 1/5, 0], [0, 2/5, 0], [1/5, 2/5, 0], [2/5, 2/5, 0], [0, 0, 1/5], [1/5, 0, 1/5], [2/5, 0, 1/5], [0, 1/5, 1/5], [1/5, 1/5, 1/5], [2/5, 1/5, 1/5], [0, 2/5, 1/5], [1/5, 2/5, 1/5], [2/5, 2/5, 1/5], [0, 0, 2/5], [1/5, 0, 2/5], [2/5, 0, 2/5], [0, 1/5, 2/5], [1/5, 1/5, 2/5], [2/5, 1/5, 2/5], [0, 2/5, 2/5], [1/5, 2/5, 2/5], [2/5, 2/5, 2/5]]

This table contains some point of R^3. How can I put these point.

Hi All,

I have a stupid question.

I wanted to distinguish the fuction of map and Map and tried some test.

Sometimes, the Map is identical to the map:

But sometimes not:

Maple Help seems not to explain the difference very clearly.

map_and_Map.mw

I'd appreciate any help on this topic. Thank a lot.

Dear all,

Please help in this question.

Using   I want to plot in R^3, the set of point u[i,j]^k . This point has as cordinate  (x[i],y[j],t[k]).

x := i -> (1/5)*i;  #  x[i] the x-coordinate
y := j -> (1/5)*j; # y[j] the y-coordinate
t := k -> (1/5)*k;  #  t[k] the t-coordinate

The name of point is u[i,j]^k

How can I  plot all the point.

 with(geom3d):

point(u[i,j]^k, x(i),y(j),t(k));

Thank you.

Hi all

Assume that we have the following colomn vector:

and we want to construct the following matrix:

How can we do this with maple????

Note that r is a variable and the dimension of matrix depend on r.

regards

Mahmood   Dadkhah

Ph.D Candidate

Applied Mathematics Department

Dear all;

I need a help in this question.

u(x,y,t) my solution of PDE. x,y space, and t=time.

In the case without t. .i.e. u(x,y). Here is a visualization of the lattice u(x_i,y_j). i=1..3, anf j=1..3.  Please try this example, it's working. 

my question, if i would like to add a third variable t, i.e.  and get u(x_i,y_j,t_k)  on each points. I would like to show the lattice in (x,y) plan for each t_k. I Think I will get many parallel lattice. 

Please can some one, modify this code to get the parallel lattice. Thanks.

with(LinearAlgebra):
with(plots):
with(PDEtools):
interface(rtablesize=20):

Here,

L := 'L':
N := 'N':
g := 'g':
Z := i -> -L+2*L/(N+1)*i;
x[0] = Z(0),x[N+1] = Z(N+1),y[0] = Z(0),y[N+1] = Z(N+1);

N := 4;
L := 1;
r := L/(N+1)/4;
ngon := (n,x,y,r,phi) -> [seq([x+r*cos(2*Pi*i/n+phi), y+r*sin(2*Pi*i/n+phi)], i = 1 .. n)]:
p[1] := display([seq(polygonplot(ngon(4,Z(0),Z(j),r,Pi/2),color=magenta),j=0..N+1),

seq(polygonplot(ngon(4,Z(N+1),Z(j),r,Pi/2),color=magenta),j=0..N+1),

seq(polygonplot(ngon(4,Z(i),Z(0),r,0),color=magenta),i=1..N),

seq(polygonplot(ngon(4,Z(i),Z(N+1),r,0),color=magenta),i=1..N),

seq(seq(polygonplot(ngon(20,Z(i),Z(j),r,0),color=white),i=1..N),j=1..N),

textplot([seq(seq([Z(i+0.1),Z(j),typeset(u[i,j])],i=0..N+1),j=0..N+1)],align={above,right})],view=[Z(-1)..Z(N+2),Z(-1)..Z(N+2)],tickmarks=[[seq(Z(i)=typeset(x[i]=evalf[2](Z(i))),i=0..N+1)],[seq(Z(i)=typeset(y[i]=evalf[2](Z(i))),i=0..N+1)]],

axes=boxed,scaling=constrained,labels=[``,``]):
p[1];