Sudal's Garage

Implicit Euler Method 를 구현한다 function [T,Y] = odeIEuler(odefun,t,y0) Y = zeros(length(t),1); Y(1) = y0; T = t; for n = 1:length(t) - 1 h = T(n+1) - T(n); Y(n+1) = fzero(@(y) y - Y(n) - odefun(T(n+1),y) * h,0); end end Forward, Backward and central difference formulas 를 이용해 Local Truncation Error 를 구현하는 script를 구현%% (a) clear; close; fun = @(x) sin(sqrt(x)); x = 2; k = 0:1:5; fw = zeros(length(k)..

2차 미분방정식을 MATLAB 으로 바꾸어푸는 문제. 2차 미분방정식을 2개의 1차 미분방정식으로 변환시켰다! MATLAB 의 ode45 함수는 1차 미분방정식만 풀 수 있다. function dudt = odefun(t,u) dudt = [0;0]; dudt(1) = u(2); % u1' = u2 dudt(2) = -u(1)^n; % u2' = -u1^n end odefun 에 두개의 1차 미분방정식을 넣고 코딩을 한다. Solution:clear; close; figure(1) for n = [1 5 9 15] [T,Y] = draw_oscillator(n); plot(T,Y) hold on end legend({'n = 1','n = 5','n = 9','n = 14'},'NumColumns',2)..