QUICK

等間隔格子です。
XOJOです。

必要とあればTVD化のソースや不等間隔TVD-QUICKもあります。

'以下ソース

  dim f1 as double
   if u >= 0 then
    f1 = f-u*dt*(3*fdn1+3*f-7*fup1+fup2)/(8*dx)
  else
    f1 = f-u*dt*(3*fup1+3*f-7*fdn1+fdn2)/(8*dx)
  end if
  return f1