Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Code was computed on two different PC:
1) Win 8.1, Maple 17.02 x64
2) Win 8.1, Maple 18.01 x64

In brief:

n:=3: 
Grid:-Seq('f(i)', i = 1 .. n):

works well in Maple 17.02 but gives an error in Maple 18.01. Why?

However,

(j -> Grid:-Seq('f(i)', i = 1 .. j))(n):

works well for both Maple versions.

according to examples from Maple help page
http://www.maplesoft.com/support/help/Maple/view.aspx?path=Grid%2fSeq
we have to pass exactly 'f(i)' (not just f(i)) into Grid:-Seq(...). Why?

Why doesn't Digits:=30: work properly for Grid:-Seq(...)?

For details see attached file. Output was produced in Maple 17.02. In Maple 18.01 the output will differ.


restart; n := 3; f1 := proc (i) options operator, arrow; evalf(sqrt(i)) end proc; f2 := proc (i) options operator, arrow; evalf(sqrt(i+1)) end proc; f3 := proc (i) options operator, arrow; evalf(i*Pi) end proc; g1 := proc (f, n) local za, zb, zc, str; str := sprintf("g1(%a,n)", f); za := seq(f(i), i = 1 .. n); zb := Grid:-Seq(f(i), i = 1 .. n); zc := (proc (j) options operator, arrow; Grid:-Seq(f(i), i = 1 .. j) end proc)(n); print("------------------------------"); print(cat("(a):   ", str), za); print(cat("(b):   ", str), zb); print(cat("(c):   ", str), zc); print("------------------------------") end proc; g2 := proc (f, n) local za, zb, zc, str; str := sprintf("g2(%a,n)", f); za := seq('f(i)', i = 1 .. n); zb := Grid:-Seq('f(i)', i = 1 .. n); zc := (proc (j) options operator, arrow; Grid:-Seq('f(i)', i = 1 .. j) end proc)(n); print("------------------------------"); print(cat("(a):   ", str), za); print(cat("(b):   ", str), zb); print(cat("(c):   ", str), zc); print("------------------------------") end proc

3

(1)

g1(f1, n);

"------------------------------"

 

"(a):   g1(f1,n)", 1., 1.414213562, 1.732050808

 

"(b):   g1(f1,n)", 1, 2^(1/2), 3^(1/2)

 

"(c):   g1(f1,n)", 1, 2^(1/2), 3^(1/2)

 

"------------------------------"

 

"------------------------------"

 

"(a):   g1(f2,n)", 1.414213562, 1.732050808, 2.

 

"(b):   g1(f2,n)", 1.414213562, 1.732050808, 2.000000000

 

"(c):   g1(f2,n)", 1.414213562, 1.732050808, 2.000000000

 

"------------------------------"

 

"------------------------------"

 

"(a):   g1(f3,n)", 3.141592654, 6.283185308, 9.424777962

 

"(b):   g1(f3,n)", 3.141592654, 6.283185308, 9.424777962

 

"(c):   g1(f3,n)", 3.141592654, 6.283185308, 9.424777962

 

"------------------------------"

 

"+++++++++++++++++++++++++++++++++++++++++++"

 

"------------------------------"

 

"(a):   g2(f1,n)", 1., 1.414213562, 1.732050808

 

"(b):   g2(f1,n)", 1., 1.414213562, 1.732050808

 

"(c):   g2(f1,n)", 1., 1.414213562, 1.732050808

 

"------------------------------"

 

"------------------------------"

 

"(a):   g2(f2,n)", 1.414213562, 1.732050808, 2.

 

"(b):   g2(f2,n)", 1.414213562, 1.732050808, 2.

 

"(c):   g2(f2,n)", 1.414213562, 1.732050808, 2.

 

"------------------------------"

 

"------------------------------"

 

"(a):   g2(f3,n)", 3.141592654, 6.283185308, 9.424777962

 

"(b):   g2(f3,n)", 3.141592654, 6.283185308, 9.424777962

 

"(c):   g2(f3,n)", 3.141592654, 6.283185308, 9.424777962

 

"------------------------------"

(2)

Digits := 30;

30

 

"------------------------------"

 

"(a):   g1(f1,n)", 1., 1.41421356237309504880168872421, 1.73205080756887729352744634151

 

"(b):   g1(f1,n)", 1, 2^(1/2), 3^(1/2)

 

"(c):   g1(f1,n)", 1, 2^(1/2), 3^(1/2)

 

"------------------------------"

 

"------------------------------"

 

"(a):   g1(f2,n)", 1.41421356237309504880168872421, 1.73205080756887729352744634151, 2.

 

