%20.9.03. Finn Haugen
clear

%Løsning deloppgave 1:
q=0.2;
V=1;
A=-q/V;
B=1/V;
C=1;
D=0;

%Løsning deloppgave 2:
blandetank=ss(A,B,C,D);
tslutt=60;
[cA,t]=step(blandetank,tslutt);
figure(1)
plot(t,cA*0.8)
grid
ylabel('cA [mol]')
xlabel('t [sek]')

%Løsning deloppgave 3:
tslutt=60;
t=[0:0.1:tslutt]';
wA=0.8*ones(length(t),1);
cA=lsim(blandetank,wA,t);
figure(2)
subplot(2,1,1)
plot(t,wA)
grid
ylabel('wA [mol]')
xlabel('t [sek]')
subplot(2,1,2)
plot(t,cA)
grid
ylabel('cA [mol]')
xlabel('t [sek]')


