Run Code
|
API
|
Code Wall
|
Misc
|
Feedback
|
Login
|
Theme
|
Privacy
|
Patreon
EulerPendulo
# Movimento de pêndulo simples # # g=9.8 # gravidade l=1 # comprimnto do fio do pêndulo a=0 # tempo inicial theta0=pi/6 # arco inicial thetav=0 # velocidade inicial b=45 # tempo final em segundos n=100000 h=(b-a)/n t=seq(a,b,by=h) # vetor com invervalos entre a e b, de h em h. y=matrix(0,2,n+1) # espaço para guardar as aproximações y_0, y_1, y_2, ..., y_15. y[,1]=c(theta0,thetav) # y_0 F<-function(s,u){ # Equação y'(t)=F(t,y(t)) c( u[2] , -sin(u[1])*g/l ) } for ( i in 1:n){ # Método de Euler y[,i+1]=y[,i]+h*F(t[i],y[,i]) # t_i=a+(i-1)*h } plot(t,y[1,],'l')
run
|
edit
|
history
|
help
0
Pendulum Experiment
lab1
Linear and Log linear Model
Integracao-Interpolacaa
jcruet_crickets and temp
Scatterplot using ggplot Sample
03-08-2020-Integral-Numerica
G
R Basic
24-09-2020Def-função