Differentiation of the position vector.

I have written a code for postion of a vibrating system. I am unable to determine the velocity and acceleration from the expression u. I want to diffrentiate u(i) in the same loop and plot it with t. How this can be done? Please Help.

m1 = 40;

m = m1*1000/9.81;

k = 3500*1000;

u(1) = 0.01; v(1) = 0.1;

wn = sqrt(k/m);

i=1;

for t=0:.001:1

u(i) = u(1)*cos(wn*t) + (v(1)/wn)*sin(wn*t);

i=i+1;

end

t= 0:0.001:1;

plot(t,u);

grid on

Answer 1

[fred ssemwogerere]

Hello, you could make use of the symbolic math toolbox. The page below should get you started:

https://mathworks.com/help/symbolic/performing-symbolic-computations.html