function y=egyenes(x0,x1)
				% egyenes(x0,x1) az f fuggveny minimumat
				%keresi az x0,x1 pontok egyenesen
				%parabolikus modszerrel
f0=f(x0);f1=f(x1);if f1>f0 x=x1;x1=x0;x0=x;fx=f1;f1=f0;f0=fx;end;
%f0>f1 miatt x0->x1 lefele indulunk
iter=1;A=(3-sqrt(5))/2;B=(1-A)/A;x=x1+B*(x1-x0);fx=f(x);
while iter<1000&fx<f1 x0=x1;f0=f1;x1=x;f1=fx;x=x1+B*(x1-x0);fx=f(x);iter=iter+1;
end;x2=x;%x0,x1,x2 kozrefog egy minimumot
t0=0;t1=1;t2=1/A;%t parameterre parabolikus minimumkereses
f0=f(x0);f1=f(x1);f2=f(x2);
for i=1:30 f12=f1-f2;f10=f1-f0;t10=t1-t0;t12=t1-t2;
     nevezo=2*(f12*t10-f10*t12);
     if nevezo!=0 %a 3 pont nincs egy egyenesen
       t=t1-(f12*t10^2-f10*t12^2)/nevezo;
       ft=f(x0+t*(x1-x0));
       if t<t1&ft<f1 t2=t1;f2=f1;t1=t;f1=ft;
       elseif t<t1&ft>=f1 t0=t;f0=ft;
       elseif t>t1&ft<f1 t0=t1;f0=f1;t1=t;f1=ft;
       elseif t>t1&ft>=f1 t2=t;f2=ft;
       else t1=t1+eps;f1=f(x0+t1*(x1-x0));
       end;%ha egy egyenesen vannak:aranymetszeses lepes
     else if abs(t1-t0)>abs(t2-t1) T=t0;Y=t2;fT=f0;fY=f2;
          else T=t2;Y=t0;fT=f2;fY=f0;
          end;
          y=(1-A)*t1+A*T;fy=f(x0+y*(x1-x0));
          if fy<f1 t0=t1;f0=f1;t1=y;f1=fy;t2=T;f2=fT;
          else t0=Y;t2=y;f0=fY;f2=fy;
          end;
     end;
	  y=x0+t1*(x1-x0);%t parameterbol visszateres x-re
end;