I'm trying to solve a partial differential equation with two boundary conditions below. The general solution contains arbitrary functions of the non-differentiated variable. These functions are solved for and assigned but do not appear in the final solution return. Can anybody help me with this?

**> restart;**

**> l:=lambda;**

Hi, I want to draw the contourplot of this two family of lines x*y and x*y^2. with(plots): contourplot( {x*y, x*y^2}, x=0..10, y=0..20); I would like that the first function (x*y) is colored in red ad the secod (x*y^2) in blue. How could I do it? I tried with "colouring= [red, blue]" but the coloring is not what I want. Thanks in advance for the help. JJ (By the way..I wish you the best to all of you during 2006!!)

This is the simple code that I came up with (after some frustration) to solve the problem of separating the sequence of positive integers from 1 to 20 into one sequence of primes and one sequence of nonprimes, and eventually displaying those sequences. I'm just curious to see what other ways this same question could have been solved.
>

Hi. I'm a very new Maple user, but I don't think that this qualifies as a Newbie question. Has anybody had problems when trying to write directly to a PS file from within the Maple GUI on OS 10.4.3 for Mac? Consider the following:
currentdir(kernelopts(homedir)):
plotsetup(ps,plotoutput="testplot.ps"):
plot(sin(x),x=0..2*Pi);
Running this code from the command-line works as expected, but running this from the GUI interface usually returns:
Error, Error in device driver: plot terminated
The strange thing is that about 10% of the time, the above code actually will work from the GUI. I cannot produce failure or success consistently in the GUI. I have tried this on both a G4 and a G5. Am I the only user who has had this problem?

I try to deduce distribution fonction from a parametrized set of points with Maple.
At this stage I have a dataset define through a relationship like y=f(x), and I want to obtain the distribution of y points, given that x in [-A,A]. The theoretic formula is : F[Y](z) = P(y<z).
i.e. F[Y](z) = Int(delta[f(x)< z],x=-A..A);
and I use the **piecewise** Maple function to implement it ( Int(piecewise(f(x)< z,1,0),x=-A..A) ), but for Maple :
Int(piecewise(f(x)< z,1,0),x=-A..A) = piecewise(Int(f(x)< z, x=-A..A),1,0))
which is totally different !

I was able to create a stand-alone application from a M-file with MatLab 6.0 R12. Is it possible to do the same with Maple 9.5? I mean, can I take a Maple worksheet and convert it into a executable file which can be run on a PC that does not have Maple 9.5 installed? I would like to know since it is easier for me to create GUI applications with Maple 9.5 than with MatLab.

Using the following integrand

.71428491807021e-1*(1/2-1/2*erf(.50253255206611e-2*xi*2^(1/2)))/Pi^(1/2)*

exp(-1/2*(.14285698361405*xi-5.6995353453900)^2)*2^(1/2)*(1/2+1/2*

erf(1.1171516020037*(5.4091420603773+.44756682898115*xi)*2^(1/2)))

I get -0 for Int( %, xi= -99.496042185589..infinity): evalf(%,14);

which certainly is false. But if I cut off at a reasonable upper bound

(say exp(...) <= 1E-16) I get what I expect (up to rounding errors).

I consider that as a bug and wonder whether it is in the NAG library

or through Maple calling it - any explanation?

Edited to add: `Mapl

I was wondering if some one had the knowledge as to why the plot3d command will not work for the dirichlet elliptical wave equation, at fixed time t, yet maple is able to evaluate this function at specific values for theta and r. (the coord system for the 3d plot should work with "ellcylindrical" but does not)

the modes of the function look like (for Cosine-Elliptic)

MathieuCE(j,q,theta)*MathieuCE(j,q,I*r)*cos(lambda*t)

j=0, 1, 2...n where q is a 'zero' value of MathieuCE(j,q,I*1) 0<theta<2*Pi, 0<r<1, t>=0

the error stated is "Plotting error, MESH must be a list of lists of vertices, or an hfarray"

I would like to split a polynomial into even and odd terms. Has this capability been provided in a package? PolynomialTools seems the obvious choice, but doesn't do this. Here's one approach

SplitPolynomialEvenOdd := proc(poly::polynom(anything,v), v)
description "return the even and odd parts of a polynomial in v";
local p;
p := collect(poly,v);
if p::`+` then
return selectremove(t -> degree(t,v)::even, p);
elif degree(p,v)::even then
return (p,0);
else
return (0,p);
end if;
end proc:
SplitPolynomialEvenOdd(x^2 + 3*x + 1, x);

i want to know that is ther any procedure or builtin functin which return the exponent of a variable.For example if i have a variabl x^2 result will be 2.and also i want to know the procedure by which if i have two lists of variables we can compare there type of variables and cobine the coefficients of same variables.For example
>s:=[v1,2*x^2*y,v2,1*x*y^2,v0,3*x^3];
>p:=[a*x^3,b*x*y^2,c*x^2*y];
i want to get result
v1,2+c ,v2,1+b, v0,3+a
how can i get this result any body help me

I have a system of matrix equations and would like to solve it for a certain vector, but without stating anything else but the names of the matrices and vectors involved. This is to be used for further studies in a numerical matlab model where the matrices and vectors are specified. Example:
Let A, B and C be regular matrices where:
A*B=C
Solving for B we get:
B=inv(A)*C
How can I make maple do this for me without specifying the elements in A, B, and C.

I have maple10 and I'm noticing something that's confusing. I recently played with a couple options from the tools menu, namely (tools -> options) then under the display tab the I played with the input display choices, and under the interface tab I played with the "Default format for new worksheets choices". However for this question I have the input display on 2-D math notation and I'm using worksheet mode as the default for new worksheets. I'm not sure if I applied those choices globally or applied them to the session. I have the italics choice button chosen on the formatting row and also clicked the choice "math" next to a drop down menu full of other lists.

I Need to Generate Primitive and Irreducible Polynomials in Galois Extension Field GF(2^32). How do I do it?
Any Pointers to any code / theory is welcome.
Thanks

I want to 'revert' the product rule from differentiation:
I want to collect terms like f(x)*diff(g(x),x)+g(x)*diff(f(x),x)
into diff(f(x)*g(x),x)

I want to combine terms in a differentiation: combine(Diff(f(x)+g(x),x)-Diff(g(x),x)) It does not work!