Greetings,

Maple 15 allows the following syntax

omega := sqrt(w0^2*(1+((z-zf)/z0)^2))

But it does not allow this one 

omega^2 := sqrt(w0^2*(1+((z-zf)/z0)^2))

Why is this so for functions and variables? Is there any way around this, I am really bugged by this issue. The Maple Math told me that the syntax was invalid so I had type the commands directly.

Hi,

I've encountered a problem with Maple.

Is it possible to create a procedure which both displays some plots and also returns a value?

When I write a simple example:

It plots a circle as expected.

But when I add return:

It returns value 2, but no plot is displayed.

Is it possible to do both, display and return a value?

I used to write two copies of procedure, one for plotting and one for returning value.

But thats not "pretty" solution.

Thanks for help.

I have a formula, which is called the cube law. It requires input between 0 <= x <= 1. When creating  a procedure I run into all kinds of problems and get nothing but error codes. Trying to use assume works, it puts a tilde mark on the variable x, but doesn't restrict the input. Here is the formula the way it is writtin in a book on calculus:

f(x) := x^3/x^3+(1-x)^3 (0 <= x <= 1).

Trying to restrict the input to be restricted as required by the relational statement has been proven fruitless.

Any help would be greatly appreciated, thanks,

Rolf.

How do I setup/plot this probability

As a thought example, time waiting in a line at cash register.  Let's say 40% of the time you are waiting for ~1 minute while 35% of the time you are served immediately and the rest of the time trailing off waiting up to 15 minutes.

Is this is like a skewed normal distribution?  and how could I set it up?  And then how can I randomly generate some values from this distribution?

I have a characteristic equation. some times It has polar roots . sometimes It has real roots and sometimes both of them.

I want to extract real roots and extract polar roots if they are.

for instance:

q:=m3*r^3+m2*r^2+m1*r+m0:

rot:=solve(q=0,r);

I want to know how can I use if in this part ?

Hallo, I have a list with derrivatives. I'd like to choose elements with highest derrivative.

Is there any procedure in Maple to determine order of derrivative?

Hi there,

I used the implicitplot command to plot an implicit equation as follows:

with(plots):

implicitplot(x^3+y^3-3xy=0,x=-3..3,y=-3..3, scaling=constrained)

but maple13 takes error and specifies {xy} as an extra unknown.

Your help will be appreciated.

Regards

Hallo. I'd like to use two commands in loop like:

for i from 1 by 1 to k  do r:=Search(polyQ[i],[Q]) and T[r]:=T[r]-polyeq[i] end do

But Maple: "Error, unterminated loop". Hmm...

Can maple solve maximization problem like

q := proc (a, b, c) options operator, arrow; .2*b+.1*c end proc;
print(`output redirected...`); # input placeholder
(a, b, c) -> 0.2 b + 0.1 c
w := proc (a, b, c) options operator, arrow; .7*a+.1*c end proc;
print(`output redirected...`); # input placeholder
(a, b, c) -> 0.7 ab + 0.1 c
e := proc (a, b, c) options operator, arrow; .7*a+.2*b*c end proc;
print(`output redirected...`); # input placeholder
(a, b, c) -> 0.7 a + 0.2 b c

with(Optimization)

Maxmize(int(min(100+(.7*a+q)*(1/2), a), q)+int(min(100+(.2*b+w)*(1/2), b), w)+int(min(100+(.1*c+e)*(1/2), c), e)-a-b-c-ab-ac-bc)

Error, (in Optimization:-NLPSolve) cannot convert procedures to piecewise

Thanks alot if you can help me.Urgent! Really appreciate.

Hi,

Anyone knows how to define a ring of polynomials with variables x_1,...,x_n such that x_i*x_j = -x_j*x_i if i \neq j?

I tried using the Physics package and AntiCommutator but the problem is that in that case the variables anticommute with themselves so I have x_i^2=0.

If there is no direct way to this, I guess I could define a procedure that would look at the monomials of a polynomial and order them in lexicographic order and each time it switches two variables with different index it would multiply the monomial by -1.

I also don't know how to do that since I don't know how to look at a specific term in a monomial. For example, if my monomial is x_ix_jx_k, is there a way to find the first two variables and then switch them if LexOrder(variable1,variable2)= LexOrder(x_i,x_j)= false, i.e. if i>j?

Thanks,

M-A

Hi there,

I used the following plot command in maple13.

with(plots):

implicitplot(abs(x)+abs(y)=2, x=-2..2, y=-2..2, scaling=constrained)

But the output(square) is not complete at vertices.In fact at the vertex the graph is like a trapezoid.

Any help will be appreciated.

Yegan

hi

i have solved my equation as folllow :

pde:= diff(T(x, y), x)-1.555*10^(-7)*(diff(T(x, y), y, y))/ ...........

sol := pdsolve(pde, {T(0, y) = 0, (D[2](T))(x, 0) = 1325.754092, (D[2](T))(x, 0.25e-4) = 1970434.783}, numeric)

I wana know that maple has used which of numeric method to solve my equation ?

1.ForwardTime1Space[forward/backward]

2.CenteredTime1Space[forward/backward]

3.BackwardTime1Space[forward/backward]

4.ForwardTimeCenteredSpace or Euler

5.CenteredTimeCenteredSpace or CrankNicholson

6.BackwardTimeCenteredSpace or BackwardEuler

7.Box

8.LaxFriedrichs

or ... ?

Tahnks.

For different reasons I need to ocasionally export a number of Maple worksheet in a folder to pdf files. Is there a way to automate this? I would want that the worksheet is opened, output removed, then executed and eventually exported to pdf. It can take quite a while to do this manually for about 50 worksheets.

Dear people in Maple primes,

Please tell me why algsubs would not work well in the following case.

algsubs(-eta*M[n]/(gamma*(eta*M[n]+gamma))-1/(eta*M[n]+gamma) = factor(-eta*M[n]/(gamma*(eta*M[n]+gamma))-1/(eta*M[n]+gamma)), (-eta*M[n]/(gamma*(eta*M[n]+gamma))-1/(eta*M[n]+gamma))*p[i]+(int(p[i], i = 0 .. M[n]))*eta/(gamma*(eta*M[n]+gamma))+alpha/(eta*M[n]+gamma));

Following error message appeared when the above code is executed:

Error, (in algsubs) cannot compute degree of pattern in eta

Best wishes.

taro

HI, dear all. When I tried to use the plot option 'adaptive' to make my plot more smooth and realistic, I encountered the following erro. I cannot understand why. From the help guide, I learn that adaptive can be assigned n or true or false, but errors appeared.  Thanks for your help.

> implicitplot(-x^3+3*x+a = 0, a = -3 .. 3, x = -4.0 .. 4.0, view = [-3 .. 3, -4 .. 4], adaptive = 2, resolution = 1000, numpoints = 2000);
Error, (in plot/options2d) unexpected option: adaptive = 2
> help("adaptive");