I=6;
N=5;
B= [0 ;1];
C=eye(2);
A=[0 1 ;-I -N];
D=[0 ;0];
h=1;
s = sym('s');

WPC=[1;-1];

Aprim=expm(h*A);
Bprim=double( int( expm( s*A )*B, s, 0, h ) );
%Bprim=integral(fun,0,h);
Cprim=C;



sim('lab3JBmodel2');
figure()
hold on
plot(tout,x1);
plot(h*(1:length(x2)),x2,'o');
xlabel('t');
ylabel('x(t)');
hold off

----------------

x0=1.5;

WPC=[-1.7, -1, -0.75, 0, 0.75, 1, 1.7,];


% analiza równania różnicowego
for i=1:length(WPC)
    figure()
    hold on
    lama=WPC(i);
    sim('lab3JBmodel');
    plot(1:length(out),out,'o');
    plot(tout,out2);
    grid on;
    xlabel('k');
    ylabel('x(k)');
    title(['Układ dyskretny, lambda/alfa = ' num2str(lama)]);
    hold off
end

%analiza równania różniczkowego