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.

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;

I'm having difficulties in using the Heaviside function in Maple.
My heaviside function is (t-tau), and I'm having trouble integrating tau from zero to t. After integration, the expression still shows an integral.
Is it because when tau = t, the heaviside becomes zero? Any and all help is appreciated.
Attached is the actual sheet. The 2nd to last and last steps are the steps or concern....
**** This problem has been overcome *****

How do use an equation to create a list/array?

e.g y=x*5+1

for x from 0 to 1 in increments of 0.1

answer should look similar to:

[0,1]

[0.1,1.5]

[0.2,2]

[0.3,2.5]

Has anyone been able to use the statopts option in the Bootstrap command which is in the Statistics package in Maple 10? Note the following sequence of commands
> with(Statistics)
> Bootstrap('Mean', [1, 2, 3, 4, 5], replications = 6, output = array)
print(Matrix(1, 6, {(1, 1) = 2., (1, 2) = 3.39999999999999990, (1, 3) = 2.20000000000000016, (1, 4) = 2.79999999999999982, (1, 5) = 3.60000000000000008, (1, 6) = 2.20000000000000016}));
[2., 3.39999999999999990, 2.20000000000000016, 2.79999999999999982,
3.60000000000000008, 2.20000000000000016]
According to the documentation which is attached as a PDF file it appears that this command should work.

Does any body know if there is a function like interp1 available in Maple or at least what is the easiest way of establishing a look-up function in Maple.
I have a n*2 matrix of data (x, and f(x) in each column) and I would want to interpolate data from this matrix.
This is actually one of the applications of interp1 function in MATLAB...