lab11exe4

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

function F=F(Y,t)
 F=[Y(2);
 -9.8*sin(Y(1))]
endfunction
t0=0
y0=[0.5; //é o primeiro o primeiro porque consideramos o o como o primeiro da matriz
 0]
h=0.1
N=100
[y,t]=runge4(F,t0,y0,h,N)
plot(t,y,'o-')