Finite difference scheme and iteration what is wrong with my code? Please help me
조회 수: 3 (최근 30일)
이전 댓글 표시
function Perturbation(x,N,epsilon)
clc
close all
x=10;
N=10;
M=10;
h=1/M;
alpha=1;
beta=-1;
epsilon=0.01;
dtau=1/N; % step size
tau = 0:dtau:4; % Calculates upto y(3)
w = sym(zeros(1,length(tau)));
t= sym(zeros(1,length(tau))); % initial condition
F_0 = @(r) r; % change the function as you desire
F_10 = @(r) 2*r;
F_11 = @(r) 3*r;
for i=1:(length(h)-1) % calculation loop
tk=i*tau;
t=tk;
w(0) = F_0(t(i));
w(1) = F_10(t(i))+F_11(t(i));
w(i+1) = w(i) + w(i-1)+x % main equation
end
댓글 수: 4
답변 (0개)
참고 항목
카테고리
Help Center 및 File Exchange에서 Symbolic Math Toolbox에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!