Run Code
|
API
|
Code Wall
|
Misc
|
Feedback
|
Login
|
Theme
|
Privacy
|
Patreon
EulerPendulo
Language:
Ada
Assembly
Bash
C#
C++ (gcc)
C++ (clang)
C++ (vc++)
C (gcc)
C (clang)
C (vc)
Client Side
Clojure
Common Lisp
D
Elixir
Erlang
F#
Fortran
Go
Haskell
Java
Javascript
Kotlin
Lua
MySql
Node.js
Ocaml
Octave
Objective-C
Oracle
Pascal
Perl
Php
PostgreSQL
Prolog
Python
Python 3
R
Rust
Ruby
Scala
Scheme
Sql Server
Swift
Tcl
Visual Basic
Layout:
Vertical
Horizontal
# 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')
[
+
]
Show input
Absolute running time: 1.09 sec, cpu time: 1.57 sec, memory peak: 57 Mb, absolute service time: 1,23 sec
edit mode
|
history
|
discussion