Script Programming(스크립트 프로그래밍)

필요한데가 많을 것 같아서 올려 놓는다.

실행 결과:

1.100000 : Floating point
0 : Integer
0.000000 : Floating point
10 : Integer
Something else : Not matched.

F# 프로그래밍 중에 한글 출력을 시도하였다.
그러나 당연히 되리라고 믿었던 printfn 에 의한 한글 출력이 되지 않았다.

printfn "한글"

이 한 줄 짜리 F# 소스를 컴파일하고 실행하면

???

가 출력된다. 소스를 UTF8 인코딩으로 저장해 보기도 하고, 명령창을 chcp 65001 하여 UTF8 창으로 바꾸어 보기도 하였지만 허사였다.
인터넷 검색을 해보아도 해결책을 찾지 못했다.
아마도 F# 언어를 사용하는 사람들은 적어도 한글을 안써도 되는 듯...
(.NET에서 한글을 못 쓰면 문제가 많을텐데...)

그러나 알고 보면 쉬운 것.
fsc 로 컴파일할 때  --codepage:숫자 옵션으로 인코딩 숫자(예를 들어 949)를 지정하면 된다.

우선 F# 소스를 하나 보자

컴파일하기:
D:\test\f#>fsc --codepage:949 hello.fs
Microsoft (R) F# 2.0 Compiler build 2.0.0.0
Copyright (c) Microsoft Corporation. All Rights Reserved.

실행하기:
D:\test\f#>hello
Hello
안녕하세요?
F# 프로그래밍에서 한글 출력을 연습합니당...

* 참고: fsi 에서는 한글 입력과 출력이 잘 되지 않는다.


* Visual Studio 통합 환경에서는 한글 문제가 발생하지 않는다.

이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-변경금지 2.0 대한민국 라이센스에 따라 이용하실 수 있습니다.

목적:
http://stackoverflow.com/questions/1907709/a-simple-wrapper-for-f-to-do-matrix-operations 에 올려진 소스에서 String.spit 관련를 문제 해결하자.

테스트를 위해 필요한 자료:   dnAnalytics
                                        F# Power Pack


// 파일명: TestFSharpMatrix.fs (원래 소스)

#r @"D:\WORK\tools\dnAnalytics_windows_x86\bin\dnAnalytics.dll" 
#r @"FSharp.PowerPack.dll" 
 
open dnAnalytics.LinearAlgebra 
open Microsoft.FSharp.Math 
open dnAnalytics.LinearAlgebra.Decomposition 
 
// F# matrix -> ndAnalytics DenseMatrix 
let mat2dnmat (mat:matrix) =  
    let m = new DenseMatrix(mat.ToArray2D()) 
    m 
 
// ndAnalytics DenseMatrix -> F# matrix  
let dnmat2mat (dnmat:DenseMatrix) =  
    let n = dnmat.Rows 
    let m = dnmat.Columns 
    let mat = Matrix.create n m 0. 
    for i=0 to n-1 do 
        for j=0 to m-1 do 
            mat.[i,j] <- dnmat.Item(i,j) 
    mat 
 
// random matrix 
let randmat n m = 
    let r = new System.Random() 
    let ranlist m =  
        [ for i in 1..m do yield r.NextDouble() ] 
    matrix ([1..n] |> List.map (fun x-> ranlist m)) 
 
// is square matrix 
let issqr (m:matrix) = 
    let n, m = m.Dimensions 
    n = m 
 
// is postive definite  
let ispd m = 
    if not (issqr m) then false 
    else  
        let m1 = mat2dnmat m 
        let qrsolver = dnAnalytics.LinearAlgebra.Decomposition.Cholesky(m1) 
        qrsolver.IsPositiveDefinite() 
 
// determinant 
let det m = 
    let m1 = mat2dnmat m 
    let lusolver = dnAnalytics.LinearAlgebra.Decomposition.LU(m1) 
    lusolver.Determinant () 
 
// is full rank 
let isfull m = 
    let m1 = mat2dnmat m 
    let qrsolver = dnAnalytics.LinearAlgebra.Decomposition.GramSchmidt(m1) 
    qrsolver.IsFullRank() 
 
// rank 
let rank m = 
    let m1 = mat2dnmat m 
    let svdsolver = dnAnalytics.LinearAlgebra.Decomposition.Svd(m1, false) 
    svdsolver.Rank() 
 