"(b):   g1(f2,n)", 1.414213562, 1.732050808, 2.000000000

 

"(c):   g1(f2,n)", 1.414213562, 1.732050808, 2.000000000

 

"------------------------------"

 

"------------------------------"

 

"(a):   g1(f3,n)", 3.14159265358979323846264338328, 6.28318530717958647692528676656, 9.42477796076937971538793014984

 

"(b):   g1(f3,n)", 3.14159265358979323846264338328, 6.283185308, 9.424777962

 

"(c):   g1(f3,n)", 3.14159265358979323846264338328, 6.283185308, 9.424777962

 

"------------------------------"

 

"+++++++++++++++++++++++++++++++++++++++++++"

 

"------------------------------"

 

"(a):   g2(f1,n)", 1., 1.41421356237309504880168872421, 1.73205080756887729352744634151

 

"(b):   g2(f1,n)", 1., 1.414213562, 1.732050808

 

"(c):   g2(f1,n)", 1., 1.414213562, 1.732050808

 

"------------------------------"

 

"------------------------------"

 

"(a):   g2(f2,n)", 1.41421356237309504880168872421, 1.73205080756887729352744634151, 2.

 

"(b):   g2(f2,n)", 1.414213562, 1.732050808, 2.

 

"(c):   g2(f2,n)", 1.414213562, 1.732050808, 2.

 

"------------------------------"

 

"------------------------------"

 

"(a):   g2(f3,n)", 3.14159265358979323846264338328, 6.28318530717958647692528676656, 9.42477796076937971538793014984

 

"(b):   g2(f3,n)", 3.141592654, 6.283185308, 9.424777962

 

"(c):   g2(f3,n)", 3.141592654, 6.283185308, 9.424777962

 

"------------------------------"

(3)

``

NULL


Download Grid[Seq].mw

Hello guys...

please help me to solve these equation.I have some equation with some boundary condition.

Why can't I have good results for it?
q3:=-0.1080877976e-1:
q2:=0.04351279185:
q1:=-0.7037638761e-1:
q0:=0.05575044494:
c:=3:

eq := q3*m^6+q2*m^4+q1*m^2+q0;
rishe := solve(eq = 0, m);
a := 1/(2*c); rho := a/`#msub(mi("t"),mn("1"))`;
f0 := evalf(cosh(abs(rishe[1])*xi));
f1 := evalf(cosh(abs(Re(rishe[3]))*xi)*cos(abs(Im(rishe[3]))*xi));
f2 := evalf(sinh(abs(Re(rishe[3]))*xi)*sin(abs(Im(rishe[3]))*xi));
pin := C0*k0*f0+C1*f1+C2*f2;
ppin := diff(pin, xi); pis := C0*f0+(k1*C1+C2*k2)*f1+(C2*k1-k2*C1)*f2;
pin1 := evalf(subs(xi = a, pin));
pin2 := evalf(subs(xi = -a, pin));
pin3 := evalf(subs(xi = a, ppin));
pin4 := evalf(subs(xi = -a, ppin));
pis1 := evalf(subs(xi = a, pis));
pis2 := evalf(subs(xi = -a, pis));
hal := solve({pin1 = sigma01xx, pin2 = sigma01xx, pin3 = 0, pin4 = 0, pis1 = sigma01xy, pis2 = sigma01xy}, {C0, C1, C2, k0, k1, k2});
how can I have a good result with this method?
please help me to solve this multi equations multi variable.

I have a big expression with multiple terms. The expression is kind of :

expr = a*ln(u) + b*ln(u) + a*ln(v)+b*ln(v)

I want to factor ln(u) and ln(v) to get (a+b)*ln(u) + (a+b)*ln(v)

I tried factor(expr,ln(u)) but I get the error :
Error, (in factor) 2nd argument, ln(2*da+(4*da^2+4*d^2-4*da*Sigma[a]+Sigma[a]^2+Delta[b]^2)^(1/2)-Sigma[a]), is not a valid algebraic extension.

I tried without Delta[b] and Sigma[a] : deltab and sigmaa but it didn't work either.

Do you have an idea how to factor the many ln ?

Thanks,
Mathieu

Well, I'm having issues getting this expression to simplify.  I have a bunch of polynomial results I store in an array, and the simplify command isn't cancelling obvious terms which should cancel out.

June_18.mw

Curiously, it will simplify if I use the command by selecting it from the drop down menu when I select the expression (see lines (7) vs. (8))

i am trying to export an animated gif file but the export windows freezes and exports an empty file with zero bytes

the animation was created by: display(map(p, convert(A(() .. (), 2), list)), insequence = true)

