Gram-Schmidt orthogonalization and approximation We want to find the polynomial of degree 4 that is the best approximation to NiMtSSRjb3NHNiI2I0kieEdGJQ== on the interval NiM3JCwkSSNQaUdJKnByb3RlY3RlZEdGJiEiIkYl . To do this we first take the standard basis for polynomials, { NiM2JyIiIkkieEc2IiokRiUiIiMqJEYlIiIkKiRGJSIiJQ== } and use Gram-Schmidt to construct an orthonormal basis. While we could do this more efficiently using built-in commands of maple, we will just use maple as a calculator.
First, we define our inner product as a function on two functions. Our inner product is the usual one, that is NiMvLUkkPCw+RzYiNiRJInBHRiY2I0kicUdGJi1JJGludEdGJjYkKiYtRig2I0kieEdGJiIiIi1GKkYwRjIvRjE7LCRJI1BpR0kqcHJvdGVjdGVkR0Y4ISIiRjc= . Note that ||p-q|| is small precisely when the values of NiMtSSJwRzYiNiNJInhHRiU= are close to those of NiMtSSJxRzYiNiNJInhHRiU= for all x in the interval NiM3JCwkSSNQaUdJKnByb3RlY3RlZEdGJiEiIkYl . For simplicity, we'll operate on expressions in NiNJInhHNiI= rather than general functions. InnerProd:=(p,q) -> int(p*q, x=-Pi..Pi); NiM+SSpJbm5lclByb2RHNiJmKjYkSSJwR0YlSSJxR0YlRiU2JEkpb3BlcmF0b3JHRiVJJmFycm93R0YlRiUtSSRpbnRHNiRJKnByb3RlY3RlZEdGMEkoX3N5c2xpYkdGJTYkKiY5JCIiIjklRjUvSSJ4R0YlOywkSSNQaUdGMCEiIkY7RiVGJUYl We'll also want to compute the norm of our functions. However, Norm is a maple built-in function, so we'll be a little more formal and call Norm by his full name. Norman:=p->sqrt(InnerProd(p,p)); NiM+SSdOb3JtYW5HNiJmKjYjSSJwR0YlRiU2JEkpb3BlcmF0b3JHRiVJJmFycm93R0YlRiUtSSVzcXJ0RzYkSSpwcm90ZWN0ZWRHRi9JKF9zeXNsaWJHRiU2Iy1JKklubmVyUHJvZEdGJTYkOSRGNUYlRiVGJQ== Now we do Gram-Schmidt to construct an orthonormal basis (with respect to this inner product) for polynomials of degree 4. The first element of this basis is the constant function NiMvLSZJImZHNiI2IyIiITYjSSJ4R0YnIiIi . However, we need to normalize it to be of unit length in our norm. f[0] := 1/Norman(1);
NiM+JkkiZkc2IjYjIiIhLCQqJiIiIyMiIiJGK0kjUGlHSSpwcm90ZWN0ZWRHRi8jISIiRitGLA== Next, we confirm that NiMvLSZJImZHNiI2IyIiIjYjSSJ4R0YnRis= is orthogonal to NiMmSSJmRzYiNiMiIiE= . Since we want an orthonormal basis, we still want to divide NiMmSSJmRzYiNiMiIiI= by its length InnerProd(x,f[0]); NiMiIiE= f[1]:= x/Norman(x); NiM+JkkiZkc2IjYjIiIiLCQqKEkieEdGJkYoIiInI0YoIiIjSSNQaUdJKnByb3RlY3RlZEdGMCMhIiRGLkYt Notice that while NiMtJkkiZkc2IjYjIiIiNiNJInhHRiY= is orthogonal to NiMqJEkieEc2IiIiIw== , NiMmSSJmRzYiNiMiIiE= is not. InnerProd(x^2,f[0]); NiMsJComIiIjIyIiIkYlSSNQaUdJKnByb3RlY3RlZEdGKSMiIiZGJSNGJyIiJA== InnerProd(x^2,f[1]); NiMiIiE= This means we have to use Gram-Schmidt to find a function in the span of NiM8JSZJImZHNiI2IyIiISZGJTYjIiIiKiRJInhHRiYiIiM= which is orthogonal to both NiMmSSJmRzYiNiMiIiE= and NiMmSSJmRzYiNiMiIiI= . Note that since NiMmSSJmRzYiNiMiIiE= and NiMmSSJmRzYiNiMiIiI= are already normalized, we don't have to divide by their lengths. ftemp:= x^2 - InnerProd(x^2,f[0])*f[0] - InnerProd(x^2,f[1])*f[1];
f[2] := collect(ftemp/Norman(ftemp),x); NiM+SSZmdGVtcEc2IiwmKiRJInhHRiUiIiMiIiIqJEkjUGlHSSpwcm90ZWN0ZWRHRi1GKSMhIiIiIiQ= NiM+JkkiZkc2IjYjIiIjLCYqKCIjNSMiIiJGKEkjUGlHSSpwcm90ZWN0ZWRHRi8jISImRihJInhHRiZGKCMiIiQiIiUqJkYuIyEiIkYoRitGLCNGOEY1 Now we play the same game with NiMqJEkieEc2IiIiJA== and NiMqJEkieEc2IiIiJQ== . Note that if you were working these integrals out by hand, you'd quickly realize that the even powers are orthogonal to the odd ones. ftemp:= x^3 - InnerProd(x^3,f[0])*f[0] - InnerProd(x^3,f[1])*f[1] - InnerProd(x^3,f[2])*f[2]:
f[3] := collect(ftemp/Norman(ftemp),x); NiM+JkkiZkc2IjYjIiIkLCYqKCIjOSMiIiIiIiNJI1BpR0kqcHJvdGVjdGVkR0YwIyEiKEYuSSJ4R0YmRigjIiImIiIlKihGLyMhIiRGLkYrRixGM0YtI0Y5RjY= ftemp:= x^4 - InnerProd(x^4,f[0])*f[0] - InnerProd(x^4,f[1])*f[1] - InnerProd(x^4,f[2])*f[2] - InnerProd(x^4,f[3])*f[3]:
f[4] := collect(ftemp/Norman(ftemp),x); NiM+JkkiZkc2IjYjIiIlLCgqKCIiIyMiIiJGK0kjUGlHSSpwcm90ZWN0ZWRHRi8jISIqRitJInhHRiZGKCMiJDAiIiM7KihGLiMhIiZGK0YrRixGMkYrIyEjWCIiKSomRitGLEYuIyEiIkYrIyIiKkY1 So, we now have an orthonormal basis for fourth degree polynomials on NiM3JCwkSSNQaUdJKnByb3RlY3RlZEdGJiEiIkYl . Just in case you are disturbed by all the NiNJI1BpR0kqcHJvdGVjdGVkR0Yk s and square roots, let's also write their approximate values, and also draw their graphs. OrthoBasis:=[f[0],f[1],f[2],f[3],f[4]]; NiM+SStPcnRob0Jhc2lzRzYiNycsJComIiIjIyIiIkYpSSNQaUdJKnByb3RlY3RlZEdGLSMhIiJGKUYqLCQqKEkieEdGJUYrIiInRipGLCMhIiRGKUYqLCYqKCIjNUYqRiwjISImRilGMkYpIyIiJCIiJSomRixGLkY4RiojRi9GPSwmKigiIzlGKkYsIyEiKEYpRjJGPCMiIiZGPSooRixGNEZCRipGMkYrI0Y1Rj0sKCooRilGKkYsIyEiKkYpRjJGPSMiJDAiIiM7KihGLEY5RilGKkYyRikjISNYIiIpRigjIiIqRk8= evalf(OrthoBasis,4); NiM3JyQiJSkpUiEiJSwkSSJ4RzYiJCIlKD4jRiYsJiokRigiIiMkIiVhOEYmJCElZVdGJiIiIiwmKiRGKCIiJCQiJSVcKSEiJkYoJCElTl1GJiwoKiRGKCIiJSQiJWZgRjlGLSQhJVNYRiYkIiUnWyVGJkYz plot(OrthoBasis,x=-Pi..Pi, color=[black,red,blue,green,brown],thickness=2); LSUlUExPVEc2Ki0lJ0NVUlZFU0c2JDdTNyQkITMqKioqKjR0ayNmVEohIzwkIjMuRlYsL0dVKilSISM9NyQkITN3IipwcikzUFkrJEYsRi03JCQhM0AzRVsqeXNhKUdGLEYtNyQkITNmWCYpKT4yZzl2I0YsRi03JCQhMyk0NTc3LmZsaCNGLEYtNyQkITN2IlFNOjoqSCNbI0YsRi03JCQhM3lcQWhjSSN5TiNGLEYtNyQkITMmZT1kLUZOKkdBRixGLTckJCEzdyUqR0JQJFJjNCNGLEYtNyQkITNtQ2omZikzeGk+RixGLTckJCEzJ29yNEFOKjRFPUYsRi03JCQhM1FRUSozKio9ZHEiRixGLTckJCEzZCE+amdAKj5xOkYsRi03JCQhM2A7bSwlKUg3TTlGLEYtNyQkITN2KDRGLEspKUhJIkYsRi03JCQhM1MqNCFSI1swUj0iRixGLTckJCEzcFlAUXFXSVU1RixGLTckJCEzV3AzdyRSKVxCIypGL0YtNyQkITNFSmM4bzM5R3lGL0YtNyQkITMnPVtCUC04SWYnRi9GLTckJCEzQ0IzPEA/KXlCJkYvRi03JCQhM2AqKVI8WW9aWlJGL0YtNyQkITM4VHFYPVcyLEVGL0YtNyQkITNxVEpZKW9jWU8iRi9GLTckJCEzOCx6N1ZLSixKISM/Ri03JCQiMzpTRmYzeEVhOEYvRi03JCQiMydSMigpSHZlLGMjRi9GLTckJCIzSWN3UjxUYmlRRi9GLTckJCIzKjQqPUpvZjAzX0YvRi03JCQiM3BdMVhJcU9DbEYvRi03JCQiMyVIb0JsbWx6eihGL0YtNyQkIjNedSdmIzQpej9AKkYvRi03JCQiMzlIdjxFQ0ZbNUYsRi03JCQiM05saiVHJTMlUj0iRixGLTckJCIzXC0pUkVmd29JIkYsRi03JCQiM2NNMnZReUZUOUYsRi03JCQiM3ooMzJzJCpReGMiRixGLTckJCIzNytLP0JzIyoqcCJGLEYtNyQkIjM8I0gkZU1hO0g9RixGLTckJCIzJT4+Q1onZVlrPkYsRi03JCQiM2tZdXlmaXglNCNGLEYtNyQkIjNUNHEpKWZ1LkdBRixGLTckJCIzej93JyoqPSY+Z0JGLEYtNyQkIjMhKj4zYT1XaiJbI0YsRi03JCQiM0E/YyopeXUiM2kjRixGLTckJCIzLGktSG5XSVhGRixGLTckJCIzdj1SV2hOLnlHRixGLTckJCIzc3MoKlJTQTIwSUYsRi03JCQiMyEpKioqXC9sI2ZUSkYsRi0tJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiJGX3VGX3UtRiY2JDdTNyRGKiQhM00wV2h2SCkpNHBGLzckRjEkITNCIW94a1NeJzNtRi83JEY0JCEzJlxFWylHOWJZakYvNyRGNyQhM0hiSSM9ZiN6XmdGLzckRjokITN2JTQiKj0hKnpddiZGLzckRj0kITNwVTJlTnV4ZmFGLzckRkAkITMqSDlrekMlKmY9JkYvNyRGQyQhM1EmR0dgaDJEIVxGLzckRkYkITM3aFdLTWBLNFlGLzckRkkkITNrcVNoI0gkMzxWRi83JEZMJCEzOyUqPWdKa1o7U0YvNyRGTyQhMyRSayl6VzxxXlBGLzckRlIkITNDIj5zImYiR09YJEYvNyRGVSQhM1VPW3dhMUxhSkYvNyRGWCQhMzMwJSp5XT8hZidHRi83JEZlbiQhM0FPY1VZNylSZyNGLzckRmhuJCEzX0xYXnc/YCNII0YvNyRGW28kITNcYVx1d1BwRz9GLzckRl5vJCEzbW8sMzdyeUA8Ri83JEZhbyQhM0Umej0wckEsWCJGLzckRmRvJCEzPWJYXjBRMV82Ri83JEZnbyQhMydITCVmSkNUI28pISM+NyRGam8kITMkKT0nUUd1QDVzJkZneTckRl1wJCEzaipSUVdrUzorJEZneTckRmBwJCEzbGo1JnkxKkdAbyEjQDckRmRwJCIzJW9qSihbLXB5SEZneTckRmdwJCIzJ1JlNCYpUUc1aiZGZ3k3JEZqcCQiM3VRdmtVaWkmXClGZ3k3JEZdcSQiM3Z2dHE3T11YNkYvNyRGYHEkIjNpbGVtTFQtTjlGLzckRmNxJCIzUkk1Y0MsOjo8Ri83JEZmcSQiMyY0NndzRiM9RT9GLzckRmlxJCIzJlxYLy1dZWNJI0YvNyRGXHIkIjNBRmU9QCFmU2cjRi83JEZfciQiM14jeUtVQ2FXKEdGLzckRmJyJCIzT0Iyd1Z3MXFKRi83JEZlciQiM2sjZURdIm9AW01GLzckRmhyJCIzTWkhUVMlRycqUVBGLzckRltzJCIzImVZTDRCP0stJUYvNyRGXnMkIjMwUCU0OVA2M0slRi83JEZhcyQiM2tGUUI0cVUyWUYvNyRGZHMkIjNFeklzJipHYCtcRi83JEZncyQiM3lMaVs8PEAiPiZGLzckRmpzJCIzPSIqPiZmTzokZWFGLzckRl10JCIzczhdcypHWVd3JkYvNyRGYHQkIjMzVTNPUE1EUWdGLzckRmN0JCIzZk0nKkhsKik9SWpGLzckRmZ0JCIzViEzcmNeMyc0bUYvNyRGaXQkIjNEJnk/RClIKSk0cEYvLUZcdTYmRl51JCIjNUZhdUZfdUZfdS1GJjYkN1M3JEYqJCIzeEghcHNfPzEjKilGLzckRjEkIjMlZkZhRGwmUXp4Ri83JEY0JCIzLF0/cDQ+eUZvRi83JEY3JCIzciIqKlFVKFFnLmVGLzckRjokIjNrcydIJDQqPj0jW0YvNyRGPSQiMylcJG9MNmtwJCpRRi83JEZAJCIzLTNZTDI8KG8yJEYvNyRGQyQiMypmVndrWHFgRiNGLzckRkYkIjNrdS5TZiNRUVwiRi83JEZJJCIzKzMjKW8rNmZGd0ZneTckRkwkIjNuaCw5amspKXBnRmJwNyRGTyQhMzpfO0MwPUJkXkZneTckRlIkITMhSFRZRE4/dzYiRi83JEZVJCEzIVwzI1EkbyMpPW4iRi83JEZYJCEzJGYjKnkyMjUmZUBGLzckRmVuJCEzXXgkcGUmcCxnREYvNyRGaG4kITNhVlZqVTxTKClIRi83JEZbbyQhM1AoKWU2TmMicEkkRi83JEZebyQhM1pZeC05aVxIT0YvNyRGYW8kITNxSyEpeWFlKTQoUUYvNyRGZG8kITMnKipweHJ4WyQpMyVGLzckRmdvJCEzaXosKDQ7WSFcVUYvNyRGam8kITNSQz1eZFZlb1ZGLzckRl1wJCEzISpSME88PjFOV0YvNyRGYHAkITNNJD51LkQoSGdXRi83JEZkcCQhM1Q9Yj0hKVtXTldGLzckRmdwJCEzIUdLVkNUWjlQJUYvNyRGanAkITNWPm1bMCRRIWVVRi83JEZdcSQhMyVld1xYJEdkIzQlRi83JEZgcSQhM3R0KlJzKFI+JClRRi83JEZjcSQhMzVkIzQ5JyopKWVqJEYvNyRGZnEkITN3c1AkZmtwKDRMRi83JEZpcSQhM3MkNEtqXCdbcUhGLzckRlxyJCEzIUhhUydSTSEqZkRGLzckRl9yJCEzKlx3KW8hKT12V0BGLzckRmJyJCEzJSkpW1giKkchKlJrIkYvNyRGZXIkITN6UjRjJj0oM0c2Ri83JEZociQhMylcLGteKFFsQ2FGZ3k3JEZbcyQiM2BjNiM+T1kkKmUoRmJwNyRGXnMkIjMpSCFSYnIiUnlyKEZneTckRmFzJCIzaGlEOy5cJCopWyJGLzckRmRzJCIzTlRObDhgJSpwQUYvNyRGZ3MkIjM/QUgpMz5YPzQkRi83JEZqcyQiM0x3YkUleUEjKilRRi83JEZddCQiMyVIJ1EtOXkwX1tGLzckRmB0JCIzMGN0cC4sdGRkRi83JEZjdCQiM3czLVEpeV4ncG5GLzckRmZ0JCIzWDxNKil5NiRIeShGLzckRml0JCIzWUN0LGEwaT8qKUYvLUZcdTYmRl51Rl91Rl91RlxfbC1GJjYkN1g3JEYqJCEzOi0/dls/XWI1Riw3JCQhMzQnXDghb1s2dElGLCQhMzcmMyM0Ij5VPkAqRi83JEYxJCEzc3E4OkVyV1V6Ri83JCQhMykqKnoqNFJcMFhIRiwkITNDZi8+aWVYJypvRi83JEY0JCEzTSRRY19bXVEhZkYvNyRGNyQhM1dkeWgxJWYyJ1FGLzckRjokITNmXTNyVVskKmU/Ri83JEY9JCEzKUd4XyIqR2g1MiZGZ3k3JEZAJCIzJVxMTylSK2Vzc0ZneTckRkMkIjNuPy9DKmU4KjM9Ri83JEZGJCIzVXNSLSFmOSlHRkYvNyRGSSQiM1UwcXJhZ2BjTUYvNyRGTCQiMzhlQVVaP2c/U0YvNyRGTyQiMyk0dWJwOE1GUCVGLzckRlIkIjNTd1FUJ1JzJj1ZRi83JEZVJCIzYDVeN012RjxaRi83JEZYJCIzMCcpXDhILSZSbyVGLzckRmVuJCIzNykpKnBNUmNVYiVGLzckRmhuJCIzL2AxX0cqcCIqRyVGLzckRltvJCIzMipwMExqTzApUkYvNyRGXm8kIjNSaywocGxnb2AkRi83JEZhbyQiM1R1bGZgV3Z5SUYvNyRGZG8kIjNtJmUmSHBEVDxERi83JEZnbyQiM2UlNFlmaVRxJD5GLzckRmpvJCIzNT5lUVlWKGVIIkYvNyRGXXAkIjM/Q28pZlZsZCZvRmd5NyRGYHAkIjNdIUd3ZjheSGMiRmJwNyRGZHAkITNEZEdkNCkpKlEhb0ZneTckRmdwJCEzTjtKX3AqXGZGIkYvNyRGanAkITMpXDtaIVsqXHYqPUYvNyRGXXEkITNfQE1OJmZlV10jRi83JEZgcSQhMyRRaFh4VylwXklGLzckRmNxJCEzPG1Sa3lpTkVORi83JEZmcSQhMyhwLm5hN2ZzKFJGLzckRmlxJCEzdWdreFkhKWYtVkYvNyRGXHIkITM/ZzhjYiEzVmIlRi83JEZfciQhM1hjakYzN2snbyVGLzckRmJyJCEzZW5COyRbeGJyJUYvNyRGZXIkITM5Kkg7Iyo0TztpJUYvNyRGaHIkITMrRm9gW0JVJ1ElRi83JEZbcyQhM0RHLystayEqNFNGLzckRl5zJCEzeSNcJ0h0QVJbTUYvNyRGYXMkITMoKXBYJz5HUVR0I0YvNyRGZHMkITMvSjpCREN4Oj1GLzckRmdzJCEzRG5TO2ZvNGJxRmd5NyRGanMkIjM9YHhFO1wsK11GZ3k3JEZddCQiM1Mlb1UmeTQuN0BGLzckRmB0JCIzeUd3PWEicEl4JEYvNyQkIjNRITRuVixwOyJHRiwkIjNRSHdmLFpzWVpGLzckRmN0JCIzZzUtcSlHN095JkYvNyQkIjN1WD0jNCFIYlRIRiwkIjM0Tl8vKXBabSRvRi83JEZmdCQiM0l5N1ZSYUddekYvNyQkIjMvJylbVVhDTHRJRiwkIjNrIm9VYUAjNDsjKkYvNyRGaXQkIjN0Lj0zYj9dYjVGLC1GXHU2JkZedUZfdUZcX2xGX3UtRiY2JDdlcDckRiokIjNHJm8hRzxvI28+IkYsNyQkITNFWzxtZFBOMkpGLCQiMysjUWFkQmgmcDVGLDckRl1pbCQiM2coM3FJOTZgWypGLzckJCEzI1JDbCR5ZigpUUlGLCQiM2ZqM1BSVndOJClGLzckRjEkIjNoKkduZ1sqUVhzRi83JCQhM2wmUjNSLFlbKEhGLCQiM3dtZDcnKkdlVmpGLzckRmVpbCQiM3FgWkNKJ29WWyZGLzckJCEzSy83SGtRRTpIRiwkIjNBUiJScjBQbm0lRi83JEY0JCIzTWs0QCwoKW8qKVFGLzckJCEzIXBkTjJWbSU9R0YsJCIzTTBFJSlmWTgmRyNGLzckRjckIjNpZyN6WFFAJClwKUZneTckJCEzR0IuZ14mNFNvI0YsJCEzTllaajQrYllQRmd5NyRGOiQhMyRSQ09PPyYzXDlGLzckJCEzT1RLUCI0SCVcREYsJCEzOEQtKVJ1aipmQkYvNyRGPSQhM0VqKCo0XTwoRzckRi83JCQhM3c6TDIvNjE/Q0YsJCEzb08sRE4oZWxxJEYvNyRGQCQhM2BUJFE4VUAielRGLzckJCEzLz1aVmoiekxII0YsJCEzMVNgZFw4RGZYRi83JEZDJCEzbSpmJWVWXUJPW0YvNyQkITMzODYrKEc2Yz4jRiwkITNfMEVSIylIZVRcRi83JCQhM0lTXXUudEdpQEYsJCEzQkJGPz94XEFdRi83JCQhM2BuKilbP0wnKkdARiwkITNxdUEsNSY+KzMmRi83JEZGJCEzZSFILmNwc142JkYvNyQkITNCX1BUQ0FVaT9GLCQhM04mKilHJGVYJipHXkYvNyQkITNxNFlmNl4/SD9GLCQhM2BiWyFwMidcQV5GLzckJCEzPW5heCkqeilmKj5GLCQhMz8mKlxJeFh3JzQmRi83JEZJJCEzMEFDUno7cl9dRi83JCQhMykzLSQzPl5WJSo9RiwkITNsdyI0OGNJJTNcRi83JEZMJCEzZ1FTTjo5XClwJUYvNyRGTyQhMzMnM2J0JylbQD4lRi83JEZSJCEzbU9YK2E/LWNNRi83JEZVJCEzb2g/Jz0xYDJmI0YvNyRGWCQhMz5cbTJEPiVIbyJGLzckRmVuJCEzWyFvKnAxb1UnSClGZ3k3JEZobiQiMypbKj4oeiYpPUQjPUZneTckRltvJCIzcUtqQCQpZWAzNUYvNyRGXm8kIjNDM2EtMkdLLj5GLzckRmFvJCIzKnkob2FvPCw4RUYvNyRGZG8kIjNNImZLMFhrNEckRi83JEZnbyQiMyE9diEpKWU6YiN6JEYvNyRGam8kIjNWUnlHNEkhSD0lRi83JCQhM1QiNGdNYmxHKT5GLyQiMyVwRW1nM1IsSiVGLzckRl1wJCIzRio0XDtZLE9TJUYvNyQkITNiLkBaLyJcJHlwRmd5JCIzPCVmZUB5b2ZZJUYvNyRGYHAkIjNAai0wPHAwKVslRi83JCQiM3AwdCEzKUdGO21GZ3kkIjMkXGltY3AvI29XRi83JEZkcCQiM0oqZillWC4pW1MlRi83JCQiMzEyKip5SUtAZD5GLyQiMyZ6cyU+bElwOVZGLzckRmdwJCIza09eYyZ5YUI+JUYvNyRGanAkIjMyIVsrIm9ZaUBRRi83JEZdcSQiM0N3ayc+aj1VSCRGLzckRmBxJCIzKEdaLGZwJXpcRUYvNyRGY3EkIjNvV2xjei1uQD5GLzckRmZxJCIzOF4nMy5nJj07NUYvNyRGaXEkIjN4Kls2TWJAPFMiRmd5NyRGXHIkITNnU3dTXUQoKilIKUZneTckRl9yJCEzWCR5NjJoNTByIkYvNyRGYnIkITNAWmU8MSoqZVFFRi83JEZlciQhMyZRZFFiZ1k4VyRGLzckRmhyJCEzXSdSOSpcMyVSOyVGLzckRltzJCEzeGYnZjsmKTQjNFpGLzckJCIzL1VQbFxjIm8qPUYsJCEzZUZkZ01kZzlcRi83JEZecyQhM2Nlc2QmUSlRYl1GLzckJCIzaTArXGpNLygqPkYsJCEza15aMjchcHk0JkYvNyQkIjNIPmVEaTVpSD9GLCQhM0ZwYlg3IilwQV5GLzckJCIzJ0hqQDVtKT5pP0YsJCEzPXNlaV13KCpHXkYvNyRGYXMkITNbRCFHJWZdejpeRi83JCQiM0xQQiIpZjo0R0BGLCQhM1hkbytDJz43MyZGLzckJCIzLUdzJClmb1NoQEYsJCEzNlUvSzVRSkNdRi83JCQiM3M9QCcpZkBzJT4jRiwkITNKQDkwJFxhUyVcRi83JEZkcyQhM1pDRENhRVNSW0YvNyQkIjM1OnQjXEs7VEgjRiwkITM9J0hdeVIsYmIlRi83JEZncyQhMyNmMiFvSHYwalRGLzckJCIzTT9VRC9bIjRVI0YsJCEzWyRvNzNPQCQqcCRGLzckRmpzJCEzKy9GUlclPihISkYvNyQkIjMxPyM9KFtmQV5ERiwkITNRcTYnKSk0YXZMI0YvNyRGXXQkITNzYV4mUkp6Z1EiRi83JCQiMzZUSDR0NDEkbyNGLCQhMyU0QzxZRyI+NFJGZ3k3JEZgdCQiMy9AJzNjKVFEKlsoRmd5NyRGXmJtJCIzJT0qW2cocDxKOCNGLzckRmN0JCIzZy46WytlJz1xJEYvNyQkIjNEIylHPUpLejRIRiwkIjNJUGlJWUQuQFhGLzckRmZibSQiM1p5TUQjKmY2J1EmRi83JCQiM0I0M21xREp0SEYsJCIzazwiKVx0X0spSCdGLzckRmZ0JCIzIzNhIyo0bCMpKWVzRi83JCQiMztIQiJITS0jUklGLCQiM1huIT1rKCpSa00pRi83JEZeY20kIjMtIyk+akooM0dcKkYvNyQkIjMjSFdQemFpdTUkRiwkIjMyblMwW2MmKnA1Riw3JEZpdCQiMz8mKUdDSG8jbz4iRiwtRlx1NiZGXnUkIjFfTW1YJSllcWshIzskIjJ3bW9WKCllcWsiRixGYGluLSUrQVhFU0xBQkVMU0c2JFEieDYiUSFGZmluLSUlVklFV0c2JDskISthRWZUSiEiKiQiK2FFZlRKRl5qbjskITIpKio9SkUnWzA1IkZfaW4kIjJHeS1vU3Q9QyJGX2luLSUqVEhJQ0tORVNTRzYjIiIj The orthonormal basis we have constructed is closely related to the Legendre polynomials , one of several families of orthogonal polynomials . Now we can turn to our minimization problem, which becomes easy with an orthonormal basis. Recall that the minimum with respect to our norm is just the projection of the target function onto the subspace spanned by the basis. With an orthogonal basis, this is trivial since it is just a sum of inner products times each basis element.
While we know that the inner product of NiMtSSRjb3NHNiI2I0kieEdGJQ== with NiMmSSJmRzYiNiMiIiE= , NiMmSSJmRzYiNiMiIiI= , and NiMmSSJmRzYiNiMiIiQ= is zero, we write them anyway. approx := collect(
InnerProd(cos(x),f[0])*f[0] + InnerProd(cos(x),f[1])*f[1] + InnerProd(cos(x),f[2])*f[2] + InnerProd(cos(x),f[3])*f[3] + InnerProd(cos(x),f[4])*f[4], x);
evalf(approx); NiM+SSdhcHByb3hHNiIsKiooLCYqJEkjUGlHSSpwcm90ZWN0ZWRHRisiIiNGLCEjQCIiIkYuRiohIilJInhHRiUiIiUjISV2OiIiKSomLCYqJEYqISIlIyEjWEYsKiZGKEYuRiohIicjIiR2J0YxRi5GMEYsRi4qJEYqISIjIyIjOkYsKiZGKEYuRipGOCMhJE4iRjQ= NiMsKCokSSJ4RzYiIiIlJCIrPz8pZmgjISM2KiRGJSIiIyQhKyI0eUdfJSEjNSQiKyMqUUUkeSpGLyIiIg== The fit is pretty good, differing by less than .07 on the whole interval. plot([cos(x),approx], x=-Pi..Pi, color=[red, green], thickness=2); LSUlUExPVEc2Jy0lJ0NVUlZFU0c2JDdlbjckJCEzKioqKio0dGsjZlRKISM8JCEiIiIiITckJCEzNCdcOCFvWzZ0SUYsJCEzRGRSI1w3amwoKiohIz03JCQhM3ciKnByKTNQWSskRiwkITNuZ0BwbEJPMSoqRjU3JCQhMykqKnoqNFJcMFhIRiwkITN5PyRHPV8mWzIpKkY1NyQkITNAM0VbKnlzYSlHRiwkITMuLz8oUjMtUW4qRjU3JCQhM2ZYJikpPjJnOXYjRiwkITNrcmpSZGtlWyMqRjU3JCQhMyk0NTc3LmZsaCNGLCQhM20iKSkzOyozMmAnKUY1NyQkITN2IlFNOjoqSCNbI0YsJCEzWHQlKik0K19VIXpGNTckJCEzeVxBaGNJI3lOI0YsJCEzWmRWS0xUZCMzKEY1NyQkITMmZT1kLUZOKkdBRiwkITNnel5TM2RYO2hGNTckJCEzdyUqR0JQJFJjNCNGLCQhM0hCUnI7Ong1XUY1NyQkITNtQ2omZikzeGk+RiwkITNXQi9SLydSLCNRRjU3JCQhMydvcjRBTio0RT1GLCQhM0U1TD9PaVFEREY1NyQkITNRUVEqMyoqPWRxIkYsJCEzSlhqZCR5T15NIkY1NyQkITNkIT5qZ0AqPnE6RiwkIjMpPjpPQSdwNXJmISNANyQkITNgO20sJSlIN005RiwkIjNnKillKVFSI1tpOEY1NyQkITN2KDRGLEspKUhJIkYsJCIzWyI0Pz89Iz1ZRUY1NyQkITNTKjQhUiNbMFI9IkYsJCIzJVJmdCdmKjNKeCRGNTckJCEzcFlAUXFXSVU1RiwkIjMjW1oyJClIOkIvJkY1NyQkITNXcDN3JFIpXEIjKkY1JCIzSS0yZShcKltSZ0Y1NyQkITNFSmM4bzM5R3lGNSQiM0ZWPSVSTzskKjMoRjU3JCQhMyc9W0JQLThJZidGNSQiMzpQIVIjPlU/L3pGNTckJCEzQ0IzPEA/KXlCJkY1JCIzYmUyKSlcbklmJylGNTckJCEzYCopUjxZb1paUkY1JCIzYV0zMiFHTzRCKkY1NyQkITM4VHFYPVcyLEVGNSQiMzpVcFo4TmlqJypGNTckJCEzVCI0Z01ibEcpPkY1JCIzdzVwIyp5YTAvKSpGNTckJCEzcVRKWSlvY1lPIkY1JCIzeF90SEArLjIqKkY1NyQkITNiLkBaLyJcJHlwISM+JCIzWU5JRCw3bXYqKkY1NyQkITM4LHo3VktKLEohIz8kIjNnJlE+JDQ+JioqKioqRjU3JCQiM3AwdCEzKUdGO21GYXQkIjNdai4jKVwvN3kqKkY1NyQkIjM6U0ZmM3hFYThGNSQiMyR5O2J1LFElMyoqRjU3JCQiMzEyKip5SUtAZD5GNSQiM1BxN2xnazI0KSpGNTckJCIzJ1IyKClIdmUsYyNGNSQiM24sS1BsYTF1JypGNTckJCIzSWN3UjxUYmlRRjUkIjNwOmEuNkFFaiMqRjU3JCQiMyo0Kj1Kb2YwM19GNSQiMzFYMmJJaD11JylGNTckJCIzcF0xWElxT0NsRjUkIjMmSDglcHF4MVl6RjU3JCQiMyVIb0JsbWx6eihGNSQiMzliQih5JGZjNXJGNTckJCIzXnUnZiM0KXo/QCpGNSQiMzFkdyRHVidlW2dGNTckJCIzOUh2PEVDRls1RiwkIjNrLiUzSVUpbyEqXEY1NyQkIjNObGolRyUzJVI9IkYsJCIzMHUhXCkpWyJ5c1BGNTckJCIzXC0pUkVmd29JIkYsJCIzKUhWZDFmbCczRUY1NyQkIjNjTTJ2UXlGVDlGLCQiM19mRy4oeW06SCJGNTckJCIzeigzMnMkKlF4YyJGLCQiMyMpZmdlQ1tRZElGZ3Q3JCQiMzcrSz9CcyMqKnAiRiwkITMqXHVmd0JCeEciRjU3JCQiMzwjSCRlTWE7SD1GLCQhMzg5YzRHOC9iREY1NyQkIjMlPj5DWidlWWs+RiwkITNfLT1xbSQpek5RRjU3JCQiM2tZdXlmaXglNCNGLCQhMyNvSCFlYzBJLl1GNTckJCIzVDRxKSlmdS5HQUYsJCEzJD40Qip6Lk40aEY1NyQkIjN6P3cnKio9Jj5nQkYsJCEzIXlzWHZPKyQqNChGNTckJCIzISo+M2E9V2oiWyNGLCQhM1l5bC53JHksIXpGNTckJCIzQT9jKil5dSIzaSNGLCQhMzMzKkcyXVBWbilGNTckJCIzLGktSG5XSVhGRiwkITNHOkRVUzQrRCMqRjU3JCQiM3Y9UldoTi55R0YsJCEzdjhbLWcpKW9hJypGNTckJCIzdVg9IzQhSGJUSEYsJCEzbiE9NFZ6J2UrKSpGNTckJCIzc3MoKlJTQTIwSUYsJCEzKFFueTNgYnAhKipGNTckJCIzLycpW1VYQ0x0SUYsJCEzXXVlYms8cncqKkY1NyQkIjMhKSoqKlwvbCNmVEpGLEYtLSUmQ09MT1JHNiYlJFJHQkckIiM1Ri4kRi9GLkZqXWwtRiY2JDdcbzckRiokITNIKEh1SGE1UFAqRjU3JEYxJCEzVSdlIUdHWDEqZipGNTckRjckITMnUnBUMCNlIXlzKkY1NyQkITNsJlIzUixZWyhIRiwkITNhcSlSK0c4X3YqRjU3JEY8JCEzTTF3OzVHKmZ3KkY1NyQkITNLLzdIa1FFOkhGLCQhMyVRRmknZWRqZygqRjU3JEZBJCEzKVtbJCl6X0cnUigqRjU3JCQhMyFwZE4yVm0lPUdGLCQhM29sIzNFenZ3aipGNTckRkYkITM9Kj1OQCJRT2slKkY1NyQkITNHQi5nXiY0U28jRiwkITNPUWBNYHA6QiMqRjU3JEZLJCEzUyVHYyNmej8/KilGNTckRlAkITNsbEwmXE46TjopRjU3JEZVJCEzJypmYVQsUyJmRihGNTckRlokITNTdUBdOWk8SWlGNTckRmluJCEzJz5oXSJ5clRNXUY1NyRGXm8kITMhKTQueHozWmVQRjU3JEZjbyQhMzpzMixtcigqKlEjRjU3JEZobyQhMyFwYHNSPi46OyJGNTckRl1wJCIzKTRLYDokcCo+QSNGYXQ3JEZjcCQiM0xiIjQtcSVmKGUiRjU3JEZocCQiMzA2IjNicGklZUdGNTckRl1xJCIzdT43Iz0yInlkUkY1NyRGYnEkIjMqcGApM0YpcCR5XkY1NyRGZ3EkIjNbWVgjb1hbWzcnRjU3JEZcciQiM0JQN1I5JSopKTRyRjU3JEZhciQiMywiSHE5MiZwbXlGNTckRmZyJCIzcWlFdiNRJTNpJilGNTckRltzJCIzL2RXQFdeJFszKkY1NyRGYHMkIjN5by93QT9ZeSUqRjU3JEZlcyQiMyg0c3FMeVJlZypGNTckRmpzJCIzWUk/XkJjNypwKkY1NyRGX3QkIjMjMzs3PipbQ2goKkY1NyRGZXQkIjMpcChmRDIvQSR5KkY1NyRGW3UkIjNtdkFuIyoqcE13KkY1NyRGYHUkIjMxN20hR1krL3EqRjU3JEZldSQiMygqW1BTSTFSNScqRjU3JEZqdSQiM0BdMyQqUidSelsqRjU3JEZfdiQiM05bXiE+KFFJOSIqRjU3JEZkdiQiM2c+U29LN3R2JilGNTckRml2JCIzJyo+I0dnJlxSMHpGNTckRl53JCIzd1MhW2p5NChIckY1NyRGY3ckIjNvNV8ob25MTTgnRjU3JEZodyQiM18lKnpyO0UySF5GNTckRl14JCIzKnlDLT15anUmUkY1NyRGYngkIjNEKDRBYDgxOyNHRjU3JEZneCQiMyZveFl5I1F6OzpGNTckRlx5JCIzJ3pdMylRNnVyQ0ZhdDckRmF5JCEzSXJBSmVXQC02RjU3JEZmeSQhM2Q5Im9yaSEzQENGNTckRlt6JCEzYmFLci8jXF54JEY1NyRGYHokITNML25EWGhPRV1GNTckRmV6JCEzXiNlKWZ1W1pBaUY1NyRGanokITMpXDBUT2laUkgoRjU3JEZfW2wkITN3V3Q2cmdBXCIpRjU3JEZkW2wkITMpKjN3WCFHaTUlKilGNTckRmlbbCQhM056VFhXJVtdVypGNTckJCIzUSE0blYscDsiR0YsJCEzViJmdz1jZEtpKkY1NyRGXlxsJCEzcCMqWzk2cyo+dCpGNTckJCIzRCMpRz1KS3o0SEYsJCEzdCZlS0UhcCR6dipGNTckRmNcbCQhM1k0KG8kKWUhPm0oKkY1NyQkIjNCNDNtcURKdEhGLCQhM18lXChvSilwaHYqRjU3JEZoXGwkITNuLFE5Myl6c3MqRjU3JEZdXWwkITNHemxOIXkvJilmKkY1NyRGYl1sJCEzNSFlSS5gNVBQKkY1LUZlXWw2JkZnXWxGal1sRmhdbEZqXWwtJStBWEVTTEFCRUxTRzYkUSJ4NiJRIUZiW20tJSVWSUVXRzYkOyQhK2FFZlRKISIqJCIrYUVmVEpGalttOyQhMilRJz1RISoqKipSNSEjOyQiMnVkXVo0JioqUjVGYFxtLSUqVEhJQ0tORVNTRzYjIiIj For comparison, let's look at the Taylor polynomial for the cosine on the same interval. Note that while the Taylor polynomial fits much better for small values, over the whole interval it is much worse. tailor:=convert(taylor(cos(x),x,5),polynom); plot([cos(x), tailor], x=-Pi..Pi, color=[red, blue], thickness=2); NiM+SSd0YWlsb3JHNiIsKCIiIkYnKiRJInhHRiUiIiMjISIiRioqJEYpIiIlI0YnIiND LSUlUExPVEc2KC0lJ0NVUlZFU0c2JDdlbjckJCEzKioqKio0dGsjZlRKISM8JCEiIiIiITckJCEzNCdcOCFvWzZ0SUYsJCEzRGRSI1w3amwoKiohIz03JCQhM3ciKnByKTNQWSskRiwkITNuZ0BwbEJPMSoqRjU3JCQhMykqKnoqNFJcMFhIRiwkITN5PyRHPV8mWzIpKkY1NyQkITNAM0VbKnlzYSlHRiwkITMuLz8oUjMtUW4qRjU3JCQhM2ZYJikpPjJnOXYjRiwkITNrcmpSZGtlWyMqRjU3JCQhMyk0NTc3LmZsaCNGLCQhM20iKSkzOyozMmAnKUY1NyQkITN2IlFNOjoqSCNbI0YsJCEzWHQlKik0K19VIXpGNTckJCEzeVxBaGNJI3lOI0YsJCEzWmRWS0xUZCMzKEY1NyQkITMmZT1kLUZOKkdBRiwkITNnel5TM2RYO2hGNTckJCEzdyUqR0JQJFJjNCNGLCQhM0hCUnI7Ong1XUY1NyQkITNtQ2omZikzeGk+RiwkITNXQi9SLydSLCNRRjU3JCQhMydvcjRBTio0RT1GLCQhM0U1TD9PaVFEREY1NyQkITNRUVEqMyoqPWRxIkYsJCEzSlhqZCR5T15NIkY1NyQkITNkIT5qZ0AqPnE6RiwkIjMpPjpPQSdwNXJmISNANyQkITNgO20sJSlIN005RiwkIjNnKillKVFSI1tpOEY1NyQkITN2KDRGLEspKUhJIkYsJCIzWyI0Pz89Iz1ZRUY1NyQkITNTKjQhUiNbMFI9IkYsJCIzJVJmdCdmKjNKeCRGNTckJCEzcFlAUXFXSVU1RiwkIjMjW1oyJClIOkIvJkY1NyQkITNXcDN3JFIpXEIjKkY1JCIzSS0yZShcKltSZ0Y1NyQkITNFSmM4bzM5R3lGNSQiM0ZWPSVSTzskKjMoRjU3JCQhMyc9W0JQLThJZidGNSQiMzpQIVIjPlU/L3pGNTckJCEzQ0IzPEA/KXlCJkY1JCIzYmUyKSlcbklmJylGNTckJCEzYCopUjxZb1paUkY1JCIzYV0zMiFHTzRCKkY1NyQkITM4VHFYPVcyLEVGNSQiMzpVcFo4TmlqJypGNTckJCEzVCI0Z01ibEcpPkY1JCIzdzVwIyp5YTAvKSpGNTckJCEzcVRKWSlvY1lPIkY1JCIzeF90SEArLjIqKkY1NyQkITNiLkBaLyJcJHlwISM+JCIzWU5JRCw3bXYqKkY1NyQkITM4LHo3VktKLEohIz8kIjNnJlE+JDQ+JioqKioqRjU3JCQiM3AwdCEzKUdGO21GYXQkIjNdai4jKVwvN3kqKkY1NyQkIjM6U0ZmM3hFYThGNSQiMyR5O2J1LFElMyoqRjU3JCQiMzEyKip5SUtAZD5GNSQiM1BxN2xnazI0KSpGNTckJCIzJ1IyKClIdmUsYyNGNSQiM24sS1BsYTF1JypGNTckJCIzSWN3UjxUYmlRRjUkIjNwOmEuNkFFaiMqRjU3JCQiMyo0Kj1Kb2YwM19GNSQiMzFYMmJJaD11JylGNTckJCIzcF0xWElxT0NsRjUkIjMmSDglcHF4MVl6RjU3JCQiMyVIb0JsbWx6eihGNSQiMzliQih5JGZjNXJGNTckJCIzXnUnZiM0KXo/QCpGNSQiMzFkdyRHVidlW2dGNTckJCIzOUh2PEVDRls1RiwkIjNrLiUzSVUpbyEqXEY1NyQkIjNObGolRyUzJVI9IkYsJCIzMHUhXCkpWyJ5c1BGNTckJCIzXC0pUkVmd29JIkYsJCIzKUhWZDFmbCczRUY1NyQkIjNjTTJ2UXlGVDlGLCQiM19mRy4oeW06SCJGNTckJCIzeigzMnMkKlF4YyJGLCQiMyMpZmdlQ1tRZElGZ3Q3JCQiMzcrSz9CcyMqKnAiRiwkITMqXHVmd0JCeEciRjU3JCQiMzwjSCRlTWE7SD1GLCQhMzg5YzRHOC9iREY1NyQkIjMlPj5DWidlWWs+RiwkITNfLT1xbSQpek5RRjU3JCQiM2tZdXlmaXglNCNGLCQhMyNvSCFlYzBJLl1GNTckJCIzVDRxKSlmdS5HQUYsJCEzJD40Qip6Lk40aEY1NyQkIjN6P3cnKio9Jj5nQkYsJCEzIXlzWHZPKyQqNChGNTckJCIzISo+M2E9V2oiWyNGLCQhM1l5bC53JHksIXpGNTckJCIzQT9jKil5dSIzaSNGLCQhMzMzKkcyXVBWbilGNTckJCIzLGktSG5XSVhGRiwkITNHOkRVUzQrRCMqRjU3JCQiM3Y9UldoTi55R0YsJCEzdjhbLWcpKW9hJypGNTckJCIzdVg9IzQhSGJUSEYsJCEzbiE9NFZ6J2UrKSpGNTckJCIzc3MoKlJTQTIwSUYsJCEzKFFueTNgYnAhKipGNTckJCIzLycpW1VYQ0x0SUYsJCEzXXVlYms8cncqKkY1NyQkIjMhKSoqKlwvbCNmVEpGLEYtLSUmQ09MT1JHNiYlJFJHQkckIiM1Ri4kRi9GLkZqXWwtRiY2JDdnbzckRiokIjM+T1pbSiIqNFI3RjU3JEYxJCEzTVRCUFZzJmZ4JkZndDckRjckITN6JVFcXihvMSE9IkY1NyRGPCQhMzRFJEdicXBALSNGNTckRkEkITNfaFcqKTQhNGV1I0Y1NyQkITMhcGROMlZtJT1HRiwkITNkR3BFITNnZFUkRjU3JEZGJCEzKCp6KnA4XGhBKFJGNTckJCEzR0IuZ14mNFNvI0YsJCEzMXonZWRPJzMnUiVGNTckRkskITNDOEcoKj1pXixaRjU3JCQhMzxyRUhoUypIZSNGLCQhM0ZUI2Y/TD4+IltGNTckJCEzT1RLUCI0SCVcREYsJCEzdSdHLTdrL2cqW0Y1NyQkITNiNlFYQFQnZV4jRiwkITNWOXNVWkR1YVxGNTckRlAkITMmektAJGYyNCopXEY1NyQkITN2W1EheTchPV5DRiwkITNJJGVpaVNyKioqXEY1NyQkITN3OkwyLzYxP0NGLCQhM2d6LSspR1Y5KlxGNTckJCEzeCN5VS4zVSopUSNGLCQhMyVbI0g3IVJTVSdcRjU3JEZVJCEzbUFtTj90Mz5cRjU3JCQhMy89WlZqInpMSCNGLCQhM0whXCozKT16OnglRjU3JEZaJCEzbDdZJnkzcmpiJUY1NyRGaW4kITNtK2g1TU9DQVJGNTckRl5vJCEzKDRpISpSUyFReUlGNTckRmNvJCEzV2lfLGl1JCpSP0Y1NyRGaG8kITM6RDxlcil5LS0iRjU3JEZdcCQiMyZHU0E/bUFAMCNGYXQ3JEZjcCQiM3lucSplXXgqeTlGNTckRmhwJCIzIlJUclUoKUdAciNGNTckRl1xJCIzYFhjaVY1VDVRRjU3JEZicSQiM0Mta1ZpRnlmXUY1NyRGZ3EkIjMtRmBaI0g3ei8nRjU3JEZcciQiMyVlNSkqKT14WiM0KEY1NyRGYXIkIjNMOGx1MWhMMHpGNTckRmZyJCIzeSs6JVI7I2ZmJylGNTckRltzJCIzVlVmZSZvKSk0QipGNTckRmBzJCIzSTsuVjR5aWonKkY1NyRGZXMkIjNPUCEqXEFqMC8pKkY1NyRGanMkIjNdSjIoNDZJcSEqKkY1NyRGX3QkIjNNP28mRz9oYygqKkY1RmR0NyRGW3UkIjNfSWApNFg/InkqKkY1NyRGYHUkIjNrMiw2LiJRJTMqKkY1NyRGZXUkIjNTViRSM0N4ITQpKkY1NyRGanUkIjNgbydHO1BwU24qRjU3JEZfdiQiMzUlSCczNiMzTEUqRjU3JEZkdiQiM0d0S0pZPll1JylGNTckRml2JCIzOlEuREw0OFp6RjU3JEZedyQiM1U7JT1tOGJPNihGNTckRmN3JCIzNEgiPi08WnAwJ0Y1NyRGaHckIjMhR3R2JD4xdzNdRjU3JEZdeCQiM0ZtQWgtLTQ1UUY1NyRGYngkIjMzTngvNSkqeXZFRjU3JEZneCQiM1JFcikzPFw6VCJGNTckRlx5JCIzPS9DISkpXGgoekFGYXQ3JEZheSQhM0ctTHZtLjwkcCpGYXQ3JEZmeSQhM1leXGJTInlaMSNGNTckRlt6JCEzL0N4PT9ZRiE0JEY1NyRGYHokITM8TCdbdCg0UjxSRjU3JEZleiQhMygqb2JTP08jSGIlRjU3JCQiMzU6dCNcSztUSCNGLCQhM15QVEQiKTRtdFpGNTckRmp6JCEzbWFccUpdOUJcRjU3JCQiM3k/NDYoKlxiIVIjRiwkITNXJipbQzN3NG1cRjU3JCQiM00/VUQvWyI0VSNGLCQhMzNjKVFXcUg+KlxGNTckJCIzISo+dlI2WUZeQ0YsJCEzQ0xDQT4ibyoqKlxGNTckRl9bbCQhMyRmYCMqZTpKJiopXEY1NyQkIjM/PyZITz1Ja14jRiwkITM0LkgzPm8mUiZcRjU3JCQiMzE/Iz0oW2ZBXkRGLCQhM1djSU5xUjojKltGNTckJCIzIyo+cCFRckBnZSNGLCQhM295TzlRXjAuW0Y1NyRGZFtsJCEzL0FXJj0+eWJvJUY1NyQkIjM2VEg0dDQxJG8jRiwkITNuVjx1SUs9LFdGNTckRmlbbCQhM3p6RmdFIm9nLCVGNTckJCIzUSE0blYscDsiR0YsJCEzcE1aJj0mKlJyWyRGNTckRl5cbCQhM2UwMTA2JVwiR0dGNTckRmNcbCQhMz1ZZ180KT16MSNGNTckRmhcbCQhMyU9U1p5XWxNPCJGNTckRl1dbCQhMyFlZ2E7bWw/UiZGZ3Q3JEZiXWwkIjM+K1w1Jj4qNFI3RjUtRmVdbDYmRmddbEZqXWxGal1sRmhdbC0lK0FYRVNMQUJFTFNHNiRRIng2IlEhRmZebS0lJVZJRVdHNiQ7JCErYUVmVEohIiokIithRWZUSkZeX207JCEyKVEnPVEhKioqKlI1ISM7JCIydWRdWjQmKipSNUZkX20tJSpUSElDS05FU1NHNiMiIiMtJShTQ0FMSU5HRzYjJS5VTkNPTlNUUkFJTkVERw== More explictly, we can compare the norms of the differences, and see the dramatic difference. We also show a plot of the two differences on the same axes. evalf(Norman(cos(x)-approx));
evalf(Norman(cos(x)-tailor)); NiMkIisnKVJoLFkhIzY= NiMkIisnUiUzJip6ISM1 plot([cos(x)-approx, cos(x)-tailor], x=-Pi..Pi, color=[green,blue], thickness=2); LSUlUExPVEc2Jy0lJ0NVUlZFU0c2JDdhcjckJCEzKioqKio0dGsjZlRKISM8JCEzUVVxRHFYKkdFJyEjPjckJCEzTnVlWy1LWkNKRiwkITMpRzwpKlxbejNmJkYvNyQkITNFWzxtZFBOMkpGLCQhMydHYiZwYGsnRyZcRi83JCQhMz1BdyRHSk0tNCRGLCQhM3Q0OWVObyl5TSVGLzckJCEzNCdcOCFvWzZ0SUYsJCEzMThQVm1mKVx4JEYvNyQkITMscSQqPUJhKmYwJEYsJCEzMkldR19KQUxLRi83JCQhMyNSQ2wkeWYoKVFJRiwkITNbb3poZztuQEZGLzckJCEzJXk2VE5gYzwtJEYsJCEzKkcyXCE0KD0lUkFGLzckJCEzdyIqcHIpM1BZKyRGLCQhMywkby86WGxieSJGLzckJCEzbCZSM1IsWVsoSEYsJCEzQycqZTo+TCczMSJGLzckJCEzKSoqeio0UlwwWEhGLCQhMy4uOzJtNkZcVCEjPzckJCEzSy83SGtRRTpIRiwkIjN5RSN6U2shem06RlxvNyQkITNAM0VbKnlzYSlHRiwkIjMoeS5bNlNXRWUnRlxvNyQkITNiIzQ0LGhwPiZHRiwkIjMkb05YLSJ6KlE5IkYvNyQkITMhcGROMlZtJT1HRiwkIjNgM0wvSi4tXzpGLzckJCEzRGg/T15LJ1x5I0YsJCIzcWIwP0A8PikpPUYvNyQkITNmWCYpKT4yZzl2I0YsJCIzYXEiKVFaTnhkQEYvNyQkITNHQi5nXiY0U28jRiwkIjNHUnpYbFFHPkRGLzckJCEzKTQ1NzcuZmxoI0YsJCIzakFTWncxUHJFRi83JCQhMydlUV9pYXcoKmYjRiwkIjNTP1xBJkgwNW8jRi83JCQhMzxyRUhoUypIZSNGLCQiM0MrRUNFISpbIW8jRi83JCQhM1xjSEx3OkBtREYsJCIzYWFnPm1QTXFFRi83JCQhM09US1AiNEglXERGLCQiM2xAK2gvKHo1bCNGLzckJCEzYjZRWEBUJ2VeI0YsJCIzNGU8Mkl5PChlI0YvNyQkITN2IlFNOjoqSCNbI0YsJCIzOTEqUSdSTmojXCNGLzckJCEzdzpMMi82MT9DRiwkIjN2JzRWRXB1JltBRi83JCQhM3lcQWhjSSN5TiNGLCQiMyoqPjUiNG8pUkw+Ri83JCQhMy89WlZqInpMSCNGLCQiM1Q7UXc+PjNgOkYvNyQkITMmZT1kLUZOKkdBRiwkIjM5VipwNDEwczgiRi83JCQhM0lTXXUudEdpQEYsJCIzYEwjKWZOX2YiKW9GXG83JCQhM3clKkdCUCRSYzQjRiwkIjNFPilvTzltWE8jRlxvNyQkITNxNFlmNl4/SD9GLCQhMzU9I2YieT9fPz9GXG83JCQhM21DaiZmKTN4aT5GLCQhM3VtOSxpQyhvOydGXG83JCQhMydvcjRBTio0RT1GLCQhMy0nUUQ+cSE0YThGLzckJCEzUVFRKjMqKj1kcSJGLCQhMyk+M1FnKmVMTz1GLzckJCEzZjkmeU01ZnpqIkYsJCEzMi1mKUdMdCpHP0YvNyQkITNkIT5qZ0AqPnE6RiwkITMjPic0JD4nZUdpQEYvNyQkITMxWjowZSw9TzpGLCQhMyoqeVhtW0pyMUFGLzckJCEzYi4qUis1aEBdIkYsJCEzcE1tJCpbJik+T0FGLzckJCEzMGcjRz8vVSJvOUYsJCEzYHEjcCl6Ij00RCNGLzckJCEzYDttLCUpSDdNOUYsJCEzc25FQmpJN15BRi83JCQhMyVvQlchPSRSOFMiRiwkITNLJDNgO28nKXlCI0YvNyQkITM5ZD0yX2NibzhGLCQhM3ZgUi12VSI9QCNGLzckJCEzWHglKjQnKT54TjhGLCQhMyoqXFFBXFdJdEBGLzckJCEzdig0RixLKSlISSJGLCQhMy8mPiEpWzgwRzcjRi83JCQhM1opZmU3IXBXVjdGLCQhM1BbdGAmb2c/KyNGLzckJCEzUyo0IVIjWzBSPSJGLCQhMydHRXc5N0BuJT1GLzckJCEzcFlAUXFXSVU1RiwkITNfRTEieUdYME8iRi83JCQhM1dwM3ckUilcQiMqISM9JCEzW21XUUNmKmVgKUZcbzckJCEzRUpjOG8zOUd5Rml6JCEzYWAmUlwvMHQwI0ZcbzckJCEzJz1bQlAtOElmJ0ZpeiQiM2RtWChveDk0diRGXG83JCQhM0NCMzxAPyl5QiZGaXokIjN4TiY0R3JPQXMqRlxvNyQkITNgKilSPFlvWlpSRml6JCIzPklSY2U4LGg5Ri83JCQhMzhUcVg9VzIsRUZpeiQiM3pHWjsyXGhePUYvNyQkITNUIjRnTWJsRyk+Rml6JCIzNi8+Y2JwOiMpPkYvNyQkITNxVEpZKW9jWU8iRml6JCIzUz1LJnkoUi96P0YvNyQkITMud15YKnpYNy4iRml6JCIzZEpNJipRZHo7QEYvNyQkITNiLkBaLyJcJHlwRi8kIjNbVCgzTTRqVDkjRi83JCQhMyRvVyNSOS1DV09GLyQiM21cT1ZNcy1oQEYvNyQkITM4LHo3VktKLEpGXG8kIjNPcVNqP11KbkBGLzckJCIzemRzQ3kyMmBKRi8kIjN3RCcpXG0qPUU7I0YvNyQkIjNwMHQhMylHRjttRi8kIjNedDNbclhdWUBGLzckJCIzT05uTClcWnorIkZpeiQiMzlyXE5Ncy8+QEYvNyQkIjM6U0ZmM3hFYThGaXokIjMjPWEmW1liUCEzI0YvNyQkIjMxMioqeUlLQGQ+Rml6JCIzOjRfWi0kZW8pPkYvNyQkIjMnUjIoKUh2ZSxjI0ZpeiQiM3E0TlVhI2U3Jz1GLzckJCIzSWN3UjxUYmlRRml6JCIzaW9FSSJSJGUqWyJGLzckJCIzKjQqPUpvZjAzX0ZpeiQiMyplUXNteSpbWCkqRlxvNyQkIjNwXTFYSXFPQ2xGaXokIjNkUzZmbTlHblNGXG83JCQiMyVIb0JsbWx6eihGaXokITN6NCdvdiVbUTk+RlxvNyQkIjNedSdmIzQpej9AKkZpeiQhMzU1YXYuV3MlWylGXG83JCQiMzlIdjxFQ0ZbNUYsJCEzZzhmNFA+JVFRIkYvNyQkIjNObGolRyUzJVI9IkYsJCEzdVA8YEhII28lPUYvNyQkIjMkUjNWeHIzYUMiRiwkITM3U256SSU0bSsjRi83JCQiM1wtKVJFZndvSSJGLCQhMzNVbWtZYVNIQEYvNyQkIjNeTnY7L3BaUzhGLCQhM3BvL04+IXomekBGLzckJCIzYG9fcDpzMnU4RiwkITM0dXR5JSl6MzxBRi83JCQiM2IsSUFGdm4yOUYsJCEzSVFmXFUiZTlDI0YvNyQkIjNjTTJ2UXlGVDlGLCQhM1VzIlIiMy9GX0FGLzckJCIzKUcjW084SipHWiJGLCQhM08mR2M5JjR1XEFGLzckJCIzPTYqeXpRM1hdIkYsJCEzWT01eCR5U1lCI0YvNyQkIjNbKipIZmlPN086RiwkITM+KSl6WGpXeDFBRi83JCQiM3ooMzJzJCpReGMiRiwkITNzLCpcamwtZzsjRi83JCQiMyZSOTAtM0xRaiJGLCQhM0Iob1JWbid6UT9GLzckJCIzNytLP0JzIyoqcCJGLCQhMyRcdHVNeigzYj1GLzckJCIzPCNIJGVNYTtIPUYsJCEzJ1wqXEY0cWdSOEYvNyQkIjMlPj5DWidlWWs+RiwkITNZKnlhKSk+O1wxJ0ZcbzckJCIzSD5lRGk1aUg/RiwkITN3J3lWMSRRcCQqPkZcbzckJCIza1l1eWZpeCU0I0YsJCIzVi4ya24pZWxJI0ZcbzckJCIzLUdzJClmb1NoQEYsJCIzXCZvVTk6YDwjb0ZcbzckJCIzVDRxKSlmdS5HQUYsJCIzKD4iXHZZXENKNkYvNyQkIjM1OnQjXEs7VEgjRiwkIjNaM3NUZCZvd2IiRi83JCQiM3o/dycqKj0mPmdCRiwkIjMmZURgNGNzayU+Ri83JCQiM00/VUQvWyI0VSNGLCQiM1dIbjN6MlhfQUYvNyQkIjMhKj4zYT1XaiJbI0YsJCIzIWZqMjMmcFohXCNGLzckJCIzPz8mSE89SWteI0YsJCIzQThHeDEsXyllI0YvNyQkIjMxPyM9KFtmQV5ERiwkIjMzWFBwP0ljYEVGLzckJCIzKipwREVKUWlvREYsJCIzISlcLm0xXVBzRUYvNyQkIjMjKj5wIVFyQGdlI0YsJCIzUltqS2UhNDhvI0YvNyQkIjNIcTdOJ2Y+TWcjRiwkIjNFJHBhI3BZenpFRi83JCQiM0E/YyopeXUiM2kjRiwkIjMvJCpwRyh6WnNtI0YvNyQkIjM2VEg0dDQxJG8jRiwkIjNtJ3pPVD8meUFERi83JCQiMyxpLUhuV0lYRkYsJCIzek5tSlNdWitBRi83JCQiM1EhNG5WLHA7IkdGLCQiM3BvVCpSN25laSJGLzckJCIzdj1SV2hOLnlHRiwkIjNwTngrN14kM3QoRlxvNyQkIjNEIylHPUpLejRIRiwkIjM/KW9Kd0EqNFJERlxvNyQkIjN1WD0jNCFIYlRIRiwkITM3IUdaU2Y/J1JNRlxvNyQkIjNCNDNtcURKdEhGLCQhM1FFeTpaMHRENUYvNyQkIjNzcygqUlNBMjBJRiwkITM7ZyhbdEFkbnoiRi83JCQiMzteZ2wiSFBALSRGLCQhM2NCbyk9X0spXEFGLzckJCIzO0hCIkhNLSNSSUYsJCEzNS1nWU5xOUpGRi83JCQiM2cyJ29UUm5pMCRGLCQhM2whem4mZThmVEtGLzckJCIzLycpW1VYQ0x0SUYsJCEzNm9IKj4lKXA/eSRGLzckJCIzW2s2bydcKFIhNCRGLCQhM1I6UVU1MV1gVkYvNyQkIjMjSFdQemFpdTUkRiwkITMiKnlrX3RiIm8mXEYvNyQkIjNPQFA+KmZGWDckRiwkITNJOyMqZU0vJ0hmJkYvNyQkIjMhKSoqKlwvbCNmVEpGLCQhM0I5VXAncCUqR0UnRi8tJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiM1RmZgbUZkYG0tRiY2JDdbbzckRiokITNidCVbSiIqNFI3IkYsNyRGNiQhM3NaWTVYckJjNUYsNyRGQCQhMyRIQiFcX04hKT0qKkZpejckRkokITM+KClHWXFjPDIkKkZpejckRlQkITMoZXhVMFwmSEUoKUZpejckRmhuJCEzcCUqKipIO2VKJnkoRml6NyRGY28kITMmPmF4UzIkKnojcEZpejckRl1wJCEzNDRxQnpPcmNnRml6NyRGZ3AkITNBI1JFZydcS3dfRml6NyRGXHEkITMrbTMvLEE5dlhGaXo3JEZhcSQhM0luZ2pzWWJeUkZpejckRmVyJCEzQ0hmOSVcWDhTJEZpejckRl9zJCEzXVgibztDaF4iSEZpejckRmlzJCEzNU54J0gib1tqQEZpejckRmN0JCEzb20wYj9ZM2c6Rml6NyRGXXUkITNcQXlnIyl5XykzIkZpejckRmd1JCEzRjchKSpSKyNmPHVGLzckRlx2JCEzI1FbISk9dShbYVtGLzckRmF2JCEzJyo0aSUqPiJ6JltLRi83JEZbdyQhMyFRLytDZjZDKj5GLzckRl94JCEzK3k8Nj82Jlw7IkYvNyRGY3kkITN3YkE4RCNwWWYnRlxvNyRGXXokITNeITMwX1IzLXQkRlxvNyRGYnokITNVUUYqR1RZbnUiRlxvNyRGZ3okITN3JXpDWSpbekElKSEjQDckRl1bbCQhMyh6LWppJlxOaEpGZ2VtNyRGYltsJCEzYiJHd3VdKCk9OCJGZ2VtNyRGZ1tsJCEzXEchPnVnU1QmRyEjQTckRlxcbCQhM0Q5RTU0OmJTXyEjQjckRmFcbCQhM3cjSDZVUGBmSCUhI0M3JEZbXWwkITNyci9TLHpMbiopISNFNyRGX15sJCIzdiF6W0FdWihlSiEjTTckRmNfbCQhMyVRISlwQClSXGwmKUZdZ203JEZdYGwkITNuWngkeVliaSFSRmlmbTckRmJgbCQhMydvRyV6KDMwK2clRmVmbTckRmdgbCQhMyk0bnhfaWQiZUZGYWZtNyRGXGFsJCEzazxcP2NEO2o1RmdlbTckRmFhbCQhMzImUWlndSk+KjMkRmdlbTckRmZhbCQhMyZcMnQ5UVAyTylGZ2VtNyRGW2JsJCEzJVIqR3RPJz5zIT1GXG83JEZgYmwkITMjZkM+alByM3QkRlxvNyRGamJsJCEzdSY9SSFSPlU3bkZcbzckRl5kbCQhM1luRWFRUSMpKj4iRi83JEZiZWwkITMlenpWaixCUyg+Ri83JEZcZmwkITNPYlQlKTQ/MSU9JEYvNyRGYWZsJCEzKnBoMWEoPWotXEYvNyRGZmZsJCEzKmZ6U15ZUF9YKEYvNyRGYGdsJCEzUGo7QnomNGYzIkZpejckRmpnbCQhM05Ddl5mblVjOkZpejckRmRobCQhM3ZyMiVlTGJoPCNGaXo3JEZeaWwkITNWVFM5P3NrNUhGaXo3JEZoaWwkITM2Y2wtZj9AOk1GaXo3JEZcW20kITMwJ1t1KTMkZigpKVJGaXo3JEZhW20kITNycVxgZjEmZWMlRml6NyRGZlttJCEzaE8oPlEiRyQqM19GaXo3JEZbXG0kITMlKlJDaSgqM2B0ZkZpejckRmBcbSQhM1I1VSgqWyVSbCNvRml6NyRGalxtJCEzKnA4JHklKXptS3hGaXo3JEZkXW0kITMpUkZKSS0hXEwoKUZpejckRl5ebSQhM14/NiUpeSdlR0oqRml6NyRGaF5tJCEzR0ckUno1InpBKipGaXo3JEZiX20kITNbUz47Xm9XYzVGLDckRlxgbSQhMy0hXDUmPio0UjciRiwtRmFgbTYmRmNgbUZkYG1GZGBtRmdgbS0lK0FYRVNMQUJFTFNHNiRRIng2IlEhRmVdbi0lJVZJRVdHNiQ7JCErYUVmVEohIiokIithRWZUSkZdXm47JCEyX083KHpWI3A5IiEjOyQiMWNRXnlddiMpXEYsLSUqVEhJQ0tORVNTRzYjIiIj Just to emphasize the ease of solving the minimization problem once the basis is in hand, here are a couple more examples. Naturally, since our basis was constructed specifically for the interval NiM3JCwkSSNQaUdJKnByb3RlY3RlZEdGJiEiIkYl , we must work there. However, it would be not hard to adapt this same basis to work on any interval. In order to avoid lots of typing, we write a Maple function that computes the inner product for us, and another that shows us the graphs. MakeApprox:=(f,basis)->collect(add(InnerProd(f,basis[i])*basis[i], i=1..nops(basis)),x):
ShowApprox:=proc(f,basis)
local a;
a:=MakeApprox(f,basis);
print(a);
print(evalf(a), " norm of difference=",evalf(Norman(a-f)));
plot([f, a], x=-Pi..Pi, color=[red, green], thickness=2);
end: First, just to check, we apply it to the same problem we just solved. ShowApprox(cos(x),OrthoBasis); NiMsKiooLCYqJEkjUGlHSSpwcm90ZWN0ZWRHRigiIiNGKSEjQCIiIkYrRichIilJInhHNiIiIiUjISV2OiIiKSomLCYqJEYnISIlIyEjWEYpKiZGJUYrRichIicjIiR2J0YvRitGLUYpRisqJEYnISIjIyIjOkYpKiZGJUYrRidGNiMhJE4iRjI= NiUsKCokSSJ4RzYiIiIlJCIrPz8pZmgjISM2KiRGJSIiIyQhKyI0eUdfJSEjNSQiKyMqUUUkeSpGLyIiIlE1fm5vcm1+b2Z+ZGlmZmVyZW5jZT1GJiQiK01ZZyxZRio= LSUlUExPVEc2Jy0lJ0NVUlZFU0c2JDdlbjckJCEzKioqKio0dGsjZlRKISM8JCEiIiIiITckJCEzNCdcOCFvWzZ0SUYsJCEzRGRSI1w3amwoKiohIz03JCQhM3ciKnByKTNQWSskRiwkITNuZ0BwbEJPMSoqRjU3JCQhMykqKnoqNFJcMFhIRiwkITN5PyRHPV8mWzIpKkY1NyQkITNAM0VbKnlzYSlHRiwkITMuLz8oUjMtUW4qRjU3JCQhM2ZYJikpPjJnOXYjRiwkITNrcmpSZGtlWyMqRjU3JCQhMyk0NTc3LmZsaCNGLCQhM20iKSkzOyozMmAnKUY1NyQkITN2IlFNOjoqSCNbI0YsJCEzWHQlKik0K19VIXpGNTckJCEzeVxBaGNJI3lOI0YsJCEzWmRWS0xUZCMzKEY1NyQkITMmZT1kLUZOKkdBRiwkITNnel5TM2RYO2hGNTckJCEzdyUqR0JQJFJjNCNGLCQhM0hCUnI7Ong1XUY1NyQkITNtQ2omZikzeGk+RiwkITNXQi9SLydSLCNRRjU3JCQhMydvcjRBTio0RT1GLCQhM0U1TD9PaVFEREY1NyQkITNRUVEqMyoqPWRxIkYsJCEzSlhqZCR5T15NIkY1NyQkITNkIT5qZ0AqPnE6RiwkIjMpPjpPQSdwNXJmISNANyQkITNgO20sJSlIN005RiwkIjNnKillKVFSI1tpOEY1NyQkITN2KDRGLEspKUhJIkYsJCIzWyI0Pz89Iz1ZRUY1NyQkITNTKjQhUiNbMFI9IkYsJCIzJVJmdCdmKjNKeCRGNTckJCEzcFlAUXFXSVU1RiwkIjMjW1oyJClIOkIvJkY1NyQkITNXcDN3JFIpXEIjKkY1JCIzSS0yZShcKltSZ0Y1NyQkITNFSmM4bzM5R3lGNSQiM0ZWPSVSTzskKjMoRjU3JCQhMyc9W0JQLThJZidGNSQiMzpQIVIjPlU/L3pGNTckJCEzQ0IzPEA/KXlCJkY1JCIzYmUyKSlcbklmJylGNTckJCEzYCopUjxZb1paUkY1JCIzYV0zMiFHTzRCKkY1NyQkITM4VHFYPVcyLEVGNSQiMzpVcFo4TmlqJypGNTckJCEzVCI0Z01ibEcpPkY1JCIzdzVwIyp5YTAvKSpGNTckJCEzcVRKWSlvY1lPIkY1JCIzeF90SEArLjIqKkY1NyQkITNiLkBaLyJcJHlwISM+JCIzWU5JRCw3bXYqKkY1NyQkITM4LHo3VktKLEohIz8kIjNnJlE+JDQ+JioqKioqRjU3JCQiM3AwdCEzKUdGO21GYXQkIjNdai4jKVwvN3kqKkY1NyQkIjM6U0ZmM3hFYThGNSQiMyR5O2J1LFElMyoqRjU3JCQiMzEyKip5SUtAZD5GNSQiM1BxN2xnazI0KSpGNTckJCIzJ1IyKClIdmUsYyNGNSQiM24sS1BsYTF1JypGNTckJCIzSWN3UjxUYmlRRjUkIjNwOmEuNkFFaiMqRjU3JCQiMyo0Kj1Kb2YwM19GNSQiMzFYMmJJaD11JylGNTckJCIzcF0xWElxT0NsRjUkIjMmSDglcHF4MVl6RjU3JCQiMyVIb0JsbWx6eihGNSQiMzliQih5JGZjNXJGNTckJCIzXnUnZiM0KXo/QCpGNSQiMzFkdyRHVidlW2dGNTckJCIzOUh2PEVDRls1RiwkIjNrLiUzSVUpbyEqXEY1NyQkIjNObGolRyUzJVI9IkYsJCIzMHUhXCkpWyJ5c1BGNTckJCIzXC0pUkVmd29JIkYsJCIzKUhWZDFmbCczRUY1NyQkIjNjTTJ2UXlGVDlGLCQiM19mRy4oeW06SCJGNTckJCIzeigzMnMkKlF4YyJGLCQiMyMpZmdlQ1tRZElGZ3Q3JCQiMzcrSz9CcyMqKnAiRiwkITMqXHVmd0JCeEciRjU3JCQiMzwjSCRlTWE7SD1GLCQhMzg5YzRHOC9iREY1NyQkIjMlPj5DWidlWWs+RiwkITNfLT1xbSQpek5RRjU3JCQiM2tZdXlmaXglNCNGLCQhMyNvSCFlYzBJLl1GNTckJCIzVDRxKSlmdS5HQUYsJCEzJD40Qip6Lk40aEY1NyQkIjN6P3cnKio9Jj5nQkYsJCEzIXlzWHZPKyQqNChGNTckJCIzISo+M2E9V2oiWyNGLCQhM1l5bC53JHksIXpGNTckJCIzQT9jKil5dSIzaSNGLCQhMzMzKkcyXVBWbilGNTckJCIzLGktSG5XSVhGRiwkITNHOkRVUzQrRCMqRjU3JCQiM3Y9UldoTi55R0YsJCEzdjhbLWcpKW9hJypGNTckJCIzdVg9IzQhSGJUSEYsJCEzbiE9NFZ6J2UrKSpGNTckJCIzc3MoKlJTQTIwSUYsJCEzKFFueTNgYnAhKipGNTckJCIzLycpW1VYQ0x0SUYsJCEzXXVlYms8cncqKkY1NyQkIjMhKSoqKlwvbCNmVEpGLEYtLSUmQ09MT1JHNiYlJFJHQkckIiM1Ri4kRi9GLkZqXWwtRiY2JDdcbzckRiokITNIKEh1SGE1UFAqRjU3JEYxJCEzVSdlIUdHWDEqZipGNTckRjckITMnUnBUMCNlIXlzKkY1NyQkITNsJlIzUixZWyhIRiwkITNhcSlSK0c4X3YqRjU3JEY8JCEzTTF3OzVHKmZ3KkY1NyQkITNLLzdIa1FFOkhGLCQhMyVRRmknZWRqZygqRjU3JEZBJCEzKVtbJCl6X0cnUigqRjU3JCQhMyFwZE4yVm0lPUdGLCQhM29sIzNFenZ3aipGNTckRkYkITM9Kj1OQCJRT2slKkY1NyQkITNHQi5nXiY0U28jRiwkITNPUWBNYHA6QiMqRjU3JEZLJCEzUyVHYyNmej8/KilGNTckRlAkITNsbEwmXE46TjopRjU3JEZVJCEzJypmYVQsUyJmRihGNTckRlokITNTdUBdOWk8SWlGNTckRmluJCEzJz5oXSJ5clRNXUY1NyRGXm8kITMhKTQueHozWmVQRjU3JEZjbyQhMzpzMixtcigqKlEjRjU3JEZobyQhMyFwYHNSPi46OyJGNTckRl1wJCIzKTRLYDokcCo+QSNGYXQ3JEZjcCQiM0xiIjQtcSVmKGUiRjU3JEZocCQiMzA2IjNicGklZUdGNTckRl1xJCIzdT43Iz0yInlkUkY1NyRGYnEkIjMqcGApM0YpcCR5XkY1NyRGZ3EkIjNbWVgjb1hbWzcnRjU3JEZcciQiM0JQN1I5JSopKTRyRjU3JEZhciQiMywiSHE5MiZwbXlGNTckRmZyJCIzcWlFdiNRJTNpJilGNTckRltzJCIzL2RXQFdeJFszKkY1NyRGYHMkIjN5by93QT9ZeSUqRjU3JEZlcyQiMyg0c3FMeVJlZypGNTckRmpzJCIzWUk/XkJjNypwKkY1NyRGX3QkIjMjMzs3PipbQ2goKkY1NyRGZXQkIjMpcChmRDIvQSR5KkY1NyRGW3UkIjNtdkFuIyoqcE13KkY1NyRGYHUkIjMxN20hR1krL3EqRjU3JEZldSQiMygqW1BTSTFSNScqRjU3JEZqdSQiM0BdMyQqUidSelsqRjU3JEZfdiQiM05bXiE+KFFJOSIqRjU3JEZkdiQiM2c+U29LN3R2JilGNTckRml2JCIzJyo+I0dnJlxSMHpGNTckRl53JCIzd1MhW2p5NChIckY1NyRGY3ckIjNvNV8ob25MTTgnRjU3JEZodyQiM18lKnpyO0UySF5GNTckRl14JCIzKnlDLT15anUmUkY1NyRGYngkIjNEKDRBYDgxOyNHRjU3JEZneCQiMyZveFl5I1F6OzpGNTckRlx5JCIzJ3pdMylRNnVyQ0ZhdDckRmF5JCEzSXJBSmVXQC02RjU3JEZmeSQhM2Q5Im9yaSEzQENGNTckRlt6JCEzYmFLci8jXF54JEY1NyRGYHokITNML25EWGhPRV1GNTckRmV6JCEzXiNlKWZ1W1pBaUY1NyRGanokITMpXDBUT2laUkgoRjU3JEZfW2wkITN3V3Q2cmdBXCIpRjU3JEZkW2wkITMpKjN3WCFHaTUlKilGNTckRmlbbCQhM056VFhXJVtdVypGNTckJCIzUSE0blYscDsiR0YsJCEzViJmdz1jZEtpKkY1NyRGXlxsJCEzcCMqWzk2cyo+dCpGNTckJCIzRCMpRz1KS3o0SEYsJCEzdCZlS0UhcCR6dipGNTckRmNcbCQhM1k0KG8kKWUhPm0oKkY1NyQkIjNCNDNtcURKdEhGLCQhM18lXChvSilwaHYqRjU3JEZoXGwkITNuLFE5Myl6c3MqRjU3JEZdXWwkITNHemxOIXkvJilmKkY1NyRGYl1sJCEzNSFlSS5gNVBQKkY1LUZlXWw2JkZnXWxGal1sRmhdbEZqXWwtJStBWEVTTEFCRUxTRzYkUSJ4NiJRIUZiW20tJSVWSUVXRzYkOyQhK2FFZlRKISIqJCIrYUVmVEpGalttOyQhMilRJz1RISoqKipSNSEjOyQiMnVkXVo0JioqUjVGYFxtLSUqVEhJQ0tORVNTRzYjIiIj Now we will try several other functions, just to show off. First, sin(x), which doesn't do as well (we'd need to add a fifth degree polynomial to our basis for a good fit.) ShowApprox(sin(x),OrthoBasis); NiMsJiooLCYhIzoiIiIqJEkjUGlHSSpwcm90ZWN0ZWRHRioiIiNGJ0YnRikhIidJInhHNiIiIiQjIiNORisqJiwmKiRGKSEiI0YvKiZGJUYnRikhIiUjISNARitGJ0YtRidGJw== NiUsJiokSSJ4RzYiIiIkJCErOyhwKFEkKiEjNkYlJCIrckskKXAmKSEjNVE1fm5vcm1+b2Z+ZGlmZmVyZW5jZT1GJiQiKywyJTNtIkYt LSUlUExPVEc2Jy0lJ0NVUlZFU0c2JDdlbjckJCEzKioqKio0dGsjZlRKISM8JCEzPTVLVF9LenppISNFNyQkITN3IipwcikzUFkrJEYsJCEzVVdrYlUjeV9PIiEjPTckJCEzQDNFWyp5c2EpR0YsJCEzNGBrSTshKkdMREY1NyQkITNmWCYpKT4yZzl2I0YsJCEzbUpeMHguNi5RRjU3JCQhMyk0NTc3LmZsaCNGLCQhM0c/KSlRJyk0VTddRjU3JCQhM3YiUU06OipII1sjRiwkITNgWnE7NmRlRGhGNTckJCEzeVxBaGNJI3lOI0YsJCEzOG5ES09GYWZxRjU3JCQhMyZlPWQtRk4qR0FGLCQhMy11TyQpeSE+OCJ6RjU3JCQhM3clKkdCUCRSYzQjRiwkITMqZkMlcGhoLWEnKUY1NyQkITNtQ2omZikzeGk+RiwkITMyUjByRmJjVCMqRjU3JCQhMydvcjRBTio0RT1GLCQhM3RmZyNcOW9lbipGNTckJCEzang8YnIiNGZ3IkYsJCEzaE16bHh6RDUpKkY1NyQkITNRUVEqMyoqPWRxIkYsJCEzQDZJL3B0NjQqKkY1NyQkITNmOSZ5TTVmemoiRiwkITN2TSdmU0dhdSgqKkY1NyQkITNkIT5qZ0AqPnE6RiwkITJTI1JIPCMpKioqKioqRiw3JCQhM2IuKlIrNWhAXSJGLCQhMyhccyNHdF9YdyoqRjU3JCQhM2A7bSwlKUg3TTlGLCQhMykzIVsmR0dabiEqKkY1NyQkITM5ZD0yX2NibzhGLCQhM18sXSVmayo9J3oqRjU3JCQhM3YoNEYsSykpSEkiRiwkITNJQU9cU0RgVicqRjU3JCQhM1MqNCFSI1swUj0iRiwkITMqKXBtO0NsJzNFKkY1NyQkITNwWUBRcVdJVTVGLCQhMyQqZV9tZV1vTicpRjU3JCQhM1dwM3ckUilcQiMqRjUkITNoa11KQStCcXpGNTckJCEzRUpjOG8zOUd5RjUkITNIYDQ1Vj94X3FGNTckJCEzJz1bQlAtOElmJ0Y1JCEzclo5bG10a0RoRjU3JCQhM0NCMzxAPyl5QiZGNSQhMzVXJGU/UlM7KyZGNTckJCEzYCopUjxZb1paUkY1JCEzJFwyL0whSHZYUUY1NyQkITM4VHFYPVcyLEVGNSQhM18jUjU4IlElPWQjRjU3JCQhM3FUSlkpb2NZTyJGNSQhMy5xJmU4KlxVZzhGNTckJCEzOCx6N1ZLSixKISM/JCEzcT0zbHIjMzg1JEZldDckJCIzOlNGZjN4RWE4RjUkIjM2QVJdbD04XThGNTckJCIzJ1IyKClIdmUsYyNGNSQiMytGIWZhLiRHS0RGNTckJCIzSWN3UjxUYmlRRjUkIjM0ajJCVE5BblBGNTckJCIzKjQqPUpvZjAzX0Y1JCIzSSs4OiopM3p2XEY1NyQkIjNwXTFYSXFPQ2xGNSQiM2hHa3pwVUNyZ0Y1NyQkIjMlSG9CbG1senooRjUkIjNFI1J5NyJ5TUpxRjU3JCQiM151J2YjNCl6P0AqRjUkIjNbJHlUZGVHTCd6RjU3JCQiMzlIdjxFQ0ZbNUYsJCIzZkxhKCkzTWlsJylGNTckJCIzTmxqJUclMyVSPSJGLCQiMzchR0pzJCoqKjRFKkY1NyQkIjNcLSlSRWZ3b0kiRiwkIjNEOjhXJHBbUGwqRjU3JCQiM2BvX3A6czJ1OEYsJCIzZyZcbyVRNjgyKSpGNTckJCIzY00ydlF5RlQ5RiwkIjNaZkRgIio+QzsqKkY1NyQkIjM9Nip5elEzWF0iRiwkIjNpXU9uKmVQIXkqKkY1NyQkIjN6KDMycyQqUXhjIkYsJCIzWVw0KT1FYCoqKioqRjU3JCQiMyZSOTAtM0xRaiJGLCQiMzN0K3o7JVEsKSoqRjU3JCQiMzcrSz9CcyMqKnAiRiwkIjNuRyllciU9dTsqKkY1NyQkIjM5WUsqKUdqYWs8RiwkIjMleid6MjM2KkciKSpGNTckJCIzPCNIJGVNYTtIPUYsJCIzJClRPjdNJ3ohbycqRjU3JCQiMyU+PkNaJ2VZaz5GLCQiMy0nZi0rP3hdQipGNTckJCIza1l1eWZpeCU0I0YsJCIzNHlVJD1pWiRlJylGNTckJCIzVDRxKSlmdS5HQUYsJCIzSmpzJXBJMm8iekY1NyQkIjN6P3cnKio9Jj5nQkYsJCIzeWNDbXo/c1VxRjU3JCQiMyEqPjNhPVdqIlsjRiwkIjNCLCY9WGRRMzgnRjU3JCQiM0E/YyopeXUiM2kjRiwkIjMkNCIpcCYqKnBfdlxGNTckJCIzLGktSG5XSVhGRiwkIjNfXkBHMidvKmZRRjU3JCQiM3Y9UldoTi55R0YsJCIzSDRFd3RRPTBFRjU3JCQiM3NzKCpSU0EyMElGLCQiMyRlTHdKTW40TyJGNTckJCIzISkqKipcL2wjZlRKRiwkIjNwYXdwT016UkpGLy0lJkNPTE9SRzYmJSRSR0JHJCIjNSEiIiQiIiFGal1sRltebC1GJjYkN2huNyRGKiQiM14kcGFtVEFKLiNGNTckJCEzNCdcOCFvWzZ0SUYsJCIzc049MyRIKWZ1dyEjPjckRjEkITNpRC02VzdkdVRGaF5sNyQkITMpKip6KjRSXDBYSEYsJCEzaHciXFxZYVRRIkY1NyRGNyQhMyVRbiV5dCVwQUgjRjU3JEY8JCEzaUE8XGUqem83JUY1NyRGQSQhMyM0JXlqKVtSU3AmRjU3JEZGJCEzXjlyJ0c2PSkpKXBGNTckRkskITNVLzxeNiwoXCd6RjU3JEZQJCEzbyczZT9QUix3KUY1NyRGVSQhM001KlxYNjVXTypGNTckRlokITNFcWQjKVtXNGYoKkY1NyQkITMpMy0kMz5eViUqPUYsJCEzW2Zfc2wpW2MpKSpGNTckRmluJCEzKD1TbFlfTEUnKipGNTckRl5vJCEzZ2lRQyUqXCEzKioqRjU3JEZjbyQhMzYsJEdOIylHSikqKkY1NyRGaG8kITNaJHo/PyhwNkwqKkY1NyRGXXAkITNKay8qZWlsNCUpKkY1NyRGZ3AkITMrOj1ZP2psTiYqRjU3JEZhcSQhM3IsJ0cxVyRbKyIqRjU3JEZmcSQhM2kjUXIhPnI+J2YpRjU3JEZbciQhM3ZTUFljQipbKHlGNTckRmByJCEzRHlsb04iKmZyckY1NyRGZXIkITNfMCkpemEhKmZnaUY1NyRGanIkITNvV0lpRyJvQ1EmRjU3JEZfcyQhMyYpUi5IU2lkYVZGNTckRmRzJCEzKUcjcHJYdFpETEY1NyRGaXMkITN6JSkqXENCVkVAI0Y1NyRGXnQkITN0Sm9hIXA5cjsiRjU3JEZjdCQhM2cxNil6JjR4ZEVGZXQ3JEZpdCQiMyo9RDwuSWwjZTZGNTckRl51JCIzQCd5P3JmVSR5QEY1NyRGY3UkIjNNU0BZTiNHakQkRjU3JEZodSQiM0dWeidwSyZISlZGNTckRl12JCIzM3YnNDVvNz5MJkY1NyRGYnYkIjNFRkUjZi8rKlJpRjU3JEZndiQiMyopKioqKT1xPGBrckY1NyRGXHckIjNTTDxdVnp3MnpGNTckRmF3JCIzJVs0O0ZIaGpmKUY1NyRGZnckIjM1I29XUkdiXzYqRjU3JEZgeCQiM11GRVZUKlFiYipGNTckRmp4JCIzKG9EYjE+XG8kKSpGNTckRl95JCIzIylwWyhcJEhwRyoqRjU3JEZkeSQiM0xFJD13YVcwKSoqRjU3JEZpeSQiM09ReHVFSy0iKioqRjU3JEZeeiQiMyRHTU9iYzstJyoqRjU3JCQiMy9VUGxcYyJvKj1GLCQiMzRrMipbJFIzIykpKkY1NyRGY3okIjNgWnROITQ1YHYqRjU3JEZoeiQiM0Y/XzRnJkd3TypGNTckRl1bbCQiM14zJyk9OG8kXHcpRjU3JEZiW2wkIjMpKmVjKSopRzokW3pGNTckRmdbbCQiMzkqR0k3diRmJSpwRjU3JEZcXGwkIjM9aSFlMEM+KFtjRjU3JEZhXGwkIjM4RHFdZlxSL1VGNTckJCIzUSE0blYscDsiR0YsJCIzWD1wI1JLIXBQTEY1NyRGZlxsJCIzeSUpKlxbISkpZixDRjU3JCQiM3VYPSM0IUhiVEhGLCQiM15fXjIjR1IiUjlGNTckRltdbCQiM2NpXEBDYnksVEZoXmw3JCQiMy8nKVtVWENMdElGLCQhM0NnUzQuWF84eEZoXmw3JEZgXWwkITNzPCZwbFpBSi4jRjUtRmVdbDYmRmddbEZbXmxGaF1sRltebC0lK0FYRVNMQUJFTFNHNiRRIng2IlEhRmlqbC0lJVZJRVdHNiQ7JCErYUVmVEohIiokIithRWZUSkZhW207JCExVXhvQygqKipSNSEjOiQiMltXWSJIXyoqUjUhIzstJSpUSElDS05FU1NHNiMiIiM= ShowApprox(cos(1+x)-x,OrthoBasis); NiMsLioqLUkkY29zRzYkSSpwcm90ZWN0ZWRHRihJKF9zeXNsaWJHNiI2IyIiIkYsSSNQaUdGKCEiKSwmKiRGLSIiI0YxISNARixGLEkieEdGKiIiJSMhJXY6IiIpKiotSSRzaW5HRidGK0YsRi0hIicsJiEjOkYsRjBGLEYsRjMiIiQjISNORjEqJiwmKiZGLSEiJUYlRiwjISNYRjEqKEYlRixGLUY7Ri9GLCMiJHYnRjRGLEYzRjFGLComLCYqJkYtISIjLCZGOUY+RjBGLEYsISIiKihGOUYsRi1GREY8RiwjIiNARjFGLEYzRixGLComRi1GTUYlRiwjIiM6RjEqKEYlRixGLUZERi9GLCMhJE4iRjc= NiUsLCokSSJ4RzYiIiIlJCIrPDZVODkhIzYqJEYlIiIkJCIrXFBJZXlGKiokRiUiIiMkIStoOXNWQyEjNUYlJCErL203QDwhIiokIitUKyNmRyZGMyIiIlE1fm5vcm1+b2Z+ZGlmZmVyZW5jZT1GJiQiK1g/XD45RjM= LSUlUExPVEc2Jy0lJ0NVUlZFU0c2JDdTNyQkITMqKioqKjR0ayNmVEohIzwkIjMiUjdFbk0hSCxFRiw3JCQhM3ciKnByKTNQWSskRiwkIjMiXD1fJlsjeVVlI0YsNyQkITNAM0VbKnlzYSlHRiwkIjMsVClwRCZSJ2ZkI0YsNyQkITNmWCYpKT4yZzl2I0YsJCIzbDBaXzt2eHJERiw3JCQhMyk0NTc3LmZsaCNGLCQiMzxpSy0pSDczZCNGLDckJCEzdiJRTTo6KkgjWyNGLCQiMz1FWlJkM29xREYsNyQkITN5XEFoY0kjeU4jRiwkIjMwNEshNCkqKj1wREYsNyQkITMmZT1kLUZOKkdBRiwkIjNUJ2UoZi5qPGtERiw3JCQhM3clKkdCUCRSYzQjRiwkIjNaPC9YKVE8SmIjRiw3JCQhM21DaiZmKTN4aT5GLCQiM0E+VSRveT1TYCNGLDckJCEzJ29yNEFOKjRFPUYsJCIzQTleQSJRW1FdI0YsNyQkITNRUVEqMyoqPWRxIkYsJCIzYm9wVDpXJ29ZI0YsNyQkITNkIT5qZ0AqPnE6RiwkIjNWcnY5ckUqPlQjRiw3JCQhM2A7bSwlKUg3TTlGLCQiM1xCTTtoQU9UQkYsNyQkITN2KDRGLEspKUhJIkYsJCIzPHlSKWZWUHVEI0YsNyQkITNTKjQhUiNbMFI9IkYsJCIzSFFwJHBaVXE7I0YsNyQkITNwWUBRcVdJVTVGLCQiM1MzPHNwKDQ5LyNGLDckJCEzV3AzdyRSKVxCIyohIz0kIjNrMCRmLGVPJD4+Riw3JCQhM0VKYzhvMzlHeUZccSQiM2EhZkwkKXlAJGY8Riw3JCQhMyc9W0JQLThJZidGXHEkIjM/NUw7SkYjPWciRiw3JCQhM0NCMzxAPyl5QiZGXHEkIjNBQjY8cGdfNzlGLDckJCEzYCopUjxZb1paUkZccSQiM09DKlsnW2k1PDdGLDckJCEzOFRxWD1XMixFRlxxJCIzZ2NjXUlXWycpKipGXHE3JCQhM3FUSlkpb2NZTyJGXHEkIjM9ZDJeNWc/aXlGXHE3JCQhMzgsejdWS0osSiEjPyQiMyQqUUMnbyJvNWdhRlxxNyQkIjM6U0ZmM3hFYThGXHEkIjMnenIuJDR1PWpHRlxxNyQkIjMnUjIoKUh2ZSxjI0ZccSQiMyQ0b0xURCU9ZmAhIz43JCQiM0ljd1I8VGJpUUZccSQhM2JnXDkpWyx3LSNGXHE3JCQiMyo0Kj1Kb2YwM19GXHEkITMpMz5RLCdvTjNaRlxxNyQkIjNwXTFYSXFPQ2xGXHEkITNPLF0yJz1rKVJ0RlxxNyQkIjMlSG9CbG1senooRlxxJCEzUXVxYWhjeXMpKkZccTckJCIzXnUnZiM0KXo/QCpGXHEkITNYdlooZlgjXGs3Riw3JCQiMzlIdjxFQ0ZbNUYsJCEzTj4memlXNnldIkYsNyQkIjNObGolRyUzJVI9IkYsJCEzJFEoKil6aUdRZjxGLDckJCIzXC0pUkVmd29JIkYsJCEzIVxcbzd0ayN5PkYsNyQkIjNjTTJ2UXlGVDlGLCQhM3o1KlwpeXIiZj8jRiw3JCQiM3ooMzJzJCpReGMiRiwkITNZRnMxeFNiMkNGLDckJCIzNytLP0JzIyoqcCJGLCQhM18nUihbZiNvUmcjRiw3JCQiMzwjSCRlTWE7SD1GLCQhM0k4PHpiZHYheSNGLDckJCIzJT4+Q1onZVlrPkYsJCEzPUd6ZSYpKT4pW0hGLDckJCIza1l1eWZpeCU0I0YsJCEzXTkrZHMwbyQ0JEYsNyQkIjNUNHEpKWZ1LkdBRiwkITMiKj41d0xNSUNLRiw3JCQiM3o/dycqKj0mPmdCRiwkITM2ZThDYW1ST0xGLDckJCIzISo+M2E9V2oiWyNGLCQhM0JxWzxiXlBDTUYsNyQkIjNBP2MqKXl1IjNpI0YsJCEzZ3pXbiU0cSIzTkYsNyQkIjMsaS1IbldJWEZGLCQhMyZSZiV6YCZRJm9ORiw3JCQiM3Y9UldoTi55R0YsJCEzIjMyInkhSCgqKT1PRiw3JCQiM3NzKCpSU0EyMElGLCQhMyhvJj0/JnlvW2wkRiw3JCQiMyEpKioqXC9sI2ZUSkYsJCEzIls9ZyplXCo9byRGLC0lJkNPTE9SRzYmJSRSR0JHJCIjNSEiIiQiIiFGYltsRmNbbC1GJjYkN1M3JEYqJCIzXzA4K2B2L2tDRiw3JEYwJCIzWCVwKSlHSXBUXiNGLDckRjUkIjMnZiIzV1JpN19ERiw3JEY6JCIzUSQpcDkoR2p0ZSNGLDckRj8kIjNHa10qKVJddDhFRiw3JEZEJCIzJCp5MyNHZ14pSEVGLDckRkkkIjMwbWV2cEckXGojRiw3JEZOJCIzKFE7UEtzZCVIRUYsNyRGUyQiMzRNMk5jbGg2RUYsNyRGWCQiMykqeXJMPiQqKjNlI0YsNyRGZ24kIjMqNCg0O1JdSE5ERiw3JEZcbyQiMyJvLWMoSFMsJFsjRiw3JEZhbyQiM0pqJilRP01INUNGLDckRmZvJCIzayY9Q0shKilIQUJGLDckRltwJCIzUngwKlwqNEBCQUYsNyRGYHAkIjMyMzomPU4hNEBARiw3JEZlcCQiM1Y3JlspR0J1JSk+Riw3JEZqcCQiMz50SGFzZ3VjPUYsNyRGYHEkIjM3akZgJ0h1UHAiRiw3JEZlcSQiMy0nKipHbnNmc2AiRiw3JEZqcSQiM1QmeWE4YkNHTiJGLDckRl9yJCIzTXUiMz9wS2E7IkYsNyRGZHIkIjMoM0x3X1ElPSVlKkZccTckRmlyJCIzO2MnKnl2KDRzZShGXHE3JEZecyQiM18nKT51W1NGUmBGXHE3JEZkcyQiMyxyI3pqUU1BIkhGXHE3JEZpcyQiM15RMmQ+I1s9TChGXXQ3JEZfdCQhMyRvWndoVyE9eTtGXHE3JEZkdCQhM0FYWzFgekE+VUZccTckRml0JCEzLmsiNDpzOChSbkZccTckRl51JCEzZ1pYYGw9WSc+KkZccTckRmN1JCEzVzQuXjA8cCM+IkYsNyRGaHUkITNaOVshejxqbFYiRiw3JEZddiQhM19hN1hBcVokcCJGLDckRmJ2JCEzJD4iKXlsJW9XQD5GLDckRmd2JCEzVm83I1sqZlJqQEYsNyRGXHckITNgZGRVIVFFQFEjRiw3JEZhdyQhM2t6QVsqKVJKKmYjRiw3JEZmdyQhM29LJ29MSSsiKXojRiw3JEZbeCQhMyhvVFFjMz4kKilIRiw3JEZgeCQhM082c2dXMmhhSkYsNyRGZXgkITNKMiZ6TEAkeSxMRiw3JEZqeCQhM1tLKFI4YDtKVSRGLDckRl95JCEzb2NfW15QXjVORiw3JEZkeSQhMywhUUQpcCVHI3pORiw3JEZpeSQhM2BCdk5cLVQ0T0YsNyRGXnokITNjbWAoZj1WZmckRiw3JEZjeiQhM2hoOllZWjpsTkYsNyRGaHokITNUMCUpZk1dKHBaJEYsLUZdW2w2JkZfW2xGY1tsRmBbbEZjW2wtJStBWEVTTEFCRUxTRzYkUSJ4NiJRIUZhZWwtJSVWSUVXRzYkOyQhK2FFZlRKISIqJCIrYUVmVEpGaWVsOyQhMTFYYDo6QjNRISM6JCIybz1JalVwN3cjISM7LSUqVEhJQ0tORVNTRzYjIiIj ShowApprox(exp(x)*cos(x),OrthoBasis); NiMsMCoqLDIqJEkjUGlHSSpwcm90ZWN0ZWRHRigiIiMiIyEqRiciJDUjKiRGJyIiJSEiJSIkMCIiIiIqJi1JJGV4cEc2JEYoSShfc3lzbGliRzYiNiMsJEYnRilGMEYnRikhIyEqKiZGMkYwRidGMEYrKiZGMkYwRidGLUYtRjIhJDAiRjAtRjM2IywkRichIiJGMEYnISIqSSJ4R0Y2Ri0jISQ6JCIkRyIqKiwuISM6RjBGJ0ZIKiRGJyIiJEYpRjIiIzpGOkZIKiZGMkYwRidGSkYpRjBGPUYwRichIihGQkZKIyEjTiIjOyomLCYqKiwmISIkRjBGJkYpRjAsJkZARjBGMkYwRjBGPUYwRichIiYjRkhGUCooRiVGMEY9RjBGJ0ZNIyIkTiIiI2tGMEZCRilGMComLCYqKCwmRjBGMEYyRjBGMEY9RjBGJyEiIyNGVUYtKihGR0YwRj1GMEYnRlcjIiNARlBGMEZCRjBGMCooRlZGMEY9RjBGJ0ZAI0ZARi0qKkZURjBGVkYwRj1GMEYnRlUjIiImRlAqKEYlRjBGPUYwRidGVyMhI0ZGRQ== NiUsLCokSSJ4RzYiIiIlJCErMiJRcDYiISM1KiRGJSIiJCQhK28rSzRdRioqJEYlIiIjJCErcStyJlIjRipGJSQiKy9bai83ISIqJCIrI2U7aDciRjUiIiJRNX5ub3Jtfm9mfmRpZmZlcmVuY2U9RiYkIitUJm9abCRGKg== LSUlUExPVEc2Jy0lJ0NVUlZFU0c2JDdZNyQkITMqKioqKjR0ayNmVEohIzwkITNEc1leYD1SQFYhIz43JCQhM3ciKnByKTNQWSskRiwkITNmcXZRQClwIzRcRi83JCQhM0AzRVsqeXNhKUdGLCQhM187cjwkeUYyUyZGLzckJCEzZlgmKSk+Mmc5diNGLCQhM25rZDdcc3ouZkYvNyQkITMpNDU3Ny5mbGgjRiwkITNoXFIyJG8oUkBqRi83JCQhM3YiUU06OipII1sjRiwkITM9d2loW1YyL21GLzckJCEzeVxBaGNJI3lOI0YsJCEzJXpCXlxAYz5xJ0YvNyQkITMmZT1kLUZOKkdBRiwkITMhcG0sVT1OUmUnRi83JCQhM3clKkdCUCRSYzQjRiwkITNkUykpeUh0I0c7J0YvNyQkITNtQ2omZikzeGk+RiwkITNYOy95YykpNG1gRi83JCQhMydvcjRBTio0RT1GLCQhM25FKVFweSUqbzElRi83JCQhM1FRUSozKio9ZHEiRiwkITMubXFAbClITFcjRi83JCQhM2QhPmpnQCo+cTpGLCQiM0Z0b0skZTc/QyIhI0A3JCQhM2A7bSwlKUg3TTlGLCQiMyJRJSk+X1dFckMkRi83JCQhM3YoNEYsSykpSEkiRiwkIjNHLi5JJTNwLD4oRi83JCQhM1MqNCFSI1swUj0iRiwkIjMxUmEmNCpwKFs6IiEjPTckJCEzcFlAUXFXSVU1RiwkIjNfJWZzRWVGInk8RmZwNyQkITNXcDN3JFIpXEIjKkZmcCQiM1Z0KFslKVstN1MjRmZwNyQkITNFSmM4bzM5R3lGZnAkIjN1X0U0OkhsU0tGZnA3JCQhMyc9W0JQLThJZidGZnAkIjNiY3VlJSlROikzJUZmcDckJCEzQ0IzPEA/KXlCJkZmcCQiMzlNXGMyMG5HXkZmcDckJCEzYCopUjxZb1paUkZmcCQiMyZHNXNreW0tQSdGZnA3JCQhMzhUcVg9VzIsRUZmcCQiM3UhKj04NzxOXXVGZnA3JCQhM3FUSlkpb2NZTyJGZnAkIjMib0g3TExfS2spRmZwNyQkITM4LHo3VktKLEohIz8kIjNaRFUlbydwKSpvKipGZnA3JCQiMzpTRmYzeEVhOEZmcCQiMycpZi5cazdhTTZGLDckJCIzJ1IyKClIdmUsYyNGZnAkIjNPTzpDJSkpcCdcN0YsNyQkIjNJY3dSPFRiaVFGZnAkIjMtWyJvIipSX0lPIkYsNyQkIjMqNCo9Sm9mMDNfRmZwJCIzQGVdO1t6Pmc5Riw3JCQiM3BdMVhJcU9DbEZmcCQiM3U+SiNvOjplXyJGLDckJCIzJUhvQmxtbHp6KEZmcCQiMyplJ2UmXHJNM2IiRiw3JCQiM151J2YjNCl6P0AqRmZwJCIzdDYmUTVfKyc+OkYsNyQkIjM5SHY8RUNGWzVGLCQiMyxEJXl1JkhxQjlGLDckJCIzTmxqJUclMyVSPSJGLCQiM0BSV0swQ2xLN0YsNyQkIjNcLSlSRWZ3b0kiRiwkIjM8R2o5P3QsUScqRmZwNyQkIjNjTTJ2UXlGVDlGLCQiMyNcJnl5JCkzR2VhRmZwNyQkIjN6KDMycyQqUXhjIkYsJCIzRyEzdjpGZWlZIkYvNyQkIjM3K0s/QnMjKipwIkYsJCEzWE8mRyN5aVRbcUZmcDckJCIzPCNIJGVNYTtIPUYsJCEzJ2YoPihwMmI5ZiJGLDckJCIzJT4+Q1onZVlrPkYsJCEzRDpqSXFyTU5GRiw3JCQiM2tZdXlmaXglNCNGLCQhM0xMI1IxJUdca1NGLDckJCIzVDRxKSlmdS5HQUYsJCEzeiJITWJkdS9uJkYsNyQkIjN6P3cnKio9Jj5nQkYsJCEzPF9ePylILS5fKEYsNyQkIjNNP1VEL1siNFUjRiwkITNHQ1FQbGpTZCUpRiw3JCQiMyEqPjNhPVdqIlsjRiwkITNWL1AnekhWI1wlKkYsNyQkIjMxPyM9KFtmQV5ERiwkITN2bSU9aDZRXzEiISM7NyQkIjNBP2MqKXl1IjNpI0YsJCEzLUd2UUt1WCM+IkZeejckJCIzNlRINHQ0MSRvI0YsJCEzKjRTbWtvZj1KIkZeejckJCIzLGktSG5XSVhGRiwkITNYSls/OElGTzlGXno3JCQiM1EhNG5WLHA7IkdGLCQhMzYnSCNHJ2VNU2QiRl56NyQkIjN2PVJXaE4ueUdGLCQhMzhjPHpvUWA7PEZeejckJCIzdVg9IzQhSGJUSEYsJCEzSmswPVJddWM9Rl56NyQkIjNzcygqUlNBMjBJRiwkITNFQWtURVIpKioqPkZeejckJCIzLycpW1VYQ0x0SUYsJCEzemA7eClbWmo6I0ZeejckJCIzISkqKipcL2wjZlRKRiwkITNgRkEsYyNwU0ojRl56LSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiIUZiXWxGY11sLUYmNiQ3VzckRiokITNtKikzVjRbI3lxJEZmcDckRjEkITNTJT1NKSo9LlpyIkZmcDckRjYkITNQc0gpXCRvcHJfRi83JEY7JCIzcWozXnY3OSYzJEYvNyRGQCQiM0swODAjKTROPXNGLzckRkUkIjNaUHFAT3FiJzMpRi83JEZKJCIzKiozQXhhR2Y0b0YvNyRGTyQiM1wlPTRycE08NiVGLzckRlQkIjNBOSk0aCs6dGImRmJzNyRGWSQhMyZ6UVBgRHVZNiRGLzckRmhuJCEzV1JzaEtJdT5rRi83JEZdbyQhM1VyL0w5TkE8JilGLzckRmJvJCEzOzQwJzNTTlZkKkYvNyRGaG8kITMjW0tPKClRXVMiKilGLzckRl1wJCEzTyVvY3lRW1dTJ0YvNyRGYnAkITNCciQpXHpQOy5DRi83JEZocCQiM11XSUopXD9kYyVGLzckRl1xJCIzb3RMJGYtM1dCIkZmcDckRmJxJCIzUzU3ITQhKSpmWUJGZnA3JEZncSQiMyZbby4kKlIhPS1ORmZwNyRGXHIkIjMmKUd1aHdCJSpIXEZmcDckRmFyJCIzZzI/NndmZjhrRmZwNyRGZnIkIjNEJCpwYTV3eFshKUZmcDckRltzJCIzd21qbTs2KVxlKkZmcDckRmBzJCIzOUxXQVgkeUI3IkYsNyRGZnMkIjNaNiwmKnkyZSRHIkYsNyRGW3QkIjNgKWYkRyg0TSo0OUYsNyRGYHQkIjNfQjhxOHhKQzpGLDckRmV0JCIzSChvIT4oKmZgNDtGLDckRmp0JCIzIWZGSil5YHNdO0YsNyRGX3UkIjNnVT9EU1woNGsiRiw3JEZkdSQiM2NVKCpIQlNbZzpGLDckRml1JCIzPHZcQnNqdDg5Riw3JEZediQiMyhlbyUzP1h1bDZGLDckRmN2JCIzb1glUlZzKVF0JSlGZnA3JEZodiQiM1BJJVtlN0olSFFGZnA3JEZddyQhM2shRyZRSi1nIXoiRmZwNyRGYnckITNWLyI0NSQpUSg9IipGZnA3JEZndyQhM2xdZEc3IjQiKXkiRiw3JEZceCQhMys2XiY+Sj9JKkdGLDckRmF4JCEzWSZwTSdSKCkqcDolRiw3JEZmeCQhMy5pTCJSZiI0c2NGLDckRlt5JCEzSyk+dEIpZj48dUYsNyRGZXkkITNqJCopPkJxJyk9RCpGLDckRmB6JCEzM0N5STZQJVw7IkZeejckRmV6JCEzdipcZGsnUSwkRyJGXno3JEZqeiQhM1wkUiYpNFZLIjM5Rl56NyRGX1tsJCEzRVYib0soemRcOkZeejckRmRbbCQhM1NSYW0tT2kqcCJGXno3JEZpW2wkITNUdkMqZj15OiY9Rl56NyRGXlxsJCEzXCJHdE4jRyg+LCNGXno3JEZjXGwkITNLN1tMKnBvUz4jRl56NyRGaFxsJCEzQEpGZl42ZidRI0Zeei1GXV1sNiZGX11sRmNdbEZgXWxGY11sLSUrQVhFU0xBQkVMU0c2JFEieDYiUSFGXWhsLSUlVklFV0c2JDskISthRWZUSiEiKiQiK2FFZlRKRmVobDskITIlW0FRRFdpUEMhIzokIjJzV0VuNmU1OyNGXnotJSpUSElDS05FU1NHNiMiIiM=