% Samples for institute day clc format compact A = [1 2 3; 3 0 1; 4 5 6] det(A) inv(A) A^2 n=[1:20] format long g x=1./2.^(n); y1=sin(x)./x; y2=(1-cos(x))./x.^2; [x;y1;y2]' x1=linspace(-2,2); y3=x1.^2; plot(x1,y3,[1],[1],'o',[0 2],[-1 3])