こんにちはゲストさん。会員登録(無料)して質問・回答してみよう!

締切り済みの質問

ニュートン法

フォートランでf(x)=0の解を求めるこのようなプログラムを作りました。
program nyu-ton
implicit none
real::x,f,dfdx


read(*,*)x
do
call FDF(x,f,dfdx)
x=x-f/dfdx
write(*,*)x,f
if(abs(f)<1e-7)exit

end do
stop
end program nyu-ton

subroutine FDF(x,f,dfdx)
implicit none
real::x,f,dfdx
f=x*(x-2.0)
dfdx=2.0*x-2.0
end subroutine FDF

このプログラムをニュートン法で複素関数の解を求めるプログラムにするには、どこをどのように変えればよいのかわかりません。

z^3+z=0 や z^3=512iのような問題を数値的に解くプログラムを作りたいんですが、教えてください願いします。

投稿日時 - 2006-11-27 01:17:05

QNo.2564790

すぐに回答ほしいです

このQ&Aは役に立ちましたか?

0人が「このQ&Aが役に立った」と投票しています

回答(1)

ANo.1

実数 (real) じゃなくって複素数 (complex) を使えばいいと思うんだけど....

投稿日時 - 2006-11-28 14:44:17

あなたにオススメの質問