%Gleichung
%Minimum einer Funktion bestimmen

clear all
tic
%Definition der Konstanten Parameter

syms B C K

lgN=[3.74 3.743 2.623 2.799];
lgs=[2.882 2.882 2.919 2.919];

ysym=sum(sum((lgN-B-K./(lgs-C)).^2));

[B,C,K] = meshgrid (-25:0.2:25);

%[B] = meshgrid (-23:1:-21)
%[C] = meshgrid (1:1:3)
%[K] = meshgrid (24:1:26);

y=eval(ysym);

%alte Funktion, neuerdings durch die Doppelsumme berechnet
%y=((lgN1-B-K./(lgs1-C)).^2)+((lgN2-B-K./(lgs1-C)).^2)+((lgN3-B-K./(lgs2-C)).^2)+((lgN4-B-K./(lgs2-C)).^2);

[miny,ind]=min(y(:));
[cmin,bmin,kmin]=ind2sub(size(y),ind);
bminwert=B(1,bmin,1)
cminwert=C(cmin,1,1)
kminwert=K(1,1,kmin)
time=toc