Guest User

Untitled

a guest
Aug 18th, 2018
88
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Scilab 0.64 KB | None | 0 0
  1. function y = fun(x)
  2.     s=poly(0,"s")
  3.     p= s - 1  //Polinomio (ej: x - 1)
  4.     y = horner(p, x)
  5. endfunction
  6.  
  7. function  x = biseccion(a,b,n,tol2,tol1)
  8.  
  9. u = fun(a)
  10. v = fun(b)
  11. paso = 1
  12.  
  13. if sign(u) == sign(v) then
  14.    printf('Error la función debe cambiar de signo en (a,b)')
  15.    x = 0
  16. else
  17.     c=0
  18.     while ((b-a) * 0.5 > tol2 & paso < (n + 1))
  19.         c = (b+a) / 2
  20.         w = fun(c)
  21.         printf("Paso: %i\n", paso)
  22.        
  23.         if sign(u) == sign(w) then
  24.             a = c
  25.             u = w
  26.         else
  27.             b = c
  28.             v = w
  29.         end
  30.             paso = paso + 1
  31.     end
  32.  x = c
  33. end
  34.  
  35. endfunction
Add Comment
Please, Sign In to add comment