% Session p. 21
format long

function y=func(t);
	z=t+i*t.*(1-t); y=real(z.^(i./z-1).*(1+i*(1-2*t)));
endfunction

quad('func',0,1)