map214ab.py

""" Programa de teste do python sobre bisseccao """
from math import *

## Defino o modulo que executa o metodo da bisseccao:

def bissec(f,a,b,n):
    xi,xf=a,b
    count=n
    while count > 0:
        raiz = 0.5*(xi+xf)
        if f(xi)*f(raiz) < 0:
            xf=raiz
        elif f(xi)*f(raiz) > 0:
            xi=raiz

        else:
            break
        count=count-1
    return raiz

## defino as funcoes
#w1,w2 = 0.8,1.0
w1=float(raw_input("Qual o lado1 em metros: "))
w2=float(raw_input("Qual o lado2 em metros: "))
a=95
arad=(a/180)*pi

L = lambda b : w1/sin(b)+w2/sin(arad+b)
dL = lambda b : -w1*cos(b)/sin(b)**2 - w2*cos(arad+b)/sin(arad+b)**2

x1,x2 = 0.1,1.6
n=100
if dL(x1)*dL(x2) < 0:
    root= bissec(dL,x1,x2,n)
    print " O comprimento maximo eh ", L(root)
else:
    print "Nao sei fazer"