Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- function z = zloty_podzial(x,g,e)
- a = -100;
- b = 100;
- z=b;
- xn=[b-(b-a)*(sqrt(5)-1)/2; a+(b-a)*(sqrt(5)-1)/2];
- f= [funkcja(x-xn(1)*g);funkcja(x-xn(2)*g)];
- while norm(b-a)>0.01*e
- if f(1)>f(2)
- a=xn(1);
- xn(1)=xn(2);
- xn(2)=a+(b-xn(1));
- else
- b=xn(2);
- xn(2)=xn(1);
- xn(1)=a+(b-xn(2));
- end
- f = [funkcja(x-xn(1)*g);funkcja(x-xn(2)*g)];
- end
- if f(1)>f(2)
- z = xn(2);
- else
- z = xn(1);
- end
- if z - e < a & z + e > a
- a = a - 500;
- end
- if z + e > b & z - e < b
- b = b + 500;
- end
- end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement