LAB11EXE1

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

function F=F(Y,t)
 F=[0.1*(10-Y(2)^2)*Y(1)-Y(2);
 Y(1)]
endfunction
t0=0
y0=[1; 
 1]
h=0.1
N=200
[y,t]=runge4(F,t0,y0,h,N)
plot(t,y,'o-')
//azul=z nao importa nesse caso
//verde y
//se eu quiser plotar apenas o y no caso o verde plot(t,y(2,:),'o-') vai sair em azul agora