freact1.m

来自「Incorporating Prior Knowledge in Cubic S」· M 代码 · 共 7 行

function dc = freact1(t,c)
global b1 b2 b3 b4
dc = zeros(3,1);
dc(1) = -b1*c(1)+b2*c(2);
dc(2) = b1*c(1)-b2*c(2)-b3*c(2)+b4*c(3);
dc(3) = b3*c(2)-b4*c(3);