ONLY DO WHAT ONLY YOU CAN DO

こけたら立ちなはれ 立ったら歩きなはれ

2013-01-01から1年間の記事一覧

Clojure で 積分(台形則)して π を求める

πの求め方 積分(台形則) (defn f[x] (/ 4 (+ 1 (* x x)))) (def a 0) (def b 1) ; 台形則で積分 (doseq [j (range 1 11)] (def n (Math/pow 2 j)) (def h (/ (- b a) n)) (def x a) (def s 0) (doseq [i (range 1 n)] (def x (+ x h)) (def s (+ s (f x)))) …

Clojure で 連分数展開を使って π を求める

πの求め方 こう書いてもおk ;自作の逆正接関数 (defn myAtan [x x2 n t] (def m (quot n 2)) (def denom (/ (* (* m m) x2) (+ n t))) (def nume (- n 2)) (if (<= nume 1) (/ x (+ 1 denom)) (myAtan x x2 nume denom))) (doseq [degree (map #(- (* % 15)…

Haskell で 連分数展開を使って π を求める

πの求め方 こう書いてもおk import Text.Printf import Control.Monad -- 自作の逆正接関数 myAtan::Double->Double->Int->Double->Double myAtan x x2 n t = let m = (fromIntegral (n `div` 2)) denom = (m * m * x2) / ((fromIntegral n) + t) num = n -…

C# で 連分数展開を使って π を求める

πの求め方 こう書いてもおk using System; public class CS0512 { public static void Main() { for (int i = 11; i <= 31; i += 2) { double radian = 1; double x2 = radian * radian; // 自作の逆正接関数 double d1 = myAtan(radian, x2, i, 0.0); // i…

Objective-C で 積分(台形則)して π を求める

πの求め方 積分(台形則) #import <Foundation/Foundation.h> #import <math.h> double f(double x) { return 4 / (1 + x * x); } int main() { const double a = 0; const double b = 1; // 台形則で積分 int n = 2; int i, j; for (j = 1; j <= 10; j++) { double h = (b - a) / n; double s </math.h></foundation/foundation.h>…

Java で 積分(台形則)して π を求める

πの求め方 積分(台形則) public class Java0601 { public static void main(String []args) { final double a = 0; final double b = 1; // 台形則で積分 int n = 2; for (int j = 1; j <= 10; j++) { double h = (b - a) / n; double s = 0; double x = a; …

C# で 積分(台形則)して π を求める

πの求め方 積分(台形則) using System; public class CS0601 { public static void Main() { const double a = 0; const double b = 1; // 台形則で積分 int n = 2; for (int j = 1; j <= 10; j++) { double h = (b - a) / n; double s = 0; double x = a; f…

D言語で 積分(台形則)して π を求める

πの求め方 積分(台形則) import std.stdio; import std.math; void main(string[] args) { const double a = 0; const double b = 1; // 台形則で積分 int n = 2; for (int j = 1; j <= 10; j++) { double h = (b - a) / n; double s = 0; double x = a; for…

Go言語で 積分(台形則)して π を求める

πの求め方 積分(台形則) package main import "fmt" import "math" func main() { const a = 0 const b = 1 // 台形則で積分 var n int = 2 for j := 1; j <= 10; j++ { var h float64 = (b - a) / float64(n) var s float64 = 0 var x float64 = a for i :=…

Dart に printf はないのんか?

"TRY DART!" というページを見つけました。 Dartの環境を作らなくても Online で、試せます。 http://try.dartlang.org/以前紹介した、CompileOnline でも、Dart を試せますが。 http://www.compileonline.com/execute_dart_online.phpで、"Dart API Referen…

Access で 画面の全項目の定義を出力

Dim fs: Set fs = CreateObject("Scripting.FileSystemObject") Dim tso: Set tso = fs.OpenTextFile("C:\xxx\AllFormsControls.txt", 2, True) Dim ct As DAO.Container For Each ct In CurrentDb.Containers If ct.Name = "Forms" Then Dim doc As DAO.Doc…

Access で 帳票の全項目の定義を出力

Dim fs: Set fs = CreateObject("Scripting.FileSystemObject") Dim tso: Set tso = fs.OpenTextFile("C:\xxx\AllReportsControls.txt", 2, True) Dim ct As DAO.Container For Each ct In CurrentDb.Containers If ct.Name = "Reports" Then Dim doc As DAO…

Amazon Top100 入り

さまざまな言語で数値計算 第1巻 級数展開・連分数展開作者: 山岡直樹出版社/メーカー: ForNext発売日: 2013/11/01メディア: Kindle版この商品を含むブログ (5件) を見る

ようやく出版

さまざまな言語で数値計算 第1巻 級数展開・連分数展開作者: 山岡直樹出版社/メーカー: ForNext発売日: 2013/11/01メディア: Kindle版この商品を含むブログ (5件) を見るプレビューしてみたらレイアウトがガタガタで悲しい...こちらが、calibre で mobi に…

calibleでePubファイルをmobiファイルに変換

表紙 内容

SigilというソフトでePubファイルを作って,calibleというソフトでmobiファイルに変換してみた

ePub ファイルを calible ってソフトで読み込んで表示すると... mobiファイルに変換すると... cssが全然効いてないわけでもないし... 何なん?また, この ePub ファイルは calible ってソフトでは普通に開けるんだけど, 紀伊国屋の kinoppy で開こうとすると…

SigilというソフトでePubファイルを作ってみた

こちらの記事とか こちらの記事を参考に SigilというソフトでePubファイルを作ってみた元の html は、こんな感じ それを Sigil に貼っつけて 保存するだけで ePub ファイルが作成できる。 (改ページは挿入したけど) その ePub ファイルを calible ってソフト…

Haskell で FizzBuzz

Compile and Execute Programs Online fizzbuzz :: Int -> String fizzbuzz x | x `mod` 15 == 0 = "FizzBuzz" | x `mod` 3 == 0 = "Fizz" | x `mod` 5 == 0 = "Buzz" | otherwise = show x main = do putStr $ unlines $ map fizzbuzz [1..100] Compiling t…

『CompileOnline』すっげー

自分のPCに、わざわざ環境を作らなくても、こんなことができてしまう... Compile and Execute Programs Online Ada with Ada.Text_IO; use Ada.Text_IO; procedure Hello is begin Put_Line("Hello, world!"); end Hello; Compiling the source code.... $…

PowerShell 文字列中の "[", "]"

PowerShellでは文字列中の "[" や "]" を何か特別扱いしているのか? PS C:\Music\Unknown Artist\Unknown Album> ls ディレクトリ: C:\Music\Unknown Artist\Unknown Album Mode LastWriteTime Length Name ---- ------------- ------ ---- -a--- 2013/04/1…

System.Numerics.BigInteger.Log10

F#

前々から、BigInteger.Log10 が使えなくて、何でだろうと思ってたんですが、 c:\>fsi Microsoft (R) F# 2.0 Interactive build 2.0.0.0 Copyright (c) Microsoft Corporation. All Rights Reserved. For help type #help;; > System.Numerics.BigInteger.Log…

PowerShell で iTunes にインポートする

指定したフォルダ配下の曲を一括インポート function ImportWav($arg) { $args = $arg.fullName.split("\") # 階層構造が 一律ではない # c:\音楽\MP3\グループ分け\アーティスト名\アルバム名\ディスク番号\曲名.mp3 # c:\音楽\MP3\グループ分け\アーティス…

F# で Project Euler Problem 25

> (999.0 + (System.Math.Log10 5.0) / 2.0) / (System.Math.Log10 ((1.0 + (sqrt 5.0)) / 2.0)) - ;; val it : float = 4781.859271 > System.Convert.ToString(System.Numerics.BigInteger.Pow(1618033989I, 4781) / System.Numerics.BigInteger.Pow(10000…

Project Euler Problem 25

フィボナッチ数列で1000桁になる最初の項は何番目か?http://odz.sakura.ne.jp/projecteuler/index.php?Problem%2025フィボナッチ数列は以下の漸化式で定義される:, ただし . 最初の12項は以下である. 12番目の項, が3桁になる最初の項である.1000桁になる最…

F# で Project Euler Problem 24

まづ、0,1,2,3 の数列で試してみる > let problem24 x = - let i = ref x - let s = ref (Set [0..3]) - [3 .. -1 .. 1] - |> List.iter - ( - fun n -> - let j = if n = 1 then 1 else List.reduce(*) [n .. -1 .. 2] - let mutable k = !i / j - i := !i …

Project Euler Problem 24

0,1,2,3,4,5,6,7,8,9からなる順列を辞書式に並べたときの100万番目はいくつか?http://odz.sakura.ne.jp/projecteuler/index.php?cmd=read&page=Problem%2024順列とはモノの順番付きの並びのことである. たとえば, 3124は数 1, 2, 3, 4 の一つの順列である. …

Scala で Project Euler Problem 23

F# 版 を そのまま移植与えられた数を素因数分解して、素因数とその指数とを Map にして返す関数 scala> def get_prime_factor(map:collection.mutable.Map[Long, Long], n: Long, factor: Long = 2) { | if (n >= factor) { | if (n % factor != 0 ) | get_…

Project Euler Problem 23

2つの過剰数の和で書き表せない正の整数の総和を求めよhttp://odz.sakura.ne.jp/projecteuler/index.php?cmd=read&page=Problem%2023完全数とは, その数の真の約数の和がそれ自身と一致する数のことである. たとえば, 28の真の約数の和は, であるので, 28は…

F# で Project Euler Problem 23

与えられた数を素因数分解して、素因数とその指数とを Map にして返す関数 > let rec get_prime_factor (map:Map<int, int> byref) (n:int) (factor: int) = - if n >= factor then - if n % factor <> 0 then - get_prime_factor &map n (factor + 1) - else - if Map</int,>…

Project Euler Problem 23

2つの過剰数の和で書き表せない正の整数の総和を求めよhttp://odz.sakura.ne.jp/projecteuler/index.php?cmd=read&page=Problem%2023完全数とは, その数の真の約数の和がそれ自身と一致する数のことである. たとえば, 28の真の約数の和は, であるので, 28は…