Next: 3 Octave Up: MATLAB $B Previous: 1 $B5/8;(B

2 $B3P$(=q$-(B

$B9TC10L$NJT=85!G=!"%R%9%H%j!
$BL>A0$NBgJ8;z!">.J8;z$O6hJL$9$k!#(B
clear $BJQ?tL>(B $B$"$k$$$O(B clear who $B$H$9$k$H;HMQ$5$l$F$$$kJQ?tL>$rI=<($9$k(B whos $B$H$9$k$H;HMQ$5$l$F$$$kJQ?tL>$rI=<($9$k(B
$BJQ?tL>(B $B$H$9$k$H!"(B``$BJQ?tL>(B= $BJQ?t$NCM(B'' $B$HI=<($5$l$k!#(B [$BJQ?tL>(B] $B$H$9$k$H!"(B``ans = $BJQ?t$NCM(B'' $B$HI=<($5$l$k!#(B (ans=$BJQ?tL>(B $B$HF1CM$H$$$&$3$H!)(B)
$B4X?t8F$S=P$7(B or $BJQ?t(B = $B4X?t8F$S=P$7(B
$BA0l9g(B ans $B$H$$$&JQ?t$K7k2L$,F~$k!#(B
$B9TKv$K(B ; $B$r$D$1$k$H7k2L$OI=<($5$l$J$$!#(B
Octave $B$O8_49%7%9%F%`(B
$B2#%Y%/%H%k(B $B$O(B [1 2 3] $B$G!"=D%Y%/%H%k(B $\left( \begin{array}{c} 1\ 2\ 3 \end{array}\right)$ $B$O(B [1; 2; 3] $B$GI=$9!#(B
$B9TNs(B $B$O(B
```
	[1 2 3; 4 5 6; 7 8 9]
```
$B$HI=$o$9!#(B
ones(m,n) $B$G(B m $B9T(B n $BNs$N!"[email protected],$,(B 1 $B$G$"$k9TNs(B
zeros(m,n) $B$G(B m $B9T(B n $BNs$N!"[email protected],$,(B 0 $B$G$"$k9TNs(B
eye(n,n) $B$G(B n $B¹
hilb(n), invhilb(n) $B$O$=$l$>$l(B ${\tt n}$ $B $B
hadamard(n) $B$O(B n $B
rand(m,n) $B$G0lMMMp?t9TNs!"(B randn(m,n) $B$G@55,Mp?t9TNs(B
sylvester_matrix(n)
vander(n)
toeplitz(v,n)
2:5 $B$O(B [2 3 4 5], 0:0.2:1 $B$O(B [0.0 0.2 0.4 0.6 0.8 1.0]
$B%Y%/%H%[email protected],$O!"%Y%/%H%k$NJQ?tL>(B($B%$%s%G%C%/%9(B)
+, -, *, / $B$OIaDL$N0UL#$N;MB'(B ($B?t!"%Y%/%H%k!"<0(B)
A' $B$O9TNs$^$?$O%Y%/%H%k(B A $B$N%(%k%_!<%H6&Lr$rI=$9!#(B $BC1$J$kE>CV$O(B A.' $B$GI=$9!#(B
$B=D%Y%/%H%k(B , $B$NFb@Q$O(B
```
	x' * y
```
$B$H$J$k!#(B
$B6e6e$NI=$r:n$k(B (1:9)'*(1:9)
det($B@5J}9TNs(B) $B$O9TNs<0(B
trace()
rank(a, tol)
kron($B9TNs(B,$B9TNs(B) $B$O(B Kronecker $B@Q(B
$BFC@-B?9`<0(B poly($B9TNs(B)
$BB?9`<0$N@Q(B conv($B78?t%Y%/%H%k(B, $B78?t%Y%/%H%k(B)
inv($B9TNs(B) $B$O5U9TNs(B ($B$G$bLGB?$J$3$H$G$O;H$o$J$$$3$H!*(B)
$B9TNs$K$D$$$F$O(B B / A $B$O(B ${\tt B} {\tt A}^{-1}$ , A \ B $B$O(B ${\tt A}^{-1} {\tt B}$ $B$rI=$9!#(B $BO"N)(B$B ${\tt A} x={\tt b}$ $B$N2r(B $x={\tt A}^{-1}{\tt b}$ $B$O(B A \ b $B$G(B $B7W;;$G$-$k!#(B
LU $BJ,2r$9$k$K$O(B lu()
```
	[L,U]=lu(A)
	y=L\b
	x=U\y
```
$B$"$k$$$O%T%\%C%H$D$-(B LU $BJ,2r$O(B
```
	[L,U,P]=lu(A);
```
$B$G$G$-$k(B ( $B$H$J$k(B² , , $B$,5a$a$i$l$k(B)$B!#(B LU $BJ,2r$K0z$-B3$-(B $B$r2r$/$K$O(B
```
       [L,U,P]=lu(A);
       y=L\(P*b);
       x=U\y
```
Cholesky $BJ,2r$r$9$k$K$O(B chol()
```
	u=chol(A)
	l=u'
```
$B$H$9$k$H(B A=l u $B$,@.$jN)$D!#(B
hess(), shur(), svd()
eig($B@5J}9TNs(B) $B$O8GM-CM(B, [v,lambda]=eig($B@5J}9TNs(B) $B$O8GM-%Y%/%H%k!"8GM-CM(B
```
	A=[1 2 3;4 5 6;8 5 2];
	[v,lambda]=eig(A);
	inv(v) * A * v
```
$B9TNsJQ?t(B($B9T;O$^$j(B:$B9T=*$j(B,$BNs;O$^$j(B:$BNs=*$j(B) $B$GItJ,9TNs(B
$B9TNs$NL>A0(B(:,j) $B$H$9$k$HBh(B j $BNs%Y%/%H%k!"(B $B9TNs$NL>A0(B(i,:) $B$H$9$k$HBh(B i $B9T%Y%/%H%k!"(B
$B9TNs$r(B1$B v=A(:)
1$B A=zeros(m,n); A(:)=v; $B$N$h$&$K$9$k!#(BOctave $B$K$O(B reshape() $B$H$$$&4X?t$,$"$j!"(B
```
   A=reshape(v,m,n);
```
$B$N$h$&$K;H$($k!#(B
diag($B%Y%/%H%k(B) $B$OBP3Q9TNs(B ([email protected],$N:n$k%Y%/%H%k$r;XDj(B), $BNc(B $B$($P(B diag([1 2 3]) $B$O(B $\left( \begin{array}{ccc} 1 & 0 & 0 \\ 0 & 2 & 0 \\ 0 & 0 & 3 \end{array}\right)$ $B$rI=$9!#(Bdiag([1 2], 1) $B$O(B $\left( \begin{array}{ccc} 0 & 1 & 0 \\ 0 & 0 & 2 \\ 0 & 0 & 0 \end{array}\right)$ $B$rI=$o$9!#(B
triu($B9TNs(B) $B$O>e;03QItJ,(B
tril($B9TNs(B) $B$O2<;03QItJ,(B
[A B] $B$H$+(B [A, B] $B$OJB$Y$F:n$k(B ($B$=$l$>$l2#!"=D(B)$B!#(B
cond() $B$O>r7o?t(B, rcond() $B$O(B LINPACK $B$N>r7o?t(B ($B@53N$K$O>r7o?t$NAF$$I>2A(B). (Octave $B$K$O(B rcond() $B$O$J$$!#(B)
norm(v,p), norm(v,inf)
norm(x) $B$O(B norm(x,2) $B$KEy$7$$!#(B
norm(x,'fro') $B$O(B Frobenius $B%N%k%`!#(B
lookfor $B$O(B Octave $B$K$O$J$$!#(B
roots($B78?t$rI=$9%Y%/%H%k(B) $B$GB?9`<0$N:,(B
sum($B%Y%/%H%k(B) [email protected],$NOB!"(Bprod($B%Y%/%H%k(B) [email protected],$N@Q!"(B mean($B%Y%/%H%k(B) [email protected],$NJ?6QCM!"(Bmedian($B%Y%/%H%k(B) $B$GCf1{CM!"(B max($B%Y%/%H%k(B) [email protected],$N:GBgCM!"(Bmin($B%Y%/%H%k(B) [email protected],$N:G>.CM!"(B var($B%Y%/%H%k(B) $B$GJ,;6(B ($B$?$@$7MWAG?t(B $B!`(B (n-1) $B$G7W;;$9$k$d$D(B)$B!"(B cov($B%Y%/%H%k(B) $B$GJ,;6(B ($B$?$@$7MWAG?t(B $B!`(B (n-1) $B$G7W;;$9$k$d$D(B)$B!"(B std($B%Y%/%H%k(B) $B$GI8=`JP:9(B ($BJ,;6$N@5$NJ?J}:,(B)$B!"(B length($B%Y%/%H%k(B) $B$G%Y%/%H%k$H$7$F$N
$B$$$o$f$kFC
plot($B%Y%/%H%k(B), closeplot

$B@)8f9=J8$O(B

  while $B>r7o<0(B
    $BJ8(B1
    $BJ8(B2
    ...
    $BJ8(Bn
  end

  if $B>r7o<0(B
    $BJ8(B1
    ..
    $BJ8(Bn
  else
    $BJ8(B'1
    ...
    $BJ8(B'n
  end

  if $B>r7o<0(B1
    $BJ8(B
  elseif $B>r7o<0(B2
    $BJ8(B
  end

  for i=1:n
   $BJ8(B1
   ...
   $BJ8(Bn
  end

for $BJQ?t(B=$B=i4|CM(B:$B=*N;CM(B $B$^$?$O(B for $BJQ?t(B=$B=i4|CM(B:$BA}J,(B:$B=*N;CM(B
break $B$O%k!<%W$rH4$1=P$9(B
$B!V$+$D!W$O(B &, $B!V$^$?$O!W(B $B$O(B |, $B!VEy$7$/$J$$!W$O(B ~= $B$^$?$O(B != $B$GI=$o$9!#(B $B$D$^$jH]Dj$O(B ~ $B$^$?$O(B ! $B$H$$$&$3$H$@$m$&$J!#(B
$B%-!<%\!<%IF~NO(B
```
	$BJQ?tL>(B = input('$B2?$+F~NO$7$F2<$5$$(B')
```
$B!z(BOctave for FreeBSD $B$OF|K\8l$,F~NO$G$-$J$$!#(B
pause $B$^$?$O(B pause($BIC?t(B)
$B1i;;;R$NA0$K(B . [email protected],$4$H$H$$$&0UL#!)(B c=[1 2 3] $B$KBP$7$F(B c*c $B$O7?$,$*$+$7$$$+$i7W;;$G$-$J$$!#(B c .* c $B$H$9$k$H(B [1 4 9]
$BJQ?tL>(B=$BCM(B $B%Y%/%H%k$N>l9g$O(B [$BJQ?tL>(B1,$BJQ?tL>(B2]=$BCM(B, [$BJQ?tL>(B1 $BJQ?tL>(B2]=$BCM(B

$B#1JQ?t4X?t$N%0%i%U$rIA$/$K$O!"%Y%/%H%k$r(B2$BK\MQ0U$9$k(B

	% y=x^2 $B$N%0%i%U(B
	x=0:0.1:10;
	y=x.^2;
	plot(x,y)

	% y=sin(x) $B$N%0%i%U(B
	x=0:0.1:10;
	y=sin(x);
	plot(x,y)

gset term postscript; gset output "foo.ps"; replot

contour(z,$BEy9b@~$N%l%Y%k?t(B,x,y)

  n=10; x=(-1:1/n:1)'; y=(-1:1/n:1); z=x*y; contour(z,10,x,y)

mesh(x,y,z) $B$O(B3$B
i, I, j, J $B$O5u?tC10L!"(Bpi $B$O1_<~N(!"(B e $B$O<+A3BP?t$NDl(B (Napier $B$N?t(B)$B!"(B inf $B$OL58BBg!"(Beps $B$O7W;;5!%$%W%7%m%s(B
$B$O(B 3+5i $B$GI=$;$k!#(B5 $B$H(B i $B$N4V$K6uGr$rCV$$$F$O(B $B$$$1$J$$!#(B
real(), imag(), conj(), arg(), angle() $B$J$I$,;H$($k!#(B

eval($B<0(B) $B$O<0$NI>2A(B

	t='1/(i+j-1)'
	for i=1:n
		for j=1:n
			a(i,j)=eval(t);
		end
	end

$B2hLL=PNO$N@:EY$O(B

format long

format short

format long e $BD9$$7e?t!";X?t7A<0(B

format short e $BC;$$7e?t!";X?t7A<0(B

format free $B<+M37A<0(B

format none $B<+M37A<0(B
quad('$B4X?tL>(B',$B=i4|CM(B,$B=*C quad8() $B$O(B Octave $B$K$O$J$$!#(B

save -ascii $B%U%!%$%kL>(B $BJQ?tL>(B1 $BJQ?tL>(B2, save -binary $B%U%!%$%kL>(B $BJQ?tL>(B1 $BJQ?tL>(B2, save -mat-binary $B%U%!%$%kL>(B $BJQ?tL>(B1 $BJQ?tL>(B2

Next: 3 Octave Up: MATLAB $B Previous: 1 $B5/8;(B Masashi Katsurada $BJ?@.(B15$BG/(B5$B7n(B8$BF|(B

`format long`
`format short`
`format long e`	$BD9$$7e?t!";X?t7A<0(B
`format short e`	$BC;$$7e?t!";X?t7A<0(B
`format free`	$B<+M37A<0(B
`format none`	$B<+M37A<0(B

2 $B3P$(=q$-(B

2 $B3P$(=q$-(B