AMATH 301 Lecture Notes - Lecture 21: Euler Method, Matlab

73 views3 pages
27 Feb 2018
School
Course

Document Summary

% x(0) = 1; v(0) = v0 x0 = 1; v0 = 0; y0 = [x0; v0]; % v" = -sin(x) - v = -sin(y1) -y2 f = @(t,y)([y(2); -sin(y(1)) - y(2)]); [t_out, y_out] = ode45(f, [0 t], y0); x = y_out(:,1); plot(t_out, x, "k") t = 0:dt:t; n = length(t); y = zeros(2,n); y(:,1) = y0; % forward euler method for k = 1:n-1 y(:,k+1) = y(:,k) + dt*f(t(k), y(:,k)); end x = y(1,:); hold on plot(t, x, "r") Lorenz equations x0 = 1; y0 = 0; z0 = 2; v0 = [x0, y0, z0]; T = 200; sigma = 10; rho = 28; beta = 8/3; f = @(t, v)([sigma * (v(2) - v(1)); v(1) * (rho - v(3)) - v(2); v(1) * v(2) - beta * v(3)]); [t_out, v_out] = ode45(f, [0 t], v0); x = v_out(:,1); y = v_out(:,2); z = v_out(:,3); plot3(x, y, z, "k") view(-4, 16)

Get access

Grade+20% off
$8 USD/m$10 USD/m
Billed $96 USD annually
Grade+
Homework Help
Study Guides
Textbook Solutions
Class Notes
Textbook Notes
Booster Class
40 Verified Answers
Class+
$8 USD/m
Billed $96 USD annually
Class+
Homework Help
Study Guides
Textbook Solutions
Class Notes
Textbook Notes
Booster Class
30 Verified Answers

Related textbook solutions

Related Documents

Related Questions