Matlab

clear
clc

% Graficar un señal senoidal a través del tiempo, que represente el
% comportamiento físico de un oscilador armónico forzado. Generando una
% solución transitoria y una soluciónestacionaria. La solución transitoria
% es idéntica a la de un oscilador amortiguado NO forzado.
% -------------------------------------------------------------------------
% Recibe: Masacolgante[Kg],Constante de resorte(k)[N/m],
% Constante de amortiguacion (Gamma) [Kg/s] , Fuerza oscilante impulsora[N]
% Devuelve: Tiempo[Seg], Amplitud[m]%--------------------------------------------------------------------------

m=input('Ingrese la masa del objeto oscilante[Kg]: ');
k=input('Ingrese la constante de resorte[N/m]: ');
ga=input('Ingrese la constante de amortiguación(Gamma) [Kg/s]: ');
F0=input('Ingrese la fuerza senoidal impulsora[N]: ');

%--------------------------------------------------------------------------
% Parametros Libres

xo=5;% Posición inicial del movimiento
vo=0; % Velocidad inicial del movimiento
faseRad=pi; % Angulo de fase de movimiento%--------------------------------------------------------------------------

% FORMULAS DE APLICACIÓN

wo=sqrt(k/m); % Frecuencia angular inicial

w=sqrt(wo^2-ga^2); % Frecunciaangular (En amortiguamiento)

A=sqrt(xo^2+((vo+ga*xo)/w)^2);

Phiforzado= pi/2;

Afor=F0;

x=linspace(0,3,10000); n=100;%--------------------------------------------------------------------------
for s=0:3;
% ESCOGENCIA DEL CASO

input = menu('Elija la opción deseada para los movimientos forzados o amortiguado:', 'Imagenes estáticas de los movimientos','Imagen dinámica MOA','Imagendinámica MOF:Solucion estacionaria','Imagen dinámica MOF: Solución transitoria');

switch input

case 1,

%Posición de un oscilador amortiguado y forzado

v=linspace(0,3.24,2001);...