Matlab - Ayuda con Octave

Hola, les escribo para preguntar una duda que tengo al computar lo siguiente con Octave:
Tengo que plotear una primitiva de f(x)=(sin(x))/x
Que código insertarian al programa.Por más que lo intento no lo consigo.

Gracias

PD: Debo seguir como patrón, el modelo de integral como division de un intervalo en "n" rectangulos, siendo el intervale de "a" a "b". Creo que me entendeis lo que os digo.

Un saludo

Hola de nuevo, os dejo el código con el que he probado a plotear la funcion. Al plotearla me sale una recta paralela al eje x (cosa que no me cuadra), pues la solución deberia ser esta (segun Wolfram Alpha): http://www.wolframalpha.com/input/?i=integral+sin%28t%29%2Ft+dt+from+0+to+10

Quieres graficar o quieres integrar?

Para graficar:

% define tus valores de x (ej., de 0 a 4pi en incrementos de .01)
x = 0 : .01 : 4*pi;

% define la función
f = sin(x)./x;

% grafica
plot(x,f)

Para integrar, usa la función quad
q = quad(fun, a, b)