Baixe o app para aproveitar ainda mais
Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original
deff ('F=f(v,t)','F=9.8-0.01*v^2') t0=0 y0=0 h=0.01 N=1000//quero no tempo 10 segunos encontrar meu v no caso foi 31.1874 [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-')
Compartilhar