Advertisement
Guest User

Untitled

a guest
Nov 28th, 2015
68
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.66 KB | None | 0 0
  1. function z = zloty_podzial(x,g,e)
  2. a = -100;
  3. b = 100;
  4. z=b;
  5. xn=[b-(b-a)*(sqrt(5)-1)/2; a+(b-a)*(sqrt(5)-1)/2];
  6. f= [funkcja(x-xn(1)*g);funkcja(x-xn(2)*g)];
  7. while norm(b-a)>0.01*e
  8. if f(1)>f(2)
  9. a=xn(1);
  10. xn(1)=xn(2);
  11. xn(2)=a+(b-xn(1));
  12. else
  13. b=xn(2);
  14. xn(2)=xn(1);
  15. xn(1)=a+(b-xn(2));
  16. end
  17. f = [funkcja(x-xn(1)*g);funkcja(x-xn(2)*g)];
  18. end
  19. if f(1)>f(2)
  20. z = xn(2);
  21. else
  22. z = xn(1);
  23. end
  24.  
  25. if z - e < a & z + e > a
  26. a = a - 500;
  27. end
  28.  
  29. if z + e > b & z - e < b
  30. b = b + 500;
  31. end
  32. end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement