Viewing All Docs
Title |
---|
HW20
% problem 3.5.5 a=[0 3 -3 0;2 0 -3 2]; syms t; T=[1 t t^2 t^3]; B1=a*T.’; B2=a*T.’-[0;1]; set(groot,’defaultAxesColorOrder’,[0 0 […] |
HW21
%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 […] |