2 years, 79 days

## complicated case...

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})},{}]);

## Many thanks...

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.

## Thanks a lot....

This new way also works fine!

Thanks a lot.

Thank you very very much, dear friend.

Thanks for spending time. Thank you so much.

## Thanks a lot....

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

Dear friend,

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.

## A point...

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.

## Thanks a lot....

Dear friend,

Thank you very much.

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.

## Thanks a lot....

Dear friend,

Thank you very much for your explanations.

Best Wishes

## @Carl Love  Yes, you are right. I...

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.

## Thanks a lot...

Dear friend,

Thank you very much.

I also used Laplace-approach but in another way.

Thanks a lot.

Yours Sincerely,

J.

## Thanks a lot...

Dear friend,

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.

## A solution...

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

## Thanks a lot....

Dear Friend,

Thanks a lot. You made me happy.

I really do not know how I can appreciate you.

Thanks again.

Best wishes.

## Thanks a lot....

Dear Friend,

Thanks a lot. You made me happy.

I really do not know how I can appreciate you.

Thanks again.

Best wishes.

## Thanks a lot....

Dear Friend,

Thanks a lot. Yes, you are right.

 1 2 3 4 Page 2 of 4
﻿