例

Euler法による $\frac{d^2}{dx^2}y=-y, y(0)=0, \frac{d}{dx}y(0)=1$ の数値解を例にとる。
まずプログラムを

        ...
      open(10,file='euler.dat')
      write(10,200) '# x, y, z'
        ...
      do i=1,N
        ...
        write(10,210) x,y,z
      end do
 200  format(a)
        ...

としてデータファイル（ここではeuler.datとした）にデータを書き出す。 gnuplotを起動して

    gnuplot> set title'd^2y/dx^2=-y'
    gnuplot> set xlabel'x'
    gnuplot> set ylabel'y(x)'
    gnuplot> plot'euler.dat' title"Euler's method",sin(x)

とすると、図のようなグラフが現れる。