Evo, to je izvorna koda

	program bezierjeva_krivulja

	real kt(100, 3), t, tbk(100,3)
	real bernstein(4), suma	
	integer i, j, koord, n

	kt(1,1) = 1.0
	kt(1,2) = 0.0
	kt(1,3) = 0.0

	kt(2,1) = 2.0
	kt(2,2) = 2.0
	kt(2,3) = 1.0

	kt(3,1) = 5.0
	kt(3,2) = 0.0
	kt(3,3) = 0.0

	kt(4,1) = 6.0
	kt(4,2) = 3.0
	kt(4,3) = 1.0

	n = 51

	do 30 koord=1,3
	  t = 0.0
	  do 10 i=1,n
	    t=(i-1.0)/(n-1.0)
	    bernstein(1) = (1-t)**3
	    bernstein(2) = 3*t*(1-t)**2
	    bernstein(3) = 3*t*t*(1-t)
	    bernstein(4) = t**3
	    t = t+0.1
	  
	    suma = 0.0
	    do 20 j=1,4
	      suma = suma+bernstein(j)*kt(j,koord)
20	    continue
	    tbk(i, koord) = suma

10	  continue
30	continue

	call pisi('tbk.vrml', tbk, n)
	call pisi('kt.vrml', kt, 4)

	end