lab11exe3

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

function F=F(Y,t)
 F=[100-0.1*Y(2);
 0.01*Y(1)-0.1*Y(2)]
endfunction
t0=0
y0=[0; 
 0]
h=0.1
N=20000
[y,t]=runge4(F,t0,y0,h,N)
plot(t,y,'o-')
//Esse faz direto porque é de primeira ordem nao precisa trocar , apenas Y(1) direto )
//NESSE EXERCICIO EXERCICIO O PRIEMRIO ´P é a potencia o segundo o verde é a temperatura 
//quando pede calcule a solução é o grafico
//O tempo que demora para atingir 900 k olhar o grafico verde t=220
//Esse exercicio nao preciso trocar as variaveis 
//T'(t)=0.01P(t)-0.01T(t)
//P'(t)=100-0.1T