using IntervalArithmetic


function f(x, y)
    return 9 * x^4 - y^4 + 2 * y^2
end


x = @interval(40545)
y = @interval(70226)
println(@sprintf "x = %s y = %s f(x,y) = %s" x y f(x, y) )