function [z] = calcg(t, T, cc)
    exec('**/place.sci', -1)

    i = place(T, t)
    dt = t-T(i)
    dt2 = dt * dt

    z = cc(i, 1)+cc(i, 2)*dt+cc(i,3)*dt2+cc(i,4)*dt2*(t-T(i+1))
endfunction