Something like this

Chebyshev polynomials have the form:

Tn(x) = cos(n*arccos(x))

orthogonality interval is (-1, 1) for n = (2, 3, 4,

In the name of God,

Hello dear experts,

I have a nonlinear multivariate system. By running fsolve and with the same initial values, 2 times, at the first time it yields the result but at the second time it doesn't solve. All of the settings in the both fsolve executions are the same.

Do you have any idea?

Best regards,

How do you parse a string containing quoted material?

If the material is read in from a file
readline("f:/test.txt")
and the line contains quotes, if you parse you get an error.

Recursive algorithm for finding the sum of the first n odd positive integers.

 procedure   odd_sum ( n: positive integer )

    if   n = 1   then   return   1

    else   return   

I think i'm missing something easy her , but I am trying to plot a circle and cannot figure out

what I am doing wrong. The equation for the circle is:

x^2+(y-3)^2=25

I define the equation

c := x^2+(y-3)^2 = 25

Solve it for y so I should be able to plot it with plot(cplot); :

cplot:=solve(c, y)

But what I get as answer is:

3+sqrt(25-x^2), 3-sqrt(25-x^2)

The equation of the circle is part of a homework assignment,...

Hi,

I'm working at the moment on percolation theory and, in order to get interesting images I use Maple to produce random spheres packing using the code below:

with(RandomTools):
with(plots):

x:=Generate(list(integer(range=0..20),100)): # random coordinate for x parameter
y:=Generate(list(integer(range=0..20),100)): # random coordinate for y parameter
z:=Generate(list(integer(range=0..20),100)): # random coordinate for z parameter
R:=Generate(list(float(range=1..3...

Hi,

I m trying to plot the gravitational potentialfield of two masses in maple. Im doing my work with inspirtion from this paper:

http://ocw.mit.edu/courses/aeronautics-and-astronautics/16-07-dynamics-fall-2009/lecture-notes/MIT16_07F09_Lec18.pdf

Page 4 and 5. Im able to plot the first graph on page 5, the contour plot of the gravitational potential in graph, by setting y=0. But when i try and plot the function of 2 variables in maple...

Hi all,

I read the help on using parameters, then I tried to modify the example

f := proc(a, b, c)

    if _params['b','c'] = NULL then
    a
    elif _params['c'] = NULL then
    a + b
    else
    a + b + c

    end if

end proc:

I can't really spot where it went wrong, but it just wouldnt execute as I would expect...

Hi all,

A := [[1,2],[3,4]];

A1:=convert(A, Matrix);

I am very happy with conversion between Matrices and Vectors.

However, I does not seem to be able to do this backwards to [row] or [column] Vectors???

Say how do I get back to A from A1???

I tried to use

convert(A1,Vector);

convert(A1,Vector[row]);

convert(A1,Vector[row],2);

And many others...

Hi, when I was learning to use the taylor-command I always get an error message that says: Error, (in series/abs) unable to compute series.

The exercise I was trying to solve was as easy as this:

> f:=x->sqrt(abs(x));

> plot(f(x),x=-5..5);

> taylor(f(x),x=1);

Error, (in series/abs) unable to compute series

Can someone tell me what am I doing wrong?

Many Thanks!

I use Hermit as example and try the following way, it can not sumtohyper,

how to get the hyper form from a differential equation?

restart;
with(gfun):
test4 := diff(P(x),x$2)-2*x*diff(P(x),x)+2*n*P(x)=0;
test4 := {diff(P(x),x$2)-2*x*diff(P(x),x)+2*n*P(x), P(0) = 1, (D(P))(0) = 1};
newtest4 := borel(test4, P(x), diffeq);
rec4 := diffeqtorec(newtest4, P(x), a(v));
rec4 := diffeqtorec(test4, P(x), a(v));
genfun := rsolve(rec4, a(v), 'genfunc'(z));

restart;
with(SumTools):
an := -n+L;
bn := -x+L;
gen := Summation(Product(an, L=0..k-1)*Product(bn, L=0..k-1)/k!*(-1/a)^k, k=0..infinity);
genfun1 := subs(a=z,simplify(gen));
genfun2 := exp(z)*(1-(z/mu))^x;

correct is genfun2

but using famous hypergeom of charlier written in many books and papers, it can not summation to a correct generating function

Hello everyone,

                       I have a par.mw file which contains a list of parameters. Now I have another reac.mw file in which I have to use some of those paramters. My question is what is the way for this without searching par.mw for those specific parametes values and then putting manually in reac.mw.

Following an earlier thread (Collect-And-Sqrt) where a solution using applyrule was proposed (and that solution works), I tried to use applyrule to implement some of the trig half-angle relationships. Immediately I ran into trouble:

Say I want to use the rule

r:=cos(th::algebraic)=2*cos(th/2)^2-1;

so I say

applyrule(r,cos(alpha));

and Maple promptly goes into an...