fconc12.m

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

function y = fconc12(x)

global b1 b2 b3 b4
global c12
global cnew

b1 = x(1);
b2 = x(2);
b3 = x(3);
b4 = x(4);
[t,cnew] = ode45('freact1',[0:0.02:2],[1 0 0]);
c12new = cnew(round(linspace(1,101,12)),:);

y = sum(sum((c12new-c12).*(c12new-c12)));