포트란 뉴턴 랩슨법과 근의 공식 소스
- 최초 등록일
- 2010.12.19
- 최종 저작일
- 2010.10
- 1페이지/ 한컴오피스
- 가격 1,000원
소개글
3차방정식에서 뉴턴랩슨법으로 하나의 근을 구한 다음
남은 두개의 근을 근의 공식으로 구하는 프로그램입니다
목차
없음
본문내용
real a,b,c,d,r,i
write(*,*) `f(x)=a*x**3+b*x**2+c*x+d`
write(*,*) `input a,b,c,d?`
read(*,*) a,b,c,d
i=0.
r=1.
E=1.E-6
do
i=r-f(r)/g(r)
if(abs(f(i)).le.E)then
answer=i
exit
endif
if(abs(f(i)).gt.E)then
r=i
endif
enddo
write(*,*) `answer=`, answer
q=(b+a*i)*(b+a*i)-4.*a*(d/i)
if(q.eq.0.)then
x1=-(b+a*i)/(2*a)
x2=-(b+a*i)/(2*a)
elseif(q.gt.0.)then
x1=(-(b+a*i)+sqrt(q))/(2*a)
x2=(-(b+a*i)-sqrt(q))/(2*a)
else
xreal=-(b+a*i)/(2*a)
ximag=sqrt(-q)/(2*a)
참고 자료
없음