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++.

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" .

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.

I want to plot the volume that is enclosed by the following four areas: x^2+3*y^2=z z=0 y^2=x x^2=y I got so far that my plot looks like this .

eqn :=-p*(l^2-a^2)^(3/2)+a^3*q; l is letter not number if want to find a .....a =??......plz use maple to find the answer i use hand to do it, so the answer is a = (p^(1/3)*l)/sqrt(p^(2/3)+q^(2/3)); but i use Maple it gave me different answer.... so plz help!

This is a program discribing the dynamic behavior of a N pendelum system. When I tried to comparae the result of different initial conditions, I found that every time i run it, the graph will be different.

I know that double pendelum system is a typical chaotic system, is this the reason why i have different result for each time?

Here I attach the file, but plz rename it to ***.mws file before run, b/c i wrote it in maple classic worksheet.

> restart:
interface(warnlevel=0):
with(SolveTools):with( plots ): with( plottools ):
with( DEtools ): with( PDEtools ):
g:=9.81:
x[0][2]:=0:
y[0][2]:=0:
#initialization function of each beam,in this function, we should give 3 variables, n=the number of this #beam, la=the length of this beam, ma=the  mass of this beam,q=the initial angle of the beam,and qt=the initial angle velocity of the beam
beams:=proc(n,la,ma,q,qt)
global l,m,x,y,theta,icq,icqt:
icq[n]:=theta[n](0)=q:
icqt[n]:=D(theta[n])(0)=qt:
l[n]:=la:           #length of  beam n
m[n]:=ma:           #mass of beam n
#initialize the positions of each beam
x[n][2]:=0:        #the x position of the end point of each beam
y[n][2]:=0:        #the x position of the end point of each beam
x[n][1]:=0:        #the y position of the center point of each beam
y[n][1]:=0:        #the y position of the center point of each beam
end proc;
#define the forces that applied on the end of  beam n
Force:=proc(n,fm,fn)
global fx,fy:
fx[n]:=fm:
fy[n]:=fn:
end proc;