Roots of a function: Difference between revisions

=={{header|FreeBASIC}}==

Simple bisection method.

const iterations=20000000

dim as double last,st=(max-min)/iterations,v

for n as double=min to max step st

if sgn(v)<>sgn(last) then

redim preserve a(1 to ubound(a)+1)

a(ubound(a))=n

end if

last=v

next

function roots(f1 as function(as double) as double,min as double,max as double, a() as double) as long

dim as double last,O,st=(max-min)/iterations,v

for n as double=min to max step st

if sgn(v)<>sgn(last) and n>min then bisect(f1,n-st,n,O,a()):n=O

last=v

print "in range -20 to 20"

print "All roots approximate"

for n as long=1 to ubound(r)

print n,CRound(r(n),6),,defn(r(n))

next n

end if

All roots approximate

number root to 6 dec places function value at root

1 0 -2.929925652002424e-009

2 1 1.477781779325033e-009

=={{header|Go}}==