clc
clear all
close all
global V S0
V=5;
S0=50;
x0=[0.2; 0.3; 0.1; 0.4];
w0=S0*x0;
opcje=odeset('Events',@fstop);
[t,w]=ode45(@fdestylacja,[0 10], w0, opcje);
subplot(2,1,1)
plot(t,w)
xlabel('czas [s]')
ylabel('ilosc moli składnika w wyparce [kmol]')
legend('A','B','C','D')
[lw,lk]=size(w)
SK=sum(w(lw,:))
w01=w(lw,:)'+x0*25;
[t1,w1]=ode45(@fdestylacja,[0 10], w0, opcje);
subplot(2,1,2)
plot(t1,w1)
xlabel('czas [s]')
ylabel('ilosc moli składnika w wyparce [kmol]')
legend('A','B','C','D')
figure
plot(t,w,t1,w1)
{"html5":"htmlmixed","css":"css","javascript":"javascript","php":"php","python":"python","ruby":"ruby","lua":"text\/x-lua","bash":"text\/x-sh","go":"go","c":"text\/x-csrc","cpp":"text\/x-c++src","diff":"diff","latex":"stex","sql":"sql","xml":"xml","apl":"apl","asterisk":"asterisk","c_loadrunner":"text\/x-csrc","c_mac":"text\/x-csrc","coffeescript":"text\/x-coffeescript","csharp":"text\/x-csharp","d":"d","ecmascript":"javascript","erlang":"erlang","groovy":"text\/x-groovy","haskell":"text\/x-haskell","haxe":"text\/x-haxe","html4strict":"htmlmixed","java":"text\/x-java","java5":"text\/x-java","jquery":"javascript","mirc":"mirc","mysql":"sql","ocaml":"text\/x-ocaml","pascal":"text\/x-pascal","perl":"perl","perl6":"perl","plsql":"sql","properties":"text\/x-properties","q":"text\/x-q","scala":"scala","scheme":"text\/x-scheme","tcl":"text\/x-tcl","vb":"text\/x-vb","verilog":"text\/x-verilog","yaml":"text\/x-yaml","z80":"text\/x-z80"}