請問有關RUNGE KUTTA 換式問題

[Jerryhk]

小弟是新手, 只是上網查了一部分ODE 相關而又適用的指令, 如果問上一些低級問題請別見怪...

現時情況是這樣, 我要用RUNGE KUTTA 去解5條 ODE, 寫成這樣:

function output = systemRHS(t,x)
vm=0.5*t
imax=30000;
i1= x(1);
i2= x(2);
i3= x(3);
v1= x(4);
v2= x(5);
if(v2>vm){
output1=(v2-(500*i1))/1000;
output2=(v1-v2-(50*i2))/0.00005;
output3=(v1-380)/100;
output4=-(i3+0)/.001;
output5=(i2-i1)/.0001;
output=[output1; output2; output3; output4; output5];

但MATLAB RUN 很久都說BUSY, 之後PLOT 圖係有WAVE, 但又不知道對不對....
想問一下各位高手我有沒有寫錯?
另外還有一個情況, 因為那5條式是用來模擬真空斷路器, 因此中途要按不同的 V 和 I 而決定替換 I2的式, 請問這個又可以如何做到? 我網上查了很久都找不到相關做法:(....

先謝謝大家:)