## repartition function in Maple...

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 !

## strange numerical integration...

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?

## Plotting Modes of Elliptical Wave Equation...

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"

## Splitting a polynomial into even and odd terms...

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);
```

## maple programming...

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

## How do I solve a system of matrix equations withou...

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.

## Primitive Polynomials In GF(2^n)...

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

## 'revert' the product rule...

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)

## Problems with Compiling Procedures...

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]);

## A simple Maplet example

Maple
Here is a simple Maplet example,

## Pattern Based Transforms

Maple
This topic was developed in a discussion with Mad Math. In Maple, pattern based transforms can be done using following Subs procedure

## OpenMaple example with Dev-C++

Maple
Here are the steps for making StartMaple example with Dev-C++.
1. Start Dev-C++.
2. Click File - New - Project, select Console Application, check C Project radio button and type the name of the project - something like SMexample.
3. Copy the source from ?StartMaple help page and paste it in main.c (instead of pretyped text there).
4. Click Alt+P (or Project - Project Options), go to Directories, select Include Directories and add C:\Program Files\Maple 10\extern\include there.
5. 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" .