@Jarekkk Thanks man, your help took me out from this frustration which I experienced over 2 days. Finally, it's done:))...Jarekk, you might have noticed that my animation is a kind of related to a work of civil engineers, indeed this is who I am and explains why I am not so familiar with the programming. Thus I would really appreciate if your expertise knowledge in maple could help me to sort out one more issue. Over these few weeks I was thinking how to make that movement looks to be more natural and for that what command I need to use. Certainly, both columns and beam itself have particular cubic kind of deformation formula (might be quite complex depending on a nature of applied forces) which need to be incorporated with one of the relevant maple commands smth like polygon,line...

But, one thing I don't know is how to input those formulas into commands (for vertical columns and horizontal beam) and more worse make each of these structure elements to be oscillating. Things get even complicated when later I need to extend this algorithm to number of stories.....As a beginner, I am completely blind and don't know how and what to do, but somehow I have to make it...so, hope that more experienced maple users could help..Jarekk, have any idea and way how?

Thanks 

@Jarekkk Thanks man, your help took me out from this frustration which I experienced over 2 days. Finally, it's done:))...Jarekk, you might have noticed that my animation is a kind of related to a work of civil engineers, indeed this is who I am and explains why I am not so familiar with the programming. Thus I would really appreciate if your expertise knowledge in maple could help me to sort out one more issue. Over these few weeks I was thinking how to make that movement looks to be more natural and for that what command I need to use. Certainly, both columns and beam itself have particular cubic kind of deformation formula (might be quite complex depending on a nature of applied forces) which need to be incorporated with one of the relevant maple commands smth like polygon,line...

But, one thing I don't know is how to input those formulas into commands (for vertical columns and horizontal beam) and more worse make each of these structure elements to be oscillating. Things get even complicated when later I need to extend this algorithm to number of stories.....As a beginner, I am completely blind and don't know how and what to do, but somehow I have to make it...so, hope that more experienced maple users could help..Jarekk, have any idea and way how?

Thanks 

@Jarekkk Sorry, I missed out the "ball" procedure. I attached the amended version, could take a look...I would really appreciate if you could point out where my errors are, and give some explanations to your changes...Thanks



Polygon_Problem.mw

@Jarekkk Thank you so much, u finally helped me to sort out the problem on which I was struggling for almost two days..Also, by just simply increasing the number of Digits to 15, it also helped me to solve those 2-additional problems that I posted later. 
Thanks

Could you explain me, how come even by using default Digits=10, the answer still can't be correctly obtained. I thought that Digits=10 is already too high  

@Jarekkk Thank you so much, u finally helped me to sort out the problem on which I was struggling for almost two days..Also, by just simply increasing the number of Digits to 15, it also helped me to solve those 2-additional problems that I posted later. 
Thanks

Could you explain me, how come even by using default Digits=10, the answer still can't be correctly obtained. I thought that Digits=10 is already too high  

Maple Experts please heeeelp
1-prob
I am experiencing another problem with the same question I posted above.
Why does the maximize command not perform the way it supposed to do?
The maximum point of "sol(t)" in a range "t=0..td" is smth around 0.314159009425
(hand calc)
while

>maximize(sol(t),t=0..dr);  

returns

>0.4003673052

why on earth the answer is different...aaahhh:(((( 

2-prob

also, why 

>sol(td);
>0.6359172952

it has to be 0.314159009425

can someone help to get the right value

thx

 

I tried before, however, the graphs are still different, besides now they look
as a kind of pulsation, not a solid line as before were

I tried before, however, the graphs are still different, besides now they look
as a kind of pulsation, not a solid line as before were

to gkokovidis
the issue u mentioned above about 'pi'  being wrong, as it has to be 'Pi', actually due to the way an internet converts and shows 'Pi'...The coding I included above is the same as I wrote into a worksheet, nothing esle to add...the main problem, is that why animate command does not recognise points...

to gkokovidis
the issue u mentioned above about 'pi'  being wrong, as it has to be 'Pi', actually due to the way an internet converts and shows 'Pi'...The coding I included above is the same as I wrote into a worksheet, nothing esle to add...the main problem, is that why animate command does not recognise points...

Thank u,I made those changes, but the problem still arise, sorry it previous post, I missed out the value of dur

k:=20:
m:=4:

p[0] := 1:

tr := 4:

dur:=evalf(3*pi):

omega[n]:=sqrt(k/m):

T[n]:=(2 pi/omega[n]):

> maxpoint := proc (x) 
local tr, p1;

tr := evalf(x*T[n]):
p1 := piecewise(t1 <= tr, p[0]*t1/tr, p[0]):
maximize(int(p1*sin(omega[n]*(t-t1))/(m*omega[n]), t1 = 0 .. t), t = 0 .. dur) :

end proc;
> 
y:=ceil(tr/T[n]):
plot(maxpoint(tr1), tr1 = 0 .. y, sample = ['i' $ 'i'=0..y], numpoints = 20, adaptive = false);

Error, (in testeq) too many levels of recursion

Thank u,I made those changes, but the problem still arise, sorry it previous post, I missed out the value of dur

k:=20:
m:=4:

p[0] := 1:

tr := 4:

dur:=evalf(3*pi):

omega[n]:=sqrt(k/m):

T[n]:=(2 pi/omega[n]):

> maxpoint := proc (x) 
local tr, p1;

tr := evalf(x*T[n]):
p1 := piecewise(t1 <= tr, p[0]*t1/tr, p[0]):
maximize(int(p1*sin(omega[n]*(t-t1))/(m*omega[n]), t1 = 0 .. t), t = 0 .. dur) :

end proc;
> 
y:=ceil(tr/T[n]):
plot(maxpoint(tr1), tr1 = 0 .. y, sample = ['i' $ 'i'=0..y], numpoints = 20, adaptive = false);

Error, (in testeq) too many levels of recursion

uau, it works, thx u, I am just starting my learning process of maple, so most of the things I still don't know, but have to pick up fast as my final project is purely based on this software language..thx, and hope for your future help 

uau, it works, thx u, I am just starting my learning process of maple, so most of the things I still don't know, but have to pick up fast as my final project is purely based on this software language..thx, and hope for your future help 

Thank you for your respond, but my problem is that my function can be quite complex, it is defined by the user how extensive it can  go and what the period it will hold, then in that case I don't really know hot to manipulate each function element..
Let's say I have the following function...

pforce:=p*sin(omega1*t)+3*cos(3*omega1*t) ;
        here, omega1 is the natural frequency(defined by user), which can yield us the period (T=(2*pi/omega1));
                p any user defined value;
                t time range 
so if the period and function are defined by user, can be any value, then how to plot this function with that period in x-axis 

Thank u