# respuesta al examen parcial 1, por euler

clear all
clc
close all

A = 1; B = 1; Tp = 1;

N = 10000;
deltat = 0.0001;
o = zeros(N, 1);
n = 1;
while n<=N-1
  in(n) = rand;
  m = -((B+Tp)/A)*o(n)+Tp*in(n)/A;

  o(n+1) = m*deltat+o(n);
  n = n+1;
end


figure(1, "visible", "off");
plot(o);
grid;
print -deps "../media/examen_euler.eps";
