himmelss.m

来自「蒙特卡罗方法及其例子」· M 代码 · 共 12 行

function ss = himmelss(k,data)
% sum-of-squares for Himmelblau 9.9

time = data.ydata(:,1);
Aobs = data.ydata(:,2);
y0   = data.y0;

[t,y] = ode45(@himmelode,time,y0,[],k);
Amodel = y(:,1);

ss = sum((Aobs-Amodel).^2);