Resolution of an example of a twostory building using the od

Solution

one of the method is here: first we have to consider a function and next we have to initialise the values and by simple steps we get the ans, have a look: function Second_Oder_Ode % SOLVE d2x/dt2+5 dx/dt - 4 x = sin(10 t) % initial conditions: x(0) = 0, x\'(0)=0 t=0:0.001:3; % time scale initial_x = 0; initial_dxdt = 0; [t,x]=ode45( @rhs, t, [initial_x initial_dxdt] ); plot(t,x(:,1)); xlabel(\'t\'); ylabel(\'x\'); function dxdt=rhs(t,x) dxdt_1 = x(2); dxdt_2 = -5*x(2) + 4*x(1) + sin(10*t); dxdt=[dxdt_1; dxdt_2]; end end