Starting diary file "diary-2016-01-04.txt" on Mon Jan 04 2016  1:05 PM.

1+2
ans = 3

rand(5)
ans =
    0.8147    0.0975    0.1576    0.1419    0.6557
    0.9058    0.2785    0.9706    0.4218    0.0357
    0.1270    0.5469    0.9572    0.9157    0.8491
    0.9134    0.9575    0.4854    0.7922    0.9340
    0.6324    0.9649    0.8003    0.9595    0.6787

pi
ans = 3.1416

for i=1:10, sqrt(i), end
ans = 1
ans = 1.4142
ans = 1.7321
ans = 2
ans = 2.2361
ans = 2.4495
ans = 2.6458
ans = 2.8284
ans = 3
ans = 3.1623

myfern
fastfern
edit myfern

p  = [ .85  .92  .99  1.00];
A1 = [ .85  .04; -.04  .85];  b1 = [0; 1.6];
A2 = [ .20 -.26;  .23  .22];  b2 = [0; 1.6];
A3 = [-.15  .28;  .26  .24];  b3 = [0; .44];
A4 = [  0    0 ;   0   .16];

x  = [.5 ; .5];  

p
p =
    0.8500    0.9200    0.9900    1.0000
A1
A1 =
    0.8500    0.0400
   -0.0400    0.8500

b1
b1 =
         0
    1.6000
clc

A1
A1 =
    0.8500    0.0400
   -0.0400    0.8500

b1
b1 =
         0
    1.6000

x = [.5;.5]
x =
    0.5000
    0.5000

x = A1*x+b1
x =
    0.4450
    2.0050

x = A1*x+b1
x =
    0.4584
    3.2864

x = A1*x+b1
x =
    0.5211
    4.3751

x = A1*x+b1
x =
    0.6180
    5.2980

x = A1*x+b1
x =
    0.7372
    6.0786

x = A1*x+b1
x =
    0.8698
    6.7373

x = A1*x+b1
x =
    1.0088
    7.2919

x = A1*x+b1
x =
    1.1492
    7.7578

I = eye(2)
I =
     1     0
     0     1

A = I - A1
A =
    0.1500   -0.0400
    0.0400    0.1500

b1
b1 =
         0
    1.6000

x = A \ b1;
x
x =
    2.6556
    9.9585

clc
clear
load temperature
whos
  Name          Size                Bytes  Class     Attributes

  A100      10000x10000            873608  double    sparse    
  A4           16x16                 2048  double              
  b100      10000x1                 80000  double              
  b4           16x1                   128  double              

b100
b100 = ...

b100(50)
ans = 100

t = A100 \ b100;

t
t =
    0.0588
    0.1181
    0.1784
    0.2402
    ...
    0.0188
    0.0125
    0.0063

T = reshape(t,100,100);
size(T)
ans =
   100   100

surfc(T)
shg

size(A100)
ans =
       10000       10000

nnz(A100)
ans =
       49600

spy(A100) ; shg

size(b100)
ans =
       10000           1

plot(b100);shg

diary off