10/01/10 04:48:39
ｽﾚﾘﾝｸ(tech板:636番)
# [1] 授業単元：数値計算
# [2] 問題文(含コード&リンク)： URLﾘﾝｸ(ime.nu)
#
# ３元連立非線形方程式
#
# x_1^2 + x_2^2 + x_3^2 - 4 = 0
# x_1^2 - 2x_1 + x_2^2 + x_3 - 2 = 0
# x_1^2 + sin^2x_2 + x_3^2 - 3 = 0
#
# を解くNewton法のC++のプログラムをつくり、領域 0 < x_1,x_2,x_3 < 2 にある解を求めよ
# ただし、Jacobianの計算に必要となる偏導関数は数値微分により計算し、連立１次方程式の解法には
# Gaussの消去法あるいはLU分解法を用いること