function value = poly_int ( a, b )
  ya = 2*a.^5 - 3*a.^4 + a.^2;
  yb = 2*b.^5 - 3*b.^4 + b.^2;
  value = yb - ya;
  return
end

