{\rtf1\mac\ansicpg10001\cocoartf824\cocoasubrtf330
{\fonttbl\f0\fswiss\fcharset77 Helvetica;\f1\fnil\fcharset78 HiraKakuPro-W3;\f2\fnil\fcharset77 Monaco;
}
{\colortbl;\red255\green255\blue255;\red106\green117\blue255;}
\paperw11900\paperh16840\margl1440\margr1440\vieww14600\viewh10420\viewkind0
\pard\tx560\tx1120\tx1680\tx2240\tx2800\tx3360\tx3920\tx4480\tx5040\tx5600\tx6160\tx6720\ql\qnatural\pardirnatural

\f0\fs24 \cf0 /*  
\f1 \'98\'41\'97\'a7\'95\'fb\'92\'f6\'8e\'ae\'82\'c9\'82\'c2\'82\'a2\'82\'c4\'82\'cc\'83\'76\'83\'8d\'83\'4f\'83\'89\'83\'80\'82\'f0\'8f\'91\'82\'a2\'82\'c4\'83\'65\'83\'58\'83\'67\'82\'c9\'82\'b7\'82\'e9
\f0 . \
*/\
\pard\tx560\tx1120\tx1680\tx2240\tx2800\tx3360\tx3920\tx4480\tx5040\tx5600\tx6160\tx6720\ql\qnatural\pardirnatural

\f2\fs20 \cf0 def gaussE(A,Y) \{\
  N = size(A)[0];\
  for (K=1; K<=N-1; K++) \{\
    for (I=K+1; I<=N; I++) \{\
      P = A[I-1][K-1]/A[K-1][K-1];\
      for (J=K; J<=N; J++) \{\
        A[I-1][J-1] = A[I-1][J-1]-P*A[K-1][J-1];\
      \}\
      Y[I-1] = Y[I-1]-P*Y[K-1];\
    \} \
    print(A); print(Y); print("----------------");\
  \}\
\}\
\
def gtest1() \{\
  A=newmat(3,3,\cf2 [[1,-1,2],[-1,2,-3],[3,1,1]]\cf0 );\
  Y=newvect(3,[5,-6,8]);\
  print(A); print(Y); print("----------------");\
  gaussE(A,Y);\
  return [A,Y];\
\}\
\
def laplacian(N) \{\
   A = newmat(N,N);\
   for (I=0; I<N; I++) \{\
	 if (I > 0) A[I][I-1] = 1;\
     	A[I][I] = -2;\
          if  (l < N-1) A[I][I+1] = 1;\
   \}\
   return A;\
\}\
[gtest1(), laplacian(3)];\
}