LoginSignup
0
0

積分方程式(定数型)「2009 奈良教育大学前期【3】」をChatGPTとWolframAlphaとsympyでやってみたい。

Last updated at Posted at 2023-09-30

・次のリンクで積分方程式をまとめています。ガンバリマス。
 (本ページは,モトネタです。)
https://qiita.com/mrrclb48z/items/d05574b3c60327e96ab2#fxgx%CE%BBintegratehtfttab
WolframAlpha先生が、
「Solution as a Fredholm integral equation」Fredholm? とおっしゃています。
 matahematicaを調べる必要がありそうです。

オリジナル

上と同じです。大学入試数学問題集成 様>テキスト【3】

ChatGPT-3.5先生へ(???私の方が、わからなくなってきました。???)

入力文(一括です。)

次の等式を満たす関数f(x)を求めよ.
f(x)=cos(x)+2*integrate(t*f(t)*sin(t),(t,0,pi/2))  

??? 返信抜粋

...詳細を教えていただければ、もっと具体的なサポートを提供できるかもしれません。
sympyで
...# f(t)を解析的に求める試み(一般的な解を求めることは難しいかもしれません)

私の方が、わからなくなってきました。私の質問が変なのかもしれません。
私のChatGPTは、中断です。

WolframAlphaで

・これが一番でした。

Solution as a Fredholm integral equation
f(x) = cos(x) - π/4

sympyで(T氏の数学日記 様の方法を参考に)

・次のリンクで積分方程式をまとめています。ガンバリマス。
https://qiita.com/mrrclb48z/items/d05574b3c60327e96ab2#fxgx%CE%BBintegratehtfttab
・以下が、上のモトネタです。
・同値変形?でやってみました。

from sympy import *
var('x,t,A',real=True)
f=Function('f')
eq1=Eq(f(x),cos(x)+2*integrate(t*f(t)*sin(t),(t,0,pi/2))  )  ;print("#1",eq1)
eq2=Eq(              integrate(t*f(t)*sin(t),(t,0,pi/2)),A)  ;print("#2",eq2)
eq3=eq1.subs({eq2.lhs:eq2.rhs})                              ;print("#3",eq3)
eq4=eq3.subs({x:t})                                          ;print("#4",eq4)
eq5=Eq(eq2.rhs,eq2.lhs)                                      ;print("#5",eq5)
eq6=eq5.subs({eq4.lhs:eq4.rhs})                              ;print("#6",eq6)
eq7=Eq(eq6.lhs,eq6.rhs.doit())                               ;print("#7",eq7)
eq8=Eq(A,solve(eq7)[0])                                      ;print("#8",eq8)
eq9=eq3.subs({eq8.lhs:eq8.rhs})                              ;print("#9",eq9)
print()
eqA=eq1                                                      ;print("#A",eqA)
eqB=eq9.subs({x:t})                                          ;print("#B",eqB)
eqC=eqA.subs({eqB.lhs:eqB.rhs})                              ;print("#C",eqC)
eqD=Eq(eqC.lhs,eqC.rhs.doit())                               ;print("#D",eqD)
eqE=Eq(eq9.rhs,eqD.rhs)                                      ;print("#E",eqE)
#1 Eq(f(x), cos(x) + 2*Integral(t*f(t)*sin(t), (t, 0, pi/2)))
#2 Eq(Integral(t*f(t)*sin(t), (t, 0, pi/2)), A)
#3 Eq(f(x), 2*A + cos(x))
#4 Eq(f(t), 2*A + cos(t))
#5 Eq(A, Integral(t*f(t)*sin(t), (t, 0, pi/2)))
#6 Eq(A, Integral(t*(2*A + cos(t))*sin(t), (t, 0, pi/2)))
#7 Eq(A, 2*A + pi/8)
#8 Eq(A, -pi/8)
#9 Eq(f(x), cos(x) - pi/4)

#A Eq(f(x), cos(x) + 2*Integral(t*f(t)*sin(t), (t, 0, pi/2)))
#B Eq(f(t), cos(t) - pi/4)
#C Eq(f(x), cos(x) + 2*Integral(t*(cos(t) - pi/4)*sin(t), (t, 0, pi/2)))
#D Eq(f(x), cos(x) - pi/4)
#E True

sympyの実行環境

①私の環境は,pycharmです。
②よく聞くのは、Jupyterです。
③web上で、上記のソースを「SymPy Live shell」に、コピー貼り付けでもできました。
黒背景の右上に、マウスを移動すると、コピーマークが発生します。
??? タブレット環境で、コピー貼り付けが実行できませんでした。???

参考

以下、いつもの?おすすめです。

0
0
0

Register as a new user and use Qiita more conveniently

  1. You get articles that match your needs
  2. You can efficiently read back useful information
  3. You can use dark theme
What you can do with signing up
0
0