codigo area 3 numérico UFRGS

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

deff('dy=f(y,t)','dy=(t-y)/t')
h=0.001
N=1000
t0=2
y0=3
[y,t]=euler(f,t0,y0,h,N)
plot(t,y,'-o')
[y2,t2]=eulermelhorado(f,t0,y0,h,N)
plot(t2,y2,'-ko')
//t1=0:0.001:10
//y1=exp(t1)./(1+exp(t1))
//y1=exp(-2*t1)
//plot(t1,y1,'r')