Hi,

I used the dsolve command mentioned by Maple experts in the earlier blog to solve the following ode with bcs.

ode:= diff(y(x),x$2) + 1/x *diff(y(x),x) + 0.6*y(x) + y(x)^3 - y(x)^5  = 0

bcs:= D(y)(0) = 0,   y(d/2) = 0.8 * y(0)

The ode is a spherical coordinate system. The radius is defined by ‘x’ and the diameter by 'd'. 

The...

Hello!

I have a question about array/Array, as the picture illustrated, if applying a function to the entries of an array or Array, for example, map, then the results which maple yields lack of commas, as we can see, each entry in array/Array is so close and if the entries are very complex, then a confusion will be easily caused.

Is there any way to conquer it?

Hey peeps

I'm quite new to maple, and to this site, so I hope I do this properly.

I'm trying to make a procedure, that througout the procedure divides 2 numbers. But I need it to round down to the nearest integers, at all time, as I need an integer as output. I have searched quite a lot, and haven't found anything that helps me to do this.

The procedure is the following:

walla := proc (x) local y, z, w;

y := x;

Ques.mw In the WS attached,I want to specify that x1,x2,y1,y2 are real.I don't want the complex conjugate sign over the products in the hermitian conjugate matrix.

Also there are 2 examples of "Upper Shift" matrix- 3by3 & 4by4.I need a general index notation for this kind of matrix.(e.g. kronicker delta_i,j represents identity matrices).I tried with kronicker delta_i+1,j;but it doesn't give...

Maple 14.01 is now available. It provides updates in many areas, including:

• Mathematics: Updates to the VectorCalculus, DifferentialAlgebra, MathematicalFunctions, and Student packages, the convert command, and tools for solving differential equations
• Interface: Enhancements to context menus, tables, embedded components, document blocks, and the start-up code region
• Plotting:

Hej!

I've some problems in factorising terms in an equation.

As one can see in the maple file below, the terms "cosh((1/2)*a+(1/2)*y)", "cosh((1/2)*a-(1/2)*y)", "sinh((1/2)*a+(1/2)*y)" and "sinh((1/2)*a-(1/2)*y)" recur in the equation. So want to factorise those terms, but it doesn't work with factor().

Does anybody have a solution, how I can simplify that equation?

test.mw

For a eqation such as y(t)=2*sin(3*t)

I want to solve its fourier series under the different step,

and plot corresponding figure like the picture below

How should i do?

help me please.

Hello,

I have spent all morning trying to create a surface plot from a set of survey data that I have. The data I have is in the format xi, yi, zi. It is in Excel in columns. I read some previous questions on this topic, but my plot just comes out all jumbled.  Here is what I coded.

As soon as I figure out how to post images, i will post the plot...

Hello!!

I have a polynomial where each element is of the form  L * x1^a * x2^b * X3^c *x4^d. I want to replace X3^c by a sum of other elements , the number of this elements depends on c.

Is this possible? if so how can I do it on maple?

Hope you can help, thanks!

I built the shock model in MapleSim

then I got the response accleration

For this curve, I want to solve its fourier series in Maple.

Do it work?

thinks in advance.

this is the msim file

test.msim

Perhaps you have heard the terms "ordering difference" or "session dependent" applied to results of some Maple computation. It used to get heard more often back before Maple 12, when elements of sets in Maple were ordered according to address.

What would typically happen would be something like this. A computation gets done by looping either over the variables (symbol indets) of an expression, or over named elements of some set. And the computation happens to be very sensitive to which variables get chosen and used first. But the order in which they get chosen is according to the rule by which the elements of Maple sets are ordered. So if the address in memory of some of the variable names just happens to be different, then they get selected in an different order, and an entirely different computation path might get taken.

Now, merely using a name at some apparently unrelated earlier point in the Maple session can cause that name to get a different memory address than if the name is first used later on. The upshot was that apparently innocuous and unrelated use of variables at some earlier stage in a session could actually completey change a later result of an algorithm which happened to be sensitive to the order of set member access.

Since the difference in the final result in such cases would depend on ostensibly harmless earlier activity in the session, the whole effect was labelled "session dependence". And it caused a lot of grief.

As a wild (ficticious) example, you merely enter the name `x` early in a worksheet, and somehow the later result of an integral has a totally different (valid, but more complicated or unsimplified) form.

But since Maple 12 sets' entries are "ordered" (internally, by Maple's kernel) and a reduction in the appearance of session dependence has come about. But I was reminded recently (on usenet) that ordering differences can also apply to adding and multiplying expressions. And by coincidence this also ties in to something else that I've been thinking about: roundoff error.

