-->A=[3 3 0 0; 2 2 0 1; 3.3 5.5 7 %pi; 5 8 0 14]
 A  =
 
    3.     3.     0.    0.         
    2.     2.     0.    1.         
    3.3    5.5    7.    3.1415927  
    5.     8.     0.    14.        
 
-->det(A)
 ans  =
 
  - 63.  
 
-->E=[
-->    2.1756236   -0.6122359
-->   -0.6122359    1.6621008
-->]
 E  =
 
    2.1756236  - 0.6122359  
  - 0.6122359    1.6621008  
 
-->A=E; X=[1;0]
 X  =
 
    1.  
    0.  
 
-->exec iterate2d.in;
N = 2.260126, X = [0.962612; -0.270886], i = 1
N = 2.487751, X = [0.908502; -0.417881], i = 2
N = 2.558916, X = [0.872400; -0.488792], i = 3
N = 2.577045, X = [0.852632; -0.522512], i = 4
N = 2.581405, X = [0.842528; -0.538652], i = 5
N = 2.582439, X = [0.837505; -0.546429], i = 6
N = 2.582684, X = [0.835038; -0.550192], i = 7
N = 2.582741, X = [0.833834; -0.552016], i = 8
N = 2.582755, X = [0.833247; -0.552901], i = 9
N = 2.582758, X = [0.832961; -0.553332], i = 10
 
-->exec iterate2d.in;
N = 2.582759, X = [0.832822; -0.553540], i = 1
N = 2.582759, X = [0.832755; -0.553642], i = 2
N = 2.582759, X = [0.832722; -0.553691], i = 3
N = 2.582759, X = [0.832706; -0.553715], i = 4
N = 2.582759, X = [0.832698; -0.553727], i = 5
N = 2.582759, X = [0.832695; -0.553732], i = 6
N = 2.582759, X = [0.832693; -0.553735], i = 7
N = 2.582759, X = [0.832692; -0.553737], i = 8
N = 2.582759, X = [0.832691; -0.553737], i = 9
N = 2.582759, X = [0.832691; -0.553738], i = 10
 
-->exec iterate2d.in;
N = 2.582759, X = [0.832691; -0.553738], i = 1
N = 2.582759, X = [0.832691; -0.553738], i = 2
N = 2.582759, X = [0.832691; -0.553738], i = 3
N = 2.582759, X = [0.832691; -0.553738], i = 4
N = 2.582759, X = [0.832691; -0.553738], i = 5
N = 2.582759, X = [0.832691; -0.553738], i = 6
N = 2.582759, X = [0.832691; -0.553738], i = 7
N = 2.582759, X = [0.832691; -0.553738], i = 8
N = 2.582759, X = [0.832691; -0.553738], i = 9
N = 2.582759, X = [0.832691; -0.553738], i = 10
 
-->Emax=X; Lmax=N;
 
-->E*X ./ X
 ans  =
 
    2.5827592  
    2.5827592  
 
-->A=inv(E); X=[1;0]
 X  =
 
    1.  
    0.  
 
-->exec iterate2d.in;
N = 0.546475, X = [0.938365; 0.345647], i = 1
N = 0.682731, X = [0.800424; 0.599434], i = 2
N = 0.762004, X = [0.687237; 0.726433], i = 3
N = 0.787990, X = [0.621359; 0.783526], i = 4
N = 0.794708, X = [0.587167; 0.809466], i = 5
N = 0.796331, X = [0.570106; 0.821571], i = 6
N = 0.796716, X = [0.561719; 0.827328], i = 7
N = 0.796807, X = [0.557623; 0.830095], i = 8
N = 0.796828, X = [0.555627; 0.831432], i = 9
N = 0.796833, X = [0.554656; 0.832080], i = 10
 
-->exec iterate2d.in;
N = 0.796834, X = [0.554184; 0.832394], i = 1
N = 0.796835, X = [0.553955; 0.832547], i = 2
N = 0.796835, X = [0.553843; 0.832621], i = 3
N = 0.796835, X = [0.553789; 0.832657], i = 4
N = 0.796835, X = [0.553763; 0.832675], i = 5
N = 0.796835, X = [0.553750; 0.832683], i = 6
N = 0.796835, X = [0.553744; 0.832687], i = 7
N = 0.796835, X = [0.553741; 0.832689], i = 8
N = 0.796835, X = [0.553739; 0.832690], i = 9
N = 0.796835, X = [0.553738; 0.832691], i = 10
 
-->exec iterate2d.in;
N = 0.796835, X = [0.553738; 0.832691], i = 1
N = 0.796835, X = [0.553738; 0.832691], i = 2
N = 0.796835, X = [0.553738; 0.832691], i = 3
N = 0.796835, X = [0.553738; 0.832691], i = 4
N = 0.796835, X = [0.553738; 0.832691], i = 5
N = 0.796835, X = [0.553738; 0.832691], i = 6
N = 0.796835, X = [0.553738; 0.832691], i = 7
N = 0.796835, X = [0.553738; 0.832691], i = 8
N = 0.796835, X = [0.553738; 0.832691], i = 9
N = 0.796835, X = [0.553738; 0.832691], i = 10
 
-->Emin=X; Lmin=1/N
 Lmin  =
 
    1.2549652  
 
-->E*X ./ X
 ans  =
 
    1.2549652  
    1.2549652  
 
-->Emax, Emin
 Emax  =
 
    0.8326911  
  - 0.5537378  
 Emin  =
 
    0.5537378  
    0.8326911  
 
-->Emax*Emin
       !--error 10 
Inconsistent multiplication.

 
 
-->Emax'*Emin
 ans  =
 
    8.563D-10  
 
-->help xarrows
 
-->exec arrowdemo.in
 
-->// Some xarrow stuff
 
-->// and graphics demo
 
-->scf(0); // open graphics window 0 (if it is not already open)
 
-->clf(); // clear it
 
-->t=0:0.1:1;
 
-->y=t.*(1-t);
 
-->plot(t,y,'b*-')
 
-->xtitle("y = t(1-t) in blue and matching arrow","t","y");
 
-->xarrows([0.1,0.2],[0.2,0.1],1,2) // the 2 is blue, the 1 is the size of arrowhead
 
-->//
 
-->scf(1); // open graphics window 1 (if it is not already open)
 
-->clf(); // clear it
 
-->plot(t,y,'rs-')
 
-->xtitle("y = t(1-t) in red and matching arrow","t","y");
 
-->xarrows([0,0.5],[0,0.25],1,5) // the 5 is red, the 1 is the size of arrowhead
 
-->//
 
-->scf(2); // open graphics window 2 (if it is not already open)
 
-->clf(); // clear it
 
-->plot(t,y,'gd-')
 
-->xtitle("y = t(1-t) in green and matching arrow","t","y");
 
-->xarrows([0.5,0],[0.25,0],1,3) // the 3 is green? the 1 is the size of arrowhead
 
 
-->E=[
-->    2.1756236   -0.6122359
-->   -0.6122359    1.6621008
-->]
 E  =
 
    2.1756236  - 0.6122359  
  - 0.6122359    1.6621008  
 
-->function y = squareroot(x)
-->y = sqrt(x);
-->end
-->endfunction
end
 !--error 34 
Incorrect control instruction syntax.

at line       3 of function squareroot called by :  
endfunction
 
 
-->function y = squareroot(x)
-->y = sqrt(x);
-->endfunction