Dies ist der Test für die Gleichung Differential Runge-kutta45 Verwendung: f (x, y) = (-5 * x - y/5)^1/8 + 10Warum ist das numerische Ergebnis anders (RK45)?
warum die numerischen Ergebnis ist anders? Früher habe ich:
function Rk_JL()
f(x,y)= (-5*x - y/5)^1/8 + 10
tspan = 0:0.001:n
y0 = [0.0, 1.0]
return ODE.ode45(f, y0,tspan);
end
und
function [X1,Y1] = RK_M()
f = @(x,y) (-5*x - y/5)^1/8 + 10;
tspan = 0:0.001:n;
y0 = 1
[X1,Y1]= ode45(f,tspan,1);
end
Wissen Sie, wie viel Toleranz in beiden Fällen angenommen wird? – Jagte
cross-gepostet zu https://groups.google.com/forum/#!topic/julia-users/tNbuR8xOGXo. –
tatsächlich gibt es eine Beziehung mit der Toleranz zwischen den Methoden. –