I like the way you solved it by changing variables! Thanks a lot!

I think the first plot in my file and the second plot in your file are correct.

I like the way you solved it by changing variables! Thanks a lot!

I think the first plot in my file and the second plot in your file are correct.

@Markiyan Hirnyk I wasn't trying to criticize you.If I made you feel that way, I am sorry. I do appreciate your help.

@Markiyan Hirnyk I wasn't trying to criticize you.If I made you feel that way, I am sorry. I do appreciate your help.

Markiyan, I did the same thing yesterday as you did above. Both maple and mathematica didn't work. I think it's probably a bug for maple and hopefully maple can fix it in the near future.

Anyway, thanks for your help!

Markiyan, I did the same thing yesterday as you did above. Both maple and mathematica didn't work. I think it's probably a bug for maple and hopefully maple can fix it in the near future.

Anyway, thanks for your help!

Thanks for your help. Two questions:

1. You mentioned that "#This will indicate that the result res2, where the substitution is done first, is the correct one:", I don't understand how to tell which one is the correct result from

res1:=eval(value(u),[cl=6300,ct=3140]):
res2:=value(eval(u,[cl=6300,ct=3140])):
simplify(res1-res2):
evalf(eval(%,{omega=1,t=1,tau=1,H=1}));

It just tell us res1-res2 is nonzero, how do you know res1 is wrong and res2 is correct?

2. Is there a way that I can carry out the pure symbolic result of u? Actually, that is what I really want, but "Int" is an inert representation of int and you mentioned value(u) did not give the right result.

Thanks again!

Thanks for your help. Two questions:

1. You mentioned that "#This will indicate that the result res2, where the substitution is done first, is the correct one:", I don't understand how to tell which one is the correct result from

res1:=eval(value(u),[cl=6300,ct=3140]):
res2:=value(eval(u,[cl=6300,ct=3140])):
simplify(res1-res2):
evalf(eval(%,{omega=1,t=1,tau=1,H=1}));

It just tell us res1-res2 is nonzero, how do you know res1 is wrong and res2 is correct?

2. Is there a way that I can carry out the pure symbolic result of u? Actually, that is what I really want, but "Int" is an inert representation of int and you mentioned value(u) did not give the right result.

Thanks again!

@Markiyan Hirnyk 

Hi, Marklyan, please take a look at this new file test.mw. Do you have any idea why u1test != u2test? I plot these two variables by assigning them the same values of parameters, but the results look so different. Thanks!

@Markiyan Hirnyk 

Hi, Marklyan, please take a look at this new file test.mw. Do you have any idea why u1test != u2test? I plot these two variables by assigning them the same values of parameters, but the results look so different. Thanks!

@Markiyan Hirnyk Thank you for your help!

@Markiyan Hirnyk Thank you for your help!

Thank you so much for the update.

Just want to make sure I understand you correctly. Do you mean that the value given by u1test is always correct while the u2test is correct only when (497700/157)*tau[l]+4720*t < H is satisfied?

Does that mean if I need a correct symbolic solution for u1 when (497700/157)*tau[l]+4720*t > H, I can never do that by using maple 12? Do I need maple 16?

Thanks again.

Thank you so much for the update.

Just want to make sure I understand you correctly. Do you mean that the value given by u1test is always correct while the u2test is correct only when (497700/157)*tau[l]+4720*t < H is satisfied?

Does that mean if I need a correct symbolic solution for u1 when (497700/157)*tau[l]+4720*t > H, I can never do that by using maple 12? Do I need maple 16?

Thanks again.

Thank you for your reply. I tried your method, however, when the program goes to "u1 := int(int(convert(sin(omega*(ts+ys/ct))*Heaviside(t-ts-ys/ct), piecewise, t),
 ys = H-cl*ts .. ct*ts), ts = H/(cl+ct) .. (H+cl*tau[l])/(cl+ct));", it just kept evaluating without giving me any result. Is that because I am using a different version of maple? I am using Maple 12.

Here is my orignal file "test" and modified file "test.m". Really appreciate your help.

test.mw   test_m.mw