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

2 vues (au cours des 30 derniers jours)
OZGUR YILDIRIM
OZGUR YILDIRIM le 28 Déc 2020
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 commentaires
Asvin Kumar
Asvin Kumar le 5 Jan 2021
https://www.mathworks.com/help/symbolic/subs.html Should do just that.

Connectez-vous pour commenter.

Réponses (0)

Catégories

En savoir plus sur Symbolic Math Toolbox dans Help Center et File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by