% Matlab diary from cs110a class #9, Jan 24 2005.
%
% This is what I typed in class.  
%
% Type it into Matlab yourself to
% reproduce the output I showed.

b = [ones(4,1)  ; zeros(12,1) ]
size(b)
load temperature
A4
b
b = 100*b
t = A4 \ b
T = reshape(t,4,4)
T = T'
b = [ 100 * ones(100,1) ; zeros(9900,1) ]
size(b)
load temperature
size(A100)
t = A100 \ b
size(t)
T = reshape(t,100,100)'
max(t)
min(t)
size(T)
surf(T)
clc
help surf
help surfc
surfc(T)
size(A100)
A=A100;
prod(size(A))
nnz(A)
whos
spy(A4)
spy(A)
R = chol(A)
clc
A4
A100
clc
A4
R=chol(A4)
clc
R' * R - A4
spy(R)
clc
spy(A)
R = chol(A);
spy(R)
nnz(A)
nnz(R)
10000 * 10000 / 2
clc