euler_convergence_and_bias.sci
clc clear a=10.0 y0 = 1.0 h = 1.1*2.0/a //0.01, 0.2, 1.1 L = 1.0 function y=an(x) y = y0*exp(-a*x) endfunction function yn=f(x) yn = zeros(1,length(x)) for i=0:length(x)-1 yn(i+1)=(1.0-a*h)^i*y0 end endfunction x = 0:0.01:L; plot(x, an(x),'r') x = 0:h:L; plot(x, f(x),'b')