load temperature whos Name Size Bytes Class Attributes A100 10000x10000 635204 double sparse A4 16x16 2048 double ans 1x56 112 char b100 10000x1 80000 double b4 16x1 128 double A4 A4 = Columns 1 through 14 4 -1 0 0 -1 0 0 0 0 0 0 0 0 0 -1 4 -1 0 0 -1 0 0 0 0 0 0 0 0 0 -1 4 -1 0 0 -1 0 0 0 0 0 0 0 0 0 -1 4 0 0 0 -1 0 0 0 0 0 0 -1 0 0 0 4 -1 0 0 -1 0 0 0 0 0 0 -1 0 0 -1 4 -1 0 0 -1 0 0 0 0 0 0 -1 0 0 -1 4 -1 0 0 -1 0 0 0 0 0 0 -1 0 0 -1 4 0 0 0 -1 0 0 0 0 0 0 -1 0 0 0 4 -1 0 0 -1 0 0 0 0 0 0 -1 0 0 -1 4 -1 0 0 -1 0 0 0 0 0 0 -1 0 0 -1 4 -1 0 0 0 0 0 0 0 0 0 -1 0 0 -1 4 0 0 0 0 0 0 0 0 0 0 -1 0 0 0 4 -1 0 0 0 0 0 0 0 0 0 -1 0 0 -1 4 0 0 0 0 0 0 0 0 0 0 -1 0 0 -1 0 0 0 0 0 0 0 0 0 0 0 -1 0 0 Columns 15 through 16 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -1 0 0 -1 0 0 -1 0 4 -1 -1 4 b4 b4 = 100 100 100 100 0 0 0 0 0 0 0 0 0 0 0 0 spy(A4) A4(7,:) ans = Columns 1 through 14 0 0 -1 0 0 -1 4 -1 0 0 -1 0 0 0 Columns 15 through 16 0 0 t4 = A4 \ b4 t4 = 45.4545 59.4697 59.4697 45.4545 22.3485 32.9545 32.9545 22.3485 10.9848 17.0455 17.0455 10.9848 4.5455 7.1970 7.1970 4.5455 T = reshape(t4,4,4) T = 45.4545 22.3485 10.9848 4.5455 59.4697 32.9545 17.0455 7.1970 59.4697 32.9545 17.0455 7.1970 45.4545 22.3485 10.9848 4.5455 surfc(T) size(A100) ans = 10000 10000 spy(A100) A100(777,:) ans = (1,677) -1 (1,776) -1 (1,777) 4 (1,778) -1 (1,877) -1 t100 = A100 \ b100; T = reshape(t100,100,100); surfc(T) edit jacobi.m jacobi jacobi(1000); jacobi step = 1 step = 2 step = 3 step = 4 step = 5 step = 6 step = 7 step = 8 step = 9 step = 10 step = 11 step = 12 step = 13 step = 14 step = 15 step = 16 step = 17 step = 18 step = 19 step = 20 step = 21 step = 22 step = 23 step = 24 step = 25 step = 26 step = 27 step = 28 step = 29 jacobi(1000) figure,surfc(T) diary off