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 !

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 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'm having problems using Compiler:-Compile() on the following procedure:
findallroots:=proc(eqs,x,rng::range(numeric))
local roots,pts,i;
roots:={fsolve}(eqs,x,rng,'avoid'={x=lhs(rng),x=rhs(rng)});
if roots={} or not roots::set(numeric) then
NULL
else
pts:=sort([op(rng),op(roots)]);
op(roots),seq(procname(eqs,x,pts[i-1]..pts[i]),i=2..nops(pts))
end if;
Error, (in IssueError) only named functions are supported
If you could help me with what this error is referring to and what I could do to over come it that would be most appreciated.
Thank you for your help

I am a novice in using maple-matlab link.
I followed the help manuel instruction in maple, as below.
The command "with(Matlab);" executed successfully.
And I found new Matlab command window activated.
Even openlink() executed without error.
But I can't use any matlab command in maple due to link lost error below.
My maple version is 8.00.
And my matlab version is 6.5.0.
I will be very grateful for anyone's help.
> restart;
> with(Matlab);
[chol, closelink, defined, det, dimensions, eig, evalM, fft, getvar,inv, lu, ode45, openlink, qr, setvar, size, square, transpose]
>
> A:=matrix(2,2,[1, 1, 1, 2]);

Here is a simple Maplet example,

This topic was developed in a

discussion with Mad Math.
In Maple, pattern based transforms can be done using following Subs procedure

Here are the steps for making StartMaple example with Dev-C++.

- Start Dev-C++.
- Click File - New - Project, select Console Application, check C Project radio button and type the name of the project - something like SMexample.
- Copy the source from ?StartMaple help page and paste it in main.c (instead of pretyped text there).
- Click Alt+P (or Project - Project Options), go to Directories, select Include Directories and add C:\Program Files\Maple 10\extern\include there.
- In Parameters (in Project Options) in the Compiler window type -D_MSC_VER and in the Linker window click Add Library or Object and add "../Program Files/Maple 10/bin.win/maplec.lib" .

I am trying to invoke Maple from a C program, using OpenMaple, copying the Maple help screen example (see StartMaple). I've of course changed the path to the header file.
However, I get a linker error: 'StartMaple is an unrecognized function'. I take this to mean that the .lib and or .dll maplec files are not being linked in, but I can't seem to figure out how to.
I'm using Dev-C++ as the compiler. Anyone have experience with this

Is there a straightforward procedure to have a radical appear on the plot of a function? E.g, how can one display (2, sqrt(2)) in the plot but with the actual square-root symbol appearing. I've read in the knowledge base about the use of textplot. However, the routine is limited to exponents.