lab10exe5

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

deff ('F=f(h,t)','F=0.005*(sin(t)+1)-0.025*h')
t0=0
y0=0.5
h=0.01
N=1000//função tende a ficar constante a partir de y=o.22 altura da piscina 
[y1,t1]=euler(f,t0,y0,h,N)
plot(t1,y1,'o-')
[y2,t2]=eulermelhorado(f,t0,y0,h,N)
plot(t2,y2,'ro-')
//[y3,t3]=predicao(f,t0,y0,h,N)
//plot(t3,y3,'go-')
//[y4,t4]=runge4(f,t0,y0,h,N)
//plot(t4,y4,'ko-')
//[y5,t5]=adams2(f,t0,y0,h,N)
//plot(t5,y5,'mo-')
//[y6,t6]=adams3(f,t0,y0,h,N)
//plot(t6,y6,'co-')
//[y7,t7]=adams4(f,t0,y0,h,N)
//plot(t7,y7,'yo-')