View Options

	diff --git a/darksusy-6.2.2/src/rd/dsrdquad.f b/darksusy-6.2.2/src/rd/dsrdquad.f
	index 408a35c..e72553e 100644
	--- a/darksusy-6.2.2/src/rd/dsrdquad.f
	+++ b/darksusy-6.2.2/src/rd/dsrdquad.f
	@@ -1,37 +1,37 @@
	real*8 function dsrdquad(x0,y0,x1,y1,x2,y2,x,reslin,err)
	**********************************************************************
	*** Function dsrdquad performs a quadratic interpolation between the
	*** points (x0,y0), (x1,y1) and (x2,y2). It will return value from
	-*** the interpolation at point x. x is assumed to be x0 ≤ x ≤ x2.
	+*** the interpolation at point x. x is assumed to be x0 <= x <= x2.
	*** The points are assumed to be ordered in x, x0<x1<x2.
	*** Upon return, err is the relative error compared to a linear
	*** interpolation. The setup is done with Newton divided differences.
	***
	*** Author: Joakim Edsjo, edsjo@fysik.su.se
	*** Date: 2018-10-27
	**********************************************************************

	implicit none

	real*8 x0,y0,x1,y1,x2,y2,x,err
	real*8 f0,f1a,f1b,f2
	real*8 res,reslin

	f0=y0
	f1a=(y1-y0)/(x1-x0)
	f1b=(y2-y1)/(x2-x1)
	f2=(f1b-f1a)/(x2-x0)
	res=f0+f1a(x-x0)+f2(x-x0)*(x-x1) ! quad interpolation

	c..Compare with linear
	if (x.le.x1) then
	reslin=y0+(y1-y0)/(x1-x0)*(x-x0)
	else
	reslin=y1+(y2-y1)/(x2-x1)*(x-x1)
	endif
	err=(res-reslin)/reslin
	dsrdquad=res

	return
	end

File Metadata