# 2010-09-30

clc
clear all
close all

A = 1; B = 1;
deltat = 0.001;

entrada = 1;

N = 1000;
n = 1;

f = zeros(N,1);
o = zeros(N,1);


while n <= N -1;
  f(n+1) = deltat*entrada + f(n);
   o(n+1) = -(B/A)*deltat*o(n) + (deltat/A)*f(n) + o(n);
  n += 1;
end


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