function [x,fit,param]=sinfit(x,y)
param=fminsearch(@fun,[1 1 1],[],x,y);  % Löse Optimierungsproblem, Gesucht sind die Koeffizieten c1,c2,c3
fit=param(1)*sin(x*param(2)+param(3));  % Lösung + Plot
figure(1); 
plot(x,y); hold on; plot(x,fit); hold off;
legend('Original','Fit');

function ret=fun(c,x,y)     % Optimierungsfunktion
yy=c(1)*sin(x*c(2)+c(3));   % c1*sin(x*c2+c3) entspricht A*sin(w0*x+phi), A: Amplitude, w0: Frequenz, phi: Phasenverschiebung
ret=sum(abs(y-yy));         % Summe über alle absoluten Differenezen von Testdaten und momentaner Lösung