I=quadl(@(t) (t-3*t.^2+2*t.^3).^(-1/3),eps,0.5)