# Exemplo de método de reta tangente ou de Euler
F<-function(s,u){   # s é o parâmtro do tempo e u é um vetor.
2                   # expressão da função
}

a=0                 # início do intervalo de tempo
b=4                 # final do intervalo de tempo
n=4                 # número de divisões do intervalo

h=(b-a)/n
t=seq(a,b,by=h)
y=0*t;
y[1]=1              # Foi dado que y(0)=1

y[2]=y[1]+h*F(t[1],y[1])
y[2]


n=length(y)

for ( i in 1:(n-1)){
y[i+1]=y[i]+h*F(t[i])
}


plot(t,y,col="red")