# Answer to Question #38062 in MatLAB | Mathematica | MathCAD | Maple for Sondos

Question #38062

k=input(' Enter the value of the spring stiffness k: ');

c=input(' Enter the value of the damping coefficient c: ');

m=input(' Enter the value of the mass: ');

t0=input(' Enter the initial time t0: ');

tf=input(' Enter the final time tf: ');

x0=input(' Enter the initial displacement x0: ');

v0=input(' Enter the initial velocity v0: ');

x=[x0 v0];

t=[t0 tf];

cr=2*sqrt(k*m);

z=c/cr;

disp(['the result ' num2str(z)]);

if z<1

disp('It is an Under-Damped Response')

[t1,x1] = Pro_1(k,c,m,t,x);

plot(t1,x1(:,1),t,x1(:,2),'Linewidth',2.5)

end

The Function is :

function [t1,x1]=Pro_1(k,c,m,t,x)

[t1,x1] = ode45(rif,t,x);

function xdot = rif(t1,x1)

xdot= zeros(2,1);

xdot(1)=x1(2);

xdot(2)=(-k/m)*x1(1)-(c/m)*x1(2);

end

end

Can anyone help me please ???

c=input(' Enter the value of the damping coefficient c: ');

m=input(' Enter the value of the mass: ');

t0=input(' Enter the initial time t0: ');

tf=input(' Enter the final time tf: ');

x0=input(' Enter the initial displacement x0: ');

v0=input(' Enter the initial velocity v0: ');

x=[x0 v0];

t=[t0 tf];

cr=2*sqrt(k*m);

z=c/cr;

disp(['the result ' num2str(z)]);

if z<1

disp('It is an Under-Damped Response')

[t1,x1] = Pro_1(k,c,m,t,x);

plot(t1,x1(:,1),t,x1(:,2),'Linewidth',2.5)

end

The Function is :

function [t1,x1]=Pro_1(k,c,m,t,x)

[t1,x1] = ode45(rif,t,x);

function xdot = rif(t1,x1)

xdot= zeros(2,1);

xdot(1)=x1(2);

xdot(2)=(-k/m)*x1(1)-(c/m)*x1(2);

end

end

Can anyone help me please ???

Expert's answer

Dear visitor,

Unfortunately, the statement of your question is a bit confusing.

What do you need to be done? Please, give us the details so we could help you.

Unfortunately, the statement of your question is a bit confusing.

What do you need to be done? Please, give us the details so we could help you.

## Comments

## Leave a comment