LAB11.2

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

function F=F(Y,t)
 F=[(Y(2)^2-Y(1))/10;
 Y(3);
 0.1*(10-Y(1))*Y(3)-Y(2)]
endfunction
t0=0
y0=[1; 
 1;
 1]
h=0.1
N=1000
[y,t]=runge4(F,t0,y0,h,N)
plot(t,y,'o-')
//igual ao exercicio anterior