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

1 view (last 30 days)
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 Comments

Sign in to comment.

Answers (0)

Categories

Find more on Get Started with MATLAB in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!