John2020

200 Reputation

5 Badges

2 years, 79 days

MaplePrimes Activity


These are replies submitted by John2020

@acer 

Hello,

A complicated example that the codes not work for it is as follows:

restart;

with(Physics):

r := x*(diff(theta(t), t))^2+y*(diff(varphi(t), t))^2+z*(diff(theta(t), t$2))+w*diff(varphi(t), t$2);
g := (4*(f+T))*(diff(theta(t), t))^2+u*(diff(varphi(t), t))^2+(f+9)*(diff(theta(t), t$2))+4*s*diff(varphi(t), t$2);

solve(identity(subsindets(r=g,specfunc(:-diff),freeze),
               freeze(:-diff(theta(t),t))),
      [x,y,z,w]);
frontend(solve,[identity(r=g,diff(theta(t),t)),[x,y,z,w]],[{`+`,`*`,`=`,list,specfunc({identity})},{}]);

 

Could you please help me to extend the previous codes for this problem?

@acer 

WoW! Fantastic!

Thanks a lot for your great help.

I really do not how I can appreciate you.

You helped me so much.

Thank you very very very very .... very much.

I wish all the best to you and your family.

Sincerely yours,

J.

@acer 

This new way also works fine!

Thanks a lot.

Thank you very very much, dear friend.

Thanks for spending time. Thank you so much.

@acer 

Dear friend,

Thanks a lot. I really appreciate you. You solved the problem. Many thanks for your response and effort in helping me.

Do you have any idea about Talyor expansion I referred to it as well??

How I can freeze theta(t) in

taylor(cos(theta(t)),theta(t)=0 , 7);

to get the same result with 

taylor(cos(theta),theta=0 , 7)

?? 

The latter one works but the former one gives error.

However, perhaps freezing theta(t) is needless and I must try with other codes. I do not know. 

Sorry for bothering you.

@acer 

Dear friend,

Thanks a lot for your reply.

I did not present further examples because it is needless as the problem occurs in the first simple example as well.

You can see that when we wrote with(Physics) at the top of our codes, the answer changes.

As you stated, It seems that I must change the frozen part, but "Physic:-diff" does not work.

I really wonder about "freeze" in general; I think that freeze is one of the maple's Bugs, I mean that maple should identify these automatically because the issue does not just occur about that; For example, in Taylor's expansion of

taylor(cos(theta(t)),theta(t)=0 , 7);
maple returns error while when I remove the t-dependence as

taylor(cos(theta),theta=0 , 7)
I can get the result. Indeed, maple is incapable of recognizing such things !!

I feel that for this problem, I have to use freeze as well to reach a conclusion. However, I failed for this problem as well.

Anyway, if I find a solution for both, I shall write here.

Thanks a lot for your evergreen helps. 

@acer 

Dear friend, as I promised, I tried complicated cases and for one of them still, it works fine.

But there is an important point about using this way:

We have to not use the "physics package".

I am using two packages: with(LinearAlgebra): and with(Physics)

when I disabled the physics package this approach works, otherwise it returns the problem that I had with other approaches.

@acer 

Dear friend,

Thank you very much. 

You made me happy.

Thanks a lot.

For the complicated examples, I will try using your proposed approach. If there is a problem, I will write here.

Thank you so much.

Best Wishes.

@Kitonum 

Dear friend, 

Thank you very much for your explanations.

Best Wishes

@Carl Love 

Yes, you are right.

I think the reason is that in fracdiff, the form of the function must be specified when we use Laplace--option.

In my case, the form of first element, i.e. function, is unknown, hence we ecounter with error.

@MapleMathMatt 

Dear friend,

Thank you very much.

You made me happy.

I also used Laplace-approach but in another way. 

Your way is so interesting.

Thanks a lot.

Yours Sincerely,

J.

@Carl Love 

Dear friend,

Thanks a lot for your reply.

You made me happy.

Yes you are right, but since this is a part of codes that I am writing and this one is an example, hence I should wirte it in general form.

Yours Sincerely,

J.

I found an alternative way to get a solution:

restart;
with(inttrans):
eq := fracdiff(u(x, t), t, 1/2) = t^(1/2)*sin(x)/GAMMA(3/2):
eval(invlaplace((solve(laplace(eq, t, s), laplace(u(x, t), t, s))), s, t),u(x,0)=0);

sin(x)*t

@Kitonum 

Dear Friend,

Thanks a lot. You made me happy. 

I really do not know how I can appreciate you.

Thanks again.

Best wishes.

@tomleslie 

Dear Friend,

Thanks a lot. You made me happy. 

I really do not know how I can appreciate you.

Thanks again.

Best wishes.

@vv 

Dear Friend,

Thanks a lot. Yes, you are right.

1 2 3 4 Page 2 of 4