function [J]=my_optim(x)
    global fit;
    global name;
    global SimOpt;
    global Result;
    tic
    sim(name, [0 10], SimOpt);
    origD=origData(:,2:end);
    fitD=fitData(:,2:end);
    J=sum(sqrt(power(origD-fitD,2)));
    fprintf('J = %.3f\tdt=%.1fs\t', J, toc);
for Cnt = 1:length(x)
    fprintf('x(%i)=%.2e\t', Cnt, x(Cnt));
end
fprintf('\n');
Fid = fopen('result.txt', 'a');
fprintf(Fid, 'J = %.3f \t', J);
for Cnt = 1:length(x)
    fprintf(Fid, 'x(%i)=%.2e \t', Cnt, x(Cnt));
end
fprintf(Fid, '\n');
fclose(Fid);