Prévia do material em texto
<p>[13]: #########################################################################</p><p>#O PROGRAMA ABAIXO RESOLVE O PVI:</p><p>#y'=y</p><p>#y(0)=1</p><p>#y(1)=?</p><p>#PELO MÉTODO DE EULER, COM UM PASSO h DADO PELO USUÁRIO</p><p>#########################################################################</p><p>def f(x,y): #DEFINE A FUNÇÃO QUE REPRESENTA A DERIVADA</p><p>return y</p><p>x0=0 #VALOR INICIAL DE x</p><p>xf=1 #VALOR FINAL DE x</p><p>y0=1 #VALOR INICIAL DE y</p><p>h=float(input('Diga o valor do passo h: ')) #PEDE AO USUÁRIO O VALOR DO PASSO h</p><p>n=round((xf-x0)/h) #CALCULA O NÚMERO DE ITERAÇÕES NECESSÁRIAS#</p><p>x=x0 #DEFINE O VALOR INICIAL DE x PARA ITERAÇÃO</p><p>y=y0 #DEFINE O VALOR INICIAL DE y PARA ITERAÇÃO</p><p>for i in range(n): #FAZ AS ITERAÇÕES SEGUNDO O MÉTODO DE EULER</p><p>x=x+h</p><p>y=y+f(x,y)*h</p><p>print('O valor de y({}) é {}'.format(1,y)) #IMPRIME UMA MENSAGEM COM O VALOR DE␣</p><p>↪→y(1) NA TELA</p><p>Diga o valor do passo h: 0.0001</p><p>O valor de y(1) é 2.7181459268252266</p><p>145</p>