freact1.m

Incorporating Prior Knowledge in Cubic Spline Approximation - Application to the Identification of R

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);