codigo eq1

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

clear
clc
m=10;
F0=10;
wn=60;
wa=90;
wb=62;
k=36000;
Sst=10/36000;
t = 0:0.1:10;
plot(t, ((F0/(k-m*wa^2))*cos(wn*t)+(Sst/((wa/wn)^2-1))*cos(wa*t)), '-b;x(t) da Eq a;');