Starting diary file diary-29-Sep-2010.txt

edit fern
finitefern(100000)
x=[.5;.5]
hold on
plot(x(1),x(2),'ro')
A1 = [ .85  .04; -.04  .85];  b1 = [0; 1.6];
A1
b1
plot(x(1),x(2),'ro')

for i=1:1000;x = A1*x + b1;plot(x(1),x(2),'ro');pause;end

for i=1:1000;x = A1*x + b1;plot(x(1),x(2),'ro');end
xtip = x
A1 * xtip + b1   ==    xtip
A1 * xtip + b1     -    xtip

eye(2)
eye(10)
I = eye(2)
A = I - A1
myx = A \ b1

myx  == xtip
myx - xtip
A1 * xtip + b1     -    xtip
A1 * myx + b1 - myx

residual1 = A*xtip - b1
residual2 = A*myx - b1

A*myx - b1
A1 * myx + b1 - myx

plot(myx(1),myx(2),'k*')
clear
clc




load temperature
whos
A = A100;
size(A)
nnz(A)
b= b100;
size(b)
nnz(b)
t = A \ b;
t
max(t)
min(t)
size(t)

T = reshape(t,100,100);
size(T)
clf
surfc(T)

diary off