Ayuda con SCILAB

[Katerina Dobreva]

Tengo que hacer este ejercicio:
Representar la siguiente función ya definida en el dominio de la frecuencia (se trata de la integral de Fourier de un pulso rectangular), en el rango de frecuencias de -4pi a 4pi rad/s .

Realizar la representación para valores de V = 1 y 5, y para valores de τ = 1 y 3. Observar cómo cambia el espectro en frecuencia en cada caso. Nota: No es necesario utilizar la función espectro puesto que ya se trata de una integral de Fourier.
Este es mi código, pero no me funciona. ¿Alguien sabe por qué? Gracias.

Código [Seleccionar] Expandir

w = -4*%pi:4*%pi
V1 = 1
V2 = 5
t1 = 1
t2 = 3
F1 = V1 * t1 * sin(w * t1 / 2) / (w * t1 / 2)
F2 = V2 * t2 * sin(w * t2 / 2) / (w * t2 / 2)
plot(w, F1, "r")
plot(w, F2, "b")


Mod: Tema corregido, la Imagen era muy grande, los códigos deben ir en etiquetas de código y este tema no es una "duda general" sino sobre programación