Skip to main content

integrate_euler

Function integrate_euler 

Source
pub fn integrate_euler<Field, U, T>(
    rate: impl FnMut(Quantity<T>, &Field::Point) -> Result<Derivative<Field::Increment, T>, String>,
    time: &[Quantity<T>],
    initial_condition: Field::Point,
) -> Result<(Times<T>, U), IntegrationError>
where Field: Integrable, Field::Point: Clone, Field::Increment: Differentiable<T>, for<'a> &'a Derivative<Field::Increment, T>: Mul<Quantity<T>, Output = Field::Increment>, U: TensorVec<Item = Field::Point>,
Expand description

Explicit Euler for a single Integrable, one step per interval of time.

\mathbf{x}_{n+1} = \mathrm{reconstruct}\!\left(\mathbf{x}_n,\ h\,\mathbf{f}(t_n, \mathbf{x}_n)\right)