Lista 04 Questão 06

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

# -*- coding: utf-8 -*-
"""
Created on Mon Oct 09 02:32:50 2017
@author: Vinícius
"""
from raiz import ipf
def f(d):
 k1 = 40000.
 k2 = 40.
 m = 95.
 g = 9.81
 h = 0.43
 return (((2/k1)*((m*g*(d+h)) - (0.4*k2*d**2.5))))**0.5
 
d = ipf(f, 0.0001, 1E-7, 1000, True)
print d
"""
1 0.141568304185
2 0.16317972131
3 0.166229019648
4 0.166654639893
5 0.166713958728
6 0.166722224285
7 0.166723375984
8 0.166723536458
9 0.166723558818
0.166723558818
"""