load Statischen_Kennlinie_M1.mat;  %Daten aufnehmen
sdata= M(:,1);			%Werte Lage	
idata=M(:,2);			%Werte Strom
%xdata= [sdata, idata];
ydata=M(:,3);			%Werte Kraft
x= [1,1];			%Lösungsvektor: x(1) entspricht a, x(2) entspricht theta_0
%Berechnung mit einer Eingangsmatrix xdata
% f = @(x, xdata)x(1)*(((xdata+x(2)).*(xdata+x(2))./((k2-2*xdata)...
%     .*(k2-2*xdata))) - ((-xdata+x(2)).*(-xdata+x(2))./((k1+2*xdata)...
%     .*(k1+2.*xdata))));
%Berechnung mit zwei Eingangsvektoren sdata und idata
f = @(x, sdata, idata)x(1).*(((idata+x(2)).*(idata+x(2))./((k2-2*sdata)...
    .*(k2-2*sdata))) - ((-idata+x(2)).*(-idata+x(2))./((k1+2*sdata)...
    .*(k1+2.*sdata)))); 
x = lsqcurvefit(f,0.3, sdata, idata, ydata);
% x = lsqcurvefit(f,0.3, xdata, ydata);