Finite difference scheme and iteration what is wrong with my code? Please help me

2 次查看（过去 30 天）

OZGUR YILDIRIM 2020-12-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/703777-finite-difference-scheme-and-iteration-what-is-wrong-with-my-code-please-help-me

评论： OZGUR YILDIRIM 2021-1-5

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