Hello people in mapleprimes,

I hope that you will give me an answer about the calculation surrounding D.
D((x+y)^2) can be calculated with maple, but D((x+y)^a), where a is a constant, cannot at least directly.
And, I have the code, which Mr. Carl Love kindly gave me at
http://www.mapleprimes.com/questions/218045-Logarithmic-Derivative
, which is

h:= F-> expand(evalindets(D(F)/F, specfunc(D), d-> op(d)*'h'(op(d))));

This works for (x+y)^2, but does not work for (x+y)^a.
And, I think this reason is that the behavior of D to (x+y)^2 is not the same
 as that to (x+y)^a.
And, I want to modify this code.
But, I can't .
So, I am asking this question. How should I modify the above code so that it works for (x+y)^a.
I'm so sorry for a lot of impoliteness.
I will be very glad if you give me an answer.

Take care.

taro

Hello people in mapleprimes,
I want to collect
aa:=sigma*omega*alpha = mu(H, H)*s*omega*gamma*L+omega^(-sigma+1)*tau*mu(F, F)*s*L-omega^(-sigma+1)*tau*mu(F, F)*s*L*gamma;
with omega^(-sigma+1);

I know that the above expression could be collected with tau, which
produces the same expression as the one collected with omega^(-sigma).
And, in collecting with omega^(-sigma+1), I know that freeze-thaw could be used.

But, I want to know another way to do the same thing, using frontend or 
subsindets. Is this possible?

Thanks in advance.

taro
 

Hello people in mapleprimes,
I have a question.
I appended two pictures where from the same code, two different orders of
expression appear.
How can I do for this so as not to get error messages?
The cause of this is simplify(%,symbolic) brings different order of term a__0^(-k)*F__D ahead of a parenthesis in a jpg.file and F__D*a__0^(-k) after
that parenthesis in another jpg.fine both in the line above that of  "dairihensu1."

In this case, What I can do?
Please help me.
Best wishes.

taro

my_code.mw

Original code is

e7_4:=F__D*(Omega+1)*beta/(beta-1) = F__I*a__D^(-k)*a__0^k+T^((sigma-k-1)/(-1+sigma))*F__D*phi^(k/(-1+sigma))+F__D;

a1:=beta=k/(sigma-1);
subs_free:=
  proc(a,b,c)
    local b1;
    b1:=isolate(b,c);
    subs(b1,a);
  end proc;
isolate(e7_4,a__D^(-k));simplify(%,symbolic);dairihensu1:=subs_free(%,a1,sigma);e7_5:=applyop(simplify,[2,4,1,3,2],dairihensu1);

A case without error.

A case with a error.

Hello people in mapleprimes,

I want to simplify the next expression which has 1/k as its exponent,

especially, I want to collect for T.　I hope you will teach me how to do it.

(F__X*(Omega+1)/(F__I*(beta-1)*T))^(1/k)*(T/phi)^(beta/k)

If I do as

simplify(%)assuming(symbolic);

the output is

F__X^(1/k)*(Omega+1)^(1/k)*F__I^(-1/k)*T^((beta-1)/k)*(1/(beta-1))^(1/k)*phi^(-beta/k)

But, as all variables has 1/k as its exponent, I want to collect it to (...)^(1/k).

Is this possible?

taro

Hello people in mapleprimes,

I put a function h into my original package named z5.

And, I used it as

> h(x+y);

then, the output is

x*`z5:-h`(x)/(x+y)+y*`z5:-h`(y)/(x+y)

I want to have

x*h(x)/(x+y)+y*h(y)/(x+y),

not z5:-h for h.

How should I do for "h" to be kept from changing  to "z5:-h"?

Best wishes.

taro

part other than z5 needs designating a directory for library.**********

>new_lib_dir3 :=Please put an appropriate directory

>z5:=module() option package;
>export numer_expand,h;
>numer_expand := a->expand(numer(a))/denom(a); h:= F-> expand
>(evalindets(D(F)/F, specfunc(D), d-> op(d)*'h'(op(d)))); end module;

>LibraryTools:-Save('z5',new_lib_dir3);
>restart;

>libname:=new_lib_dir3,libname;

>with(z5);

>h(x+y);

               x*`z5:-h`(x)/(x+y)+y*`z5:-h`(y)/(x+y)