Project 2, Application 9.3
unit := (t,a,b) -> Heaviside(t-a) - Heaviside(t-b);
Zio2JUkidEc2IkkiYUdGJUkiYkdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCYtSSpIZWF2aXNpZGVHRiU2IywmOSQiIiI5JSEiIkYxLUYtNiMsJkYwRjE5JkYzRjNGJUYlRiU=
Define the Fourier coefficients for even (a(n)) or odd (b(n)) extensions and the corresponding Fourier series
a := (f,n) -> (2/Pi)*int(f(t)*cos(n*t),t=0..Pi);
b := (f,n) -> (2/Pi)*int(f(t)*sin(n*t),t=0..Pi);
cossum := (f,N,t) -> a(f,0)/2 + sum(a(f,n)*cos(n*t),n=1..N);
sinsum := (f,N,t) -> sum(b(f,n)*sin(n*t),n=1..N);
Zio2JEkiZkc2IkkibkdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCQqJkkjUGlHJSpwcm90ZWN0ZWRHISIiLUkkaW50R0YlNiQqJi05JDYjSSJ0R0YlIiIiLUkkY29zRzYkRi1JKF9zeXNsaWJHRiU2IyomOSVGN0Y2RjdGNy9GNjsiIiFGLEY3IiIjRiVGJUYl
Zio2JEkiZkc2IkkibkdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCQqJkkjUGlHJSpwcm90ZWN0ZWRHISIiLUkkaW50R0YlNiQqJi05JDYjSSJ0R0YlIiIiLUkkc2luR0YlNiMqJjklRjdGNkY3RjcvRjY7IiIhRixGNyIiI0YlRiVGJQ==
Zio2JUkiZkc2IkkiTkdGJUkidEdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCYtSSJhR0YlNiQ5JCIiISMiIiIiIiMtSSRzdW1HRiU2JComLUYtNiRGL0kibkdGJUYyLUkkY29zRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliR0YlNiMqJkY6RjI5JkYyRjIvRjo7RjI5JUYyRiVGJUYl
Zio2JUkiZkc2IkkiTkdGJUkidEdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUkkc3VtR0YlNiQqJi1JImJHRiU2JDkkSSJuR0YlIiIiLUkkc2luR0YlNiMqJkYzRjQ5JkY0RjQvRjM7RjQ5JUYlRiVGJQ==
Function from fig 9.3.9
f1 := t -> unit(t,0,Pi/2) - unit(t,Pi/2,Pi);
Zio2I0kidEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCYtSSV1bml0R0YlNiU5JCIiISwkSSNQaUclKnByb3RlY3RlZEcjIiIiIiIjRjMtRis2JUYtRi9GMCEiIkYlRiVGJQ==
plot(f1(t),t=0..Pi);
NictJSdDVVJWRVNHNiQ3Ym83JCQiIiEhIiIlKnVuZGVmaW5lZEc3JCQiMlk8JjMvMSQqUkAhIz4kIiM1ISIiNyQkIjFcLjwzNycpelUhIz0kIiM1ISIiNyQkIjFDYkQ3PXo+ayEjPSQiIzUhIiI3JCQiMSpwU2pUQShmJikhIz0kIiM1ISIiNyQkIjJcNV5DT2VSRyIhIz0kIiM1ISIiNyQkIjIoUiJvS1tXPnIiISM9JCIjNSEiIjckJCIyKDRBIVxzO3pjIyEjPSQiIzUhIiI3JCQiMiV6aWBtKikpUVUkISM9JCIjNSEiIjckJCIyJj5XISlcTSRlOCYhIz0kIiM1ISIiNyQkIjFmRDJMenhabyEjPCQiIzUhIiI3JCQiMWI0VTphKW8jKSohIzwkIiM1ISIiNyQkIjFOcHgqRypmIUciISM7JCIjNSEiIjckJCIwNiNleEdtXT4hIzokIiM1ISIiNyQkIjJYVCw9M29eaSMhIzwkIiM1ISIiNyQkIjBcRDBbbmtIJCEjOiQiIzUhIiI3JCQiMSI9WmEmeiUpPVIhIzskIiM1ISIiNyQkIjJ3Yld2KW9HalghIzwkIiM1ISIiNyQkIjElMyoqSGJtKEhfISM7JCIjNSEiIjckJCIxUFJyNCkzVCplISM7JCIjNSEiIjckJCIxI1speXlrWXhsISM7JCIjNSEiIjckJCIxcydvY0dvJHpyISM7JCIjNSEiIjckJCIxUTA7Z3IncCZ5ISM7JCIjNSEiIjckJCIxR010PyRbdGApISM7JCIjNSEiIjckJCIxI3AzMGtAST4qISM7JCIjNSEiIjckJCIwQyJcSGVWKXkqISM6JCIjNSEiIjckJCIyR1spKSopM1cnXDUhIzskIiM1ISIiNyQkIjEnW0BYU0AnNDYhIzokIiM1ISIiNyQkIjJWaFAzRypRejYhIzskIiM1ISIiNyQkIjJVMyozdGM5VDchIzskIiM1ISIiNyQkIjIkeTBEQkEhKjM4ISM7JCIjNSEiIjckJCIyUG1LP1tBTVAiISM7JCIjNSEiIjckJCIxYUBYLkV1UzkhIzokIiM1ISIiNyQkIjJRdSMpKipbakRdIiEjOyQiIzUhIiI3JCQiMkM/Kyp5WCFmYCIhIzskIiM1ISIiNyQkIjI3bTx5bVgjcDohIzskIiM1ISIiNyQkIjFWRzhwLFRyOiEjOiQhIzUhIiI3JCQiMlotWy9udU5kIiEjOyQhIzUhIiI3JCQiMmo/ajw8UmRkIiEjOyQhIzUhIiI3JCQiMiMpUXlJbi56ZCIhIzskISM1ISIiNyQkIjI6djNkbktBZSIhIzskISM1ISIiNyQkIjE6IlIkeTtjJ2UiISM6JCEjNSEiIjckJCIyPSUpKmYkbz5fZiIhIzskISM1ISIiNyQkIjIpbzAnKSlveFFnIiEjOyQhIzUhIiI3JCQiMkctI1EqcCQ+QDshIzskISM1ISIiNyQkIjJuWi4qNCg0JlE7ISM7JCEjNSEiIjckJCIxLCVHZylwbG87ISM6JCEjNSEiIjckJCIyYEtgQEUvKSlwIiEjOyQhIzUhIiI3JCQiMnhSMS8uQ1J3IiEjOyQhIzUhIiI3JCQiMHQmKUg3Kj5KPSEjOSQhIzUhIiI3JCQiMiRmYjd3WSwoKj0hIzskISM1ISIiNyQkIjJZLmh6ZyVwZz4hIzskISM1ISIiNyQkIjFwSjg6LlNKPyEjOiQhIzUhIiI3JCQiLzJ6UEQkXDQjISM4JCEjNSEiIjckJCIyeCxmaHVtRjsjISM7JCEjNSEiIjckJCIyW1gnM0BZQkNBISM7JCEjNSEiIjckJCIxL2I8V19WIkgjISM6JCEjNSEiIjckJCIwR05NemxZTiMhIzkkISM1ISIiNyQkIjJtT21rJCpmMlUjISM7JCEjNSEiIjckJCIxMSgpPVUhemBbIyEjOiQhIzUhIiI3JCQiMWtISGQjSEliIyEjOiQhIzUhIiI3JCQiMV9yJltYJT09RSEjOiQhIzUhIiI3JCQiMUstJVwwOltvIyEjOiQhIzUhIiI3JCQiMjBjOCsjUiozdiMhIzskISM1ISIiNyQkIjEvMExNTmg2RyEjOiQhIzUhIiI3JCQiMkdvVVgxMDcpRyEjOyQhIzUhIiI3JCQiMiU0NnhlJltNJUghIzskISM1ISIiNyQkIjJMNSIpZTU4KTRJISM7JCEjNSEiIjckJCIyTHghUlhDTHRJISM7JCEjNSEiIjckJCItWF1FZlRKISM2JCEjNSEiIi0lJkNPTE9SRzYmJSRSR0JHJCIjNSEiIiQiIiEhIiIkIiIhISIiLSUlVklFV0c2JDskIiIhISIiJCIrYUVmVEohIiolKERFRkFVTFRHLSUqQVhFU1NUWUxFRzYjJSdOT1JNQUxHLSUoU0NBTElOR0c2IyUuVU5DT05TVFJBSU5FREctJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiQ/JSEiIi0lKUJPVU5EU19ZRzYjJCIkPyIhIiItJS1CT1VORFNfV0lEVEhHNiMkIiVJTiEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiVdUCEiIi0lKUNISUxEUkVORzYi
cossum(f1,50,t);
LFQqJkkjUGlHJSpwcm90ZWN0ZWRHISIiLUkkY29zRzYkRiVJKF9zeXNsaWJHNiI2IywkSSJ0R0YrIiNIIiIiIyIiJUYvKiZGJEYmLUYoNiMsJEYuIiNKRjAjISIlRjcqJkYkRiYtRig2IywkRi4iIiZGMCNGMkY+KiZGJEYmLUYoNiMsJEYuIiNERjAjRjJGRComRiRGJi1GKDYjLCRGLiIjNkYwI0Y5RkoqJkYkRiYtRig2IywkRi4iIzhGMCNGMkZQKiZGJEYmLUYoNiMsJEYuIiIqRjAjRjJGViomRiRGJi1GKDYjLCRGLiIjWEYwI0YyRmZuKiZGJEYmLUYoNiMsJEYuIiIoRjAjRjlGXG8qJkYkRiYtRig2IywkRi4iI1JGMCNGOUZibyomRiRGJi1GKDYjLCRGLiIjVEYwI0YyRmhvKiZGJEYmLUYoNiMsJEYuIiNWRjAjRjlGXnAqJkYkRiYtRig2IywkRi4iI05GMCNGOUZkcComRiRGJi1GKDYjRi5GMEYyKiZGJEYmLUYoNiMsJEYuIiIkRjAjRjlGXXEqJkYkRiYtRig2IywkRi4iIzpGMCNGOUZjcSomRiRGJi1GKDYjLCRGLiIjPEYwI0YyRmlxKiZGJEYmLUYoNiMsJEYuIiM+RjAjRjlGX3IqJkYkRiYtRig2IywkRi4iI0BGMCNGMkZlciomRiRGJi1GKDYjLCRGLiIjQkYwI0Y5RltzKiZGJEYmLUYoNiMsJEYuIiNGRjAjRjlGYXMqJkYkRiYtRig2IywkRi4iI0xGMCNGMkZncyomRiRGJi1GKDYjLCRGLiIjUEYwI0YyRl10KiZGJEYmLUYoNiMsJEYuIiNaRjAjRjlGY3QqJkYkRiYtRig2IywkRi4iI1xGMCNGMkZpdA==
The pattern here is:LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUkjbWlHRiQ2I1EhRic=
f1bis := (t,N) -> (4/Pi)*sum((-1)^k * cos((2*k+1)*t)/(2*k+1),k=0..N);
Zio2JEkidEc2IkkiTkdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCQqJkkjUGlHJSpwcm90ZWN0ZWRHISIiLUkkc3VtRzYkRi1JKF9zeXNsaWJHRiU2JCooKUYuSSJrR0YlIiIiLUkkY29zR0YxNiMqJiwmRjYiIiNGN0Y3Rjc5JEY3RjdGPEYuL0Y2OyIiITklRjciIiVGJUYlRiU=
Just checking:
simplify(f1bis(t,24)-cossum(f1,50,t));
IiIh
plot({cossum(f1,50,t),f1(t)},t=0..Pi,color=black);
NiktJSdDVVJWRVNHNiM3Ym83JCQiIiEhIiIlKnVuZGVmaW5lZEc3JCQiMlk8JjMvMSQqUkAhIz4kIiM1ISIiNyQkIjFcLjwzNycpelUhIz0kIiM1ISIiNyQkIjFDYkQ3PXo+ayEjPSQiIzUhIiI3JCQiMSpwU2pUQShmJikhIz0kIiM1ISIiNyQkIjJcNV5DT2VSRyIhIz0kIiM1ISIiNyQkIjIoUiJvS1tXPnIiISM9JCIjNSEiIjckJCIyKDRBIVxzO3pjIyEjPSQiIzUhIiI3JCQiMiV6aWBtKikpUVUkISM9JCIjNSEiIjckJCIyJj5XISlcTSRlOCYhIz0kIiM1ISIiNyQkIjFmRDJMenhabyEjPCQiIzUhIiI3JCQiMWI0VTphKW8jKSohIzwkIiM1ISIiNyQkIjFOcHgqRypmIUciISM7JCIjNSEiIjckJCIwNiNleEdtXT4hIzokIiM1ISIiNyQkIjJYVCw9M29eaSMhIzwkIiM1ISIiNyQkIjBcRDBbbmtIJCEjOiQiIzUhIiI3JCQiMSI9WmEmeiUpPVIhIzskIiM1ISIiNyQkIjJ3Yld2KW9HalghIzwkIiM1ISIiNyQkIjElMyoqSGJtKEhfISM7JCIjNSEiIjckJCIxUFJyNCkzVCplISM7JCIjNSEiIjckJCIxI1speXlrWXhsISM7JCIjNSEiIjckJCIxcydvY0dvJHpyISM7JCIjNSEiIjckJCIxUTA7Z3IncCZ5ISM7JCIjNSEiIjckJCIxR010PyRbdGApISM7JCIjNSEiIjckJCIxI3AzMGtAST4qISM7JCIjNSEiIjckJCIwQyJcSGVWKXkqISM6JCIjNSEiIjckJCIyR1spKSopM1cnXDUhIzskIiM1ISIiNyQkIjEnW0BYU0AnNDYhIzokIiM1ISIiNyQkIjJWaFAzRypRejYhIzskIiM1ISIiNyQkIjJVMyozdGM5VDchIzskIiM1ISIiNyQkIjIkeTBEQkEhKjM4ISM7JCIjNSEiIjckJCIyUG1LP1tBTVAiISM7JCIjNSEiIjckJCIxYUBYLkV1UzkhIzokIiM1ISIiNyQkIjJRdSMpKipbakRdIiEjOyQiIzUhIiI3JCQiMkM/Kyp5WCFmYCIhIzskIiM1ISIiNyQkIjI3bTx5bVgjcDohIzskIiM1ISIiNyQkIjFWRzhwLFRyOiEjOiQhIzUhIiI3JCQiMlotWy9udU5kIiEjOyQhIzUhIiI3JCQiMmo/ajw8UmRkIiEjOyQhIzUhIiI3JCQiMiMpUXlJbi56ZCIhIzskISM1ISIiNyQkIjI6djNkbktBZSIhIzskISM1ISIiNyQkIjE6IlIkeTtjJ2UiISM6JCEjNSEiIjckJCIyPSUpKmYkbz5fZiIhIzskISM1ISIiNyQkIjIpbzAnKSlveFFnIiEjOyQhIzUhIiI3JCQiMkctI1EqcCQ+QDshIzskISM1ISIiNyQkIjJuWi4qNCg0JlE7ISM7JCEjNSEiIjckJCIxLCVHZylwbG87ISM6JCEjNSEiIjckJCIyYEtgQEUvKSlwIiEjOyQhIzUhIiI3JCQiMnhSMS8uQ1J3IiEjOyQhIzUhIiI3JCQiMHQmKUg3Kj5KPSEjOSQhIzUhIiI3JCQiMiRmYjd3WSwoKj0hIzskISM1ISIiNyQkIjJZLmh6ZyVwZz4hIzskISM1ISIiNyQkIjFwSjg6LlNKPyEjOiQhIzUhIiI3JCQiLzJ6UEQkXDQjISM4JCEjNSEiIjckJCIyeCxmaHVtRjsjISM7JCEjNSEiIjckJCIyW1gnM0BZQkNBISM7JCEjNSEiIjckJCIxL2I8V19WIkgjISM6JCEjNSEiIjckJCIwR05NemxZTiMhIzkkISM1ISIiNyQkIjJtT21rJCpmMlUjISM7JCEjNSEiIjckJCIxMSgpPVUhemBbIyEjOiQhIzUhIiI3JCQiMWtISGQjSEliIyEjOiQhIzUhIiI3JCQiMV9yJltYJT09RSEjOiQhIzUhIiI3JCQiMUstJVwwOltvIyEjOiQhIzUhIiI3JCQiMjBjOCsjUiozdiMhIzskISM1ISIiNyQkIjEvMExNTmg2RyEjOiQhIzUhIiI3JCQiMkdvVVgxMDcpRyEjOyQhIzUhIiI3JCQiMiU0NnhlJltNJUghIzskISM1ISIiNyQkIjJMNSIpZTU4KTRJISM7JCEjNSEiIjckJCIyTHghUlhDTHRJISM7JCEjNSEiIjckJCItWF1FZlRKISM2JCEjNSEiIi0lJ0NVUlZFU0c2IzdqeTckJCIiISEiIiQiMUZlaTd0czc1ISM6NyQkIjIlemlgbSopKVFVJCEjPSQiMCJ5YTEmKSo+KSoqISM6NyQkIjFmRDJMenhabyEjPCQiMUxoTEpsX3gpKiEjOzckJCIxYjRVOmEpbyMpKiEjPCQiMiMqKilbJ1IneUQrIiEjOzckJCIxTnB4KkcqZiFHIiEjOyQiMXdnPidHT0YsIiEjOjckJCIyQ196TzNKY2giISM8JCIxJ1IkcDRfIzQoKiohIzs3JCQiMDYjZXhHbV0+ISM6JCIxO1o6QlRQeCkqISM7NyQkIjFPcGp5IipHPkAhIzskIjFoOyg+KFwjKlwqKiEjOzckJCIyQXcicHphInpHIyEjPCQiMGlCJz4hb2MrIiEjOTckJCIyYDw+LWpHQVAjISM8JCIyMS8uSTg2KywiISM7NyQkIjIlKWVZMnlUbFgjISM8JCIyKSpwOzppOUUsIiEjOzckJCIyOSt1NyRcJjNhIyEjPCQiMmg9dnRlOUksIiEjOzckJCIyWFQsPTNvXmkjISM8JCIyaDltbU9KNiwiISM7NyQkIjFMQ1tKSCpIeiMhIzskIjJkJ28ibyRSPi01ISM7NyQkIjJCWGo2eTwzJ0ghIzwkIjE9WUkiXGdxIioqISM7NyQkIjI7J1IrMS10V0khIzwkIjF2PVQxSWIiKSkqISM7NyQkIjI6WlczalUnR0ohIzwkIjFiYy9AKW9rJykqISM7NyQkIjEiKVxvYl1iN0shIzskIjFETnYmPVNYKCkqISM7NyQkIjBcRDBbbmtIJCEjOiQiMVI+WDI6XC8qKiEjOzckJCIxOGZEKmZpP1gkISM7JCIxPUtbIj53JSkqKiohIzs3JCQiMmNMJyl6cmR3ZyQhIzwkIjJOckF2JypSJTQ1ISM7NyQkIjFaOk54X1gmbyQhIzskIjJrZS5vZSdbNzUhIzs3JCQiMiVlbnJPR0RqUCEjPCQiMjF3VCp5I3pPLCIhIzs3JCQiMCg+MydSXTUlUSEjOiQiMkxIQncwR0csIiEjOzckJCIxIj1aYSZ6JSk9UiEjOyQiMjQkM2UxeC81NSEjOzckJCIxRDpaKW9kKnpTISM7JCIwPnBkWWwsKyIhIzk3JCQiMiRwZVxAdTFUVSEjPCQiMS9rNDBIVSsqKiEjOzckJCIxVCEzISlHQTtLJSEjOyQiMV9FKnlgUCZwKSohIzs3JCQiMTktX2FyPC1XISM7JCIxJz5gWysjUmYpKiEjOzckJCIyZFFLNS1LRlslISM8JCIxYGc4UD54cikqISM7NyQkIjJ3Yld2KW9HalghIzwkIjBGcltaaVshKiohIzo3JCQiMiQqPTNSIW8hKkhaISM8JCIydzlCMSVlUiw1ISM7NyQkIjFAPUY/bl8nKlshIzskIjElWyg+eS5kNjUhIzo3JCQiMVBPWHltJCl6XCEjOyQiMiNlNzotbDc5NSEjOzckJCIxYGFqT205al0hIzskIjInej1NL3VGOTUhIzs3JCQiMW9zIltmY2s5JiEjOyQiMik9Ij0xMXY+LCIhIzs3JCQiMSUzKipIYm0oSF8hIzskIjJrQ09DIz5mMjUhIzs3JCQiMSh6eHI2X2VSJiEjOyQiMTdSOjcnPXomKiohIzs3JCQiMF5jOG9QPmMmISM6JCIxX3k0IipbQmwpKiEjOzckJCIxbmVXai8pXGsmISM7JCIxamU0dExcXCkqISM7NyQkIjFCX2BYSy1HZCEjOyQiMXQpKj1QbEpmKSohIzs3JCQiMSJlQ3ctbTUiZSEjOyQiMDxXJUcmKkgkKikqISM6NyQkIjFQUnI0KTNUKmUhIzskIjFETDpUMCFmJSoqISM7NyQkIjF1RClwQVtcMSchIzskIjImKVxLNypcTDI1ISM7NyQkIjBAXlVrKHlOaSEjOiQiMnRqemZNJ0c6NSEjOzckJCIxRmIpR04yN0snISM7JCIyQykpKXlmNlc6NSEjOzckJCIxWCk+OjFGbVMnISM7JCIyeCUpM3UtL0csIiEjOzckJCIxalQ6cW4vI1wnISM7JCIxeGwoKSk0O3krIiEjOjckJCIxI1speXlrWXhsISM7JCIyO1AuRGNbOCsiISM7NyQkIjFzJ29jR28kenIhIzskIjF2ZktBKD5UJyoqISM7NyQkIjFRMDtncidwJnkhIzskIjFZaEVlaHMkKioqISM7NyQkIjFcQEIwWyxVeiEjOyQiMWpILlVpTj4qKiEjOzckJCIxaFBJXUMxRiEpISM7JCIxPC0qUiNwPmUpKiEjOzckJCIxc2BQJjQ1QDYpISM7JCIxbExuJm9tOSMpKiEjOzckJCIxJClwV1N4Oig+KSEjOyQiMXAjZSd5Pj87KSohIzs3JCQiMTAtZklJRG4kKSEjOyQiMWpeTnUhZS0hKiohIzs3JCQiMUdNdD8kW3RgKSEjOyQiMmMwcHJfJG8wNSEjOzckJCIxT2BxJlsyJD4nKSEjOyQiMiR6UUsiPjNGLCIhIzs3JCQiMVdzbl1tRSwoKSEjOyQiMjF3MXRMU3gsIiEjOzckJCIxXyJcYyJlQSR5KSEjOyQiMWAvLG0rISo+NSEjOjckJCIxZjVpISlcPWwpKSEjOyQiMkxndShwJXAoPTUhIzs3JCQiMXZbYzVMNUghKiEjOyQiMXgtX3onZncrIiEjOjckJCIxI3AzMGtAST4qISM7JCIxamRKKXpYNSIqKiEjOzckJCIwXkpUVFt1RSohIzokIjAyS1MiKVJsJSkqISM6NyQkIjFIVnYoPXY9TSohIzskIjFWeT1UMistKSohIzs3JCQiMVpyUGg+STslKiEjOyQiMW5ZWSQpeilSeSohIzs3JCQiMW0qKipcdEcyXCohIzskIjEoUj0sZk5ieiohIzs3JCQiMS1jQyNHI2VSJyohIzskIjEmXGldU0MlKiopKiEjOzckJCIwQyJcSGVWKXkqISM6JCIyM2VlVE5KaSsiISM7NyQkIjFSLy9payRwKCkqISM7JCIyZUU7eVRCXywiISM7NyQkIjFRJyplJTRQYScqKiEjOyQiMjIhZlhBMFtANSEjOzckJCIyUClRcnNQUjA1ISM7JCIyJlwkPSFbZ3FCNSEjOzckJCIyTiEpb2YkUUM5NSEjOyQiMjhLdUhaXzgtIiEjOzckJCIxVid5QydSJT4uIiEjOiQiMiY0UVtPNDEwNSEjOzckJCIyR1spKSopM1cnXDUhIzskIjBrKGUkelUlUikqISM6NyQkIjIiMyxWYDc5ZDUhIzskIjFQTEBGeidceCohIzs3JCQiMk90cnlUUVkxIiEjOyQiMS0pKj5qWl9TKCohIzs3JCQiMiJmTEojZU5AMiIhIzskIjFIYjIoUXc8dSohIzs3JCQiMlcpXHZZRmp6NSEjOyQiMTYqKjMzVmB6KCohIzs3JCQiMl9CUWMyRlk0IiEjOyQiMWpJPzpweFUqKiEjOzckJCIxJ1tAWFNAJzQ2ISM6JCIyWF1tSSFmdTk1ISM7NyQkIjJCK2ghKlFVJD02ISM7JCIyRkQyciQpPVQtIiEjOzckJCIyIz0wZ3RMMUY2ISM7JCIyZSNmJTMvPyNINSEjOzckJCIyVS5TImVWeU42ISM7JCIyTWkwLGk+Kkc1ISM7NyQkIjItYnpFTTBYOSIhIzskIjJLSSdRdjAyQjUhIzs3JCQiMkBlZTxKWj47IiEjOyQiMTNhalVMWiMqKiohIzs3JCQiMlZoUDNHKlF6NiEjOyQiMSI9eDlFJykqUigqISM7NyQkIjFbIXApSCkzcj0iISM6JCIxdiRHbDdsXG4qISM7NyQkIjI+WysqeSRHWz4iISM7JCIxQVVCR2BVYycqISM7NyQkIjJjIj4keiN6YS03ISM7JCIxJ2ZKXS0iKSopbyohIzs3JCQiMiNcTCdwWm4tQCIhIzskIjEpW0xZNzQpcCgqISM7NyQkIjJtQEVdZDFkQSIhIzskIjIwVVJnamdIKyIhIzs3JCQiMlUzKjN0YzlUNyEjOyQiMjFROHR4IVxINSEjOzckJCIyM0ZmPUM6J1w3ISM7JCIyUS8sWmo1JVE1ISM7NyQkIjJ2WEgxIlszZTchIzskIjJEKFswTS55UzUhIzs3JCQiMV5YLCZmPkJFIiEjOiQiMlsiMzYrckBSNSEjOzckJCIyVmsqUnpWYm03ISM7JCIyJyp6I1s9SCNlLiIhIzs3JCQiMnh0JXlqIip5cTchIzskIjF4XysvPHFJNSEjOjckJCIyNiQpcCJbUi12NyEjOyQiMWlPTClvUFMtIiEjOjckJCIyWj81ZDNqPkgiISM7JCIxOmViZ0U/ISkpKiEjOzckJCIyJHkwREJBISozOCEjOyQiMTd0eWxQKT5jKiEjOzckJCIxUiQpZmJzJ3BKIiEjOiQiMS5xVm9oKnpcKiEjOzckJCIyJyo0WXpHS11LIiEjOyQiMW8vWilIIipSXiohIzs3JCQiMi4nUUg/dDRMOCEjOyQiLyQ+NFJWRGgqISM5NyQkIjFAO2tfQjtUOCEjOiQiMXgoekRRJypHeSohIzs3JCQiMkI5UHRUI0hkOCEjOyQiMik0IkhCMlVPLSIhIzs3JCQiMlBtSz9bQU1QIiEjOyQiMjBhJltLKik+ZzUhIzs3JCQiMigzQ3VVRDMyOSEjOyQiMiMzKCplKEgnWzs1ISM7NyQkIjFhQFguRXVTOSEjOiQiMDorLV8oUVwhKiEjOjckJCIxXHVyWUlscjkhIzokIjFFXHNPISkzPioqISM7NyQkIjJRdSMpKipbakRdIiEjOyQiMiJwIzRdemVAPCIhIzs3JCQiMUV1NEU2dDE6ISM6JCIyelE6TmlLKHk2ISM7NyQkIjInM0BAaSgpKjNeIiEjOyQiMlFcdj1Ya288IiEjOzckJCIxInpFJClSbV1eIiEjOiQiMiNHcj5MU1VsNiEjOzckJCIyS1pUVy5NIz46ISM7JCIydFNzJXB1WFY2ISM7NyQkIjJ5JDNuMSRwdl8iISM7JCIyaCdwODdGQGw1ISM7NyQkIjJDPysqeVghZmAiISM7JCIxayIqW3EhZT1SKiEjOzckJCIxbiZINiYpUlVhIiEjOiQiMXU8dDNcKHptKCEjOzckJCIyPSQqZUw3dkRiIiEjOyQiMXknRyQ0Lz9SYiEjOzckJCIxOU9aZkZ1YzohIzokIjE5U29CJEhETiUhIzs3JCQiMm1IKWUmUjU0YyIhIzskIjJqUzRoMy9XNSQhIzw3JCQiMip5SHFKIXldYyIhIzskIjJOW0ZvJ3AlPiI9ISM8NyQkIjI3bTx5bVgjcDohIzskIjEtLylvT2BVJFwhIzw3JCQiMlotWy9udU5kIiEjOyQhMUVhOkUuRk0pKSEjPDckJCIyIylReUluLnpkIiEjOyQhMkcpMzRgRGBZQSEjPDckJCIyOnYzZG5LQWUiISM7JCExdWd2R2slW2QkISM7NyQkIjE6IlIkeTtjJ2UiISM6JCEyJGYyV2phSFtbISM8NyQkIjI9JSkqZiRvPl9mIiEjOyQhMS56KWVBVCVlciEjOzckJCIyKW8wJykpb3hRZyIhIzskITF3TlQlcCZcYiEqISM7NyQkIjJHLSNRKnAkPkA7ISM7JCExZWpmSCZHYDgiISM6NyQkIjJuWi4qNCg0JlE7ISM7JCEyd2t0UWEnUnQ2ISM7NyQkIjEsJUdnKXBsbzshIzokITI5b3AnUU4tKzUhIzs3JCQiMmBLYEBFLykpcCIhIzskITFOJHBRJ3lXSyEqISM7NyQkIjI6Jyl6aTlrOHQiISM7JCEyIipbMmtCJ1svNSEjOzckJCIyeFIxLy5DUnciISM7JCEyRTxPJzREIVsxIiEjOzckJCIyMkJeTiFHdSF5IiEjOyQhMj5pMlspW15MNSEjOzckJCIyUTEncHc6Yyh6IiEjOyQhMS9tJ2VnTW4mKSohIzs3JCQiMi9bb0snNChmIT0hIzskITEnRy8mKXkzIWYnKiEjOzckJCIxKCozJSlcLlE5PSEjOiQhMWo9VztocEsmKiEjOzckJCIyUEo4a3QqeUE9ISM7JCExJEgkSHBgSSVcKiEjOzckJCIwdCYpSDcqPko9ISM5JCExN2RjdyNcWWEqISM7NyQkIjJ0b3E3LGB3JT0hIzskITFHWStPUnJOKSohIzs3JCQiMlprYiYqKm81az0hIzskITJmJFIoZithJj41ISM7NyQkIjJNNylwVlFMcz0hIzskITJqKj4qenRURi4iISM7NyQkIjEtMSV5eWcwKT0hIzokITIkXHpjdigpKSlSNSEjOzckJCIyMTMkKT50KHkpKT0hIzskITI3dnVaUlMsLyIhIzs3JCQiMiRmYjd3WSwoKj0hIzskITJ2QlghNG16TDUhIzs3JCQiMiRHVzNmWSRIIj4hIzskITImM0FRUmddMzUhIzs3JCQiMShIVj9rYSlHPiEjOiQhMUciKVt2TVMvKSohIzs3JCQiMjl0QU5qOW8kPiEjOyQhMVJcRWAnKW8yKCohIzs3JCQiMmU7LV1pdVolPiEjOyQhMUt4KGUqKjMqZicqISM7NyQkIjItZyJbO1l0Xz4hIzskITFsOWZyRlRtJyohIzs3JCQiMlkuaHpnJXBnPiEjOyQhMUZ0UjomcFNzKiEjOzckJCIyI29Tdk01UHk+ISM7JCExaUIwbDsvcyoqISM7NyQkIjI9NVo6WVpnKj4hIzskITJVPVpqJlwuQTUhIzs3JCQiMic9TyVcbiYpWysjISM7JCEyJWU9eTBoYkc1ISM7NyQkIjJhOFMkKSlRczg/ISM7JCEyWEdDP2NHJUg1ISM7NyQkIjFfbXQsQGNBPyEjOiQhMVJUI2Ymb3BDNSEjOjckJCIxcEo4Oi5TSj8hIzokITJLJ29ZRmhYOjUhIzs3JCQiMV92enFMR1o/ISM6JCExXmh0Qm5dUSoqISM7NyQkIjJYJD5ZRWs7aj8hIzskITE6Om9lMjZyKCohIzs3JCQiMUVUSGF6NXI/ISM6JCExV1BNbF1KUSgqISM7NyQkIjE8ajcjW1whej8hIzokITBMZmNnOm51KiEjOjckJCIyJjMmZSo0NSpwMyMhIzskITE3cHlPXyZReiohIzs3JCQiLzJ6UEQkXDQjISM4JCExV1tlUlNVcikqISM7NyQkIjIlekYpKSozIio9NiMhIzskITJWZVV2Qzt4KyIhIzs3JCQiMWZbKD5rXClHQCEjOiQhMlVuZDwycj8tIiEjOzckJCIyJikqMy09Kkd0OCMhIzskITJMeGwlM1FuQjUhIzs3JCQiMiNRcDElPjNlOSMhIzskITJsISp6Jz56M0A1ISM7NyQkIjF5SDZxdUdhQCEjOiQhL0cmKm9WKFssIiEjODckJCIyeCxmaHVtRjsjISM7JCEyRypmOys9QDE1ISM7NyQkIjF4MypbckwieUAhIzokITE6KnAseGZYKikqISM7NyQkIjJrdEFPbytOPiMhIzskITFMI0dlMztCeiohIzs3JCQiMW0nKSl6OyU9LEEhIzokITE6IyopKWYpXFx5KiEjOzckJCIxJ2ZhQmxuKTNBISM6JCExWSc0SycqPikzKSohIzs3JCQiMmJfP244XmxAIyEjOyQhMVomZiJ5eSQpZikqISM7NyQkIjJbWCczQFlCQ0EhIzskITFcKz9NPTRJKiohIzs3JCQiMTwoZW94TTVDIyEjOiQhMksoR05KXiIpNDUhIzs3JCQiMiZ6NGpLXCR5RCMhIzskITIoeVRsIW9AJj41ISM7NyQkIjIyNjwwLE5pRSMhIzskIS9ZQSt6Wz41ISM4NyQkIjFVS1MpM05ZRiMhIzokITJ5KD1gUmU0OzUhIzs3JCQiMkxQKkdtXi4kRyMhIzskITI2WlBDWiYqKjQ1ISM7NyQkIjEvYjxXX1YiSCMhIzokITJ2TmpoOypHLTUhIzs3JCQiMEdOTXpsWU4jISM5JCExaj1HW1RsIykqKiEjOzckJCIybU9tayQqZjJVIyEjOyQhMFB6UDMpSCkpKiohIzo3JCQiMTEoKT1VIXpgWyMhIzokITJQQTIqPSJwRCsiISM7NyQkIjFRbjI+YCRRXCMhIzokITJMbCgpKVtQejM1ISM7NyQkIjIweGtmZiJILUQhIzskITIwU2klUSEqUjg1ISM7NyQkIjJHIUcmRyh5dTVEISM7JCExOz5FJXAqZjo1ISM6NyQkIjFOM3VcVD8+RCEjOiQhMXcqNCpbLy86NSEjOjckJCIyKCoqb14ubjZPRCEjOyQhMjN4RyUzZmkxNSEjOzckJCIxa0hIZCNISWIjISM6JCExS01iMFMiKlIqKiEjOzckJCIydFspKXBsdDZjIyEjOyQhMWBnJGY7eScqKSkqISM7NyQkIjE2U29jIT0kcEQhIzokITEmMyY9JUdPeiYpKiEjOzckJCIyVmB6alhpdWQjISM7JCExVlUza2RsXCkqISM7NyQkIjFlXTJjb2cmZSMhIzokITE+RVxmMSVmJykqISM7NyQkIjEwaFliYyo9ZyMhIzokITFBeSNIL0ByJioqISM7NyQkIjFfciZbWCU9PUUhIzokITJgYjo4YT90KyIhIzs3JCQiMVB2Jyl6S15FRSEjOiQhMkhubypHS3o2NSEjOzckJCIxQXooWzVVW2ojISM6JCExQTZFaGlAOTUhIzo3JCQiMTIkKSkpSDQ8VkUhIzokITIkRyw/KWUnPjk1ISM7NyQkIjEjcCkqW3YqXF5FISM6JCEyLixAMyMpZjwsIiEjOzckJCIxaSU+XFNkIm9FISM6JCEyOHgheltBciw1ISM7NyQkIjFLLSVcMDpbbyMhIzokITFGdm8vZkEyKiohIzs3JCQiMndSdUkiXDIkcCMhIzskITB5UEElKipbcykqISM6NyQkIjFrJjM3eE04cSMhIzokITEmeVBzJ3lSZikqISM7NyQkIjB0VSRIWWY0RiEjOSQhMSlSKlFaQCgqcCkqISM7NyQkIjEnKm9aKFtheXIjISM6JCExLmI/QjFFLSoqISM7NyQkIjImR191LlVQTUYhIzskITIlZU8xPVNxKzUhIzs3JCQiMjBjOCsjUiozdiMhIzskITJPQUBEOyRlNTUhIzs3JCQiMXlKIT0oUVtlRiEjOiQhMkRmclg8XUksIiEjOzckJCIxJ3ojZkJRMm1GISM6JCEyWlQuQnpcTywiISM7NyQkIjJWVCNRdlBtdEYhIzskITIpKW9kYUMzQiwiISM7NyQkIjJCLnNyc2A3eSMhIzskITAmUjlic0E0NSEjOTckJCIyJm83dklPVid6IyEjOyQhMWptOiNvYyIpKioqISM7NyQkIjEvMExNTmg2RyEjOiQhMSI+MTZpb2chKiohIzs3JCQiMUVxZ3ZDSj9HISM6JCExPF02I3ooZnUpKiEjOzckJCIyJltOKW9UNiFIRyEjOyQhMSY+KEc6WmdtKSohIzs3JCQiMXIrO2UuclBHISM6JCExQmM5biVSTSkpKiEjOzckJCIyTGZPJSpINGslRyEjOyQhMTJlX3coUj0jKiohIzs3JCQiMVEnKik+PTJRJ0chIzokITIoR0QxTlM8LjUhIzs3JCQiMkdvVVgxMDcpRyEjOyQhMmtbdm1XNj0sIiEjOzckJCIxTzdLLGIpKikpRyEjOiQhMmVdZHRgN0osIiEjOzckJCIxKnkqNFFmdycqRyEjOiQhMTMxIVEhPlg3NSEjOjckJCIyQ015W1BZWCFIISM7JCExI1s5eHhOKjQ1ISM6NyQkIjEnKm9sNm9LN0ghIzokITElXChRJFJZZisiISM6NyQkIjEuU0AmbygpeSNIISM6JCExPmMmNCopbzInKiohIzs3JCQiMiU0NnhlJltNJUghIzskITFzazRtdEMmKSkqISM7NyQkIjFmKWZyN1c8JkghIzokIS9GY1hAdHEpKiEjOTckJCIxMydbYnBSKydIISM6JCExJSkqKVtcLjt5KSohIzs3JCQiMWR0JFJFTiRvSCEjOiQhMD5FImZqPzEqKiEjOjckJCIyazVFQiQzandIISM7JCExYVAhelRoKyYqKiEjOzckJCIyWWcuInA+QSQqSCEjOyQhMkJQdF9JJFIwNSEjOzckJCIyTDUiKWU1OCk0SSEjOyQhMndhQSRHRUI3NSEjOzckJCIxUWZqdkZkVEkhIzokITJsdVVJN2xPKyIhIzs3JCQiMkx4IVJYQ0x0SSEjOyQhMCZIJj0wWnIoKSohIzo3JCQiMmpRP3phaXU1JCEjOyQhMFInSEFbbyMpKiohIzo3JCQiLVhdRWZUSiEjNiQhMUZlaTd0czc1ISM6LSUlVklFV0c2JDskIiIhISIiJCIrYUVmVEohIiolKERFRkFVTFRHLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpBWEVTU1RZTEVHNiMlJ05PUk1BTEctJShTQ0FMSU5HRzYjJS5VTkNPTlNUUkFJTkVERy0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJD8lISIiLSUpQk9VTkRTX1lHNiMkIiMhKiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJUlOISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJTVRISIiLSUpQ0hJTERSRU5HNiI=
Functions from fig 9.2.4 and 9.3.10
f2 := t -> t*unit(t,0,Pi/2) + (Pi-t)*unit(t,Pi/2,Pi);
Zio2I0kidEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCYqJjkkIiIiLUkldW5pdEdGJTYlRisiIiEsJEkjUGlHJSpwcm90ZWN0ZWRHI0YsIiIjRixGLComLCZGMkYsRishIiJGLC1GLjYlRitGMUYyRixGLEYlRiVGJQ==
plot(f2(t),t=0..Pi);
NictJSdDVVJWRVNHNiQ3aW43JCQiIiEhIiIlKnVuZGVmaW5lZEc3JCQiMlk8JjMvMSQqUkAhIz4kIjJZPCYzLzEkKlJAISM+NyQkIjFcLjwzNycpelUhIz0kIjFcLjwzNycpelUhIz03JCQiMUNiRDc9ej5rISM9JCIxQ2JENz16PmshIz03JCQiMSpwU2pUQShmJikhIz0kIjEqcFNqVEEoZiYpISM9NyQkIjJcNV5DT2VSRyIhIz0kIjJcNV5DT2VSRyIhIz03JCQiMihSIm9LW1c+ciIhIz0kIjIoUiJvS1tXPnIiISM9NyQkIjIoNEEhXHM7emMjISM9JCIyKDRBIVxzO3pjIyEjPTckJCIyJXppYG0qKSlRVSQhIz0kIjIlemlgbSopKVFVJCEjPTckJCIyJj5XISlcTSRlOCYhIz0kIjImPlchKVxNJGU4JiEjPTckJCIxZkQyTHp4Wm8hIzwkIjFmRDJMenhabyEjPDckJCIxYjRVOmEpbyMpKiEjPCQiMWI0VTphKW8jKSohIzw3JCQiMU5weCpHKmYhRyIhIzskIjFOcHgqRypmIUciISM7NyQkIjA2I2V4R21dPiEjOiQiMDYjZXhHbV0+ISM6NyQkIjJYVCw9M29eaSMhIzwkIjJYVCw9M29eaSMhIzw3JCQiMFxEMFtua0gkISM6JCIwXEQwW25rSCQhIzo3JCQiMSI9WmEmeiUpPVIhIzskIjEiPVphJnolKT1SISM7NyQkIjJ3Yld2KW9HalghIzwkIjJ3Yld2KW9HalghIzw3JCQiMSUzKipIYm0oSF8hIzskIjElMyoqSGJtKEhfISM7NyQkIjFQUnI0KTNUKmUhIzskIjFQUnI0KTNUKmUhIzs3JCQiMSNbKXl5a1l4bCEjOyQiMSNbKXl5a1l4bCEjOzckJCIxcydvY0dvJHpyISM7JCIxcydvY0dvJHpyISM7NyQkIjFRMDtncidwJnkhIzskIjFRMDtncidwJnkhIzs3JCQiMUdNdD8kW3RgKSEjOyQiMUdNdD8kW3RgKSEjOzckJCIxI3AzMGtAST4qISM7JCIxI3AzMGtAST4qISM7NyQkIjBDIlxIZVYpeSohIzokIjBDIlxIZVYpeSohIzo3JCQiMkdbKSkqKTNXJ1w1ISM7JCIyR1spKSopM1cnXDUhIzs3JCQiMSdbQFhTQCc0NiEjOiQiMSdbQFhTQCc0NiEjOjckJCIyVmhQM0cqUXo2ISM7JCIyVmhQM0cqUXo2ISM7NyQkIjJVMyozdGM5VDchIzskIjJVMyozdGM5VDchIzs3JCQiMiR5MERCQSEqMzghIzskIjIkeTBEQkEhKjM4ISM7NyQkIjJQbUs/W0FNUCIhIzskIjJQbUs/W0FNUCIhIzs3JCQiMWFAWC5FdVM5ISM6JCIxYUBYLkV1UzkhIzo3JCQiMlF1IykqKltqRF0iISM7JCIyUXUjKSoqW2pEXSIhIzs3JCQiMkM/Kyp5WCFmYCIhIzskIjJDPysqeVghZmAiISM7NyQkIjI3bTx5bVgjcDohIzskIjI3bTx5bVgjcDohIzs3JCQiMilvMCcpKW94UWciISM7JCIyVkFIWidcclA6ISM7NyQkIjJuWi4qNCg0JlE7ISM7JCIybEonb1ZIMy46ISM7NyQkIjJgS2BARS8pKXAiISM7JCIyellPOVIpeVU5ISM7NyQkIjJ4UjEvLkNSdyIhIzskIjJhUiQ9Qidvd1AiISM7NyQkIjB0JilINyo+Sj0hIzkkIjFqU2dJTlI1OCEjOjckJCIyJGZiN3dZLCgqPSEjOyQiMlFCa3UoemRXNyEjOzckJCIyWS5oemclcGc+ISM7JCIyJ2UoR2MvKSozPSIhIzs3JCQiMXBKODouU0o/ISM6JCIyVGljJVFCPjU2ISM7NyQkIi8yelBEJFw0IyEjOCQiMko0Kno6LG1ZNSEjOzckJCIyeCxmaHVtRjsjISM7JCIxYnhJdSFmIyl5KiEjOzckJCIyW1gnM0BZQkNBISM7JCIxJVFMXUshZXQiKiEjOzckJCIxL2I8V19WIkgjISM6JCIwKkc5JTR1Ol0pISM6NyQkIjBHTk16bFlOIyEjOSQiMUteYSwnbyNweSEjOzckJCIybU9tayQqZjJVIyEjOyQiMW1VQnJySzNzISM7NyQkIjExKCk9VSF6YFsjISM6JCIxczMsOWg4aWwhIzs3JCQiMWtISGQjSEliIyEjOiQiMSRIb0gnUmomKWUhIzs3JCQiMV9yJltYJT09RSEjOiQiMThrSygpPjNNXyEjOzckJCIxSy0lXDA6W28jISM6JCIyRGgmXCcpZnhuWCEjPDckJCIyMGM4KyNSKjN2IyEjOyQiMmtLaWRMKClwIVIhIzw3JCQiMS8wTE1OaDZHISM6JCIxKilHZiM+InoqSCQhIzs3JCQiMkdvVVgxMDcpRyEjOyQiMk41ci8qZShRZyMhIzw3JCQiMiU0NnhlJltNJUghIzskIjJuJG89WzRXIik+ISM8NyQkIjJMNSIpZTU4KTRJISM7JCIyJCkqbzN4YXo8OCEjPDckJCIyTHghUlhDTHRJISM7JCIxej4hKj4zLUVvISM8NyQkIi1YXUVmVEohIzYkIjF2SEI5THpSSiEjQy0lJkNPTE9SRzYmJSRSR0JHJCIjNSEiIiQiIiEhIiIkIiIhISIiLSUlVklFV0c2JDskIiIhISIiJCIrYUVmVEohIiolKERFRkFVTFRHLSUqQVhFU1NUWUxFRzYjJSdOT1JNQUxHLSUoU0NBTElOR0c2IyUuVU5DT05TVFJBSU5FREctJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiRJJCEiIi0lKUJPVU5EU19ZRzYjJCIjISohIiItJS1CT1VORFNfV0lEVEhHNiMkIiU/TyEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiUhXCQhIiItJSlDSElMRFJFTkc2Ig==
f3 := t -> (Pi/2-t); plot(f3(t),t=0..Pi);
Zio2I0kidEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCZJI1BpRyUqcHJvdGVjdGVkRyMiIiIiIiM5JCEiIkYlRiVGJQ==
NictJSdDVVJWRVNHNiQ3UzckJCIiISEiIiQiMm0qW3pFanpxOiEjOzckJCIxZkQyTHp4Wm8hIzwkIjFyVFlaJj1CXSIhIzo3JCQiMU5weCpHKmYhRyIhIzskIjEucyJ5Uk9GVyIhIzo3JCQiMDYjZXhHbV0+ISM6JCIyY29PIVIrdHY4ISM7NyQkIjJYVCw9M29eaSMhIzwkIjFiWmg9JnojMzghIzo3JCQiMFxEMFtua0gkISM6JCIydU1VKHkmXDZDIiEjOzckJCIxIj1aYSZ6JSk9UiEjOyQiMiR5LERKOiIqeTYhIzs3JCQiMndiV3Ypb0dqWCEjPCQiMjNXUyFRd1k5NiEjOzckJCIxJTMqKkhibShIXyEjOyQiMiIpKVxccic+eS8iISM7NyQkIjFQUnI0KTNUKmUhIzskIjFHXUJlVyZRIikqISM7NyQkIjEjWyl5eWtZeGwhIzskIjElW2ciKnknXEkiKiEjOzckJCIxcydvY0dvJHpyISM7JCIxJEghRyMpXGZHJikhIzs3JCQiMVEwO2dyJ3AmeSEjOyQiMUYlKXkyaCo0JnkhIzs3JCQiMUdNdD8kW3RgKSEjOyQiMVFiQFpcaHFyISM7NyQkIjEjcDMwa0BJPiohIzskIjF1LVdGOyVcXichIzs3JCQiMEMiXEhlVil5KiEjOiQiMUV4WFF1Xz5mISM7NyQkIjJHWykpKikzVydcNSEjOyQiMVBUMXlCXzZfISM7NyQkIjEnW0BYU0AnNDYhIzokIjEwVHRBI1w8aCUhIzs3JCQiMlZoUDNHKlF6NiEjOyQiMkUjR2RmLzI5UiEjPDckJCIyVTMqM3RjOVQ3ISM7JCIyUDdlcWAxbEgkISM8NyQkIjIkeTBEQkEhKjM4ISM7JCIxJD1WYS5UKj1FISM7NyQkIjJQbUs/W0FNUCIhIzskIjImR0JpWiVRUCg+ISM8NyQkIjFhQFguRXVTOSEjOiQiMmRVRk1CUDBJIiEjPDckJCIyUXUjKSoqW2pEXSIhIzskIjF2X0AibyRHQm8hIzw3JCQiMjdtPHltWCNwOiEjOyQiMmpRTnMoKmUxYiIhIz43JCQiMm5aLio0KDQmUTshIzskITAsZTNKUTh4JyEjOzckJCIyYEtgQEUvKSlwIiEjOyQhMShHJWVgJHorRyIhIzs3JCQiMnhSMS8uQ1J3IiEjOyQhMTZdNk9xRko+ISM7NyQkIjB0JilINyo+Sj0hIzkkITFOJDM+J3otL0UhIzs3JCQiMiRmYjd3WSwoKj0hIzskITJ1aTFMXCQ9aUshIzw3JCQiMlkuaHpnJXBnPiEjOyQhMFFoOyJHKSopKlEhIzo3JCQiMXBKODouU0o/ISM6JCEyVnMjUSQpKVJnZyUhIzw3JCQiLzJ6UEQkXDQjISM4JCExTSFlKjRAT1RfISM7NyQkIjJ4LGZodW1GOyMhIzskITE2N2skPi8oPmYhIzs3JCQiMltYJzNAWUJDQSEjOyQhMSNlOkglSFFNbCEjOzckJCIxL2I8V19WIkgjISM6JCExd2chUTwqUTFzISM7NyQkIjBHTk16bFlOIyEjOSQhMU1RU21ZcFF5ISM7NyQkIjJtT21rJCpmMlUjISM7JCEvWnInNE8nKlwpISM5NyQkIjExKCk9VSF6YFsjISM6JCExJDRRUjpGZTkqISM7NyQkIjFrSEhkI0hJYiMhIzokITF0MSlcSUhCIykqISM7NyQkIjFfciZbWCU9PUUhIzokITJfRGkhRyIpUVo1ISM7NyQkIjFLLSVcMDpbbyMhIzokITJgTFgiRyg9UzYiISM7NyQkIjIwYzgrI1IqM3YjISM7JCExayc9S2YoNCE9IiEjOjckJCIxLzBMTU5oNkchIzokITJ3Z052PzwzQyIhIzs3JCQiMkdvVVgxMDcpRyEjOyQhMml5Wnh0My9KIiEjOzckJCIyJTQ2eGUmW00lSCEjOyQhMkhAdz5CX0VQIiEjOzckJCIyTDUiKWU1OCk0SSEjOyQhMm4/JzN6bixSOSEjOzckJCIyTHghUlhDTHRJISM7JCEybyhlZj1oYC06ISM7NyQkIi1YXUVmVEohIzYkITJLNWJPSyd6cTohIzstJSZDT0xPUkc2JiUkUkdCRyQiIzUhIiIkIiIhISIiJCIiISEiIi0lJVZJRVdHNiQ7JCIiISEiIiQiK2FFZlRKISIqJShERUZBVUxURy0lKkFYRVNTVFlMRUc2IyUnTk9STUFMRy0lKFNDQUxJTkdHNiMlLlVOQ09OU1RSQUlORURHLSUlUk9PVEc2Jy0lKUJPVU5EU19YRzYjJCIkPyUhIiItJSlCT1VORFNfWUc2IyQiIyEqISIiLSUtQk9VTkRTX1dJRFRIRzYjJCIlSU4hIiItJS5CT1VORFNfSEVJR0hURzYjJCIlNVEhIiItJSlDSElMRFJFTkc2Ig==
The EVEN extension (fig 9.3.10)
plot(cossum(f3,50,t),t=-2*Pi..3*Pi); cossum(f3,50,t);
NictJSdDVVJWRVNHNiQ3X3I3JCQhLWklSCY9JEcnISM2JCIyTXldIUdjMWU6ISM7NyQkITElSGtrJTMqPjYnISM6JCIyM104KDMzJzNTIiEjOzckJCExKWUzJClSJ3pTZiEjOiQiMiU+dSkqemNvRjchIzs3JCQhMVsrTUM1JT16JiEjOiQiMjxiY1U3SyJ6NSEjOzckJCEydl1yLmwmKUdrJiEjOyQiMXklPXVQcXJJKiEjOzckJCExV3dlYChwYFomISM6JCIvOShwOmI5aighIzk3JCQhMSN5Lm8mUSZ5SSYhIzokIjEmengkKik9ZV5mISM7NyQkITEoXDxmYkYjUl4hIzokIjJtQkl3IlEucVUhIzw3JCQhMTY3LmI3Z3FcISM6JCIwWGYmemsuJGUjISM6NyQkITFLel4wa3gtWyEjOiQiMXojeT4nNG04ISohIzw3JCQhMWBZK2M6JlxqJSEjOiQhMS85VSk+bW1zKCEjPDckJCExWihIdVZjJHpXISM6JCEydE9RInkkcCNHQiEjPDckJCEyMCVbJik9OHdCViEjOyQhMSN5Mj4nUSN5KVEhIzs3JCQhMUM6KmZlXkU7JSEjOiQhMThlLjVWVipcJiEjOzckJCExMyNHSiY9YSxTISM6JCExQTskM0Q7YzUoISM7NyQkITE3NiRwJD4jXCRRISM6JCExam94RXlVdSgpISM7NyQkITJjLE0yLS0kb08hIzskITJbYEtqWVZYLyIhIzs3JCQhMS44c2NrQC1OISM6JCEyak4qRyc9Lyc0NyEjOzckJCEyOGYzRiozOE9MISM7JCEyb0ZaLm07aFAiISM7NyQkITIoXCdmVD02MkQkISM7JCEyY2slXFZadWY5ISM7NyQkITEzMmh2OUhsSiEjOiQhMjt6WClHZVlcOiEjOzckJCExbjwxbjwoKXpJISM6JCEyKT1cWHokZjJeIiEjOzckJCExREdeZT9YJSpIISM6JCEyMz0qXCZmaltVIiEjOzckJCExd0olcGd3UiVHISM6JCEyTDckZXpsKFFGIiEjOzckJCEyeF90YDosTnAjISM7JCExQChSdjE6RDciISM6NyQkITF2KT5wVixUXyMhIzokITElNGlTVSllSSYqISM7NyQkITFBaVk9PHFhQiEjOiQhMVVITjN4alV5ISM7NyQkITIoM0lcR2tnJT0jISM7JCEwZSI0c1c1T2ghIzo3JCQhMmF6PiZRNl45PyEjOyQhMmw/L0FgXmxWJSEjPDckJCExOilSKDNHZl09ISM6JCEyVydleUZrZStHISM8NyQkITJYJClmKnlXbidvIiEjOyQhMmYwVSgpZillZDYhIzw3JCQhMmMoSCc9JDQjeWAiISM7JCIxVm1wbXchSEskISM8NyQkITJvNm1aUW4qKVEiISM7JCIxWj4oUWQ8Ij49ISM7NyQkITFKWSUpPmgnPkAiISM6JCIyZGt5NlInPSdlJCEjPDckJCEyXzlCXCZbJ1wuIiEjOyQiMDoxSW1fN08mISM6NyQkITFWXVNpY0BdKSkhIzskIjEqSFFjJmVeZW8hIzs3JCQhMU4neWJ4I3lddCEjOyQiMShRS2huPFNOKSEjOzckJCExVig9JXBlZTFjISM7JCIyRHZAXWQhRzU1ISM7NyQkITI6JillSycqKVFpUSEjPCQiMi8icDQpeSd6JT0iISM7NyQkITFLX15lInolPUIhIzskIjJzYzQzK3F6TCIhIzs3JCQhMURoclBOcFh4ISM8JCIxXidcdztPQ1wiISM6NyQkITInemAmbzo0NF4kISM9JCIydyZcKW9tRSthIiEjOzckJCIxZk8wU0F2UXMhIz0kIjJiKy0iKm9Nc2IiISM7NyQkIjI4Nm1bZ2YnZVwhIz0kIjIoM10hcDZfW18iISM7NyQkIjFkb3MmKVJXJD4qISM8JCIyYUwsOnh3bVoiISM7NyQkIjJbJFt1dTdJbTwhIzwkIjEzVCU+OmpeUiIhIzo3JCQiMSUpcCI0OmVLaCMhIzskIjJKKipcYnExKTQ4ISM7NyQkIjIoSHAmKT1ZRUVVISM8JCIyJWZyYzAmcCdbNiEjOzckJCIxd296JzNyI1JlISM7JCIxMWtmSkE/bCkqISM7NyQkIjElKT5nP1RGQXYhIzskIjFFLnMzbT8mPSkhIzs3JCQiMSM0Mlc6eF8/KiEjOyQiMUVYOSplcWJdJyEjOzckJCIyLXI3OyQqel0yIiEjOyQiMSQqR1wqM1NpJlwhIzs3JCQiMjdyJXlaQGpINyEjOyQiMjokM0Q3MSs0TSEjPDckJCIyUVEwQWZQalIiISM7JCIxaiwvdDJDWTwhIzs3JCQiMmswRW0uVkljIiEjOyQiMCVbQyFmOSZceSEjPDckJCIyR19uOzkuaXQiISM7JCEydCZ6RHkkPmtsIiEjPDckJCIyIyopKjNuQ2okND4hIzskITJWUjFsSUNKUSQhIzw3JCQiMV5IPUYnKjRnPyEjOiQhMmRJKipIRylIIipbISM8NyQkIjE4cGwyZyQzQCMhIzokITFmJnkqSG1BLWshIzs3JCQiMnMobzdHYWp0QiEjOyQhMTA5IVwwRDEuKSEjOzckJCIyOiVvZltbVk9EISM7JCExO3EqPloyQWwqISM7NyQkIjIncHEoKXpEaS9GISM7JCEyJyplc3hhZ1I4IiEjOzckJCIyd0hkNko1RyhHISM7JCEyRClvUk1Day04ISM7NyQkIjFjKSpcXyh6XSZIISM6JCEyODVwaC9xWVEiISM7NyQkIjE6QyVRPlx0LiQhIzokITJaKj4mZmx2VlkiISM7NyQkIjJXcDhYIlJbeUkhIzskITI4WyEqSDFCIjQ6ISM7NyQkIjJOKFw9Tic9Jz5KISM7JCEyTmFBKHpAal06ISM7NyQkIjJMaD9iKmY9U0ohIzskITJOUyZmU1QuZTohIzs3JCQiMkVEY2VOYDI7JCEjOyQhMTYlPWB1b0JiIiEjOjckJCIyQyo9PjsySyI9JCEjOyQhMjJNOmhoIlxOOiEjOzckJCIxS3ZfdyEpKT0/JCEjOiQhMmRnaUkoKVFDXiIhIzs3JCQiMSxyJmYhejNoTCEjOiQhMShwdUFBJz1dOCEjOjckJCIyJnBtUU54Rz9OISM7JCEyYHpsLjojNCM+IiEjOzckJCIxSi45Lj8wKHAkISM6JCEyLmFnaEMlZjo1ISM7NyQkIjJCKlIqM0Y7UShRISM7JCExYko6SCUqRyNRKSEjOzckJCIwMnp0I29rS1MhIzkkITFcWmdEaTQqeichIzs3JCQiMnY5a1FReDk+JSEjOyQhMVAqei1NdjlAJiEjOzckJCIyRFc6WSFIMWhWISM7JCEyLiZvZklebjVOISM8NyQkIjFRbk9EJVsxYCUhIzokITJPPSg9PFhDPT0hIzw3JCQiMVI8YDciPVZvJSEjOiQhMnYnZTomeSE9S0chIz03JCQiMHUncCp6KCl6JFshIzkkIjJWal1veSMqZkQiISM8NyQkIjFaOURkJCopZismISM6JCIxIm9FPCg9RVFIISM7NyQkIjFgaCFbIjQqUjwmISM6JCIyTU1aXlZ4TWglISM8NyQkIjEoZih6KEduP0wmISM6JCIxIjRoI1JxYiY+JyEjOzckJCIwLyp5Z085IVwmISM5JCIxdiNbTV4uNHkoISM7NyQkIjJiZywjPSF6YGwmISM7JCIxaThVYlRiSCUqISM7NyQkIjFyVGh2Vmg/ZSEjOiQiMjFxViRRPHQyNiEjOzckJCIxYCZlKFJAOyMpZiEjOiQiMicqWyNcVSRvLkYiISM7NyQkIjFOSCFSISo0UDknISM6JCIyQShmOURWL0s5ISM7NyQkIjFQM0kpRyIqZj0nISM6JCIyMTYjKVxQLThaIiEjOzckJCIxUSgpcHNFRkdpISM6JCIyRXZ0RXRVKD06ISM7NyQkIjEqbygqW084JVxpISM6JCIyUXUmRyZccjdhIiEjOzckJCIxUm00ZFNicWkhIzokIjJaT0djQmFiYiIhIzs3JCQiMGYmSFxacCJIJyEjOSQiMkUpeU9cKj1wYiIhIzs3JCQiMGElXFRhJEdKJyEjOSQiMm9GdDxCIipbYSIhIzs3JCQiMVUuSDUjKVIoUichIzokIjJEczNxKzldWCIhIzs3JCQiMVhoM3o0Jz5bJyEjOiQiMiQzV2o0MGlyOCEjOzckJCIwIlEkRygqW1trJyEjOSQiMmk2ZzA0KmUzNyEjOzckJCIxdTllbXB0Mm8hIzokIjIzKFsiPTZ4bS8iISM7NyQkIjEjb0FtWThWKHAhIzokIjByZEU+eG96KSEjOjckJCIxIipRbW0qKikzOSghIzokIjFmZCRRd3d5NyghIzs3JCQiMSo0Iz1IcjMxdCEjOiQiMSpSZyIqZiM+IlsmISM7NyQkIjExLnEiSCVHcnUhIzokIjEsKj5oLSh5RlEhIzs3JCQiMVc0TUZMM0J3ISM6JCIyLTw4aHpubkkjISM8NyQkIjI6ZSIpSE8jKVt4KCEjOyQiMTAoKnk6S0wiKnkhIzw3JCQiMVAiUSQpPWgpW3ohIzokITFQTDdcRFVpJSohIzw3JCQiMSNwJXA4KyVHNykhIzokITBBMlRgSDFwIyEjOjckJCIxSC02XChbJXkjKSEjOiQhMGBHaGhzbEMlISM6NyQkIjFtZF8lW2RTVikhIzokITE0Inl3IVtuKXomISM7NyQkIjFCXDhfKW8qKmYpISM6JCExeiEpeSM+VydldSEjOzckJCIwM1coPi0pZXcpISM5JCExVWdeUzg2QiIqISM7NyQkIjFSJmYkbyZ5WSMqKSEjOiQhMmp2KXBgP2ZxNSEjOzckJCIxKSpcKHAicFokMyohIzokITI4MkVSUiR5RzchIzs3JCQiMSpcXyVIdTdhIyohIzokITJzKTRlN25VLDkhIzs3JCQiLSQ+JXp4QyUqISM2JCEyKnoyMEdjMWU6ISM7LSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJSVWSUVXRzYkOyQhKzNgPSRHJyEiKiQiK2l6eEMlKiEiKiUoREVGQVVMVEctJSpBWEVTU1RZTEVHNiMlJ05PUk1BTEctJShTQ0FMSU5HRzYjJS5VTkNPTlNUUkFJTkVERy0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiI3EhIiItJSlCT1VORFNfWUc2IyQiIyEqISIiLSUtQk9VTkRTX1dJRFRIRzYjJCIlISlRISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJTVRISIiLSUpQ0hJTERSRU5HNiI=
LFQqJkkjUGlHJSpwcm90ZWN0ZWRHISIiLUkkY29zRzYkRiVJKF9zeXNsaWJHNiI2IywkSSJ0R0YrIiNIIiIiIyIiJSIkVCkqJkYkRiYtRig2IywkRi4iI0pGMCNGMiIkaCoqJkYkRiYtRig2IywkRi4iIiZGMCNGMiIjRComRiRGJi1GKDYjLCRGLkZBRjAjRjIiJEQnKiZGJEYmLUYoNiMsJEYuIiM2RjAjRjIiJEAiKiZGJEYmLUYoNiMsJEYuIiM4RjAjRjIiJHAiKiZGJEYmLUYoNiMsJEYuIiIqRjAjRjIiIyIpKiZGJEYmLUYoNiMsJEYuIiNYRjAjRjIiJUQ/KiZGJEYmLUYoNiMsJEYuIiIoRjAjRjIiI1wqJkYkRiYtRig2IywkRi4iI1JGMCNGMiIlQDoqJkYkRiYtRig2IywkRi4iI1RGMCNGMiIlIm8iKiZGJEYmLUYoNiMsJEYuIiNWRjAjRjIiJVw9KiZGJEYmLUYoNiMsJEYuIiNORjAjRjIiJUQ3KiZGJEYmLUYoNiNGLkYwRjIqJkYkRiYtRig2IywkRi4iIiRGMCNGMkZaKiZGJEYmLUYoNiMsJEYuIiM6RjAjRjIiJEQjKiZGJEYmLUYoNiMsJEYuIiM8RjAjRjIiJCpHKiZGJEYmLUYoNiMsJEYuIiM+RjAjRjIiJGgkKiZGJEYmLUYoNiMsJEYuIiNARjAjRjIiJFQlKiZGJEYmLUYoNiMsJEYuIiNCRjAjRjIiJEgmKiZGJEYmLUYoNiMsJEYuIiNGRjAjRjIiJEgoKiZGJEYmLUYoNiMsJEYuIiNMRjAjRjIiJSozIiomRiRGJi1GKDYjLCRGLiIjUEYwI0YyIiVwOComRiRGJi1GKDYjLCRGLiIjWkYwI0YyIiU0QSomRiRGJi1GKDYjLCRGLkZkb0YwI0YyIiUsQw==
Here the pattern repeats skips every other even number in the cos summation (it repeats every 4)
f3e := (t,N) -> (4/Pi)*sum(cos((2*k+1)*t)/((2*k+1)^2),k=0..N);
Zio2JEkidEc2IkkiTkdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCQqJkkjUGlHJSpwcm90ZWN0ZWRHISIiLUkkc3VtRzYkRi1JKF9zeXNsaWJHRiU2JComLUkkY29zR0YxNiMqJiwmSSJrR0YlIiIjIiIiRjxGPDkkRjxGPEY5ISIjL0Y6OyIiITklRjwiIiVGJUYlRiU=
Checking
simplify(f3e(t,24)-cossum(f3,50,t));
IiIh
The ODD extension (fig 9.2.4)
plot(sinsum(f2,50,t),t=-Pi..Pi); sinsum(f2,50,t);
NictJSdDVVJWRVNHNiQ3VzckJCEtSlpFZlRKISM2JCExTl0xYCI9KGZqISNDNyQkITI9KnByKTNQWSskISM7JCEyL05wS0hUNFAiISM8NyQkITEzRVsqeXNhKUchIzokITJqXDpJQFw9YyMhIzw3JCQhMmNhKSk+Mmc5diMhIzskITFFR3YiZVNJIVIhIzs3JCQhMSxAQEohZmxoIyEjOiQhME0+Km8hM0lEJiEjOjckJCEyPFFNOjoqSCNbIyEjOyQhMUEtJjRsXGhmJyEjOzckJCEwRDdtMEJ5TiMhIzkkITB4LVBNeTcleSEjOjckJCExJz1kLUZOKkdBISM6JCExRFtkKm8mb0kiKiEjOzckJCEyWipHQlAkUmM0IyEjOyQhMjxUP0ZLKlFZNSEjOzckJCEyWktjZikzeGk+ISM7JCEyaytpR0tCI3o2ISM7NyQkITJwcjRBTio0RT0hIzskITIkNGdVJzNXYkoiISM7NyQkITIlUVEqMyoqPWRxIiEjOyQhMkV6cExxQ2lWIiEjOzckJCEyWV55TTVmemoiISM7JCEyO14mZjgmXFZdIiEjOzckJCEyMT5qZ0AqPnE6ISM7JCEyZFBoTyYqZiFlOiEjOzckJCEyTyEqUis1aEBdIiEjOyQhMSlIS1BWQUVdIiEjOjckJCEybGg7UylIN005ISM7JCEyRCUqUihcTGlNOSEjOzckJCEyeTRGLEspKUhJIiEjOyQhMmsmR3V1T2UuOCEjOzckJCEyJSo0IVIjWzBSPSIhIzskITI+NENNeGtUPSIhIzs3JCQhMm45I1FxV0lVNSEjOyQhMltsKj49bnhVNSEjOzckJCExcDN3JFIpXEIjKiEjOyQhMUtvWjtXI3BBKiEjOzckJCExSmM4bzM5R3khIzskITFuJCl6LSczPCR5ISM7NyQkITEjW0JQLThJZichIzskITE9cjlpd0EnZichIzs3JCQhMUIzPEA/KXlCJiEjOyQhMTRhJCpIMlpTXyEjOzckJCEyJiopUjxZb1paUiEjPCQhMSZcVDE7VSdcUiEjOzckJCExVHFYPVcyLEUhIzskITFKKDRHPj1BZyMhIzs3JCQhMjw5aiUpb2NZTyIhIzwkITJCT0VyJSkpKWZPIiEjPDckJCExLHo3VktKLEohIz0kITFWZjhwbGlTSiEjPTckJCIyLHUjZjN4RWE4ISM8JCIyJGVAKillUVtiOCEjPDckJCIxdXEpSHZlLGMjISM7JCIxOUkqUS8nemdEISM7NyQkIjJqbChSPFRiaVEhIzwkIjIubFs/Sj9RJ1EhIzw3JCQiMSIqPUpvZjAzXyEjOyQiMSNIXEglZlM1XyEjOzckJCIxXjFYSXFPQ2whIzskIjEjNFxxdShRRmwhIzs3JCQiMSRvQmxtbHp6KCEjOyQiMTFgelVUVix5ISM7NyQkIjF2J2YjNCl6P0AqISM7JCIxdHNzTFFqOiMqISM7NyQkIjIiSHY8RUNGWzUhIzskIjJDRzBDWXgnWzUhIzs3JCQiMmFPWUclMyVSPSIhIzskIjFOXkdoISo+JT0iISM6NyQkIjJEISlSRWZ3b0kiISM7JCIyeTJRMCNSSzI4ISM7NyQkIjJZdF0oUXlGVDkhIzskIjIxM00rKD44VDkhIzs3JCQiMjciKnl6UTNYXSIhIzskIjIwIzRBOnlQMDohIzs3JCQiMnkzMnMkKlF4YyIhIzskIjI6eHJCJnAiemIiISM7NyQkIjFXXj8hM0xRaiIhIzokIjExJG82JXo/NDohIzo3JCQiL0s/QnMjKipwIiEjOCQiMjslZkRSQ1tUOSEjOzckJCIyQUgkZU1hO0g9ISM7JCIyVSh6cWxyaTc4ISM7NyQkIjEjPkNaJ2VZaz4hIzokIjIwbTxQbXB2PCIhIzs3JCQiMm1XKHlmaXglNCMhIzskIjJVPy8vZFRzLyIhIzs3JCQiMiU0cSkpZnUuR0EhIzskIjEyandjKlInUiIqISM7NyQkIjFAdycqKj0mPmdCISM6JCIxbTdgRWRdPHkhIzs3JCQiMCMzYT1XaiJbIyEjOSQiMGRyaSFmei1tISM6NyQkIjItaSYqKXl1IjNpIyEjOyQiMSpwPzssKjQ1XyEjOzckJCIxaS1IbldJWEYhIzokIjFyVTFiTDxsUiEjOzckJCIyKT1SV2hOLnlHISM7JCIyTCJcdFtrN1BFISM8NyQkIjJGeCpSU0EyMEkhIzskIjJCTyZSMkJhbTghIzw3JCQiLVhdRWZUSiEjNiQiMW9GQ1xXdnpKISNDLSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJSVWSUVXRzYkOyQhK2FFZlRKISIqJCIrYUVmVEohIiolKERFRkFVTFRHLSUqQVhFU1NUWUxFRzYjJSdOT1JNQUxHLSUoU0NBTElOR0c2IyUuVU5DT05TVFJBSU5FREctJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiNxISIiLSUpQk9VTkRTX1lHNiMkIiMhKiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJSEpUSEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiU1USEiIi0lKUNISUxEUkVORzYi
LFQqKCMiIiUiJEQjIiIiSSNQaUclKnByb3RlY3RlZEchIiItSSRzaW5HNiRGKUkoX3N5c2xpYkc2IjYjLCQqJiIjOkYnSSJ0R0YvRidGJ0YnRioqKCNGJSIlNEFGJ0YoRiotRiw2IywkKiYiI1pGJ0Y0RidGJ0YnRioqKCNGJSIlRD9GJ0YoRiotRiw2IywkKiYiI1hGJ0Y0RidGJ0YnRicqKCNGJSIlQDpGJ0YoRiotRiw2IywkKiYiI1JGJ0Y0RidGJ0YnRioqKCNGJSIlIm8iRidGKEYqLUYsNiMsJComIiNURidGNEYnRidGJ0YnKigjRiUiJVw9RidGKEYqLUYsNiMsJComIiNWRidGNEYnRidGJ0YqKigjRiUiIyIpRidGKEYqLUYsNiMsJComIiIqRidGNEYnRidGJ0YnKigjRiUiJSozIkYnRihGKi1GLDYjLCQqJiIjTEYnRjRGJ0YnRidGJyooI0YlIiUsQ0YnRihGKi1GLDYjLCQqJiIjXEYnRjRGJ0YnRidGJyooI0YlIiRUJUYnRihGKi1GLDYjLCQqJiIjQEYnRjRGJ0YnRidGJyooI0YlIiRAIkYnRihGKi1GLDYjLCQqJiIjNkYnRjRGJ0YnRidGKiooI0YlIiRIJkYnRihGKi1GLDYjLCQqJiIjQkYnRjRGJ0YnRidGKiooI0YlIiQqR0YnRihGKi1GLDYjLCQqJiIjPEYnRjRGJ0YnRidGJyooI0YlRl5vRidGKEYqLUYsNiMsJComIiIkRidGNEYnRidGJ0YqKigjRiUiJGgqRidGKEYqLUYsNiMsJComIiNKRidGNEYnRidGJ0YqKigjRiUiJUQ3RidGKEYqLUYsNiMsJComIiNORidGNEYnRidGJ0YqKigjRiUiJXA4RidGKEYqLUYsNiMsJComIiNQRidGNEYnRidGJ0YnKigjRiUiJEQnRidGKEYqLUYsNiMsJComIiNERidGNEYnRidGJ0YnKigjRiUiJEgoRidGKEYqLUYsNiMsJComIiNGRidGNEYnRidGJ0YqKigjRiUiJFQpRidGKEYqLUYsNiMsJComIiNIRidGNEYnRidGJ0YnKihGJUYnRihGKi1GLDYjRjRGJ0YnKigjRiVGZXRGJ0YoRiotRiw2IywkKiYiIiZGJ0Y0RidGJ0YnRicqKCNGJUZecEYnRihGKi1GLDYjLCQqJiIiKEYnRjRGJ0YnRidGKiooI0YlIiRwIkYnRihGKi1GLDYjLCQqJiIjOEYnRjRGJ0YnRidGJyooI0YlIiRoJEYnRihGKi1GLDYjLCQqJiIjPkYnRjRGJ0YnRidGKg==
Here the pattern is alternating signs for all odd numbers in teh sines summation.
f2o := (t,N) -> (4/Pi)*sum((-1)^k*sin((2*k+1)*t)/(2*k+1)^2,k=0..N);
Zio2JEkidEc2IkkiTkdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCQqKCIiJSIiIkkjUGlHJSpwcm90ZWN0ZWRHISIiLUkkc3VtRzYkRi9JKF9zeXNsaWJHRiU2JCooKUYwSSJrR0YlRi0tSSRzaW5HRjM2IyomLCYqJiIiI0YtRjhGLUYtRi1GLUYtRiRGLUYtKUY9Rj9GMC9GODsiIiFGJkYtRi1GJUYlRiU=
simplify(f2o(t,24)-sinsum(f2,50,t));
IiIh
Function from fig 9.2.5
f3 := t-> unit(t,0,Pi/3)*t + unit(t,Pi/3,2*Pi/3)*(Pi/3) + (Pi-t)*unit(t,2*Pi/3,Pi);
Zio2I0kidEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCgqJi1JJXVuaXRHRiU2JUYkIiIhLCQqJiMiIiIiIiRGMkkjUGlHJSpwcm90ZWN0ZWRHRjJGMkYyRiRGMkYyKihGMUYyLUYsNiVGJEYvLCQqJiMiIiNGM0YyRjRGMkYyRjJGNEYyRjIqJiwmRjRGMkYkISIiRjItRiw2JUYkRjlGNEYyRjJGJUYlRiU=
plot(f3(t),t=0..Pi);
NictJSdDVVJWRVNHNiQ3W283JCQiIiEhIiIlKnVuZGVmaW5lZEc3JCQiMlk8JjMvMSQqUkAhIz4kIjJZPCYzLzEkKlJAISM+NyQkIjFcLjwzNycpelUhIz0kIjFcLjwzNycpelUhIz03JCQiMUNiRDc9ej5rISM9JCIxQ2JENz16PmshIz03JCQiMSpwU2pUQShmJikhIz0kIjEqcFNqVEEoZiYpISM9NyQkIjJcNV5DT2VSRyIhIz0kIjJcNV5DT2VSRyIhIz03JCQiMihSIm9LW1c+ciIhIz0kIjIoUiJvS1tXPnIiISM9NyQkIjIoNEEhXHM7emMjISM9JCIyKDRBIVxzO3pjIyEjPTckJCIyJXppYG0qKSlRVSQhIz0kIjIlemlgbSopKVFVJCEjPTckJCIyJj5XISlcTSRlOCYhIz0kIjImPlchKVxNJGU4JiEjPTckJCIxZkQyTHp4Wm8hIzwkIjFmRDJMenhabyEjPDckJCIxYjRVOmEpbyMpKiEjPCQiMWI0VTphKW8jKSohIzw3JCQiMU5weCpHKmYhRyIhIzskIjFOcHgqRypmIUciISM7NyQkIjA2I2V4R21dPiEjOiQiMDYjZXhHbV0+ISM6NyQkIjJYVCw9M29eaSMhIzwkIjJYVCw9M29eaSMhIzw3JCQiMFxEMFtua0gkISM6JCIwXEQwW25rSCQhIzo3JCQiMSI9WmEmeiUpPVIhIzskIjEiPVphJnolKT1SISM7NyQkIjJ3Yld2KW9HalghIzwkIjJ3Yld2KW9HalghIzw3JCQiMSUzKipIYm0oSF8hIzskIjElMyoqSGJtKEhfISM7NyQkIjFQUnI0KTNUKmUhIzskIjFQUnI0KTNUKmUhIzs3JCQiMSNbKXl5a1l4bCEjOyQiMSNbKXl5a1l4bCEjOzckJCIxcydvY0dvJHpyISM7JCIxcydvY0dvJHpyISM7NyQkIjFRMDtncidwJnkhIzskIjFRMDtncidwJnkhIzs3JCQiMUdNdD8kW3RgKSEjOyQiMUdNdD8kW3RgKSEjOzckJCIxI3AzMGtAST4qISM7JCIxI3AzMGtAST4qISM7NyQkIjBDIlxIZVYpeSohIzokIjBDIlxIZVYpeSohIzo3JCQiMk4hKW9mJFFDOTUhIzskIjJOISlvZiRRQzk1ISM7NyQkIjJHWykpKikzVydcNSEjOyQiMndmJz5edj5aNSEjOzckJCIyVylcdllGano1ISM7JCIyd2YnPl52Plo1ISM7NyQkIjEnW0BYU0AnNDYhIzokIjJ3Zic+XnY+WjUhIzs3JCQiMlZoUDNHKlF6NiEjOyQiMndmJz5edj5aNSEjOzckJCIyVTMqM3RjOVQ3ISM7JCIyd2YnPl52Plo1ISM7NyQkIjIkeTBEQkEhKjM4ISM7JCIyd2YnPl52Plo1ISM7NyQkIjJQbUs/W0FNUCIhIzskIjJ3Zic+XnY+WjUhIzs3JCQiMWFAWC5FdVM5ISM6JCIyd2YnPl52Plo1ISM7NyQkIjJRdSMpKipbakRdIiEjOyQiMndmJz5edj5aNSEjOzckJCIyN208eW1YI3A6ISM7JCIyd2YnPl52Plo1ISM7NyQkIjJuWi4qNCg0JlE7ISM7JCIyd2YnPl52Plo1ISM7NyQkIjJgS2BARS8pKXAiISM7JCIyd2YnPl52Plo1ISM7NyQkIjJ4UjEvLkNSdyIhIzskIjJ3Zic+XnY+WjUhIzs3JCQiMHQmKUg3Kj5KPSEjOSQiMndmJz5edj5aNSEjOzckJCIyJGZiN3dZLCgqPSEjOyQiMndmJz5edj5aNSEjOzckJCIyWS5oemclcGc+ISM7JCIyd2YnPl52Plo1ISM7NyQkIjFwSjg6LlNKPyEjOiQiMndmJz5edj5aNSEjOzckJCIyWCQ+WUVrO2o/ISM7JCIyd2YnPl52Plo1ISM7NyQkIi8yelBEJFw0IyEjOCQiMko0Kno6LG1ZNSEjOzckJCIxZlsoPmtcKUdAISM6JCIyViRcaDZJdTc1ISM7NyQkIjJ4LGZodW1GOyMhIzskIjFieEl1IWYjKXkqISM7NyQkIjJbWCczQFlCQ0EhIzskIjElUUxdSyFldCIqISM7NyQkIjEvYjxXX1YiSCMhIzokIjAqRzklNHU6XSkhIzo3JCQiMEdOTXpsWU4jISM5JCIxS15hLCdvI3B5ISM7NyQkIjJtT21rJCpmMlUjISM7JCIxbVVCcnJLM3MhIzs3JCQiMTEoKT1VIXpgWyMhIzokIjFzMyw5aDhpbCEjOzckJCIxa0hIZCNISWIjISM6JCIxJEhvSCdSaiYpZSEjOzckJCIxX3ImW1glPT1FISM6JCIxOGtLKCk+M01fISM7NyQkIjFLLSVcMDpbbyMhIzokIjJEaCZcJylmeG5YISM8NyQkIjIwYzgrI1IqM3YjISM7JCIya0tpZEwoKXAhUiEjPDckJCIxLzBMTU5oNkchIzokIjEqKUdmIz4ieipIJCEjOzckJCIyR29VWDEwNylHISM7JCIyTjVyLyplKFFnIyEjPDckJCIyJTQ2eGUmW00lSCEjOyQiMm4kbz1bNFciKT4hIzw3JCQiMkw1IillNTgpNEkhIzskIjIkKSpvM3hhejw4ISM8NyQkIjJMeCFSWENMdEkhIzskIjF6PiEqPjMtRW8hIzw3JCQiLVhdRWZUSiEjNiQiMXZIQjlMelJKISNDLSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJSVWSUVXRzYkOyQiIiEhIiIkIithRWZUSiEiKiUoREVGQVVMVEctJSpBWEVTU1RZTEVHNiMlJ05PUk1BTEctJShTQ0FMSU5HRzYjJS5VTkNPTlNUUkFJTkVERy0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJEkkISIiLSUpQk9VTkRTX1lHNiMkIiMhKiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJT9PISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJVNPISIiLSUpQ0hJTERSRU5HNiI=
plot(sinsum(f3,48,t),t=-Pi..Pi); sinsum(f3,48,t);
NictJSdDVVJWRVNHNiQ3YW83JCQhLUpaRWZUSiEjNiQhMSlcKnBGJSpvJD0nISNDNyQkITI9KnByKTNQWSskISM7JCEyKEhHLzMiKmZvOCEjPDckJCExM0VbKnlzYSlHISM6JCEya1hqLD9XP2MjISM8NyQkITJjYSkpPjJnOXYjISM7JCEybE1HLEZoOyFSISM8NyQkITEsQEBKIWZsaCMhIzokITEoUS45K2MpXF8hIzs3JCQhMjxRTTo6KkgjWyMhIzskITEyJDQoSCFSOWYnISM7NyQkITBEN20wQnlOIyEjOSQhMSVbJTMoPUd1JHkhIzs3JCQhMSc9ZC1GTipHQSEjOiQhMUxZYS4qcGU3KiEjOzckJCEyLi9YUEkoR2lAISM7JCExMyE0eS1kJSl6KiEjOzckJCEyWipHQlAkUmM0IyEjOyQhMnUpcFc2KVEqUjUhIzs3JCQhMlBLQjN5SSF6PyEjOyQhMktNYGQiPU9ZNSEjOzckJCEyQXY4V0FBQzEjISM7JCEyaiZwMm8hcCFcNSEjOzckJCEyMj0vIW9PImUvIyEjOyQhMic+MmQtejpcNSEjOzckJCEyKDRZZjZeP0g/ISM7JCEyVSMpZWV1Jyp6LyIhIzs3JCQhMnNZdigpKnopZio+ISM7JCEyJFwuSkF5PVk1ISM7NyQkITJaS2NmKTN4aT4hIzskITInPSp6Qyk0JnAvIiEjOzckJCEycHI0QU4qNEU9ISM7JCEyLUtIYCNSQFo1ISM7NyQkITIlUVEqMyoqPWRxIiEjOyQhMmoyWCQqekNyLyIhIzs3JCQhMjE+amdAKj5xOiEjOyQhMkJIOjBWZXIvIiEjOzckJCEybGg7UylIN005ISM7JCEyKlIpPTIpUjZaNSEjOzckJCEyeTRGLEspKUhJIiEjOyQhMWsleidSdSxaNSEjOjckJCEyJSo0IVIjWzBSPSIhIzskITBrZUU9YHIvIiEjOTckJCEyODYpUUhfXVs2ISM7JCEyVkV5JSllbWgvIiEjOzckJCExQmhRd1w1ODYhIzokITI4Y3VAQFF6LyIhIzs3JCQhMilHXikpXFtTJjQiISM7JCEyPE82Qnl1Ilw1ISM7NyQkITJbOCVRQlpxeDUhIzskITIwbCVmYVomKls1ISM7NyQkITI0OSQpb2YvKzEiISM7JCEyR3JnTyhbalg1ISM7NyQkITJuOSNRcVdJVTUhIzskITF0UGlyNih6LiIhIzo3JCQhMW82els6RkIpKiEjOyQhMXB3Jj5xKWZKKSohIzs3JCQhMXAzdyRSKVxCIyohIzskITEnNDQqeWtWPCMqISM7NyQkITFKYzhvMzlHeSEjOyQhMSZIOjIreiJHeSEjOzckJCExI1tCUC04SWYnISM7JCExbTxqO19eImYnISM7NyQkITFCMzxAPyl5QiYhIzskITFVJmZNIilbd0ImISM7NyQkITImKilSPFlvWlpSISM8JCEyTWpLeG4ucSVSISM8NyQkITFUcVg9VzIsRSEjOyQhMigqND5cJFJHLEUhIzw3JCQhMjw5aiUpb2NZTyIhIzwkITFcJzNpRHRQTyIhIzs3JCQhMSx6N1ZLSixKISM9JCExXlRePnEtYUkhIz03JCQiMix1I2YzeEVhOCEjPCQiMmleWydHM2FgOCEjPDckJCIxdXEpSHZlLGMjISM7JCIxSjpXYDAtaEQhIzs3JCQiMmpsKFI8VGJpUSEjPCQiMmNSbChcTmJqUSEjPDckJCIxIio9Sm9mMDNfISM7JCIxYG06L0tUM18hIzs3JCQiMV4xWElxT0NsISM7JCIxOzInPk4rWF8nISM7NyQkIjEkb0JsbWx6eighIzskIjE2LCYzVmIqKXooISM7NyQkIjF2J2YjNCl6P0AqISM7JCIxZiNIJEdHajEjKiEjOzckJCIxJFs8Yi4tdSUpKiEjOyQiMWNEakZ0MGUpKiEjOzckJCIyIkh2PEVDRls1ISM7JCIybUY3LWokM1Q1ISM7NyQkIjJzJT5eTW01OzYhIzskIjJOKVFiTlpwWjUhIzs3JCQiMmFPWUclMyVSPSIhIzskIjEzW0phWDpaNSEjOjckJCIyRCEpUkVmd29JIiEjOyQiMkApXCJlO3lxLyIhIzs3JCQiMll0XShReUZUOSEjOyQiMnd1SjAtZnIvIiEjOzckJCIyeTMycyQqUXhjIiEjOyQiMkEnUUdKJWVyLyIhIzs3JCQiL0s/QnMjKipwIiEjOCQiMk9pIioqejk7WjUhIzs3JCQiMkFIJGVNYTtIPSEjOyQiMnJpbDtVbXIvIiEjOzckJCIxIz5DWidlWWs+ISM6JCIyOWU0a2IjKW8vIiEjOzckJCIyYyshXGpNLygqPiEjOyQiMmQkSCdRVyxpLyIhIzs3JCQiMiQ+ZURpNWlIPyEjOyQiMikpRyUzdTEuWzUhIzs3JCQiMUUoUTsnKTRmLyMhIzokIjIlNHkiSElpIlw1ISM7NyQkIjFMOy1oJyk+aT8hIzokIjFYWSx0WTNcNSEjOjckJCIwYS8vWShbeT8hIzkkIjIxdyVHUFldWTUhIzs3JCQiMm1XKHlmaXglNCMhIzskIjEjcChHbldQUzUhIzo3JCQiMUdzJClmb1NoQCEjOiQiMShlOmJtNSIzKSohIzs3JCQiMiU0cSkpZnUuR0EhIzskIjFaN2o9QE5NIiohIzs3JCQiMUB3JyoqPSY+Z0IhIzokIjB1TisiSFk5eSEjOjckJCIwIzNhPVdqIlsjISM5JCIxVjxVKD1LemYnISM7NyQkIjItaSYqKXl1IjNpIyEjOyQiMSUqelBuXDYzXyEjOzckJCIxaS1IbldJWEYhIzokIjE9eTZZMjlpUiEjOzckJCIyKT1SV2hOLnlHISM7JCIxPCsoKlI0Q05FISM7NyQkIjJGeCpSU0EyMEkhIzskIjFrMXkwREprOCEjOzckJCItWF1FZlRKISM2JCIyKVtpViM9Vjw0JCEjRC0lJkNPTE9SRzYmJSRSR0JHJCIjNSEiIiQiIiEhIiIkIiIhISIiLSUlVklFV0c2JDskISthRWZUSiEiKiQiK2FFZlRKISIqJShERUZBVUxURy0lKkFYRVNTVFlMRUc2IyUnTk9STUFMRy0lKFNDQUxJTkdHNiMlLlVOQ09OU1RSQUlORURHLSUlUk9PVEc2Jy0lKUJPVU5EU19YRzYjJCIjcSEiIi0lKUJPVU5EU19ZRzYjJCIjISohIiItJS1CT1VORFNfV0lEVEhHNiMkIiUhKVEhIiItJS5CT1VORFNfSEVJR0hURzYjJCIlNVEhIiItJSlDSElMRFJFTkc2Ig==
LEIqKiIiIyIiIkkjUGlHJSpwcm90ZWN0ZWRHISIiKSIiJCNGJUYkRiUtSSRzaW5HNiRGJ0koX3N5c2xpYkc2IjYjSSJ0R0YwRiVGJSoqI0YkIiNERiVGJkYoRilGJS1GLTYjLCQqJiIiJkYlRjJGJUYlRiVGKCoqI0YkIiNcRiVGJkYoRilGJS1GLTYjLCQqJiIiKEYlRjJGJUYlRiVGJSoqI0YkIiRAIkYlRiZGKEYpRiUtRi02IywkKiYiIzZGJUYyRiVGJUYlRigqKiNGJCIkcCJGJUYmRihGKUYlLUYtNiMsJComIiM4RiVGMkYlRiVGJUYlKiojRiQiJCpHRiVGJkYoRilGJS1GLTYjLCQqJiIjPEYlRjJGJUYlRiVGKCoqI0YkIiRoJEYlRiZGKEYpRiUtRi02IywkKiYiIz5GJUYyRiVGJUYlRiUqKiNGJCIkSCZGJUYmRihGKUYlLUYtNiMsJComIiNCRiVGMkYlRiVGJUYoKiojRiQiJEQnRiVGJkYoRilGJS1GLTYjLCQqJkY1RiVGMkYlRiVGJUYlKiojRiQiJFQpRiVGJkYoRilGJS1GLTYjLCQqJiIjSEYlRjJGJUYlRiVGKCoqI0YkIiRoKkYlRiZGKEYpRiUtRi02IywkKiYiI0pGJUYyRiVGJUYlRiUqKiNGJCIlRDdGJUYmRihGKUYlLUYtNiMsJComIiNORiVGMkYlRiVGJUYoKiojRiQiJXA4RiVGJkYoRilGJS1GLTYjLCQqJiIjUEYlRjJGJUYlRiVGJSoqI0YkIiUibyJGJUYmRihGKUYlLUYtNiMsJComIiNURiVGMkYlRiVGJUYoKiojRiQiJVw9RiVGJkYoRilGJS1GLTYjLCQqJiIjVkYlRjJGJUYlRiVGJSoqI0YkIiU0QUYlRiZGKEYpRiUtRi02IywkKiYiI1pGJUYyRiVGJUYlRig=
Here the pattern repeats every 6 terms
f3o := (t,N) -> (2*sqrt(3)/Pi)*sum(sin((6*k+1)*t)/(6*k+1)^2 - sin((6*k+5)*t)/(6*k+5)^2,k=0..N);
Zio2JEkidEc2IkkiTkdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCQqKiIiIyIiIi1JJXNxcnRHNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHRiU2IyIiJEYtSSNQaUdGMSEiIi1JJHN1bUdGMDYkLCYqJi1JJHNpbkdGMDYjKiYsJiomIiInRi1JImtHRiVGLUYtRi1GLUYtRiRGLUYtKUZARixGNkYtKiYtRj02IyomLCZGQUYtIiImRi1GLUYkRi1GLSlGSUYsRjZGNi9GQzsiIiFGJkYtRi1GJUYlRiU=
Checking:
simplify(f3o(t,7)-sinsum(f3,48,t));
IiIh