added error bands calculation (using parabolic approx.)