National Public Radio in the USA carries Car Talk, a humorous phone-in program in which Tom and Ray Magliozzi (Click and Clack, the Tappet Brothers) diagnose and offer solutions for mysterious auto-related maladies. It's an amusing hour on Saturday mornings.

I am trying to solve a system of nonlinear 36 differential-algebraic equations numerically in Maple. The last line of my code is the dsolve(.) command which is used to solve the DAE system. My problem is that it takes forever to evaluate the dsove(.) command. Maple does not give me any errors/warning, however, it keeps evaluating, and evaluating...nothing happens...

I have used stiff solvers as well hoping to make the process faster but that dowsn't work either.

Hey there, this is my first post here, so hello to everyone.

I have (so far) one question. I have a circle equation like

x²+y²-6y+6=0

I know, that on a paper, I can convert it to the "normal" circle equation, from which I can read the middle point and radius of the circle. And now my question is, how to convert in in maple? I am using the 11th version.

I have tried solving this equation and (x-a)

How can I plot standing waves with extra parameters. That is how can I plot something like 2A*sin(kx)*cos(wt) and vary the parameters to see what would happen if different types of waves traveling in different directions collided.

I would like to export-to-file a list containing equations, variables, and constants, e.g.

xyz := [dx = x*y-(1/3)*(y-z), dy = -0.01*x, dz = z/10];

I have used fprintf to write to a file. The closest I got to what I want is by doing:

fd := fopen("xyz.m", WRITE) :
fprintf( fd, %s,  CodeGeneration[Matlab](xyz,'output'=string) ):
fclose(fd) :

Using CodeGeneration[Matlab] was the only way I found to export a list like xyz above, simply writing fprintf...

Dear Maple Team,

     I found a very strange thing in my Maple 14: When I typed the two command line as follows:sum(n^4, n); expand(sum(i^4, i = 1 .. n)), the results are not the same. Now I copy them here.

                   1  5   1  4   1  3   1  
       ...

Here's an example exhibited by Nusc, which I have tweaked slightly to make it look more like your mathematica example.

### Reference: http://www.mapleprimes.com/questions/36580-Bifurcation-Diagram

### xexpr is the logistic function to be iterated (we always start off at x=1/2, which will eventually attract).
### [ra,rb] is the range of the parameter.
### acc is the number of points sampled in [ra,rb]

Bifurcation := proc(initialpoint,xexpr,ra,rb,acc)

I am having trouble understanding how to put different tickmark spacing on the x and y axes.  Would someone be kind enough to show me the syntax of a plot command which plots the function y=x from x=-1..1 with tickmarks spaced 0.25 apart on the x-axis and 0.2 apart on the y-axis?

Thank you

The question is to use a program of modified newton raphson , incorporating the Romberg intergal procedure which i have already created, to create a new program which evaluates the integral f(alpha) = 1-10*int(tan(x)^alpha),from 0 to Pi/4.
The following needs to be incorporated in your program:

-Let ci be the approximation of alpha* on iterate i of your "modified newton raphson" method, then the program should run until:
             |ci-c(i-1)|<10^(-6)

I have defined a function as follows:

f:=alpha->1-RombInt(10*tan(x)^alpha,0,Pi/4,1,2);

and now i need to plot this function to see where the root lies, i have tried to do this using the following command:

plot(10*tan(x)^alpha);

and i have tried putting in boundaries but it doesnt seem to work. Can anyone help me please. 

Hi,

As a piece of work, I am required to define a function, and plot it in 3d.

The problem is regarding the 'gamma' function. We have to define it as:

myGamma:= x->int(exp(-t)*t^(x-1),t=0..infinity);

Then a later part of the questions asks you to define phi as the following function:

How may I program the following with Maple?

(1) Define the function:

H(p) := p * c1 + p^2 * c2 + p^3 * c3  + ......

(2) Now construct the following expression:

(1-p) [ v''' + 1 ]  = H(p) [ v''' - 25 * v' + 1]

Here, v is some funciton of x.

(3) Now assume: v = u0 + u1 * p + u2 * p^2 + ...

and substitute this into the expression defined at the point (2).

This worksheet is created with Maple 14, which works fine. Running it in Maple 15 I got a message "Connection to kernel is lost" or something the like (the message was in Chinese). I am under Mac OSX and the worksheet makes use of the grtensor ii package freely available on the net.  Any assist please?

C-metric-4d.mwC-metric-4d.mw

Hello, I'm Tam and I had a problem with Conditional of Jacobi Iteration.

So, the proplem is:

>restart;

>with(LinearAlgebra):

A:=Matrix([[5,1,-1],[0,4,2],[-1,2,-7]]);

And I want to check the condition

with the maple program.

Fisrt, I typed

b:=sum(abs('A'[1,j]),j=1..ColumnDimension(A)); and I got b:=7