// inversion by lu 
let inv m = 
    let m1 = mat2dnmat m 
    let lusolver = dnAnalytics.LinearAlgebra.Decomposition.LU(m1) 
    lusolver.Inverse() 
 
// lu 
let lu m = 
    let m1 = mat2dnmat m 
    let lusolver = dnAnalytics.LinearAlgebra.Decomposition.LU(m1) 
    let l = dnmat2mat (DenseMatrix (lusolver.LowerFactor ())) 
    let u = dnmat2mat (DenseMatrix (lusolver.UpperFactor ())) 
    (l,u) 
 
// qr  
let qr m = 
    let m1 = mat2dnmat m 
    let qrsolver = dnAnalytics.LinearAlgebra.Decomposition.GramSchmidt(m1) 
    let q = dnmat2mat (DenseMatrix (qrsolver.Q())) 
    let r = dnmat2mat (DenseMatrix (qrsolver.R())) 
    (q, r) 
 
// svd  
let svd m = 
    let m1 = mat2dnmat m 
    let svdsolver = dnAnalytics.LinearAlgebra.Decomposition.Svd(m1, true) 
    let u = dnmat2mat (DenseMatrix (svdsolver.U())) 
    let w = dnmat2mat (DenseMatrix  (svdsolver.W())) 
    let vt = dnmat2mat (DenseMatrix (svdsolver.VT())) 
    (u, w, vt.Transpose) 


