function cost=lscost(par,tdata,xdata)

m=par(1);
c=par(2);
k=par(3);
x0=[1; 0];
options=odeset('RelTol',1e-5,'AbsTol',1.e-5);
[t,xmodel]=ode45('fcnspring',tdata,x0,options,m,c,k);
cost=sum(abs(xmodel(:,1)-xdata(:,1)).^2);