Estude fácil! Tem muito documento disponível na Docsity
Ganhe pontos ajudando outros esrudantes ou compre um plano Premium
Prepare-se para as provas
Estude fácil! Tem muito documento disponível na Docsity
Prepare-se para as provas com trabalhos de outros alunos como você, aqui na Docsity
Encontra documentos específicos para os exames da tua universidade
Prepare-se com as videoaulas e exercícios resolvidos criados a partir da grade da sua Universidade
Responda perguntas de provas passadas e avalie sua preparação.
Ganhe pontos para baixar
Ganhe pontos ajudando outros esrudantes ou compre um plano Premium
Trabalho computacional dois da disciplina de EDC
Tipologia: Trabalhos
1 / 9
Esta página não é visível na pré-visualização
Não perca as partes importantes!
Trabalho Computacional 2 – Equações Diferenciais C Fernanda Lima Abreu 1- Letra A (t) = var('t') y = function('y')(t) de = diff(y,t) == -(0.000120968094y) ys = desolve(de, y, ivar=t, ics=[0,1]) desolve_rk4(de,y,ivar=t,ics=[0,1],output='plot') show(ys) Letra B- (t) = var('t') y = function('y')(t) de = diff(y,t) == -(0.000120968094y) ys = desolve(de, y, ivar=t, ics=[0,1])
desolve_rk4(de,y,ivar=t,ics=[0,1],end_points=22000,xmax=50000,xmin=0,output='plot' ) show(ys) show(P) Letra C- Decaimento para essa amostra chegar a 0: (t) = var('t') y = function('y')(t) de = diff(y,t) == -(0.000120968094*y) ys = desolve(de, y, ivar=t, ics=[0,1]) P = desolve_rk4(de,y,ivar=t,ics=[0,0.015],end_points=20000,output='plot') show(ys) show(P)
y = function('y')(t) de = diff(y,t) == (0.06(30-70)) ys = desolve(de, y, ivar=t, ics=[0,70]) P = desolve_rk4(de,y,ivar=t,ics=[0,70],xmin=0,xmax=15,output='plot') show(ys) show(P) Letra A- (t) = var('t') y = function('y')(t) de = diff(y,t) == (0.06(70-78)) ys = desolve(de, y, ivar=t, ics=[0,78]) P = desolve_rk4(de,y,ivar=t,ics=[0,78],xmin=0,xmax=15,output='plot') show(ys) show(P)
Letra B - (t) = var('t') y = function('y')(t) de = diff(y,t) == (0.06*(70-72)) ys = desolve(de, y, ivar=t, ics=[0,72]) P=desolve_rk4(de,y,ivar=t,ics=[0,72],ymax=82,ymin=68,xmin=0,xmax=20,output='plo t',end_points=20) show(ys) show(P) L1 = line([(0,79.9),(20,79.9)],color='grey',linestyle='--') L2 = line([(0,70),(20,70)],color='grey',linestyle='--') L3 = line([(10,0),(10,82)],color='green',linestyle='--') L4 = line([(15,0),(15,82)],color='green',linestyle='--') P+L1+L2+L3+L
Letra C- (t) = var('t') y = function('y')(t) de = diff(y,t) == (0.06*(71-72)) ys = desolve(de, y, ivar=t, ics=[0,72]) P = desolve_rk4(de,y,ivar=t,ics=[0,72],ymax=82,ymin=68,xmin=0,xmax=20,output='plot',e nd_points=20) show(ys) show(P) L1 = line([(0,79.9),(20,79.9)],color='grey',linestyle='--') L2 = line([(0,70),(20,70)],color='grey',linestyle='--') L3 = line([(10,0),(10,82)],color='green',linestyle='--') L4 = line([(15,0),(15,82)],color='green',linestyle='--') P+L1+L2+L3+L
Letra D- Pode-se adicionar o queijo no tempo zero, ou seja, no instante inicial. A temperatura será ,portanto, de 72 graus, mantendo as propriedades do queijo. Caso a temperatura abaixe para menos de 70 graus, podemos dizer que em aproximadamente 17 minutos, o queijo ira talhar ou coagular, perdendo suas características. 3 - Letra A (J,kb,km,ka,i,Ta,Tb,Tm) = var('J kb km ka i Ta Tb Tm') J = 0.01; kb = 0.025; km = 0.025; ka = 0.015; i(t)=1; Ta = -(kadiff(x^2,t,2)); Tb = kbi ; Tm = -(km90) F(t)= kbi x(t) = function('x')(t) eqdif = Jdiff(90,t,2) == F(t) F(t) =Ta+Tb+Tm Delta = 0 - 4J show(F) show('Delta = ',Delta)