// and test code
(* todo: read matrix market format   ref: http://math.nist.gov/MatrixMarket/formats.html *) 
let readmat (filename:string) =  
    System.IO.File.ReadAllLines(filename) |> Array.map (fun x-> (x |> String.split [' '] |> List.toArray |> Array.map float)) 
    |> matrix 
 
let timeit f str=  
    let watch = new System.Diagnostics.Stopwatch() 
    watch.Start() 
    let res = f() 
    watch.Stop() 
    printfn "%s Needed %f ms" str watch.Elapsed.TotalMilliseconds 
    res 
 
let test() =  
    let testlu() =  
        for i=1 to 10 do 
            let a,b = lu (randmat 1000 1000) 
            () 
        () 

    let testsvd() = 
        for i=1 to 10 do 
            let u,w,v = svd (randmat 300 300) 
            () 
        () 

    let testdet() = 
        for i=1 to 10 do 
            let d = det (randmat 650 650) 
            () 
        () 

    timeit testlu "lu"  
    timeit testsvd "svd" 
    timeit testdet "det" 

1. 우선 위의 소스를 그대로  컴파일하였다.

컴파일하니 아래와 같이 에러가 뜬다.
 
D:\test\f#>fsc TestFSharpMatrix.fsx
Microsoft (R) F# 2.0 Compiler build 2.0.0.0
Copyright (c) Microsoft Corporation. All Rights Reserved.

TestFSharpMatrix.fsx(106,85): error FS0039: The field, constructor or member 'sp
lit' is not defined

2. 파일 확장명 바꾸고, #r 디렉티브 수정하기

dll 파일을 현재 디렉토리에 놓고 컴파일하기 위해 소스 선두의 #r 디렉티브를 수정하였다.
파일 확장명도 fsx로 바꾸었다.


// 1차 수정된 파일명: TestFSharpMatrix.fsx (소스의 선두 #r 디렉티브 부분 수정)

#r @".\dnAnalytics.dll" 
#r @".\FSharp.PowerPack.dll" 
 
open dnAnalytics.LinearAlgebra 
open Microsoft.FSharp.Math 
open dnAnalytics.LinearAlgebra.Decomposition 
 

// 이하 생략

컴파일하면 이전과 동일한 에러메세지가 뜬다,

D:\test\f#>fsc TestFSharpMatrix.fsx
Microsoft (R) F# 2.0 Compiler build 2.0.0.0
Copyright (c) Microsoft Corporation. All Rights Reserved.

TestFSharpMatrix.fsx(106,85): error FS0039: The field, constructor or member 'sp
lit' is not defined

3. FSharp.PowerPack.Compatibility.dll 도 사용하자.

F# Power Pack 의  FSharp.PowerPack.Compatibility.dll 을 현재 디렉토리에 넣어 준다.


// 2차 수정된 파일명: testFShapMatrix.fsx (FSharp.PowerPack.Compatibility.dll 을 #r 리렉티브에 추가)

2차 수정된 소스는 컴파일은 되지만 실행 결과가 없다.

D:\test\f#>fsc TestFSharpMatrix.fsx
Microsoft (R) F# 2.0 Compiler build 2.0.0.0
Copyright (c) Microsoft Corporation. All Rights Reserved.

TestFSharpMatrix.fsx(136,1): warning FS0988: Main module of program is empty: no
thing will happen when it is run


4, 실행위한 Main 모듈을 추가하자.

// 3차 수정된 파일명: testString.fsx (소스의 7째 줄 Split을 String.split로 수정)
// 3차 수정된 파일명: testFShapMatrix.fsx (FSharp.PowerPack.Compatibility.dll 을 #r 리렉티브에 추가)

이제는 컴파일도 잘 되고, 실행도 잘 된다.
실행 후 잠시 기다리면 실행에 소요된 시간이 밀리초 단위로 컨솔에 출력된다.

D:\test\f#>fsc TestFSharpMatrix.fsx
Microsoft (R) F# 2.0 Compiler build 2.0.0.0
Copyright (c) Microsoft Corporation. All Rights Reserved.

D:\test\f#>TestFSharpMatrix
lu Needed 21519.869700 ms
svd Needed 6146.898400 ms
det Needed 4221.604900 ms

* 이로서 이 글의 선두에 적어 놓은 목적이 해결되었다.




이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-변경금지 2.0 대한민국 라이센스에 따라 이용하실 수 있습니다.

목적:
http://geekswithblogs.net/MarkPearl/archive/2010/04/02/creating-my-first-f-program-in-my-new-ldquoexpert-f.aspx 에 올려진 (실행되지 않는) F# 예제가 실행되도록 하자.


// 파일명: testString.fs (에러가 많은 원래 소스)

#light

  let wordCount text =
    let words = Split [' '] text
    let wordset = Set.ofList words
    let nWords = words.Length
    let nDups = words.Length - wordSet.Count
    (nWords, nDups)

let showWordCount text =
    let nWords,nDups = wordCount text
    printfn "--> %d words in text" nWords
    printfn "--> %d duplicate words" nDups

1. 우선 위의 소스를 그대로  컴파일하였다.

컴파일하니 아래와 같이 에러가 2개 뜬다.
 
D:\test\f#>fsc testStringSplit.fs
Microsoft (R) F# 2.0 Compiler build 2.0.0.0
Copyright (c) Microsoft Corporation. All Rights Reserved.

testStringSplit.fs(5,17): error FS0039: The value or constructor 'Split' is not
defined

testStringSplit.fs(8,32): error FS0039: The namespace or module 'wordSet' is not
 defined



2. String.split 해결하기

우선 String.split 가 없다고 하는 문제를 해결하기 위해 #r 디렉티브를 써서 소스 첫 부분에
#r "FSharp.PowerPack.dll";;
#r "FSharp.PowerPack.Compatibility.dll";;
릉 추가하여 다시 컴파일하였다. (여기도 FSharp.PowerPack.Compatibility.dll 도 꼭 있어야 함이 중요하다. 이것을 모르면 엉뚱한 곳에서 삽질만 하게 된다.)


// 1차 수정된 파일명: testString.fs

이제 1차 수정된 소스를 컴파일하였다.
아래 처럼 #r 디렉티브에 관련된 에러가 뜬다.

D:\test\f#>fsc testStringSplit.fs
Microsoft (R) F# 2.0 Compiler build 2.0.0.0
Copyright (c) Microsoft Corporation. All Rights Reserved.

testStringSplit.fs(1,1): error FS0076: #r directives may only occur in F# script
 files (extensions .fsx or .fsscript). Either move this code to a script file, a
dd a '-r' compiler option for this reference or delimit the directive with '#if
INTERACTIVE'/'#endif'.

testStringSplit.fs(2,1): error FS0076: #r directives may only occur in F# script
 files (extensions .fsx or .fsscript). Either move this code to a script file, a
dd a '-r' compiler option for this reference or delimit the directive with '#if
INTERACTIVE'/'#endif'.

3. 파일 확장명 바꾸기

#r 디렉티브가 있는 소스의 확장명은 fsx라야 한다고 하니 그렇게 바꾸었다.

// 2차 수정된 파일명: testString.fsx (소스는 그대로이고 파일 확장명만 fsx로 변경)

2차 수정된 소스를 컴파일 하였다.
컴파일하기 위해서는 2개의 파일 FSharp.PowerPack.dll 과 Fharp.PowerPack.Compatibility.dll 이 현재 디렉토리에 있어야 한다. 이는  http://fsharppowerpack.codeplex.com/  에서 구할 수 있다.

D:\est\f#>fsc testStringSplit.fsx
Microsoft (R) F# 2.0 Compiler build 2.0.0.0
Copyright (c) Microsoft Corporation. All Rights Reserved.

testStringSplit.fsx(7,17): error FS0039: The value or constructor 'Split' is not
 defined

testStringSplit.fsx(10,32): error FS0039: The namespace or module 'wordSet' is n
ot defined

4, Split가 아니라 String.split

// 3차 수정된 파일명: testString.fsx (소스의 7째 줄 Split을 String.split로 수정)

3차 수정된 소스를 컴파일 하였다.


D:\test\f#>fsc testStringSplit.fsx
Microsoft (R) F# 2.0 Compiler build 2.0.0.0
Copyright (c) Microsoft Corporation. All Rights Reserved.

testStringSplit.fsx(10,32): error FS0039: The namespace or module 'wordSet' is n
ot defined

5. 4차 수정 (오타를 바로 잡다)

// 4차 수정된 파일명: testString.fsx (소스의 8째 줄 wordset을 wordSet으로 수정)

 

4차 수정된 소스를 컴파일 하였다.
컴파일이 잘된다. 그런데  Main 모듈이 없다는 경고(warning)가 뜬다.

D:\test\f#>fsc testStringSplit.fsx
Microsoft (R) F# 2.0 Compiler build 2.0.0.0
Copyright (c) Microsoft Corporation. All Rights Reserved.

testStringSplit.fsx(18,1): warning FS0988: Main module of program is empty: noth
ing will happen when it is run

Main 모듈이 없다는 위의 경고를 무시하고 실행해보았다.

D:\test\f#>TestFSharpMatrix

실행은 되지만 역시 아무 출력이 없다.

6. 실행을 위한 Main 모듈 추가하기

// 5차 수정된 파일명: testString.fsx (소스의 끝 부분에 Main 모듈 추가)

마지막으로 수정된 위의 소스를 컴파일한다.

D:\test\f#>fsc testStringSplit.fsx
Microsoft (R) F# 2.0 Compiler build 2.0.0.0
Copyright (c) Microsoft Corporation. All Rights Reserved.

컴파일된 실행파일(testStringSplit.exe)을 실행한다.

D:\test\f#>testStringSplit
--> 6 words in text
--> 2 duplicate words



* 이로서 이 글의 선두에 적어 놓은 목적이 해결되었다.




이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-변경금지 2.0 대한민국 라이센스에 따라 이용하실 수 있습니다.

2**2**5  은 (2**2)**5  와 2**(2**5) 중에 어느 것과 같을까?
즉, 지수 연산자 **는 왼쪽에서 오른쪽으로(left to right) 진행할까?
아니면 오른쪽에서 왼쪽으로(right to left) 진행할까?

만일 2**2**5 == (2**2)**5  라면 그 결과는 2**10 과 같을 것이고,
만일  2**2**5 == 2**(2**5)  라면 그 결과는 2**32 과 같을 것이다.

* 미리 결론을 말하면 Groovy 언어에서는 지수 계산의 순서가 다른 언어(Python. Ruby, Maxima 등)의 것과는 다르다는 것이다.


1. groovy 의 경우

groovysh 로 테스트하였다.

2. python의 경우

python으로 테스트하였다.

3. ruby 의 경우

irb로 테스트하였다.

4. wxMaxima 의 경우

wxMaxima로 테스트하였다.

이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-변경금지 2.0 대한민국 라이센스에 따라 이용하실 수 있습니다.

Java 7에 포함되기로 했다가 취소된 함수형 프로그래밍(closure) 미리 맛보기

    관련 자료
        [1] http://javac.info/
        [2] http://stronglytypedblog.blogspot.com/2010/06/lambdas-in-java-preview-part-1-basics.html )
        [3] http://www.infoq.com/articles/lambdas-java-analysis
        [4] http://wiki.java.net/bin/view/JDK/ClosuresSyntaxInJava7
        [5] http://www.ibm.com/developerworks/java/library/j-jtp04247.html
        [6] http://www.java.net/pub/pq/242
        [7] http://www.baptiste-wicht.com/2010/05/oracle-pushes-a-first-version-of-closures/
        [8] http://www.davidflanagan.com/2010/05/closures-in-jav.html


1. FuctionalJava 홈페이지
   http://functionaljava.org/examples/

2. 다운로드
    http://functionaljava.org/download/


3. 예제

JDK 1.5 이상 용 예제 파일 (Java 7 용은 별도로 있음)

// 소스 파일: Fibs.java

4. 컴파일하기

    javac -d . -cp functionaljava.jar Fibs.java


5. 실행하기 (옵션 3은 쓰레드 개수)

    java -classpath .;functionaljava.jar concurrent.Fibs 3

실행 결과:
ibs 3
Calculating Fibonacci sequence in parallel...
n=0 => 0
n=1 => 1
n=2 => 1
n=3 => 2
n=4 => 3
n=5 => 5
n=6 => 8
n=7 => 13
n=8 => 21
n=9 => 34
n=10 => 55
n=11 => 89
n=12 => 144
n=13 => 233
n=14 => 377
n=15 => 610
n=16 => 987
n=17 => 1,597
n=18 => 2,584
n=19 => 4,181
n=20 => 6,765
n=21 => 10,946
n=22 => 17,711
n=23 => 28,657
n=24 => 46,368
n=25 => 75,025
n=26 => 121,393
n=27 => 196,418
n=28 => 317,811
n=29 => 514,229
n=30 => 832,040
n=31 => 1,346,269
n=32 => 2,178,309
n=33 => 3,524,578
n=34 => 5,702,887
n=35 => 9,227,465
n=36 => 14,930,352
n=37 => 24,157,817
n=38 => 39,088,169
n=39 => 63,245,986
n=40 => 102,334,155
n=41 => 165,580,141
n=42 => 267,914,296
n=43 => 433,494,437
n=44 => 701,408,733
n=45 => 1,134,903,170




이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-변경금지 2.0 대한민국 라이센스에 따라 이용하실 수 있습니다.

1. Add.cs 파일

2. Mult.cs 파일

3. TestCode.cs 파일

4. dll 파일 만들기

    csc /target:library /out:MathLibrary.DLL Add.cs Mult.cs


5. 테스트 실행 파일 만들기

csc /out:TestCode.exe /reference:MathLibrary.DLL TestCode.cs


6. 실행하기

  TestCode 1234 5678

1. 다운로드
    http://www.alglib.net/ 에서 C#용 라이브러리 소스를 다운로드한다.


2. 다운로드한 파일의 압축을 적당한 디렉토리에 푼다.


3. 그 디렉토리로 가서 빌드한다. (옵션은 csc 와 mono 중에 택 1)
    build csc

    (빌드가 끝나면 out 디렉토리에 libalglib.dll 파일이 생성되어 있을 것이다.)


4. 빌드 후 테스트  (옵션은 csc 와 mono 중에 택 1)
   check csc all


5. 개별 테스트 (옵션 /reference 를 이용하여 dll 파일을 지정한다.)

   cd examples

   csc _demo_rcond_1.cs /reference:..\out\libalglib.dll
   _demo_rcond_1

   csc _demo_ode_example1.cs /reference:..\out\libalglib.dll
   _demo_ode_example1

   csc _demo_ode_example2.cs /reference:..\out\libalglib.dll
   _demo_ode_example2

소스파일명: hello.ml
컴파일하기: ocamlc -o hello.exe hello.ml
실행 명령> hello
실행 결과:

컴파일 없이 실행하기: ocaml hello.ml
실행 결과:

 
 

이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-변경금지 2.0 대한민국 라이센스에 따라 이용하실 수 있습니다.

다항식 p(x) 를 1차 다항식 ax - b 로 나눌 때의 몫과 나머지를 구하는 조립제법을
Python 언어로 구현해 보았다. 조립제법은 일명 Horner의 방법이라고도 불리우는데, 이는
x = a 에서 다항식 p(x)의 값 p(a)을 계산하는 가장 빠른 알고리즘이기도 하다.

         p(x) = (ax - b)q(x) + r

여기서 r은 나머지이며 r = p(b/a) 이다. 또 q(x)는 몫이다.

[참고]
    * 온라인으로 조립제법 표 만들기 손으로 계산하는 조립제법 표 
    * 온라인으로 구하는 다항식의 도함수: 조립제법을 이용한 다항식의 도함수


아래의 소스파일은 파일 testSyntheticDivision.py 를 수정한 것이다.

python 대신 jython이나 IronPython 으로도 수정 없이 그대로 실행된다.

1. # coding: MS949
3. #  Filename: testSyntheticDivision2.py
4. #
5. #  Purpose:  Find the quotient and remainder when some polynomial is
6. #            divided by a monic polynomial of the first degree.
7. #
8. #  Execute:  python testSyntheticDivision2.py 5 -4 7 8 6 8
9. #        Or  jython testSyntheticDivision2.py 5 -4 7 8 6 8
11. import sys
13. # 사용법 표시
14. def printUsage():
15.     print("사용법: python testSyntheticDivision2.py [제식의 계수] [피제식의 계수들]")
16.     print("조립제법(synthetic method)에 의한 다항식 나눗셈 결과를 보여준다.")
18. # 부동소수점수의 표현이 .0 으로 끝나는 경우 이를 잘라낸다.
19. # 전체 문자열 표시 너비는 매개변수 width 로 전달받아 처리한다.
20. def simplify(v, width=None):
21.     t = "%g" % v
22.     tlen = len(t)
23.     if t[tlen-2:tlen-1] == ".0":
24.         t = t[0, -2]
25.     if width != None:
26.         if tlen < width:
27.             t = "              "[0: width - tlen] + t
28.     return t
30. # "/수" 표시하기
31. # 부동소수점수의 표현이 .0 으로 끝나는 경우 이를 잘라낸다.
32. # 전체 문자열 표시 너비는 매개변수 width 로 전달받아 처리한다.
33. def simplifyReciprocal(v, width=None):
34.     t = "/%g" % v
35.     tlen = len(t)
36.     if t[tlen-2:tlen-1] == ".0":
37.         t = t[0, -2]
38.     if width != None:
39.         if tlen < width:
40.             t = "              "[0: width - tlen] + t
41.     return t
43. # 다항식을 내림차순의 스트링 표현으로 반환
44. def toPolyString(c):
45.     t = ""
46.     sc0 = simplify(c[0])
47.     if len(c) > 2:
48.         if sc0 == "1":
49.             t += "x^%d" % (len(c)-1)
50.         elif sc0 == "-1":
51.             t += "-x^%d" % (len(c)-1)
52.         else:
53.             t += sc0 + " x^%d" % (len(c)-1)
54.     elif len(c) == 2:
55.         if sc0 == "1":
56.             t += "x"
57.         elif sc0 == "-1":
58.             t += "-x"
59.         else:
60.             t += sc0 + " x"
61.     elif len(c) == 1:
62.         t += sc0
64.     for i in range(1, len(c)):
65.         k = len(c) - 1 - i
66.         sc = simplify(c[i])
67.         if k > 1:
68.             if c[i] > 0.0:
69.                 if sc == "1":
70.                     t += " + " + "x^%d" % k
71.                 else:
72.                     t += " + " + sc + " x^%d" % k
73.             elif c[i] < 0.0:
74.                 if sc == "-1":
75.                     t += " - " + "x^%d" % k
76.                 else:
77.                     t += " - " + simplify(abs(c[i])) + " x^%d" % k
78.         elif k == 1:
79.             if c[i] > 0.0:
80.                 if sc == "1":
81.                     t += " + " + "x"
82.                 else:
83.                     t += " + " + sc + " x"
84.             elif c[i] < 0.0:
85.                 if sc == "-1":
86.                     t += " - " + "x"
87.                 else:
88.                     t += " - " + simplify(abs(c[i])) + " x"
89.         elif k == 0:
90.             if c[i] > 0.0:
91.                 t += " + " + sc
92.             elif c[i] < 0.0:
93.                 t += " - " + simplify(abs(c[i]))
94.     return t
96. # 다항식 나눗셈 결과를
97. #     (피제식) = (제식)(몫) + (나머지)
98. # 형태로 출력
99. def printDivisionResult(a, c, d, b):
100.     strLine = "  " + toPolyString(c)
101.     print( strLine )
103.     strLine = "    = ( " + toPolyString( [a[0], -a[1]] ) + " )"
104.     tmp = [0.0] * (len(b) - 1)
105.     for i in range(0, len(tmp)):
106.         tmp[i] = b[i]
107.     strLine += "( " + toPolyString(tmp) + " )"
108.     r = d[len(d) - 1]
109.     if r > 0.0:
110.         strLine += " + " + simplify(r)
111.     elif r < 0.0:
112.         strLine += " - " + simplify(abs(r))
114.     print( strLine )
116. # 조립제법 계산표 출력 함수
117. def printSyntheticTable(a, c, s, d, q):
118.     strLine = "         | "
119.     strLine += simplify(c[0], 8)
120.     for i in range(1, len(c)):
121.         strLine += "  " + simplify(c[i], 8)
122.     print( strLine )
124.     strLine = simplify(a[1], 8) + " |"
125.     strLine += "           "
126.     strLine += simplify(s[1], 8)
127.     for i in range(2, len(s)):
128.         strLine += "  " + simplify(s[i], 8)
129.     print( strLine )
131.     strLine = "         |"
132.     for i in range(0, len(q)):
133.         strLine += "----------"
134.     print( strLine )
136.     strLine = "         | "
137.     strLine += simplify(d[0], 8)
138.     for i in range(1, len(d)):
139.         strLine += "  " + simplify(d[i], 8)
140.     print( strLine )
142.     strLine = simplifyReciprocal(a[0], 8) + " | "
143.     strLine += simplify(q[0], 8)
144.     for i in range(1, len(q) - 1):
145.         strLine += "  " + simplify(q[i], 8)
146.     print( strLine )
148. # 실행 시작 지점
149. if len(sys.argv) < 3:
150.     printUsage()
151.     sys.exit(1)
153. ######################################################
154. # 피제식은 c_0 x^n +  c_1 x^(n -1) + ... + c_n
155. # 제식은 a x -  b
156. a = [float(sys.argv[1]), - float(sys.argv[2])]
157. c = [0.0]*(len(sys.argv) - 3)
158. s = [0.0]*(len(sys.argv) - 3)
159. d = [0.0]*(len(sys.argv) - 3)
160. b = [0.0]*(len(sys.argv) - 3)
162. for i in range(0, len(c)):
163.     c[i] = float(sys.argv[i + 3])
165. ######################################################
166. # 조립제법의 주요 부분
167. s[0] = 0.0
168. d[0] = c[0]
169. b[0] = c[0] / a[0]
170. for i in range(1, len(c)):
171.     s[i] = b[i-1]*a[1]
172.     d[i] = c[i] + s[i]
173.     b[i] = d[i] / a[0]
175. ######################################################
176. # 몫의 계수와 나머지를 출력한다.
177. print("몫의 계수는"),
178. for i in range(0, len(b) - 2):
179.     print(simplify(b[i], None) + ", " ),
180. print(simplify(b[len(b) - 2], None)),
181. print("이고, 나머지는 " + simplify(d[len(d) - 1], None) + " 이다.")
182. print
184. ######################################################
185. # 조립제법 표를 출력한다.
186. printSyntheticTable(a, c, s, d, b)
187. print
189. ######################################################
190. # (피제식) = (제식) x (몫) + (나머지)
191. printDivisionResult(a, c, d, b)

실행> python testSyntheticDivision2.py 5 -4 7 8 6 8
몫의 계수는 1.4,  2.72,  3.376 이고, 나머지는 21.504 이다.

         |        7         8         6         8
       4 |                5.6     10.88    13.504
         |----------------------------------------
         |        7      13.6     16.88    21.504
      /5 |      1.4      2.72     3.376

  7 x^3 + 8 x^2 + 6 x + 8
    = ( 5 x - 4 )( 1.4 x^2 + 2.72 x + 3.376 ) + 21.504




이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-변경금지 2.0 대한민국 라이센스에 따라 이용하실 수 있습니다.