Question1(b), Fortran - rextester

Run Code | API | Code Wall | Misc | Feedback | Login | Theme | Privacy | Patreon

Question1(b)

Language: Layout:

program false

implicit none

real::a,b,f,x,tol,xv

integer::itrtn,i

print*, "Enter the initial values (a,b):"

read(*,*)a,b

if( (f(a)*f(b)) .gt.0) then

write(*,*)"Can not be done", f(a), f(b);

stop

end if

print*, "Enter the number of iteration:"

read(*,*)itrtn

print*, "Enter the tolerance:"

read(*,*)tol

do 10 i = 1, itrtn

xv = (a * f(b) - b* f(a)) / (f(b)-f(a))

if(f(xv).eq.0) then

print*,"The soln is: ",xv

stop

end if

if(f(a)*f(xv).le.0) then

b = xv

print*,a,b

else

a = xv

print*,a,b

end if

if(abs(a-b).le.tol) then

print*, "The approximate soln is: ",xv

stop

end if

10 continue

print*,"Maximum iteration failed"

end program

real function f(x)

real::x

f = x**3 - 2*x -5

return

end function

Show compiler warnings [ - ] Show input
Compilation time: 0.12 sec, absolute running time: 0.15 sec, cpu time: 0.01 sec, memory peak: 6 Mb, absolute service time: 0,39 sec

edit mode | history | discussion

 Enter the initial values (a,b):
 Enter the number of iteration:
 Enter the tolerance:
   2.05882359       3.00000000    
   2.08126354       3.00000000    
   2.08963919       3.00000000    
   2.09273982       3.00000000    
   2.09388375       3.00000000    
   2.09430552       3.00000000    
   2.09446073       3.00000000    
   2.09451795       3.00000000    
   2.09453917       3.00000000    
   2.09454703       3.00000000    
   2.09454989       3.00000000    
   2.09455109       3.00000000    
   2.09455109       2.09455156    
 The approximate soln is:    2.09455156

λ

                                      .NET NoSQL database for rapid development