function [ x ] = evalImplicitEuler( lambda, x0, n ) %evalImplicitEuler Evaluate the implicit Euler for dx/dt = lambda x % Detailed explanation goes here x = zeros(1, n+1); x(1) = x0; tau = 1.0 / n; for k = 1:n x(k+1) = x(k) / (1.0 - tau * lambda); end end