diff --git a/g_GMSK b/g_GMSK
new file mode 100644
index 0000000000000000000000000000000000000000..1c66951e6135f6595b7b31e7d01f64ddd875a244
--- /dev/null
+++ b/g_GMSK
@@ -0,0 +1,14 @@
+function [out] = g_GMSK(t,K,Tb,BTb)
+% g_GMSK produces the coefficient of a Gaussian filter such as used in GMSK
+% modulation.
+% * inputs: t is a vector used as support for the Gaussian
+% K is a parameter for the Gaussian amplitude
+% Tb is the bit period
+% BTb is the Bandwidth - Bit period product
+% * output: out is a vector containing the filter coefficients
+
+B = BTb/Tb;
+
+out = -K/(2*sqrt(log(2))) * (erf(2*B*pi* (t-Tb/2)/sqrt(log(4))) - erf(2*B*pi* (t+Tb/2)/sqrt(log(4))));
+
+end
\ No newline at end of file