t = linspace(0.5, 4);
alpha = 2;

fun = 1/2*(t+alpha./t);

plot(t,fun);
hold on;
plot(t,t,'r');
grid on;