메트렙 고수 횽들 ㅠㅠㅋ

여기 좆뉴비가 횽들에게 물어볼게 있어서 왔어...

내가 지금 메트렙 과제 중인데

1. Red + Blue + Green → Cyan with rate constant k1
2. Cyan → Red + Blue + Green with rate constant k2
3. Cyan → Green + Magenta with rate constant k3
4. Cyan → Blue + Magenta with rate constant k4

위에 4가지 화학식을 가지고 5개의 물질의 양 변화를 그래프로 나타내야 하는데...

일단 밑에 내가 발로 짜본 메트렙...ㅜㅜ

function Y = ee150_chemical2(sub, enzyme1, enzyme2, k)
dt = 0.1;
eps = 1e-7;
Y = zeros(5,1);
Y(1,1) = sub;
Y(2,1) = enzyme1;
Y(3,1) = enzyme2;

i=1;

while Y(1,i)>eps;
    Y(1,i+1)=Y(1,i) - k(1)*Y(1,i)*Y(2,i)*Y(3,i)*dt + k(2)*Y(4,i)*dt;
    Y(2,i+1)=Y(2,i) - k(1)*Y(1,i)*Y(2,i)*Y(3,i)*dt + k(2)*Y(4,i)*dt + k(4)*Y(4,i)*dt;
    Y(3,i+1)=Y(3,i) - k(1)*Y(1,i)*Y(2,i)*Y(3,i)*dt + k(2)*Y(4,i)*dt + k(3)*Y(4,i)*dt;
    Y(4,i+1)=Y(4,i) + k(1)*Y(1,i)*Y(2,i)*Y(3,i)*dt - k(2)*Y(4,i) - k(3)*Y(4,i)*dt - k(4)*Y(4,i)*dt;       
    Y(5,i+i)=Y(5,i) + k(3)*Y(4,i)*dt + k(4)*Y(4,i)*dt;
    Y(:,i)
    i=i+1;
end

t = 0:dt:dt*(i-1);
figure;
plot(t,Y(1,:),\'r+\');
hold on;  
plot(t,Y(2,:),\'b+\');
plot(t,Y(3,:),\'g+\');
plot(t,Y(4,:),\'c+\');
plot(t,Y(5,:),\'m+\');
end

요기까지 해서리 커맨드창에서

conc = ee150_chemical2(5e-5, 3e-5, 5e-5, [1e8; 1e-3; 1e-2; 1e-1]);

입력하면 알흠다운 그래프가 떠야하는데.....

??? Error using ==> plot
Vectors must be the same lengths.

Error in ==> ee150_chemical2 at 27
plot(t,Y(1,:),\'r+\');

계속 이렇게 뜨면서 그래프가 안나와 ㅠㅠ

교수는 제대로 가르쳐 주지도 않고 이거 해보라 그러고..ㅠㅠ 난 그저 머리 쥐싸메고 하다하다 그로기 상태..@_@ㅋ

이거 좀 속시원하게 고쳐줄 용자횽 안계신가요?...