Runge-Kutta

Runge-Kutta

Exato

Ponto Médio

Euler

Y_{0} = 1 (slider from -5 to 5)
a = 1 (slider from -5 to 5)
T: 1 / a
h = 1.2 (slider from 0 to 1.2)
f(x) = (Y_{0} * ℯ^(((-a) * x)))
A = Point(xAxis)
B = Point(Circle[A, h / 2])
i = Segment(A, B)
C = Point(Circle[B, h / 2])
j = Segment(B, C)
k = OrthogonalLine(A, xAxis)
l = OrthogonalLine(B, xAxis)
m = OrthogonalLine(C, xAxis)
D = Intersect(f, k, (0.1, 0.9048374180359596))
n = Tangent(D, f)
E = Intersect(n, l)
p = Line(E, xAxis)
F = Intersect(f, p, (1.016290732580936, 0.3619349669585751))
q = Tangent(F, f)
r = Line(D, q)
G = Intersect(r, l)
s = Line(G, xAxis)
H = Intersect(f, s, (0.37443684169925545, 0.6876764404597576))
t = Tangent(H, f)
b = Line(D, t)
I = Intersect(b, m)
c = Line(I, xAxis)
J = Intersect(f, c, (2.5304185059762347, 0.07962568948490358))
d = Tangent(J, f)
a_1 = Slope(n)
a_2 = Slope(t)
a_3 = Slope(q)
a_4 = Slope(d)
e: (1 / 6 * (a_1 + (2 * a_2) + (2 * a_3) + a_4))
c_1 = Circle(C, Distance[A, D] + (e * h))
K, L = Intersect(c_1, m)
f_1 = Segment(D, L)
M = Intersect(f, m, (1.3, 0.2725317930340126))
N = Intersect(m, r)
O = Intersect(n, m)