[1] 授業単元：数値計算
[2] 問題文(含コード&リンク)：
２元連立非線形方程式
x^2+y^2-2=0
x^2-y^2-1=0
を解く以下のプログラムhttp://kansai2channeler.hp.infoseek.co.jp/cgi-bin/joyful/img/3425.cppを例に、
３元連立非線型方程式
x1^2+x2^2+x3^2-4=0
x1^2-2x1+x2^2+x3-2=0
x1^2+sin^2(x2)+x3^2-3=0
を解くニュートン法のC++プログラムを作り、領域0<x1,x2,x3<2にある解を求めよ。
やだし、ヤコビアンの計算に必要となる変動関数は数値微分により計算し、
連立一次方程式の解法にはガウスの消去法あるいはLU分解法を用いる事。
[3] 環境
　[3.1] OS： WindowsXP
　[3.2] コンパイラ名とバージョン： cygwin
　[3.3] 言語： C++
[4] 期限： ([07年1月17日10:00まで]
[5] その他の制限：特になし

さっぱり分からないのでどなたかご教示願います。