function [y] = f_int(x)

y = 2/3*x^(3/2);