Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

I have written a program which plays the game of Multicube, a commercially made boardgame. It works OK except that when the game is replayed, the previous game output stays on the screen.  This is rather distracting and I'd like the screen cleared when a new game is run.   I'd like to know how the screen previous ouput can be cleared.   ...or what is the likely cause.

   At the start of my program I have:

restart;
interface(echo=0, verboseproc=0, warnlevel=0, prettyprint=1):

#I thought restart; would automatically have cleared any previous output.

The program reads necessary input (eg no of players, ..)  using:

x:=readline(terminal):

 

Thanks,  David

 

 

Hello,

I am trying to solve the boundary value problem (1-x^2)*y'' - 2*x*y' +12*y = 0 with y(-1) = -1 and y(1) = 1.  I have not used Maple much, but from some web surfing, it seems like the following inputs should work:

de := (1-x^2)*(diff(y(x), `$`(x, 2)))-2*x*(diff(y(x), x))+12*y(x) = 0

Y := dsolve(de, y(-1) = -1, y(1) = 1)

However, when I input these lines, I get the error: Error(in dsolve), found wrong extra argument(s): y(-1) = -1, y(1) = 1

Does this mean that Maple can't solve this problem?  Is my syntax wrong?  I would appreciate any help.

 

Thanks!

Tim

 

I have the following situation:

HB:=Bend(L,a,n);

Bend is a proc, that returns a Record with info based on its parameters. I would like to get access to the name I assign to (i.e. HB) in the proc. Any chance?

Mac Dude

PS: I can of course kludge things by adding an argument to the proc Bend. But I'd like to avoid that.

Hi,

 

I'm trying to solve the following differential equation numerically with dsolve:

but dsolve gives me this error:

> res := dsolve(DGL, numeric, parameters = [y0, A, B, C, E]);
Error, (in DEtools/convertsys) unable to convert to an explicit first-order system

I think the problem is that I use the wrong solver. Does Maple provide a solver which is capable of solving this kind of equations (nonlinear ODE)?

 

Thanks in advance!

 

Hi,

 

For my thesis I would like to illustrate the inclination of the solar system objects in a plane. Imagine the Solar system as a circular plane which is the average of the motions of all objects orbiting the Sun. Each individual planet/object is inclined towards this averaged plane - some more, some less. E.g. this image. Behind Neptune there lies the so-called Kuiper belt with many thousands of dwarf planets (Pluto is one of them and there are over 1000 objects known out there already). Now imagine that the big heavy Jupiter and other big planets perturb those small objects out there: so they are also inclined towards the average plane, see this image. Their inclination depends on their radial distance to the Sun (measured usually in Astronomical Units...1 unit is the distance Sun-Earth...we are talking about 40 to 50 units here). For one, I want to visualise this: imagine a circular plane and each orbit out there has another angle to the average plane. This is the first. But now: imagine the circle with it's 360°. Each object reaches it's highest point on its orbit around the Sun on another angle on this 360° circle. Neptune e.b. at 170°, Pluto at 250° etc. You get the picture. So not only are the objects in the Kuiper belt inclined differently, but their maximum orbit positions are also scattered across an imaginary 360° circle. I want to show this with a 3d-plane like this image but not with two peaks for one orbit. Is there a way to do/plot/visualise this in Maple? I am just interested in a visualisation of the principle without any empirical data behind this.

Thanks.

hello , to solve a differential euqations , maple gave me this error. please help me.

 

TG.mw

Hi all,

 

I'm trying to initialize a 4-vector and promote it as a tensor in Maple 16 using the Physics package. My attempt so far has been unsuccessful :) 

This 4-vector needs to be a function of the coordinates X. What I'm trying to do is to allocate terms by terms my tensor as indicated in my attempt below. Everything looks fine (the query is telling me F2 is a tensor) but when I'm trying to compute contravariant quantities, I can't obtain an answer for a given coordinate (here F2[~n](X) is not equal du v_1(X))... Finally, the last line, computing the sum over all the indexes is not giving me an answer... 

 

restart; with(Physics):
Setup(mathematicalnotation=true):
Setup(coordinatesystems = cartesian):
Define();
g_[]:

F2[1]:=(x,y,z,t)->v_1(X);F2[2]:=(x,y,z,t)->v_2(X);F2[3]:=(x,y,z,t)->v_3(X);F2[4]:=(x,y,z,t)->1;

Define(F2[mu](X));

Define(F2,query);

F2[1](X);F2[2](X);

F2[~1](X);F2[~2](X);

F2[~nu].F2[nu]

 

Surely, something is wrong in my way of defining a spacetime tensor. What would be the solution to obtain a decent contravariant F2?

 

Another more or less related question would be to know how to compute a taylor expansion of a function itself i.e. 1/(1+f(X))~ 1-f(X) since f<<1? Would there be a built-in fonction to do such a thing?

 

Many thanks!

 

Vincent

Hello people in mapleprimes,

I hope you will give an answer to my following question.

This is a question of nonunit(algebraic) in patmatch.

In help page of maple this examples are on.

 

Case i)

>patmatch(a, A::(nonunit(algebraic))+B::(nonunit(algebraic)), 'la');

false

 

Case ii)

>patmatch(a,A::(nonunit(algebraic))+B::algebraic,'la'); la

ture

[A=a,B=0]

Is unit in the case of sum zero, so that in Case i) false is shown, though in Case ii) B=0 is returned?

 

Best wishes.

 

taro

 

 

 

Dear people in mapleprimes,

 

