Jump to content
Fórum Script Brasil
  • 0

Metodo da Secante


zondak
 Share

Question

def pi_secante (f, x0, x1, epsilon, iterMax=50):
    if abs(f(x0)) <= epsilon:
        return(False, x0)
        print("k\t x.\t\t f(x)\t\t")
    for k in range (1,iterMax+1):
        x1 = x0-f(x0)/flin(x0)
        print("%d\t%f\t%f"%(k,x1,f(x1)))
        if abs(f(x1)) <= epsilon:
            return(False, x1)
        x0 = x1
    print ("O numero Máx de interações foi atingido")
    return(True,x1)

Eu sei que esse método é bem parecido com o do newton mas eu não encontro nenhum exemplo na internet da forma que eu quero, so resta alterar esse corpo do jeito certo para o método da secante. Se alguém poder me ajudar ficarei muito grato.

Link to comment
Share on other sites

1 answer to this question

Recommended Posts

  • 0

Eu fiz assim:

def Secante(f,xat,xant,*args):
    """
    Primeiro argumento livre é o contador e o segundo é a tolerância. Se não for preenchido, tolerância é de 10⁻5
    e o contador máximo é 200
    """
    if not args:
        tol = 1e-5
        contador = 200
    elif len(args) == 1:
        tol = 1e-5
        contador = args[0]
    else:
        tol = args[1]
        contador = args[0]
    cont = 0
    tol = 1e-5
    while True:
        xsuc = xat - ((xat - xant)/(f(xat) - f(xant)))*f(xat)
        erro = abs((xsuc - xat)/xsuc)
        xant = xat
        xat = xsuc
        cont += 1
        if cont>contador:
            break
        if erro<tol:
            break
    return float('%.2f'%xat)

Link to comment
Share on other sites

Join the conversation

You can post now and register later. If you have an account, sign in now to post with your account.

Guest
Answer this question...

×   Pasted as rich text.   Paste as plain text instead

  Only 75 emoji are allowed.

×   Your link has been automatically embedded.   Display as a link instead

×   Your previous content has been restored.   Clear editor

×   You cannot paste images directly. Upload or insert images from URL.

 Share

  • Forum Statistics

    • Total Topics
      149.6k
    • Total Posts
      646.2k
×
×
  • Create New...