Backward Euler with LU Decomposition
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
The code below solves forward Euler for a pendelum problem. How would I modify this code to solve with backward Euler and LU Decomposition?
length= 5;
g=9.8 ;
npoints = 5000;
dt = 0.01;
omega = zeros(npoints,1);
theta = zeros(npoints,1);
time = zeros(npoints,1);
theta(1)=0.5;
for step = 1:npoints-1
omega(step+1) = omega(step) - (g/length)*theta(step)*dt;
theta(step+1) = theta(step)+omega(step)*dt;
time(step+1) = time(step) + dt;
end
forward=[theta(end);omega(end)];
0 Kommentare
Antworten (0)
Siehe auch
Kategorien
Mehr zu Matrix Decomposition finden Sie in Help Center und File Exchange
Produkte
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!