i have tried to give maple some time (2-3 hours) and i even tried to restart the computer and lauch maple and try to export it but it didnt worked

how can i export it or how can i plot the using commands

 

I want to solve the equation sqrt(x) + sqrt(1 - x^2) = sqrt(2 - 3*x - 4*x^2) in RealDomain. I tried

RealDomain:-solve(sqrt(x) + sqrt(1 - x^2) = sqrt(2 - 3*x - 4*x^2),x);

And I got one solution. But, at here 

At here http://mathematica.stackexchange.com/questions/51316/how-can-i-get-the-exact-real-solution-of-this-equation 

they said the given equation has two real solutions. How must I understand?

Hi everyone

Title might be misleading but couldn't come up with a short version of my problem.

The idea is to show an animation of several plots of an array, which all have different values in the positive and in the negative range. For distinguishing purposes I want to show positive values in a different colour than the negative values.

For better understanding imagine a chess board with the rows and columns and there are places with values, black chess pieces as positive values and white chess pieces as negative values. Each array in the sequence from A1 to An shows a movement.

For the plots I use sparsematrixplot which unfortunately just shows all non-zero values. So I seperated the original array and made two arrays, one with the positive values and one with the negative values. I am able to animate the sequence of several plots and can play the animation of two sequences consecutively with following:

A:=seq(sparsematrixplot(convert(H||g,matrix),color=green,view=[0..10,0..10]),g=1..G):
B:=seq(sparsematrixplot(convert(F||g,matrix),color=red,view=[0..10,0..10]),g=1..G):
display([A,B], insequence=true);

Hence my questions
1. Is it possible to plot an array with different colours for different values and can I animate these plots?

2. If not, is there a way to display the two plot sequences on top of each other?

 

Thanks for your help in advance.

Hi all,

i have the following problem:

i had create this tiny procedure to compute the berstein base elements:

Berstein base

And so i would like print space, between B(i,n,t) and the formula corresponing, without quotes :D

In other word i would an output like this:

B(i,n,t),      ,(formula)

 

Thank you in advance :)

     

 

 

Hy all.

I want to solve this equation, with„dd” as numerical result. What do I do wrong? Thanks. Nico

restart;
TTot := 70;
TC := 17;
GM := .26;
QMax := 870;
V := 3600*GM*QMax*TTot;
eq := V = int(QMax*exp((-t+TC)/dd)*(1+(t-TC)/TC)^(TC/dd), t = 0 .. TTot);
fsolve(eq, dd);

 

Greetings everyone. I tried to get the roots of this polynomial by using allvalues command. But I got the roots in indeces instead of the value of the roots.

what should i do?

Hi all,

my problem concerns operation with complex number in Maple 18. The issue is the following:

i define this complex: c:=a+i*b

then i compute the square: sort(evalc(c^2))

and the output is: a^2+2*i*ab-b^2

So, how can i obtain an output like the following?  a^2-b^2b^2 +a*i*ab

In other word i want an output where the real part precede the complex part.

thank you in advance :)

Hi all. I am a new maple 18 user and Mapleprimes member.

1. I tried to plot y=x^2+2, but in the graph the value of y at x=0 is not 2(the bottom of the graph is near y=0). Same for the plot of just x^2+2(no y variable).

2. Why do I need to specify a y range for y=x^2+2? Why can't I just specify the x range only, since y is based on x?

Thanks in advance. :)

 

 

 

Maple has a powerful package for manipulation with logical expressions (Logic package). But it lacks for transfomation a boolean expression to algebraic normal form (ANF). In fact, this transfomation is very simple and implements by applying simple rules

to disjunctive normal form.

I suppose that this can be reached by introducing some custom operators and transfomations, if it is possible in Maple.

 

The problem is when initializing a Matrix with a list of strings. The worksheet excerpt below shows the normal behavior using a list of integers to initialize a square matrix: the successive list elements fill the matrix by rows.

Then trying the exact same thing with a list of strings instead of integers gives an error message!
This is not right. While it is an odd and likely rare problem, it would be better fixed.

x := [i $ i=1..25];
[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25]
Matrix(5,5,x);
Matrix(5,5,[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25]) 

y := map(convert,x,string);
["1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25"]
Matrix(5,5,y);
Error, (in Matrix) initializer defines more columns (25) than column dimension parameter specifies (5)

How to identify whether biased game is changed based on the latest toss from a sequence

can HMM identify this?

for example

if

fair fair fair biased biased fair fair biased

First 1335 1336 1337 1338 1339 1340 1341 Last Page 1337 of 2223