taro

435 Reputation

12 Badges

5 years, 242 days
Maple is to me difficult. The first version I bought was Maple9, and it was more than 15 years ago. But, I couldn't use it, feeling it too difficult. But, three years ago, I thought Maple might be helpful to my study, and since then, I have continued to learn Maple. As I got able to read the Maple help, I think that I could get to use maple better now than before. But, I feel that I am a beginner yet.

MaplePrimes Activity


These are replies submitted by taro

@tomleslie 

N := [7.5000000000, 2.5000000000, 4.0000000000]

is different from

M:={7.5, 2.5, 4}

printf("%g %g %g\n", op(N) );
outputs 7.5 2.5 4, but this cannot be used for the continuing calculation.
How can I use the outputs for the later calculation?

I'm not the questioner, but I will be very glad if you will give me an answer to this question.

 

@Kitonum 

Thank you for your detailed explanation.
I could understand the limit of inequality in maple and how to deal with inequalities owing to you.

 taro

@ThU 
Surely, in this case, it is easier to rewrite the expression manually.
But, I think I could improve my skill much in using Maple, through asking that question here and reading the code by Kitonum. I owe a lot to people who gave me answers to my questions.



 

@Kitonum 

Thank you. I owe you a lot.
I found from your answer that the basic idea, when a variable overlaps between the part already collected and the part not yet collected, is to replace the part already collected with another variable and to collect the part not yet collected.

And, thank you for showing me the use of freeze and thaw, and how to chooze the already
collected part. Your code is very helpful.
Thanks a lot.

@vv 

I'm sorry. I couldn't understand what you wrote. What is some subexpression which you wrote? and what can I do if I had defined that subexpression. You posted your
comment while I was modifying what I wrote, that is, in the part of mentioning that algsubs could not be used for appropriately getting the answer. And, if my
remove(has, c1, delta) means your writing that I have to obtain 'that subexpression,'
what can I do after doing it?

@vv 

Thank you for your comment. But, this is not for typesetting, though looking so.
My question is how to do additional collect in the case that the variable with which I want to collect the rest of expression has already been used for another collected part.


More writing to explain my problem, as op([5,2..3],a1) is not -epsilon+z, I cannot use applyop for collecting -epsilon*z + z to 
z*(1-epsilon).
   On the other hand, I can identify -epsilon*z+z with

But, this identification does not look helpful, as algsubs(d1=collect(d1,z), c1) does not play the role which I expected.

@nidojan 

I'm sorry. I think I might have lacked of enough explanation about my question.
What I want to do is to collect the part of -epsilon*z+z in a1 to z*(1-epsilon) in b1.

 

 

@Markiyan Hirnyk 

I think that your advice is to divide both of numerator and denominator with a specific expression. In that meaning, through doing so, I could simplify the expression which I wanted to simplify. Thank you.

 

@Kitonum 
Thank you for your teaching.
From your answer, I found substitution plays well.

Though a similar one, dividing first works well.

e1:=n = (theta-1)*(z-1)*(-k*psi*L*(-(epsilon-1)*psi/epsilon)^epsilon*k^epsilon+K*(-(epsilon-1)*psi*k/epsilon)^epsilon)/(gamma*((theta-1)*(z-1)*(-(epsilon-1)*psi*k/epsilon)^epsilon-(-(epsilon-1)*psi/epsilon)^epsilon*k^epsilon*theta*z*psi));
b:=(-(epsilon-1)*psi*k/epsilon)^epsilon;
c:=simplify(numer(rhs(e1))/b)/simplify(denom(rhs(e1))/b) assuming (k>0);

The problem in this case to me migh have been the use of assuming(k>0).

Thank you.

 

 

 @sand15 

Doesn't the way to ascertain whether all eigenvalues are positive become an answer?

Vector[row]([r,s]).Matrix([[1,-1/2],[-1/2,1]]).Vector[column]([r,s]);

expand(%)

a:=Matrix([[1,-1/2],[-1/2,1]]);

LinearAlgebra[Eigenvalues](a);

 

@Kitonum 

In my environment: maple2016, osx EL Capitan, the result of assignment is different from yours. A and B are reverse in my results. Could I ask you to show the way to correct this?

 

@Zeineb 

I think that what you want to obtain is the limit of g(t) when t goes to infinity.
Setting t=50, the following graph is drawn.

plot(abs(x^(2^(-50))-1),x=-0.1..2);

where it increases between -0.1 and 0, decreases between 0 and 1, and increases above 1.
I think you have to do is to show the part of below where.

By the way, how can I make that graph smooth?

 

@Kitonum 

Would you please tell me why "," that is comma, at the last of "evalindets(Expr, 'function', t->op(0, t)), " in yor code,
cannot be replaced with ";" that is semicolon?

 

 

@Scot Gould 

Sorry. As the description of two commands are the same other than the outcome of
examples there, I thought they were the same. But, they are different.
See below.

plot(-(x+2)*(x+1)*(x-1),x=-2.1..2.1);
evalf(ExtremePoints(-(x+2)*(x+1)*(x-1),x=-3..3));
evalf(CriticalPoints(-(x+2)*(x+1)*(x-1),x=-3..3));

The output of ExtremePoints contains the end points x=-3..3, and
that of CriticalPoints does not.
So, in the case of the original question, CriticalPoints, which is Markiyan's answer, is more appropriate.
 

3 4 5 6 7 8 9 Last Page 5 of 22