Archives
Trending
Support
Login
clear text
XML
Django
JavaScript
MATLAB
C
C++
C#
Python
SQL
Shell
Bash
Markdown
YAML
JSON
HTML
CSS
PHP
Java
Ruby
Go
Rust
Swift
Kotlin
Arduino
TypeScript
Perl
Autohotkey
Lua
SQF
R
Scala
Haskell
Groovy
Dart
Clojure
VB.NET
Objective-C
PowerShell
Bash
CoffeeScript
Verilog
function [X] = RK4 (a, t0, T, N, f) h = T / N; p = size(a, 1); X = zeros(p, N); t = t0; X(:, 1) = a; k = h/6; l = h/2; for i=2:N K0 = f(t, X(:, i-1)); K1 = f(t + l, X(:, i-1) + l * K0); K2 = f(t + l, X(:, i-1) + l * K1); K3 = f(t + h, X(:, i-1) + h * K2); X(:, i) = X(:, i-1) + k * ( K0 + 2 * K1 + 2 * K2 + K3 ) t = t + h; end endfunction
Mark as private (unlisted)
for 30 minutes
for 6 hours
for 1 day
for 1 week
for 1 month
for 1 year
everlasting (like CCCP)
