Low-Pass Filter Coefficients

How can I write a function that computes low-pass filter coefficients using the rectangular window? Using Thetac as the cutoff frequency in radians and nc as the number of coefficients. So far I have:

function h = lowpassfiltercoefs ( thetac, nc)

, I do not know how to continue.