MAT2630 Numerical Methods, FA2015

You are currently viewing a revision titled "HW21", saved on November 18, 2015 at 11:10 am by Ezra Halleck
Title
HW21
Content
%prob 5.1.1 df=inline('log(1+h)/h') out=[]; for i=1:16 out=[out;i 10^-i df(10^-i) df(10^-i)-1]; end format long disp(' i h=10^-i df(h) error') disp(out) %prob 5.1.2 df=inline('(exp(h)-exp(-h))/(2*h)') out=[]; for i=1:12 out=[out;i 10^-i df(10^-i) df(10^-i)-1]; end format long disp(' i h=10^-i df(h) error') disp(out) %prob 5.1.3 df=inline('(sin(pi()/3+h)-sqrt(3)/2)/h') out=[]; for i=1:16 h=10^-i; out=[out;h df(h) -h*sqrt(3)/4 df(h)-1/2 -h*sin(pi()/3+h)/2]; end format long disp(' h=10^-i df(h) -h*sqrt(3)/4 error -h*sin(pi()/3+h)/2') disp(out)
Excerpt
Footnotes


Old New Date Created Author Actions
November 18, 2015 at 4:10 pm Ezra Halleck