taro

445 Reputation

12 Badges

6 years, 100 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 
Thank you for your answer.
Especially, your use of freeze and thaw was very helpful to me.

I owe you a lot.

On the other hand, later, I found the following codes works;

Y := -L*epsilon^epsilon*(1-epsilon)^(1-epsilon)*(-omega^(-1+epsilon)*k*theta+(1-theta)*omega^epsilon)/(1-epsilon-theta);
gg:=proc(x)
expand(x/omega^epsilon)*omega^epsilon;
end proc;
applyop(gg,5,Y);

taro

@Kitonum

Thank you.

I could obtain a lot from your direct substitution and the use of empty symbol" `` ".

taro

 

@Rouben Rostamian
Thank you. I could resolve my question with you code.

taro

@tomleslie

munus 1 of the top of the expression is multiplied to the denominator of (theta -1) so that it changes to 1-theta.

So, two are equivalent. The reason I prefer the second is that the expressions in my note following the expression I wrote here have the same (theta -1) so that the terms have minus sign at the top, so I though I would change the first place
where it appeared, so that the following expressions are easy to deal with. I am supposing theta is a constant between 0 and 1. 

 

@Will 

I found a typo. So, I want to report it to maplesoft.
Though I searched for such places with type and found this thread, I am not certain this is one.

In the term of evalindets in help page of maple 2016, there is the next example contained.

evalindents['nocache'](expr,symbol,proc(s) global count; count:=count+1; cat(s,count) end);

evalindents is incorrectly spelled. evalindets is the right one.

And, ex of

evalindets(expr,specfunc(f),x->op(1,ex));

, which is written below the last one,

, I think, should have been written, intended to write expr.

@tomleslie 

The first way is

x:={b[1]=alpha*a[1]^6, b[2]=alpha*a[2]^6, c[1]=beta*a[1]^4+gamma*d[1], c[2]=beta*a[2]^4+gamma*d[2]};

y:=a[1]^2+b[1]^2+c[1]^2+d[1]^2+a[2]^2+b[2]^2+c[2]^2+d[2]^2=0;

z:=subs(x,y);

solve(z,a[1]);

which brings a solution, though including RootOf.

 

The second way is

k:= solve({op(x),y},a[1]);

But, this does not bring any answer, which I can't understand why.

 

@tomleslie 

Hello

Two ways to solve

P:=(a,b,c,d,e,f)->a^2+b^2+c^2+d^2+e^2+f^2;

that is
b:=x->x^2;
c:=x->x^2;
d:=x->x^2;
e:=x->x^2;
P(a,b(a),c(a),d(f),e(f),f)=0;
solve(%,a);

 

and

solve(P(a,b,c,d,e,f),[a,b,c,d]);

bring different solutions.

So, I couldn't understand your answer.

 

@vv 

evalf(151/11);
convert(%,rational,3);

which brings 41/3 of OP.

@tomleslie 

Thank you for your reply.

with the change of preferences of maple in Precision, Round screen display, to 3
decimal places, I can get the view of an only 3 digit below a decimal point like this.

The meaning of my question is how to do the same without changing the preference
so that without exerting this modification on other calculation on the same worksheet, that is to limit this effect on only a part of various calculation.

I thought of this question seeing the difference of the expression of the original questioner's expression of the numbers, and your result of exact expression of the numbers.

 

@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.

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