オイラー法 値が出力されない if文 数学 微分方程式

Question

前提・実現したいこと

課題でオイラー法のプログラムを組んでいて、x=0.1,0.2,....1.0の値のみを出力したい。

発生している問題・エラーメッセージ

if文の条件文をx==0.1 or 0.2 or....としたところうまくいかなかったので、いろんな値を入れていくとなぜかx==0.001~0.008まではうまく出力されるがx==0.009以降はエラー等はないが何も出力されない状態になった。

該当のソースコード

import math
import numpy as np
import sys
sys.setrecursionlimit(10**3)
def f(y):
    return y

def t(x):
    return (math.e)**x
dx = 0.001
y = 1.0
x = 0

with open('euler.txt','w') as u:
    while x < 1.0:
        y += f(y)*dx
        x += dx
        error = abs(y-t(x))
        if x == 0.009:
            print("xが" "{0:.3f}"  "近似値は" "{1:.15f}" "真の値は" "{2:.15f}" "真の値との差は" "{3:.15f}". format(x,y,t(x),error), file=u)

試したこと

if文の条件文をx==0.1 or 0.2 or....としたところうまくいかなかったので、いろんな値を入れていくとなぜかx==0.001~0.008まではうまく出力されるがx==0.009以降はエラー等はないが何も出力されない状態になった。

補足情報（FW/ツールのバージョンなど）

if文の条件文でx==0.1 or 0.2 or....ではなくもっと良い方法などがあれば教えていただけると幸いです。

Answer 1

Python に限った話ではないですが、浮動小数点の比較は注意が必要です。

いくつか、参考リンクを記載します。

• Pythonで浮動小数点数floatの誤差を考慮して比較（math.isclose） | note.nkmk.me
• 15. 浮動小数点演算、その問題と制限 — Python 3.8.3 ドキュメント
• [Python] 浮動小数点数floatの比較は要注意！！ - Qiita

2つの値が近似しているか判定: math.isclose()

Python3.5で追加された関数なのでそれより前のバージョンでは使えない。
以下のように、第一引数aと第二引数bに指定した値が厳密に等価でなくても近似値であればTrueを返す。

Comments

[Etori]

説明ありがとうございました！