>define(INT, linear, conditional(INT(a::algebraic, X::name) = a*X, _type(a, freeof(X))), INT(X::name, X::name) = (1/2)*X^2);

>INT(2*x+4, x);

x^2 + 4*x

>INT(z+x, z);

(1/2)* z ^2+ x*z

is written in the help page of maple.

This code is to make a function INT.

What it means is that if INT(a::algebraic, x::name) then this is transformed into a*X, where a is independent from x.

And, INT(X::name,X::name) is transformed into (1/2)*X^2.

 

The first of what I want to ask you is the way of use of "_type", which is not in the help page, and as for names with "_" ahead, 

it is written that "_" means internal command for Maple, so not to use it. Is it right to use a name with "_" ahead of it here?

The second of what I want to know is the way to use "conditional."

Please teach me the meaning of this.

Thanks in advance.

 

I hope you will teach me the above questions.

 

Best wishes.

 

taro

 

 

seq(i,i=1..5);  gives output 1,2,3,4,5

seq(i,i=5..1);  seems to get ignored in Maple 7.  It does not return an error message.

I was hoping to step back from a large number to a smaller one in stepsof minus one.  The for loop structure allows negative "increments" - so I experimented by putting the decrement value of -1 between the two dots.  I was expecting it to givean error message, so was surprisd when it gavesome output!

eg  seq(i, i=5.-1.2);  gives output of 38, -1

I'm curious as to what is happening here?  Also, has the seq command been upgraded in later versions of Maple?  eg is it possible to do negative decrements?

Hi there,

I would like to know whether there is some way to print non-integer axis numbers for a given plot.

Let's say I have an array for an independent variable such as [0.5 1 1.5 2]. I would like Maple to use these values in the corresponding axis instead of printing the array indices [1, 2, 3, 4].

I would like to know how this is done for a surface plot.

 

Attached is the worksheet where "b" is such a parameter/independent variable: MapleSimulation_test11.mw

 

Thank you,

jon

I am just starting with Maple 18 on Windows 8.1. I did not like the default color palette for 2D graphs and tried another one. So I edited my maple.ini file, which was in 

Hi,

 

  Excuse me, I tried to substract in a series expansion. It works for Taylor-type expansion, but not Laurent-type ones

*********

ffff:=1/(t-3);
ggg:=taylor(ffff,t=infinity,3);
coeff(ggg,t,-1);


ffff:=exp(t);
ggg:=taylor(ffff,t=0,3);
coeff(ggg,t,1);

*********

Maple gives me


Error, unable to compute coeff


1

 

. My question is, why "coeff" does not work in the first case? how to make it work? Though I can copy the expansion part, assign it to another variable, and "coeff" will work

 

Excuse me,

 

  I have the following input

***

fff:=x1^k;
int(fff,x1=0..x2);

***

I would like to obtain

1/(k+1)* (x2^(k+1)), k<>-1

 

but maple gives me just the same integral. Is there any way to get the expression?

 

Thank you!

 

 

 

Hello,

Im solving 4 ODE equations with BC. im trying to shoot the initial value but im having this error:

""Error, (in isolate) cannot isolate for a function when it appears with different arguments""

anyone could help me???

shooting92.mw

``

restart

Shootlib := "E:\\shooting/":

libname := Shootlib, libname:

with(Shoot):

with(plots):

n := 2:

FNS := {F(eta), H(eta), f(eta), g(eta), u(eta), v(eta)}:

ODE := {g(eta)*(diff(g(eta), eta))+B*(f(eta)+g(eta)) = 0, g(eta)*(diff(F(eta), eta))+F(eta)^2+B*(F(eta)-u(eta)) = 0, g(eta)*(diff(H(eta), eta))+H(eta)*(diff(g(eta), eta))+F(eta)*H(eta) = 0, diff(v(eta), eta)+f(eta)*v(eta)-u(eta)^2+B*H(eta)*(F(eta)-u(eta))-M*u(eta) = 0, diff(f(eta), eta) = u(eta), diff(u(eta), eta) = v(eta)};

{g(eta)*(diff(H(eta), eta))+H(eta)*(diff(g(eta), eta))+F(eta)*H(eta) = 0, g(eta)*(diff(g(eta), eta))+0.2e-1*f(eta)+0.2e-1*g(eta) = 0, g(eta)*(diff(F(eta), eta))+F(eta)^2+0.2e-1*F(eta)-0.2e-1*u(eta) = 0, diff(v(eta), eta)+f(eta)*v(eta)-u(eta)^2+0.2e-1*H(eta)*(F(eta)-u(eta))-3*u(eta) = 0, diff(f(eta), eta) = u(eta), diff(u(eta), eta) = v(eta)}

(1)

IC := {F(0) = gamma, H(0) = Q, f(0) = 0, g(0) = z, u(0) = 1, v(0) = alpha};

{F(0) = gamma, H(0) = Q, f(0) = 0, g(0) = z, u(0) = 1, v(0) = alpha}

(2)

BC := {F(L) = 0, H(L) = n, g(L) = -f(L), u(L) = 0};

{F(6) = 0, H(6) = 2, g(6) = -f(6), u(6) = 0}

(3)

infolevel[shoot] := 1:

S := shoot(ODE, IC, BC, FNS, [alpha = 0, gamma = 0, z = -.2, Q = 0])

Error, (in isolate) cannot isolate for a function when it appears with different arguments

 

``

``


Download shooting92.mw

First 1287 1288 1289 1290 1291 1292 1293 Last Page 1289 of 2224