Here's the simple example that came up. When you enter x/y you form a Maple expression whose internal representation (DAG) stores it as the product of x and of y^(-1). But which of those two multiplicands is stored first? It matters! When x and y get evaluated at floating-point values then it matters which is evaluated first. If y^(-1) gets evaluated first then roundoff error may occur. But if x is evaluated first then the float division might happen to occur without roundoff error.

Let's look closely at the example that came up.

> restart:

> ee := x/y;
                               x
                               -
                               y

> eval(ee, [x=2.1, y=2.4]);
                          0.8750000000

> restart:

> invee:=y/x:

> ee:=1/invee;
                               x
                               -
                               y

> eval(ee, [x=2.1, y=2.4]);
                          0.8750000001

In both sessions above, ee equals x/y. But only the second incurs roundoff error. The second result differs from the first, in the last decimal digit. It turns out that the order in which x and y^(-1) get stored in the product DAG is different in these two examples. If y^(-1) gets evaluated at y=2.4 before x gets evaluated, then Maple computes the intermediary term 1/2.4 = 0.4166666667 which of course has roundoff error. But if x gets evaluated first them Maple sees 2.1/2.4 as the division to compute, and that can get done with no roundoff due to fortuitous cancelation of factors of 3, ie. it is the same as 0.7/0.8=0.875.

You can see the internal storage distinction, as follows. I'll repeat the above pair of examples, but this time also dismantling `ee` to show its structure. Notice the order of `x` and y^(-1) insides the printed DAG.

> restart:
> ee := x/y;
                               x
                               -
                               y
> dismantle(ee);

PROD(5)
   NAME(4): x
   INTPOS(2): 1
   NAME(4): y
   INTNEG(2): -1

> eval(ee, [x=2.1,y=2.4]);
                          0.8750000000

> restart:
> invee:=y/x:
> ee:=1/invee;
                               x
                               -
                               y

> dismantle(ee);

PROD(5)
   NAME(4): y
   INTNEG(2): -1
   NAME(4): x
   INTPOS(2): 1

> eval(ee, [x=2.1,y=2.4]);
                          0.8750000001

Ok, so perhaps you are thinking: who cares, the difference is tiny. Well, small differences are easily isolated and magnified if they occur inside compound expressions. And so the ensuing results might be entirely different.

Take just a slightly more involved example,

> restart:
> ee := 1/(sqrt(0.8750001 - x/y)) - 3162;

                           1                 
                  -------------------- - 3162
                                 (1/2)       
                  /            x\            
                  |0.8750001 - -|            
                  \            y/            

> eval(ee, [x=2.1,y=2.4]);

                            0.277660

> restart:
> invee:=y/x:
> ee := 1/(sqrt(0.8750001 - 1/invee)) - 3162;

                           1                 
                  -------------------- - 3162
                                 (1/2)       
                  /            x\            
                  |0.8750001 - -|            
                  \            y/            

> eval(ee, [x=2.1,y=2.4]);
                            1.859986

Look at what caused the difference. All I did was assign y/x to invee and then re-use that to form `ee`. The expression `ee` looks the same on the surface. But inside it, the term x/y is stored with a different order of its multiplicands. And the final results are really different.

I suspect that an improvement might be had if evaluation at a point (2-argument eval) were instead done by substituting floats into numerator terms (positive integer exponent) before denominator terms (negative integer exponent). But investigation would have to precede that.

Finally, you might wonder what can be done with the errant version of `ee` from the last example. Simplification may or may not help,

> eval(simplify(ee), [x=2.1,y=2.4]);

                          0.2776606276

> eval(fnormal(ee), [x=2.1,y=2.4]);
                            0.277660

Are either of those last two results correct to ten places, you might ask? Well, no.

> Digits:=5000:

> eval(ee, [x=2.1,y=2.4]): evalf[10](%);
                          0.2776601684

acer

Have anyone toldl me how to find the family of curve at Maple 14

ax^2 + B*x*y+cy^2= D

with A=1, B=1, C =1 and D =2010

Can someone send me a clear step by step instruction to install Maple 13 for 32 bit Linux?

My OS is Ubuntu 10.04 LTS netbook edition.

I don't want to install it in home folder,instead I want to create a folder(say,MAPLE) in a partition with a lot of space and then install Maple13 inside this folder.

I have a worksheet that contains various pieces of a robot's dynamic and kinematic equations.  I want to perform different design tasks, taking the base model and working off of it.  I'd like to have only 1 model of the robot and not duplicate it across all of my different design worksheets because things change in the robot model as I switch up tooling, find better means of performing a system identification, etc.

Is there any functionality that would allow...