In this worksheet we examine some of the potential 'problems' with fitting curves to data from section 4.2
restart:
with(plots):
with(linalg):
We pick four common trigonometric functions as our set of functions
F:=A*sin(x)+B*cos(x)+C*sin(2*x)+E*cos(2*x);
LCoqJkkiQUc2IiIiIi1JJHNpbkc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkdGJTYjSSJ4R0YlRiZGJiomSSJCR0YlRiYtSSRjb3NHRilGLEYmRiYqJkkiQ0dGJUYmLUYoNiMsJEYtIiIjRiZGJiomSSJFR0YlRiYtRjFGNUYmRiY=
We pick four points, only one of which we know an exact value of sine and cosine for
pts:=[[0,Pi],[4,2],[7,9],[-1,3]];
NyY3JCIiIUkjUGlHJSpwcm90ZWN0ZWRHNyQiIiUiIiM3JCIiKCIiKjckISIiIiIk
simplify(subs({x=pts[1][1],y=pts[1][2]},F=y),trig);
LywmSSJCRzYiIiIiSSJFR0YlRiZJI1BpRyUqcHJvdGVjdGVkRw==
Sys:=seq(simplify(subs({x=pts[j][1],y=pts[j][2]},F=y),trig),j=1..4);
NiYvLCZJIkJHNiIiIiJJIkVHRiZGJ0kjUGlHJSpwcm90ZWN0ZWRHLywqKiZJIkFHRiZGJy1JJHNpbkc2JEYqSShfc3lzbGliR0YmNiMiIiVGJ0YnKiZGJUYnLUkkY29zR0YxRjNGJ0YnKiZJIkNHRiZGJy1GMDYjIiIpRidGJyomRihGJy1GN0Y7RidGJyIiIy8sKiomRi5GJy1GMDYjIiIoRidGJyomRiVGJy1GN0ZERidGJyomRjlGJy1GMDYjIiM5RidGJyomRihGJy1GN0ZKRidGJyIiKi8sKiomRi5GJy1GMDYjISIiRidGJyomRiVGJy1GN0ZTRidGJyomRjlGJy1GMDYjISIjRidGJyomRihGJy1GN0ZZRidGJyIiJA==
Solving the system is straighforward....
Soln:=evalf(solve({Sys},{A,B,C,E}));
PCYvSSJBRzYiJCErM0N3Oz0hIiovSSJCR0YlJCIrNCJwJTR4RigvSSJDR0YlJCIrWGxZYF1GKC9JIkVHRiUkIStgayh5YyVGKA==
FSoln:=subs(Soln,F);
LCotSSRzaW5HNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHNiI2I0kieEdGKCQhKzNDdzs9ISIqLUkkY29zR0YlRikkIis0InAlNHhGLS1GJDYjLCRGKiIiIyQiK1hsWWBdRi0tRi9GMyQhK2BrKHljJUYt
plot1:=pointplot(pts,symbol=circle,symbolsize=30):
plot2:=plot(FSoln,x=-2..8):
Plot the data and the curve, clearly it works.
display({plot1,plot2});
NictJSdQT0lOVFNHNic3JCQiIiEhIiIkIithRWZUSiEiKjckJCIjUyEiIiQiIz8hIiI3JCQiI3EhIiIkIiMhKiEiIjckJCEjNSEiIiQiI0khIiItJSdTWU1CT0xHNiQlJ0NJUkNMRUciI0ktJSdDVVJWRVNHNiQ3W3M3JCQhIz8hIiIkIjBYKD0oZlhSRCYhIzk3JCQhMm5tbSJIVSwiKj0hIzskIjBKZSM0TS4hKWYhIzk3JCQhMk5MTCRlJUc/eSIhIzskIjFHP2VpMUEpUichIzo3JCQhMHY9VVpAJGU8ISM5JCIxPDpNMEw+XGshIzo3JCQhMm47LyxcOVl0IiEjOyQiMSMqRyN5TmNpWychIzo3JCQhMkxlKilmXTI0ciIhIzskIjEyTig0YVUnNGwhIzo3JCQhLnY9XytzbyIhIzckIjF1Jip5NFtoPmwhIzo3JCQhMmxUZ3hgJFxqOyEjOyQiMS88J1xlbGteJyEjOjckJCEySiRla2BseVI7ISM7JCIxX2duPWZeK2whIzo3JCQhMilcN2BwJnpnaCIhIzskIjJsWlhsZjlAWichIzs3JCQhMm1tO2Flc0JmIiEjOyQiMXAlKXl4cmpKayEjOjckJCEyKSoqXDdgJ0dkWyIhIzskIjFLc1ZLQF82aCEjOjckJCEyS0xMM3MlM3o4ISM7JCIxUjZjRSZbIilmJiEjOjckJCExTEwkZS8kUWs2ISM6JCIyMCQpNCMqUVttPSUhIzs3JCQhMW5tOy8iPXFdKiEjOyQiMic9N3ckRy87bSMhIzs3JCQhMUxMJDNfPmZfKCEjOyQiMVc/S3dELUY6ISM6NyQkITFubTsvSkUrbCEjOyQiMXclPW1bR2Q5IiEjOjckJCExLCtdKG8xWVomISM7JCIxYm9GIlwielAlKiEjOzckJCExb207YUdVNF8hIzskIjElMyx3aUs5QiohIzs3JCQhMlBMTDMtUlUlXCEjPCQiMUBbZ0c5QGYiKiEjOzckJCEyKCoqKipcKD1iIXpZISM8JCIxKT5vIkgjUkRBKiEjOzckJCEyam1tVE5yUVQlISM8JCIwUDMnUkZMQSUqISM6NyQkITEqKioqXChvLk4pUSEjOyQiMnNbOi0/XEstIiEjOzckJCEyO0xMMy1PSk4kISM8JCIyeSJHJ1stLyhlNiEjOzckJCEyYW1tInpDIWVIIyEjPCQiMj1OIj1XUD4nZSIhIzs3JCQhMiIqKioqXFAqbyVRNyEjPCQiMVQjUXU7RnU/IyEjOjckJCIxb0xMTDNFbiQqISM8JCIxOlIqZiwvIWZSISM6NyQkIjFwbW1UIVJFJkchIzskIjEpR1QqPUdecmQhIzo3JCQiMS4rK11LXTRdISM7JCIxJXkxR0Bta28oISM6NyQkIjEtKysrTk8jNCchIzskIjElXFhyXGomXCUpISM6NyQkIjEtKytdUEF2ciEjOyQiMC9YJlJLI0crKiEjOTckJCIxLCsrKz8qcHAoISM7JCIwSl11IT1leiIqISM5NyQkIjEsKytdLXc9IykhIzskIjFqNyMqRyQ+OEgqISM6NyQkIjEtKyt2Vmt6JSkhIzskIjEjb1EnZiMpXEAkKiEjOjckJCIxLCsrKyZHMHUpISM7JCIxM0xPOFwkUkwqISM6NyQkIitERVQsISohIzUkIjFWYDEsNUFHJCohIzo3JCQiMSwrK11uSGkjKiEjOyQiMUMreClwIypSSSohIzo3JCQiMW87SDJCNk8oKiEjOyQiMSIqNCgzXFA7QCohIzo3JCQiMk5MZWt5Iyo0LSIhIzskIjFhLVMpcCVvYiEqISM6NyQkIi92PVVWUG81ISM4JCIwWFd3LGJeJCkpISM5NyQkIjJwbTt6KmV2OjYhIzskIjF5IipmbzZoXCYpISM6NyQkIjItK0RKU1ElRzchIzskIjFsV046Qzs2dyEjOjckJCIyTkxMJDM0N1Q4ISM7JCIxPlY7c1ErQGohIzo3JCQiMk9MTDN4eGxWIiEjOyQiMThsYDpyImUoXCEjOjckJCIyUExMTGpNP2AiISM7JCIxQkVxJD0nKSlHTSEjOjckJCIycG0iSGRPMlY7ISM7JCIxMUclb2g+TlUiISM6NyQkIi1EIm83VHYiISM2JCExeFU3M1dMKlEoISM7NyQkIjJtbSJISzVTXz0hIzskITJWelt2UWBMcyMhIzs3JCQiMktMTExRKm9dPiEjOyQhMSgqKjNxJ0c0OlohIzo3JCQiMW47SCNHRiZlPyEjOiQhMWkkM2RQSS0lbyEjOjckJCIyMitdNz1sajsjISM7JCExTyVbJlIuZEgpKSEjOjckJCIyMCsrRE9fIXBBISM7JCEyKCl6JipcdG9LMCIhIzo3JCQiMi8rXVBhUjxQIyEjOyQhMlxYI1tXT0kqPiIhIzo3JCQiMnRtInpXRykpeUMhIzskITJkYWEwTTQySyIhIzo3JCQiMlVMTGU5RWdlIyEjOyQhMiUpZlUoUktPMTkhIzo3JCQiMU0zRlc7QU5FISM6JCEyQ1xNdjUhKUdWIiEjOjckJCIxTSQzRjk8V28jISM6JCEyTis+TC5tNVgiISM6NyQkIjElM0Y+KlssNEYhIzokITFEKHplOSoqcFgiISM5NyQkIjJQJGU5VEVoTEYhIzskITFAIUdjNjozWSIhIzk3JCQiMk1lay5SNSNlRiEjOyQhMnI0c2dNNkRZIiEjOjckJCIyTkwkZVIiM0d5IyEjOyQhMidHZnc4JiozaTkhIzo3JCQiMi9EY144UyQ0RyEjOyQhMSZlbC4rbiNmOSEjOTckJCIybztIMjhzZSRHISM7JCExNEw0Jz4rU1giISM5NyQkIjEkMy1qNy9DJ0chIzokITIsMy1zIXlJWTkhIzo3JCQiLnY9N08qKSlHISM3JCEyeG1LVz45aVYiISM6NyQkIjJMJDMtOCwrVUghIzskITFDMChbZFsqMzkhIzk3JCQiMm1tbVQ1a10qSCEjOyQhMlswKVJrIzNEUCIhIzo3JCQiMW47L143STBKISM6JCEybzF2Y0tqI3A3ISM6NyQkIjFubSJ6UlFiQCQhIzokITJIRCE0eF52SzYhIzo3JCQiMlBMTCRlLF02TCEjOyQhMT1nKD5aVE4iKiohIzo3JCQiLXY9PlkyTSEjNiQhMXchNCc9XzlPJCkhIzo3JCQiMU0kZSpbSzU2TiEjOiQhMVlkTjk5UTBsISM6NyQkIjJ2bW0ielh1OU8hIzskITEnZkNNYShvMFkhIzo3JCQiMU1MZTlpIj1zJCEjOiQhMSl6JylwYSU0U0UhIzo3JCQiKiZ5KSlHUSEiKSQhMUBtXVhHLzB1ISM7NyQkIi1EY2xqTFIhIzYkIi8tKioqW3IiPSoqISM5NyQkIixERSZRUVMhIzUkIjJiRWwySzFkYSMhIzs3JCQiLnY9LU4oUlQhIzckIjEoZV0jUnN5UFEhIzo3JCQiMjArXTd5JTNUVSEjOyQiMTJtO3FcUiIqWyEjOjckJCIudiQ0a2hgViEjNyQiMTghW3JeeSVlZCEjOjckJCIsdi5baFklISM1JCIxIz4pKUdEIjQoSCchIzo3JCQiMUwkZVJQMG5eJSEjOiQiMjA7PmpKS0tWJyEjOzckJCIxbW1UNUZFblghIzokIjFwJVJxSEJlXSchIzo3JCQiMU1la3k4YSNmJSEjOiQiMUZ6NEUiKSkqPWwhIzo3JCQiLnZvLz95aCUhIzckIjE6PmU6Uz08bCEjOjckJCIxbVQ1OigpNFZZISM6JCIxdWBAVGp6K2whIzo3JCQiMUxMTCRReCRvWSEjOiQiMTxJakJrQ3FrISM6NyQkIjJtbW0ieilRanglISM7JCIwUnE5Iz1zIj4nISM5NyQkIjIwKytdUCtWKVshIzskIjEoejxoeHZlcSYhIzo3JCQiMW1tInpwZSp6XSEjOiQiMCRlI0cxWllZJSEjOTckJCIxLCsrRFwnUUgmISM6JCIyakpJdWdpYSNIISM7NyQkIjFMTGU5UzgmXCYhIzokIjFfKCpvU3EjZXAiISM6NyQkIjFubVQ1aEsrYyEjOiQiMilRSjUkKSoqUlw3ISM7NyQkIjEsK0QxIz1icSYhIzokIjEjPjRkXEFRJCkqISM7NyQkIjE8L15KJEg3dCYhIzokIjFWbmpsLSVmWyohIzs3JCQiMU0zeGMvJXB2JiEjOiQiMSM9XzxCPz9FKiEjOzckJCIwREo/ZV5FeSYhIzkkIjFATi12cShROyohIzs3JCQiMW47SDJGTzNlISM6JCIxakQyUE8qRz4qISM7NyQkIjEkM19EJFEyTWUhIzokIjFuKj09XzsrTiohIzs3JCQiL0QieSZceWZlISM4JCIxVClSbCJIdk4nKiEjOzckJCIxO0gyJDMnXCYpZSEjOiQiMkJiPlQnbywwNSEjOzckJCIxTExMM3M/NmYhIzokIjJ4PmVKdSlHZjUhIzs3JCQiMW07enBlKCk9ZyEjOiQiMi92YVpKMUxVIiEjOzckJCIyJioqKlw3YFdsNychIzskIjJlcWg3MmtiKj4hIzs3JCQiMW9tbW0qUlJMJyEjOiQiMnRxTF5hLV5kJCEjOzckJCIxb207YTwuWWwhIzokIjFqInpOLj51YiYhIzo3JCQiMlhMJGU5dE9jbiEjOyQiMiZSWiNHIzN1anUhIzs3JCQiMW87SCNlMEkmbyEjOiQiMWsiXCMqNEZLPikhIzo3JCQiMSwrK11Ra1xwISM6JCIxLnEoWzYzOXgpISM6NyQkIjFOJDMtLkJdKyghIzokIjE7MWhaIyl6PiEqISM6NyQkIjFvbVQ1QVNncSEjOiQiMWRUJGZzJHoqPiohIzo3JCQiMU0zXys9NCkzKCEjOiQiMXUvUGd3UGkjKiEjOjckJCIxLF1pIVIieTpyISM6JCIxOW1sSSdSZUkqISM6NyQkIjFvIkgyKTRaVnIhIzokIjE6JUdFISozJ0gkKiEjOjckJCIxTkwkM2RnNjwoISM6JCIxJik0RVckb0pMKiEjOjckJCIxLURKcWsjZj4oISM6JCIxUihlcyE9Jik9JCohIzo3JCQiMW87enBCcD9zISM6JCIxTlU1N2txKEcqISM6NyQkIjFNM0ZwI2VhQyghIzokIjEjMyQqKkcjeSVSIyohIzo3JCQiMSwrdm9UQXFzISM6JCIxVF85TjwmUjwqISM6NyQkIjFOJDN4J2Z2PnQhIzokIjFCc2hsIVEuKiopISM6NyQkIjFvbW1teEdwdCEjOiQiMUt3b25bJGZ0KSEjOjckJCIxTSRlUkE1XFooISM6JCIxXj5yWV5mY3ohIzo3JCQiLUQib0swZSghIzYkIjFBW1ZiL0xpbyEjOjckJCIxLCsrXW9pIm8oISM6JCIyYlgsOlwqPlJiISM7NyQkIjEsK3Y9NXMjeSghIzokIjB0RjJfPWgoUiEjOTckJCIyMSt2JDQwTyIqeSEjOyQiMm1GJlxLbDN1PyEjOzckJCIjISkhIiIkIjFbVFZRb0NOUyEjPi0lJkNPTE9SRzYmJSRSR0JHJCIjNSEiIiQiIiEhIiIkIiIhISIiLSUlVklFV0c2JDskISM/ISIiJCIjISkhIiIlKERFRkFVTFRHLSUrQVhFU0xBQkVMU0c2Jy1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRIng2Ii8lJ2ZhbWlseUdRITYiLyUlc2l6ZUdRIzEwNiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2IlEhNiItJSVGT05URzYlJShERUZBVUxURyUoREVGQVVMVEciIzUlK0hPUklaT05UQUxHJStIT1JJWk9OVEFMRy0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJEkiISIiLSUpQk9VTkRTX1lHNiMkIiRTIiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJSF6JCEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiUheiQhIiItJSlDSElMRFJFTkc2Ig==
Now we choose points that we normally see more frequently when using trig functions, e.g. multiples of pi.
pts2:=[[0,1],[Pi/2,2],[Pi,-2],[3*Pi/2,1]];
NyY3JCIiISIiIjckLCRJI1BpRyUqcHJvdGVjdGVkRyNGJSIiI0YrNyRGKCEiIzckLCRGKCMiIiRGK0Yl
We plot our four functions and the points in question. Notice that the x values of these points correspond to zeros of many of these functions. In particular, sin(2x) is the blue curve, and note that sin(2x) is zero for all the x-coordinates of the points found in pts2 data set.
plot3:=plot({sin(x),cos(x),sin(2*x),cos(2*x)},x=0..2*Pi):
plot4:=pointplot(pts2,symbol=circle,symbolsize=30):
display({plot3,plot4});
NiotJSdQT0lOVFNHNic3JCQiIiEhIiIkIiM1ISIiNyQkIitGanpxOiEiKiQiIz8hIiI3JCQiK2FFZlRKISIqJCEjPyEiIjckJCIrIikqKVE3WiEiKiQiIzUhIiItJSdTWU1CT0xHNiQlJ0NJUkNMRUciI0ktJSdDVVJWRVNHNiQ3ZW43JCQiIiEhIiIkIiM1ISIiNyQkIjFgI0dpI3p4Wm8hIzwkIjF3VUNISmN3KiohIzs3JCQiMjFsWF9lYiZwOCEjPCQiMEZmVlBpaiEqKiEjOjckJCIyMiIpPjYzeGAnPiEjPCQiMXZIR01iWzIpKiEjOzckJCIyMChSKnBkKT5oRCEjPCQiMSN5Ly01LVFuKiEjOzckJCIyRChbRV5kSyxSISM8JCIxLDEtI1snZVsjKiEjOzckJCIxY15OZWhMXV8hIzskIjFENVNDNDJgJykhIzs3JCQiMSpRaFcmXCRIZichIzskIjFRXlpUP0QveiEjOzckJCIxVDExLmZwUHkhIzskIjFfPlUhPXVEMyghIzs3JCQiMDxuZnd0bDcqISM6JCIxZlxwaGRYO2ghIzs3JCQiMlFSYSY0TCZmLyIhIzskIjF5SGV2Ong1XSEjOzckJCIyTmNrMndAKXk2ISM7JCIwJylweW1SLCNRISM6NyQkIjAoR1clSCRcOjghIzkkIjEub0sualFERCEjOzckJCIxKmUpcGJPKGVWIiEjOiQiMmAxOUgmbzhYOCEjPDckJCIyNl5oL1YkUnI6ISM7JCExOU5MNmoucmYhIz43JCQiMj0zU0NtcHVxIiEjOyQhMUAsQUJCW2k4ISM7NyQkIjJqMWtpSy8nUT0hIzskITF4KSpmN0A9WUUhIzs3JCQiMk1iVFI7KG9kPiEjOyQhMUNwWCMqKTNKeCQhIzs3JCQiMVUoeWU8KUcqNCMhIzokITFpdCxOX0pVXSEjOzckJCIyVyZbI28hR0M+QSEjOyQhMV8pKm9RJSpbUmchIzs3JCQiMVl3SmYmeShlQiEjOiQhMDpROEo7JCozKCEjOjckJCIyZDgocFY4SCNbIyEjOyQhMXRBPHRUPy96ISM7NyQkIjIkelkpUVcveWgjISM7JCExXyJHPnIxJGYnKSEjOzckJCIxdik+OCdcJW91IyEjOiQhMWA2a11pJDRCKiEjOzckJCIyUUdDUz8mWyIpRyEjOyQhMTNHaSRcQk9tKiEjOzckJCIxJTMkXCE0MUwlSCEjOiQhMUQpZU9ZYlMhKSohIzs3JCQiMlgpPSdwKHA3MEkhIzskITF3NXI1Ky4yKiohIzs3JCQiMk86SVo6ND0yJCEjOyQhMVFSIWU+aGMoKiohIzs3JCQiMkZVKVxLOFxRSiEjOyQhMSgpZjI0PiYqKioqKiEjOzckJCIxNF1edWB2MkshIzokITFzPy5iLzd5KiohIzs3JCQiMmNmSmxUPnFGJCEjOyQhMSl6VyJHIVElMyoqISM7NyQkIjIxPUAob1JKUEwhIzskITEudDV3azI0KSohIzs3JCQiMWwyIjRfM3dSJCEjOiQhMSYzdmJbbFNuKiEjOzckJCIyanQnR2QhW3lfJCEjOyQhMUVYTFRBRWojKiEjOzckJCIxPjRKVSMpUmlPISM6JCExQ1ohNDwnPXUnKSEjOzckJCIxOCFmJVskSFN6JCEjOiQhMSlvSi4jeTFZeiEjOzckJCIyJG9FKzcjKlFAUiEjOyQhMSg0MmUqZmM1ciEjOzckJCIxKWYwaWkrRzElISM6JCExOlUsKlwnZVtnISM7NyQkIjFmT1JyXScpKj0lISM6JCEyYVRyYlwpbyEqXCEjPDckJCIxJnAlKnpbTGJLJSEjOiQhMWthLW46eXNQISM7NyQkIjEncEV4QnAlW1chIzokITJqR1VGbmwnM0UhIzw3JCQiMWkveCRbcUdlJSEjOiQhMiM0ZCxzb2MiSCIhIzw3JCQiMTxPOyNlSiQ0WiEjOiQhMWUsInlYJFJkSSEjPTckJCIxcE80byk+OiVbISM6JCIyKCp5IVJeSnMoRyIhIzw3JCQiMXUiNCV6IWUyKFwhIzokIjFDYU9WNy9iRCEjOzckJCIxZUNbNCZlZzUmISM6JCIxbj45Jkcpek5RISM7NyQkIjBmIVsvKm9qQiYhIzkkIjEyK2J6L0kuXSEjOzckJCIxZE5eLyxqcGAhIzokIjEtRycqMy5ONGghIzs3JCQiMUEiR1gkeXksYiEjOiQiMU9cJlJJKyQqNCghIzs3JCQiMUwxL2pxQUJjISM6JCIxeXRJPyR5LCF6ISM7NyQkIjF4ZUtCLFRpZCEjOiQiMSVRYl5YUFZuKSEjOzckJCIxMiRlOzYoKm8pZSEjOiQiMWNvI1shNCtEIyohIzs3JCQiMUxjdTBpaT5nISM6JCIwbUZlJCkpb2EnKiEjOjckJCIxOj0+WGI5JDMnISM6JCIxPioqeXZuZSspKiEjOzckJCIxKip6aiUpW21ZaCEjOiQiMSdcXiI9YiZwISoqISM7NyQkIjEqKipHJyozRFxAJyEjOiQiMWpUOmU8cncqKiEjOzckJCItaSVIJj0kRychIzYkIjEqKioqKioqKioqKioqKioqISM7LSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJSdDVVJWRVNHNiQ3XXM3JCQiIiEhIiIkIiM1ISIiNyQkIjJMMWQ6W1c+ciIhIz0kIjEqKUh5ayFSVCoqKiEjOzckJCIybTc5SicqKSlRVSQhIz0kIjF3VUNISmN3KiohIzs3JCQiMiYqPXJZV0xlOCYhIz0kIjF6eG0meiNIWioqISM7NyQkIjFgI0dpI3p4Wm8hIzwkIjBGZlZQaWohKiohIzo3JCQiMnpCTSopb21yLSIhIzwkIjAjKVtEekUoKnkqISM6NyQkIjIxbFhfZWImcDghIzwkIjEwVVUpNC5zaSohIzs3JCQiMjIiKT42M3hgJz4hIzwkIjFrUi4sWE5QIyohIzs3JCQiMjAoUipwZCk+aEQhIzwkIjEjSCZ5OSUqWzsoKSEjOzckJCIyOlVIVDtpN0IkISM8JCIxJHA5dnpoTSl6ISM7NyQkIjJEKFtFXmRLLFIhIzwkIjFQRzMkUXFzNSghIzs3JCQiMTkrIlsmNCRlZCUhIzskIjEwRXQtJillJzQnISM7NyQkIjFjXk5laExdXyEjOyQiMkIpR2ElR0ZeKFwhIzw3JCQiMXMjM2tiTjsjZiEjOyQiMlZMRmIoUj9wUCEjPDckJCIxKlFoVyZcJEhmJyEjOyQiMiZ6Ri5uK1cmXCMhIzw3JCQiMTs1d0dhSjpzISM7JCIyOC52cHpoUUYiISM8NyQkIjFUMTEuZnBQeSEjOyQiMmRfJGZCLzlkSyEjPjckJCIxMVJeTVs4I1spISM7JCExT2wmPkcvSUQiISM7NyQkIjA8bmZ3dGw3KiEjOiQhMi5yWEl5JHo8RCEjPDckJCIxYWJ2SU0wJHoqISM7JCEyZCE9OikqKWY7eSQhIzw3JCQiMlFSYSY0TCZmLyIhIzskITJsJGZZVk9WeVwhIzw3JCQiMiZ5JWZeYChRNzYhIzskITEzNXMvLVwkMychIzs3JCQiMk5jazJ3QCl5NiEjOyQhMUpNZiRlMTgzKCEjOzckJCIyb3Iud19kckMiISM7JCExaScpKilvU0d4eiEjOzckJCIwKEdXJUgkXDo4ISM5JCExd1F6LltbQygpISM7NyQkIjIlSDIydk1vdjghIzskITE6TmprQUJbIyohIzs3JCQiMSplKXBiTyhlViIhIzokITEqWydRcDg3UScqISM7NyQkIjIlPiQqUSpmYChwOSEjOyQhMUM6OzU4XSd6KiEjOzckJCIwMCEzVk5qLjohIzkkITF0c0V6KD0qNCoqISM7NyQkIjJhVERcXnQwXyIhIzskITEneltNeSZmXCoqISM7NyQkIjIyeXFuWzh2YCIhIzskITFKam9JSyZ5KCoqISM7NyQkIjFZaGhlTVhhOiEjOiQhMTRddilvZVkqKiohIzs3JCQiMjZeaC9WJFJyOiEjOyQhMXRdTXBHKioqKioqISM7NyQkIjJEJCllJWZIUyllIiEjOyQhMUMiUS54K1EqKiohIzs3JCQiMlE6YyUpWzdhZyIhIzskITBwJil6KlIvdyoqISM6NyQkIjJeWmB1LEFDaSIhIzskITAxYCw1Vm4lKiohIzo3JCQiMmx6XWthSiVSOyEjOyQhLzRSJCk+JGYhKiohIzk3JCQiMiJSYVcvMVh0OyEjOyQhMU93KD1FIikqKnkqISM7NyQkIjI9M1NDbXB1cSIhIzskITBhI3kkUUcoRycqISM6NyQkIjF1P04lKnAudDwhIzokITFuYEVzajEkPiohIzs3JCQiMmoxa2lLLydRPSEjOyQhMSJIIlJYU2EqZikhIzs3JCQiMio0RzVYZDkpKj0hIzskITEqR09mWzhCJHohIzs3JCQiMk1iVFI7KG9kPiEjOyQhMSdcS3ImKUhGOighIzs3JCQiMnY5NSpwd1tHPyEjOyQhMT1eKip6QShbNCchIzs3JCQiMVUoeWU8KUcqNCMhIzokITFpQjk5OSw6XCEjOzckJCIxKXpeOFxsI2ZAISM6JCEyJnotJCpRKFF2JFEhIzw3JCQiMlcmWyNvIUdDPkEhIzskITJqdSEqUlk4XHEjISM8NyQkIjIvRHJJbzUhKkcjISM7JCEyYVV4Q19PJ1I4ISM8NyQkIjFZd0pmJnkoZUIhIzokIjFMSTVXXzZvXiEjPTckJCIxIlIyOiZcYD9DISM6JCIybS5aTyQzRSRHIiEjPDckJCIyZDgocFY4SCNbIyEjOyQiMmR6UUFzKUcmXCMhIzw3JCQiMnMhNHokKnkvXUQhIzskIjE4WEdEVW0heSQhIzs3JCQiMiR6WSlRVy95aCMhIzskIjEqeVtpYT1uKlwhIzs3JCQiMXhBZy1aSyNvIyEjOiQiMSJ5XW5YVilwZyEjOzckJCIxdik+OCdcJW91IyEjOiQiMSo9IW83by5VcSEjOzckJCIyJXo/biMzbFQiRyEjOyQiMXZscTk7T0p6ISM7NyQkIjJRR0NTPyZbIilHISM7JCIxUCgqSDBRN3gnKSEjOzckJCIxJTMkXCE0MUwlSCEjOiQiMSRwKipvcStSQSohIzs3JCQiMlgpPSdwKHA3MEkhIzskIjEmKWZpRShbKUgnKiEjOzckJCIyJD5nJWUxbyVRSSEjOyQiMSpbb2ZIZiEpeSohIzs3JCQiMk86SVo6ND0yJCEjOyQiMHAkcGVLdy0qKiEjOjckJCIxQEE8KnB6JSkzJCEjOiQiMWhyZTlMalYqKiEjOzckJCIxKUc5T0NdXjUkISM6JCIxVC43cjJYdCoqISM7NyQkIjJgTmMhKXk/PTckISM7JCIxVXlWJlsjPSMqKiohIzs3JCQiMkZVKVxLOFxRSiEjOyQiMSZcbW5qMikqKioqISM7NyQkIjIkcEQrVnQhZTokISM7JCIxMT5GQyJmZioqKiEjOzckJCIxO25dYEw3dEohIzokIjFDM3UoKUc3ISkqKiEjOzckJCIyQyczLGskUi8+JCEjOyQiMSkqMyopPXpKXyoqISM7NyQkIjE0XV51YHYySyEjOiQiMSopZS1qdmQ3KiohIzs3JCQiMS1MXyZSKFFVSyEjOiQiMWFJSCNlJlwoeiohIzs3JCQiMmNmSmxUPnFGJCEjOyQiMSg+PDokKUdhaiohIzs3JCQiMjE9QChvUkpQTCEjOyQiMU5WI0dFJ2ZWIyohIzs3JCQiMWwyIjRfM3dSJCEjOiQiMUIhUlFnM3ZyKSEjOzckJCIyMHYpNCpHR0ZZJCEjOyQiMU4hKnl1SFAyISkhIzs3JCQiMmp0J0dkIVt5XyQhIzskIjB6WS5aMDs7KCEjOjckJCIydyMpKXpcSjcmZiQhIzskIjE4QkpIWGpnaCEjOzckJCIxPjRKVSMpUmlPISM6JCIxYUhdWDZJW10hIzs3JCQiMW1cUSZ6OCNHUCEjOiQiMnZXM0YvVjsoUSEjPDckJCIxOCFmJVskSFN6JCEjOiQiMms6aTV3KSp6aSMhIzw3JCQiMjElM0IhRzR4JlEhIzskIjEqUiZ6TTU/IlEiISM7NyQkIjIkb0UrNyMqUUBSISM7JCIySE5jOWMoSD82ISM9NyQkIjFMVDU+XDQjKlIhIzokITJBeEY+KGVTKUgiISM8NyQkIjEpZjBpaStHMSUhIzokITEmZil5cy0jSG8jISM7NyQkIjImRycqeltHTEVUISM7JCExYTE8eSVmPylRISM7NyQkIjFmT1JyXScpKj0lISM6JCExRFUub2NnPV0hIzs3JCQiMXlUcHojKnBkVSEjOiQhMTYxKT4nZk5VaCEjOzckJCIxJnAlKnpbTGJLJSEjOiQhMXhEXVxRQWByISM7NyQkIjEncGdHTyxxUSUhIzokITFbIlttKFI5Y3ohIzs3JCQiMSdwRXhCcCVbVyEjOiQhMTE3ZCJvcypRJykhIzs3JCQiMXomWzInKXBjXiUhIzokITElcHJ6O2tmQiohIzs3JCQiMWkveCRbcUdlJSEjOiQhMTdCLS41UG0nKiEjOzckJCIwdm8kZWRbOVkhIzkkITEkKW8obykzIiozKSohIzs3JCQiMVJxJ0guLGhrJSEjOiQhMWR0VzBvQzcqKiEjOzckJCIxJD1tLW4zPm0lISM6JCEwJyk9Py95IVwqKiEjOjckJCIxR2BjMmpyeFkhIzokITFub3lbZCdmKCoqISM7NyQkIjFzVydbJVJfJHAlISM6JCExLSpRSTAkKUcqKiohIzs3JCQiMTxPOyNlSiQ0WiEjOiQhMSQpWyFwLzgpKioqKiEjOzckJCIxdFshejZiZXMlISM6JCExRjMtI1t0aioqKiEjOzckJCIwOFlPbHlCdSUhIzkkITEyYSg9eTw/KSoqISM7NyQkIjEoUShRKj0tKmVaISM6JCExMy0uQjt3YyoqISM7NyQkIjFWJ0dec0RheCUhIzokITFzPFgmZUsxIyoqISM7NyQkIjFjNmgnenMlM1shIzokITE6W0o8YiNmIikqISM7NyQkIjFwTzRvKT46JVshIzokITFNVzI8UU5vJyohIzs3JCQiMUA5dnQqUWghXCEjOiQhMUckKXlYa2NlIyohIzs3JCQiMXUiNCV6IWUyKFwhIzokITFMLiUqW0dOJXApISM7NyQkIjE7ZVclSDMlUV0hIzokITFrSDleUltbeiEjOzckJCIxZUNbNCZlZzUmISM6JCExKlFTSklJdDAoISM7NyQkIjFDOilwcTg3PCYhIzokITF5eHAsZSlvMichIzs3JCQiMGYhWy8qb2pCJiEjOSQhMU1wRWlvUiQqXCEjOzckJCIxdHFcLyYqKkhJJiEjOiQhMS0hZTlfYnp6JCEjOzckJCIxZE5eLyxqcGAhIzokITFCX1pneDtORCEjOzckJCIwJTNfcCozZFYmISM5JCEyMjQjKTM/dyRRNyEjPDckJCIxQSJHWCR5eSxiISM6JCIxIipIcDloSCwhKSEjPTckJCIxeVZ5W3VdaWIhIzokIjIlXCdIY18leiFIIiEjPDckJCIxTDEvanFBQmMhIzokIjAlNGApXGpEWyMhIzo3JCQiMldEJD0kZj1HcCYhIzskIjI4K0ovIWZZLVEhIzw3JCQiMXhlS0IsVGlkISM6JCIxPUI9ZGcjKVtdISM7NyQkIjEjNCNcPE9sQ2UhIzokIjFtKyc+QFI6MychIzs3JCQiMTIkZTs2KCpvKWUhIzokIjBFNilRJEcsLSghIzo3JCQiMCg+cWU7RWBmISM5JCIxIipbJzRaZTMheiEjOzckJCIxTGN1MGlpPmchIzokIjFcNUwuTGdVJykhIzs3JCQiMTo9PlhiOSQzJyEjOiQiMVd3JipILUk1IyohIzs3JCQiMSoqemolKVttWWghIzokIjF1a3NGTmJIJyohIzs3JCQiMSpcTHIpXHohPSchIzokIjFgKHBWOmM1eiohIzs3JCQiMSoqKkcnKjNEXEAnISM6JCIxMyRlNV1icCEqKiEjOzckJCIwdnczOSEqPkInISM5JCIxRy1qeW9pWioqISM7NyQkIi9YNyM+YiFcaSEjOCQiMT13ZVw8cncqKiEjOzckJCIxXEFQVi03bWkhIzokIjFMSUlQaTwlKioqISM7NyQkIi1pJUgmPSRHJyEjNiQiMSgqKioqKioqKioqKioqKiohIzstJSZDT0xPUkc2JiUkUkdCRyQiMidcXydlTSV5Zz4hIzwkIjFpajkuZUBSISkhIzskIjInXF8nZU0leWc+ISM8LSUnQ1VSVkVTRzYkN2VuNyQkIiIhISIiJCIiISEiIjckJCIyMWxYX2ViJnA4ISM8JCIxOSZvJ3oieV9PIiEjOzckJCIyMChSKnBkKT5oRCEjPCQiMi55PVYmKilHTEQhIzw3JCQiMkQoW0VeZEssUiEjPCQiMi92c3JKNUohUSEjPDckJCIxY15OZWhMXV8hIzskIjEsKnkoSDRVN10hIzs3JCQiMSpRaFcmXCRIZichIzskIjEjcEQqZWNlRGghIzs3JCQiMVQxMS5mcFB5ISM7JCIwS3IhKm9VJmZxISM6NyQkIjA8bmZ3dGw3KiEjOiQiMSkqUWpQIT44InohIzs3JCQiMlFSYSY0TCZmLyIhIzskIjFbJTN3N0VTbCkhIzs3JCQiMk5jazJ3QCl5NiEjOyQiMWkuWixiY1QjKiEjOzckJCIwKEdXJUgkXDo4ISM5JCIxT3BTRiJvZW4qISM7NyQkIjIlSDIydk1vdjghIzskIjFfKlxUJ3pENSkqISM7NyQkIjEqZSlwYk8oZVYiISM6JCIxJ2ZJJ2Z0NjQqKiEjOzckJCIwMCEzVk5qLjohIzkkIjE6M016VVh4KiohIzs3JCQiMjZeaC9WJFJyOiEjOyQiMDZPdEApKioqKioqISM6NyQkIjJsel1rYUolUjshIzskIjF5ODp5X1h3KiohIzs3JCQiMj0zU0NtcHVxIiEjOyQiMTxOZCNIWm4hKiohIzs3JCQiMXU/TiUqcC50PCEjOiQiMSpwUS5tKj0neiohIzs3JCQiMmoxa2lLLydRPSEjOyQiMTJEYWZEYFYnKiEjOzckJCIyTWJUUjsob2Q+ISM7JCIxQ0NiXmwnM0UqISM7NyQkIjFVKHllPClHKjQjISM6JCIxMio+YzQmb04nKSEjOzckJCIyVyZbI28hR0M+QSEjOyQiMWUqUnAxSS0oeiEjOzckJCIxWXdKZiZ5KGVCISM6JCIxO3MoZjRzRjAoISM7NyQkIjJkOChwVjhII1sjISM7JCIxM040RXVrRGghIzs3JCQiMiR6WSlRVy95aCMhIzskIjFQYXdkL2ssXSEjOzckJCIxdik+OCdcJW91IyEjOiQiMjsqUiVSKEh2WFEhIzw3JCQiMlFHQ1M/JlsiKUchIzskIjEmbzZmKVElPWQjISM7NyQkIjJYKT0ncChwNzBJISM7JCIyODFeJW9dVWc4ISM8NyQkIjJGVSlcSzhcUUohIzskIjImKXBiKT5oSixKISM+NyQkIjJjZkpsVD5xRiQhIzskITFYVDAoeUosTiIhIzs3JCQiMWwyIjRfM3dSJCEjOiQhMSMqPUdlSEdLRCEjOzckJCIyanQnR2QhW3lfJCEjOyQhMSopenNtTUFuUCEjOzckJCIxPjRKVSMpUmlPISM6JCExLEohKT0zenZcISM7NyQkIjE4IWYlWyRIU3okISM6JCExMTEkWz9XNzInISM7NyQkIjIkb0UrNyMqUUBSISM7JCExSzRwX3hNSnEhIzs3JCQiMSlmMGlpK0cxJSEjOiQhMVdxWk4mR0wneiEjOzckJCIxZk9Scl0nKSo9JSEjOiQhMSozJjNuTGlsJykhIzs3JCQiMSZwJSp6W0xiSyUhIzokITEkWyRRMCoqKjRFKiEjOzckJCIxJ3BFeEJwJVtXISM6JCExJTRnN25bUGwqISM7NyQkIjF6JlsyJylwY14lISM6JCExLVB5QDY4MikqISM7NyQkIjFpL3gkW3FHZSUhIzokITFZUFkhKT5DOyoqISM7NyQkIjFScSdILixoayUhIzokITFlZihSZVAheSoqISM7NyQkIjE8TzsjZUokNFohIzokITEnKnBoaEsmKioqKiohIzs3JCQiMVYnR15zRGF4JSEjOiQhMTElXEFVUSwpKiohIzs3JCQiMXBPNG8pPjolWyEjOiQhMV42T2U9dTsqKiEjOzckJCIxQDl2dCpRaCFcISM6JCExSykqKVs3IipHIikqISM7NyQkIjF1IjQleiFlMihcISM6JCExJDQ5bGx6IW8nKiEjOzckJCIxZUNbNCZlZzUmISM6JCExYCZ5UUJ4XUIqISM7NyQkIjBmIVsvKm9qQiYhIzkkITEjKm9NbXdNZScpISM7NyQkIjFkTl4vLGpwYCEjOiQhMXJvcWh0IW8ieiEjOzckJCIxQSJHWCR5eSxiISM6JCExNlR3VkBzVXEhIzs3JCQiMUwxL2pxQUJjISM6JCExbDRMWSdRMzgnISM7NyQkIjF4ZUtCLFRpZCEjOiQhMllzQSF6cV92XCEjPDckJCIxMiRlOzYoKm8pZSEjOiQhMigpeV9CcG8qZlEhIzw3JCQiMUxjdTBpaT5nISM6JCEyWTdPTSdSPTBFISM8NyQkIjEqKnpqJSlbbVloISM6JCEybCo9I2VWbjRPIiEjPDckJCItaSVIJj0kRychIzYkITFVRVtdJ2VmRCIhI0ItJSZDT0xPUkc2JiUkUkdCRyQiMTMkR1ZwPiFcJikhIzskIjFfTW1YJSllcWshIzskIjAleipcLi1cRCIhIzotJSdDVVJWRVNHNiQ3W3M3JCQiIiEhIiIkIiIhISIiNyQkIjFgI0dpI3p4Wm8hIzwkIjE5Jm8neiJ5X08iISM7NyQkIjIxbFhfZWImcDghIzwkIjJuUD96MikpXHEjISM8NyQkIjIyIik+NjN4YCc+ISM8JCIxYSVHLSQ0SklRISM7NyQkIjIwKFIqcGQpPmhEISM8JCIxPkhNLHNJLFwhIzs3JCQiMjpVSFQ7aTdCJCEjPCQiMSxUOmJ5KT4tJyEjOzckJCIyRChbRV5kSyxSISM8JCIxJlIyUSQqeVkuKCEjOzckJCIxOSsiWyY0JGVkJSEjOyQiMVBuNnciUm0jeiEjOzckJCIxY15OZWhMXV8hIzskIjEkemw7d21YbikhIzs3JCQiMXMjM2tiTjsjZiEjOyQiMVI6cl9sWGkjKiEjOzckJCIxKlFoVyZcJEhmJyEjOyQiMVopZSoqZU1PbyohIzs3JCQiMS03aCI+RFQhcCEjOyQiMSNRKVJybDQ/KSohIzs3JCQiMTs1d0dhSjpzISM7JCIxJ0glKVImPmA9KiohIzs3JCQiMUFmTFowIjRQKCEjOyQiMF5Fd3ZrTCYqKiEjOjckJCIxRzMiZm0wbF8oISM7JCIxJCopKilSY2YmeSoqISM7NyQkIjFOZFsleStAbyghIzskIjFmSmJ3PjQlKioqISM7NyQkIjFUMTEuZnBQeSEjOyQiMTsxL2JwJSoqKioqISM7NyQkIjFlUiNmajApKSp6ISM7JCIwKWUmXFwwZSoqKiEjOjckJCIxdXN5b2AiKmYiKSEjOyQiMUhUS1xuRyIpKiohIzs3JCQiMGZdOzVENUspISM6JCIwZiUpNHowayYqKiEjOjckJCIxMVJeTVs4I1spISM7JCIxR0dYXSUpPUAqKiEjOzckJCIxUTBDK1ZOLykpISM7JCIxXyM9N2MsKj4pKiEjOzckJCIwPG5md3RsNyohIzokIjF1Jm8/ZVl5biohIzs3JCQiMWFidklNMCR6KiEjOyQiMSNRQEgweXREKiEjOzckJCIyUVJhJjRMJmYvIiEjOyQiMXBYInBpcEVuKSEjOzckJCIyJnklZl5gKFE3NiEjOyQiMT46YkFpcE96ISM7NyQkIjJOY2syd0ApeTYhIzskIjAicDIkPjkzMSghIzo3JCQiMm9yLndfZHJDIiEjOyQiMWsmSCRmIXAsLichIzs3JCQiMChHVyVIJFw6OCEjOSQiMW8qXCRvNDEoKVshIzs3JCQiMiVIMjJ2TW92OCEjOyQiMW0hPWdnclIhUSEjOzckJCIxKmUpcGJPKGVWIiEjOiQiMm0jejApeUJlbSMhIzw3JCQiMDAhM1ZOai46ISM5JCIyUWhrSiUpPiNSOCEjPDckJCIyNl5oL1YkUnI6ISM7JCEya296JFxxPyU+IiEjPjckJCIybHpda2FKJVI7ISM7JCEyTlY7XlYoUm84ISM8NyQkIjI9M1NDbXB1cSIhIzskITFVIVtMaGAmKnAjISM7NyQkIjF1P04lKnAudDwhIzokITJjJ0hhNytVTlIhIzw3JCQiMmoxa2lLLydRPSEjOyQhMXlyRXMnM1A1JiEjOzckJCIyKjRHNVhkOSkqPSEjOyQhMUQrNSh5LiMqMychIzs3JCQiMk1iVFI7KG9kPiEjOyQhMW9RXShlXiUpKXAhIzs3JCQiMVUoeWU8KUcqNCMhIzokITFDOFwvJHAoMygpISM7NyQkIjEpel44XGwjZkAhIzokITFmc1JVVE5NIyohIzs3JCQiMlcmWyNvIUdDPkEhIzskITFSQ2tpU0FGJyohIzs3JCQiMkMwW1x1RVREIyEjOyQhMVgoKXpWZ08jeiohIzs3JCQiMi9EcklvNSEqRyMhIzskITBoImVnQycpNCoqISM6NyQkIjIjXEc4X0VYMUIhIzskITE5aUIhemIwJioqISM7NyQkIjFbVz5AWSpRSyMhIzokITE7dltDOzl6KiohIzs3JCQiMVpnRCFmTzhNIyEjOiQhMXg7MiI9JmUmKioqISM7NyQkIjFZd0pmJnkoZUIhIzokITFsQi5faycpKioqKiEjOzckJCIxIzNsdGw8VVAjISM6JCExRDJLN1NdJCoqKiEjOzckJCIyKT1EVGJubCpRIyEjOyQhMUpDVCQqUWh4KiohIzs3JCQiMWIqZk0mZTQwQyEjOiQhMXhTdVc3QF8qKiEjOzckJCIxIlIyOiZcYD9DISM6JCExTG4oW0c/dCIqKiEjOzckJCIyTUUtdzk4OVgjISM7JCExMSVcVV06Iz4pKiEjOzckJCIyZDgocFY4SCNbIyEjOyQhMDFVKCpldE9vKiEjOjckJCIycyE0eiQqeS9dRCEjOyQhMFFsXXMleWQjKiEjOjckJCIyJHpZKVFXL3loIyEjOyQhMCRmSmR4OWknKSEjOjckJCIxeEFnLVpLI28jISM6JCExRF5Sdnk4WnohIzs3JCQiMXYpPjgnXCVvdSMhIzokITFCIlF3NSEpKio0KCEjOzckJCIyJXo/biMzbFQiRyEjOyQhMSNwN1pkVi80JyEjOzckJCIyUUdDUz8mWyIpRyEjOyQhMW5fYyo+bTEoXCEjOzckJCIxJTMkXCE0MUwlSCEjOiQhMSp6UTgtKGZpUSEjOzckJCIyWCk9J3AocDcwSSEjOyQhMnRtJDQnUmFicCMhIzw3JCQiMk86SVo6ND0yJCEjOyQhMlhKOHpOVjZSIiEjPDckJCIyRlUpXEs4XFFKISM7JCExQyl6MVQtRT8nISM9NyQkIjE0XV51YHYySyEjOiQiMmNNbWI0J1I+OCEjPDckJCIyY2ZKbFQ+cUYkISM7JCIybigzUEclUmJuIyEjPDckJCIyMT1AKG9SSlBMISM7JCIxbSIpUjNLQTpRISM7NyQkIjFsMiI0XzN3UiQhIzokIjJYJ1xHTFVcKipbISM8NyQkIjIwdik0KkdHRlkkISM7JCIxRVx4aW46ISpmISM7NyQkIjJqdCdHZCFbeV8kISM7JCIxW0AhKXlkTnpwISM7NyQkIjJ3IykpelxKNyZmJCEjOyQiMVgoUjVIbHAoeSEjOzckJCIxPjRKVSMpUmlPISM6JCIxNmAmPj0oPUsnKSEjOzckJCIxbVxRJno4I0dQISM6JCIxSD94I1EyLEEqISM7NyQkIjE4IWYlWyRIU3okISM6JCIxWGYyVE9dWycqISM7NyQkIjFGXE05JHBlI1EhIzokIjF2JlJZJ1w9J3oqISM7NyQkIjIxJTNCIUc0eCZRISM7JCIxTkdqcVs6LyoqISM7NyQkIjJ3enRKRUhPKFEhIzskIjFwUUtDXDVWKiohIzs3JCQiMlV2O2hDXCYqKVEhIzskIjFwZmV1Yyg+KCoqISM7NyQkIjI4cmYhSCNwYSFSISM7JCIxc04/Ynl0ISoqKiEjOzckJCIyJG9FKzcjKlFAUiEjOyQiMWZkKnBXcyQqKioqISM7NyQkIjJWLnkoUWMxUlIhIzskIjFFPUg6VDMoKioqISM7NyQkIjIyU2BiMVVuJlIhIzskIjFwSG4nWy1CKSoqISM7NyQkIjFuKEdCXD1XKFIhIzokIjEoM2BPLlldJioqISM7NyQkIjFMVDU+XDQjKlIhIzokIjFzLlk9KVtgIioqISM7NyQkIjFsW2xzeFdGUyEjOiQiMWx0P0NnJSkpeiohIzs3JCQiMSlmMGlpK0cxJSEjOiQiMTwtUSRId0xqKiEjOzckJCIyJkcnKnpbR0xFVCEjOyQiMSNROjdZRWRAKiEjOzckJCIxZk9Scl0nKSo9JSEjOiQiMXhLLHBgW1wnKSEjOzckJCIxJnAlKnpbTGJLJSEjOiQiMXU1MFpmJXopcCEjOzckJCIxJ3BnR08scVElISM6JCIwLDIpZVAuZWchIzo3JCQiMSdwRXhCcCVbVyEjOiQiMU0mKlJhMG9PXSEjOzckJCIxeiZbMicpcGNeJSEjOiQiMShwZ0o/aU8kUSEjOzckJCIxaS94JFtxR2UlISM6JCIxJW9ydCN6XGhEISM7NyQkIjFScSdILixoayUhIzokIjJpOCl6KW95PUsiISM8NyQkIjE8TzsjZUokNFohIzokIjFLXTRPJGVaNichIz03JCQiMVYnR15zRGF4JSEjOiQhMiJmLzVeeFJkNyEjPDckJCIxcE80byk+OiVbISM6JCExYmJaOk8rYUQhIzs3JCQiMUA5dnQqUWghXCEjOiQhMVgvNlglXCh5UCEjOzckJCIxdSI0JXohZTIoXCEjOiQhMmo5MkxYby8lXCEjPDckJCIxO2VXJUgzJVFdISM6JCExVChIJmUxM29nISM7NyQkIjFlQ1s0JmVnNSYhIzokITElKkhQJG95WjMoISM7NyQkIjFDOilwcTg3PCYhIzokITFGI28kW1R2VHohIzs3JCQiMGYhWy8qb2pCJiEjOSQhMTwpR2kralNtKSEjOzckJCIxdHFcLyYqKkhJJiEjOiQhMXdPKWZdLzJEKiEjOzckJCIxZE5eLyxqcGAhIzokITFzJlstISo0TG4qISM7NyQkIjEpPjxxYHBFUyYhIzokITFFWCZ5JnlmPikqISM7NyQkIjAlM19wKjNkViYhIzkkITFyaVVoXC1CKiohIzs3JCQiMG1zZG9HQVgmISM5JCEvTzpdc15lKiohIzk3JCQiMSJbQz9TWyhvYSEjOiQhMSpcaiVwKFJKKSoqISM7NyQkIjEtakY9Im9fWyYhIzokITFjeXZWYydvKioqISM7NyQkIjFBIkdYJHl5LGIhIzokITFVKnk2Kil6JyoqKiohIzs3JCQiMSc9IzRReCdwXiYhIzokITBfMSQpKkhrIyoqKiEjOjckJCIwRGM7a1pAYCYhIzkkITEmKkcwOmtSdyoqISM7NyQkIjE5LkFYdktaYiEjOiQhMEpNXTZiNCYqKiEjOjckJCIxeVZ5W3VdaWIhIzokITEsNTpZRE07KiohIzs3JCQiMTBEImZEbkdmJiEjOiQhMVVsTkYnUig+KSohIzs3JCQiMUwxL2pxQUJjISM6JCExMCQ0PXJWcG8qISM7NyQkIjJXRCQ9JGY9R3AmISM7JCExWyE9aldeKVsjKiEjOzckJCIxeGVLQixUaWQhIzokITFQODkhPSEpPWopISM7NyQkIjEjNCNcPE9sQ2UhIzokITE2eSg0RyI+UXohIzs3JCQiMTIkZTs2KCpvKWUhIzokITFpOiJmJEhrQHIhIzs3JCQiMCg+cWU7RWBmISM5JCExVVQxMUAnKkhoISM7NyQkIjFMY3UwaWk+ZyEjOiQhMSxHVDAiWy8uJiEjOzckJCIxOj0+WGI5JDMnISM6JCEyWCE9NHZlIlwqUSEjPDckJCIxKip6aiUpW21ZaCEjOiQhMVkhbydwKDNtcCMhIzs3JCQiMSoqKkcnKjNEXEAnISM6JCExaCJcLWNuNE8iISM7NyQkIi1pJUgmPSRHJyEjNiQhMlBHbDRJPD5eIyEjQy0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiM1ISIiLSUlVklFV0c2JDskIiIhISIiJCIrM2A9JEcnISIqJShERUZBVUxURy0lK0FYRVNMQUJFTFNHNictSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1USJ4NiIvJSdmYW1pbHlHUSE2Ii8lJXNpemVHUSMxMDYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiJRITYiLSUlRk9OVEc2JSUoREVGQVVMVEclKERFRkFVTFRHIiM1JStIT1JJWk9OVEFMRyUrSE9SSVpPTlRBTEctJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiRJJSEiIi0lKUJPVU5EU19ZRzYjJCIkKyMhIiItJS1CT1VORFNfV0lEVEhHNiMkIiU/TiEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiUhZiQhIiItJSlDSElMRFJFTkc2Ig==
We expect there may be a problem with this system, let's see what happens. Notice our system is very 'sparse' to begin with.
Sys2:=seq(simplify(subs({x=pts2[j][1],y=pts2[j][2]},F=y),trig),j=1..4);
NiYvLCZJIkJHNiIiIiJJIkVHRiZGJ0YnLywmSSJBR0YmRidGKCEiIiIiIy8sJkYoRidGJUYsISIjLywmRitGLEYoRixGJw==
solve({Sys2},{A,B,C,E});
Sys2Matrix:=matrix(4,5,[0,1,0,1,1,1,0,0,-1,2,0,-1,0,1,-2,-1,0,0,-1,1]);
PTYiNiQ7IiIiIiIlO0YmIiImRVxbbDU2JEYmIiIjRiY2JEYnRiwiIiE2JEYsRixGLjYkRixGJkYmNiQiIiRGLCEiIjYkRidGKUYmNiRGJkYpRiY2JEYnRiZGMzYkRixGKUYsNiRGJ0YyRi42JEYmRiZGLjYkRjJGMkYuNiRGMkYnRiY2JEYmRidGJjYkRjJGKSEiIzYkRiZGMkYuNiRGLEYyRi42JEYyRiZGLjYkRixGJ0YzNiRGJ0YnRjM=
The rref matrix of the augmented matrix corresponding to our system. We have serious problems! This means there is NO linear combination of the four given functions which will pass through all four of the points.
rref(Sys2Matrix);
PTYiNiQ7IiIiIiIlO0YmIiImRVxbbDU2JEYmIiIjIiIhNiRGJ0YsRi02JEYsRixGJjYkRixGJkYtNiQiIiRGLEYtNiRGJ0YpRiY2JEYmRilGLTYkRidGJkYtNiRGLEYpRi02JEYnRjJGLTYkRiZGJkYmNiRGMkYyRi02JEYyRidGJjYkRiZGJ0YtNiRGMkYpRi02JEYmRjJGLTYkRixGMkYtNiRGMkYmRi02JEYsRidGLTYkRidGJ0Yt
We modify the set of points only slightly. We change the x-coordinate of the third point to pi/6 instead of pi, which will yield 2 more non-zero entries in the augmented matrix.
pts3:=[[0,1],[Pi/2,2],[Pi/6,-2],[3*Pi/2,1]];
NyY3JCIiISIiIjckLCRJI1BpRyUqcHJvdGVjdGVkRyNGJSIiI0YrNyQsJEYoI0YlIiInISIjNyQsJEYoIyIiJEYrRiU=
Sys3:=seq(simplify(subs({x=pts3[j][1],y=pts3[j][2]},F=y),trig),j=1..4);
NiYvLCZJIkJHNiIiIiJJIkVHRiZGJ0YnLywmSSJBR0YmRidGKCEiIiIiIy8sKkYrI0YnRi0qJkkiQ0dGJkYnIiIkRjBGMEYoRjAqJkYlRidGM0YwRjAhIiMvLCZGK0YsRihGLEYn
soln3:=solve({Sys3},{A,B,C,E});
PCYvSSJBRzYiIyIiIiIiIy9JIkJHRiUjIiImRigvSSJDR0YlLCYqJCIiJEYmISIiIyEiJkYoRicvSSJFR0YlIyEiJEYo
FSoln3:=subs(soln3,F);
LCotSSRzaW5HNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHNiI2I0kieEdGKCMiIiIiIiMtSSRjb3NHRiVGKSMiIiZGLSomLCYqJCIiJEYrISIiIyEiJkYtRixGLC1GJDYjLCRGKkYtRixGLC1GL0Y6IyEiJEYt
plot5:=pointplot(pts3,symbol=circle,symbolsize=30):
plot6:=plot(FSoln3,x=-Pi/2..2*Pi,color=blue):
display({plot5,plot6});
NictJSdQT0lOVFNHNic3JCQiIiEhIiIkIiM1ISIiNyQkIitGanpxOiEiKiQiIz8hIiI3JCQiK2V4KWZCJiEjNSQhIz8hIiI3JCQiKyIpKilRN1ohIiokIiM1ISIiLSUnU1lNQk9MRzYkJSdDSVJDTEVHIiNJLSUnQ1VSVkVTRzYkN2VyNyQkIS1sQmp6cTohIzYkIjFtNFtNKysrNSEjOjckJCEyOD9zJio0Kj4mWyIhIzskIjI5b0xHTDBZIj4hIzs3JCQhMkVTJVx2PWcqUiIhIzskIjJ2Si1PNi1xdyMhIzs3JCQhMD01Jik9Q15LIiEjOSQiMWNdJTRVOEJXJCEjOjckJCEydyZmXyxsa103ISM7JCIxdCQ0UFlKTS8lISM6NyQkITInelM4YCYpKW87IiEjOyQiMCg0NHV2QT1ZISM5NyQkITI5P1VaZ0lKMyIhIzskIjJtOEFnIT5xdV0hIzs3JCQhMUM2KkhhdSIpKSoqISM7JCIxKSo+RXJvOzFhISM6NyQkITFNLWNRSS9YIiohIzskIjFLOm5IcTUvYyEjOjckJCExcChvbSg+RU4qKSEjOyQiMXctXVdbT0tjISM6NyQkITEvdHg5NFtEKCkhIzskIjFgIWZWK2pBbCYhIzo3JCQhMVJlKUcmKSpwOiYpISM7JCIxRCFwYndHUW0mISM6NyQkITF0Vio0ej1mSSkhIzskIjFBViopUUw1bmMhIzo3JCQhMTNINUh4OCc0KSEjOyQiMUFbIioqeVRAbSYhIzo3JCQhMVY5QG5tTicpeSEjOyQiMk1SUWAiSCxcYyEjOzckJCExeCo+YGd2bG4oISM7JCIxMEBbJjQqekZjISM6NyQkITE2JkdNYSV6bXUhIzskIjJsbSQ+QmhmKWYmISM7NyQkITF1V2JdKj0pKW8nISM7JCIxb1ciSEosRFUmISM6NyQkITFPL29kTCUzImYhIzskIjImKj1ib3kqSFheISM7NyQkITI5SVshSGd1KkglISM8JCIxWyMzcSwiUiNIJSEjOjckJCExXkwxW2trbU0hIzskIjIxbnFSOUFDdCQhIzs3JCQhLyV5cSdvYUxFISM5JCIyQmhtZmVoVjYkISM7NyQkITFwYCxaIT5KIT0hIzskIjJFKTR0JykpKWZmQyEjOzckJCEwUUImcEEicHMqISM7JCIyXj0mR2M9Lid5IiEjOzckJCEyXGxNWVRUXD0iISM9JCIyOFw3bFExWTQiISM7NyQkIjAyYS1XSHFOKCEjOyQiMjIhZV0ib2V0RCUhIzw3JCQiMXJKKD0/ITMpWyIhIzskITFAcTInMzVkRyIhIzs3JCQiMlskNHNmdVhTQSEjPCQhMVhHSmItSG9qISM7NyQkIjJXUWNVa2FXJFIhIzwkITE6cWU7MlFiOiEjOjckJCIxND43JTRKXHklISM7JCEyJTNsMCNlalMoPSEjOzckJCIxTXUpUmEyYWomISM7JCEyJ3A8SFFZbSYzIyEjOzckJCIxRHJWb0w/WGchIzskITJjIyopWztRKnk5IyEjOzckJCIxPG8pRz4qKlxYJyEjOyQhMk8mcHVZbiJSPSMhIzs3JCQiMWo7NjByKilmbSEjOyQhMjNkImYvZTMjPiMhIzs3JCQiMTNsTDxdemtvISM7JCExJ1tTMHE0UD4jISM6NyQkIjFgOGNISHBwcSEjOyQhMk49OXAsOykpPSMhIzs3JCQiMSo+J3lUM2Z1cyEjOyQhMXRPWXBXV3hAISM6NyQkIjE9K0Z4JmUpPSEpISM7JCEyLz1gJkgrZSMzIyEjOzckJCIxUFF2N2o3aigpISM7JCEydFxhenB3byE+ISM7NyQkIjEtMk9QRThbJyohIzskITImb29fJj5RM2ciISM7NyQkIjJudic+JypRSmA1ISM7JCEyI1FpQEtIIT0/IiEjOzckJCIxLSR6Wz1kSz8iISM6JCEyJCkqW0lPMEkxTyEjPDckJCIyLj0mXHZvbng4ISM7JCIxRF5BRWhqZHchIzs3JCQiMVclcGYmeTFLOiEjOiQiMjo/JkcjZV1XdyIhIzs3JCQiMW1QP0ojZjlxIiEjOiQiMjBySChcNVs3RiEjOzckJCIxNzJnYSY0QHkiISM6JCIyYyVmLypSWSRvSSEjOzckJCIxZXcqeigpZkYnPSEjOiQiMjFhPihlQSdRTSQhIzs3JCQiMiJmRiVRJlwkWyE+ISM7JCIyMTJXdClHKj1YJCEjOzckJCIwJ3lvSCsicCU+ISM5JCIyM1tKUDpeT2AkISM7NyQkIjFoSGAwXikqKSk+ISM6JCIyJGVGYWtpKnplJCEjOzckJCIxaSF5OD1nNS4jISM6JCIyJFF5cClvJilSaCQhIzs3JCQiMidcJipSMyN6LjAjISM7JCIyM05TeCMqNGloJCEjOzckJCIydy5AYUIpcHA/ISM7JCIyZWAsYlZUQWgkISM7NyQkIjJjX1VDRTwhKjMjISM7JCIyMkhVJEdMLi1PISM7NyQkIjE4U1kqR08kM0AhIzokIjJZXTFITltiZSQhIzs3JCQiMSopcF1WVihwOSMhIzokIjIkKnp3PFBaT2AkISM7NyQkIjJXJypcdlI3Yz0jISM7JCIyTCVIU2liVmNNISM7NyQkIjJ3Q2coPl4nKm9BISM7JCIya0dITChROS9LISM7NyQkIjFKMCg+JXlKX0IhIzokIjJGTXQzNHV2JEchIzs3JCQiMTQ3XCUqeSopUUMhIzokIjIlZUZdcigpb1VCISM7NyQkIjEoKT0sWnpaREQhIzokIjJEJ1x6cS5AVTwhIzs3JCQiMCMpW3M4WTNnIyEjOSQiMmRAbyRSVF9XNiEjOzckJCIxYGRbRlZAd0UhIzokIjI5JUhdTG0+KSlbISM8NyQkIjJMbz94Ljl3diMhIzskITIkR1VXJmZ0eHAjISM8NyQkIjJOaGJ6dTghUkchIzskITI6JypbWnJ2TTEiISM7NyQkIjJQbiZmOHc1QkghIzskITI6KjNLOjdFKz4hIzs3JCQiMU1kQno5PzJJISM6JCEweGtuKUh2S0YhIzk3JCQiMjJDeTAjNFoqMyQhIzskITFbJG8kSCJbP18kISM6NyQkIjJ1dT8+T1M8PCQhIzskITFZdzQvJUdpRSUhIzo3JCQiMSJbTm1GUzhEJCEjOiQhMHJuY1czXyNcISM5NyQkIjE6LU4iPlM0TCQhIzokITEjNHFHUHR4XSYhIzo3JCQiMSopcEFEQks+TSEjOiQhMShSR1tueHUvJyEjOjckJCIyTXcuImZXcTJOISM7JCExbXFqJ3p2IWVrISM6NyQkIjJzK0QjKTQ3dWEkISM7JCExYXhUNnp6J2YnISM6NyQkIjI7RFl0dD5yZSQhIzskITIxayg+bXpzMG4hIzs3JCQiMSdcbmtQRm9pJCEjOiQhMSwyalA1MiV5JyEjOjckJCIyKFIoKWU6XWBtTyEjOyQhMiY0b1QqZTE3JG8hIzs3JCQiMVFFbzBLdChvJCEjOiQhMWk7MCJ5S00lbyEjOjckJCIybWB3ZFJKKjNQISM7JCExYDF0SiYqZVlvISM6NyQkIjFOLyhlZUgsdCQhIzokITFOJXlCYlgxJW8hIzo3JCQiMk5Ma2Z4Rjh2JCEjOyQhMChmJypwV2VEbyEjOTckJCIyLzhfaDpDUHokISM7JCEwPUp0ST8ib24hIzk3JCQiMnQjKlJqYD9oJFEhIzskITEiUU5XSWxWbichIzo3JCQiMXpCIyp6YCZIIlIhIzokITFLQU0kPVdNVCchIzo3JCQiMUpbXUIteiopUiEjOiQhMXcuSnJpPlJnISM6NyQkIjA4I0ctNXp0UyEjOSQhMSVSLiJIWlwzYiEjOjckJCIwVmY1eSJ6ZFQhIzkkITFBJUdMSDRFJ1shIzo3JCQiMSxeYm4qSG9CJSEjOiQhMUIiKjQxSjZqVCEjOjckJCIxdDIwYSJvZUolISM6JCEyRShmJkczNCYqUSQhIzs3JCQiMS5xdktlWylSJSEjOiQhMid5a1skSDIkPkQhIzs3JCQiMUxLWTZONSJbJSEjOiQhMiRHY15tLyZmZyIhIzs3JCQiMXNgYCRSeD1jJSEjOiQhMShHKHBXSVw5cCEjOzckJCIxOHZndjdsVVkhIzokIjI9Vyl6MCczNEMjISM8NyQkIjI6RS5XLzlzcyUhIzskIjImZWRDNiU+PjsiISM7NyQkIjImNCEqPjhveDZbISM7JCIyL20oKSkzPHBjPyEjOzckJCIyMHpzKyIzQSQqWyEjOyQiMkxUY05bbSJlRyEjOzckJCIxcmwlcCFbbXVcISM6JCIyMHNtXV1PSmUkISM7NyQkIjFBcidwMGB6MCYhIzokIjFbbWZFZHFIVSEjOjckJCIxdHcpcElUNzkmISM6JCIxODRKQTkqcHclISM6NyQkIjFDbkMpKSlSUUEmISM6JCIxbSEpcDxJXSI9JiEjOjckJCIxdmRdcCVRa0kmISM6JCIwKj1BQ2srcmEhIzk3JCQiMjAiZlRHIylRV2AhIzskIjFvLEc+QiM0YyYhIzo3JCQiMmIvRXQpekwjUSYhIzskIjFrJipmJ3pgTWkmISM6NyQkIjE4NnltR0osYSEjOiQiLzBXQCRcV2smISM4NyQkIjIwPU9pdShHP2EhIzskIjFablM3amhlYyEjOjckJCIxWzdwREVFUmEhIzokIi8iZUQ+ISlmbSYhIzg3JCQiMTtqOTB2QmVhISM6JCIxIylROChwdmxtJiEjOjckJCIxc2VKZVspKnphISM6JCIxb3glXEsiKiplYyEjOjckJCIxRmFbNkF0LGIhIzokIjFBNk4lR2JFayYhIzo3JCQiMSMpXGxrJnpNXyYhIzokIjFqJkdwQGt3aCYhIzo3JCQiMVJYI3kicEFYYiEjOiQiMUUiKnohPUlUZSYhIzo3JCQiMGxqVGlAKCllJiEjOSQiMXoydkxxJj5cJiEjOjckJCIxaEZdSWpAS2MhIzokIjFIKFImPTJMbmAhIzo3JCQiMFs1IylwPytyJiEjOSQiMiZHbSJcWGshb10hIzs3JCQiMSo+PWYxRHl5JiEjOiQiMS5LJj0nXFl6WSEjOjckJCIxREZzXDJ5cWUhIzokIjEmcHVtMlUrPSUhIzo3JCQiMV5zX0xrdGBmISM6JCIyJCl5XCp6WiUpM08hIzs3JCQiMUhcJHlnTkouJyEjOiQiMm5QS0BIVzosJCEjOzckJCIxM0U5I3lNRDYnISM6JCIxVTRAeSY9PVEjISM6NyQkIjEvOFFRKyd5PichIzokIjJvKio0LCZSeSlvIiEjOzckJCItaSVIJj0kRychIzYkIjJLI2UtKywrKzUhIzstJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIjNSEiIi0lJVZJRVdHNiQ7JCErRmp6cTohIiokIiszYD0kRychIiolKERFRkFVTFRHLSUrQVhFU0xBQkVMU0c2Jy1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRIng2Ii8lJ2ZhbWlseUdRITYiLyUlc2l6ZUdRIzEwNiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2IlEhNiItJSVGT05URzYlJShERUZBVUxURyUoREVGQVVMVEciIzUlK0hPUklaT05UQUxHJStIT1JJWk9OVEFMRy0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiI3EhIiItJSlCT1VORFNfWUc2IyQiI10hIiItJS1CT1VORFNfV0lEVEhHNiMkIiUhKVEhIiItJS5CT1VORFNfSEVJR0hURzYjJCIlISlRISIiLSUpQ0hJTERSRU5HNiI=
We try this again, but we remove only one zero from the matrix (replace pi/6 by pi/4 in point 3's x-coordinate). Will we still get a solution?
pts4:=[[0,1],[Pi/2,2],[Pi/4,-2],[3*Pi/2,1]];
NyY3JCIiISIiIjckLCRJI1BpRyUqcHJvdGVjdGVkRyNGJSIiI0YrNyQsJEYoI0YlIiIlISIjNyQsJEYoIyIiJEYrRiU=
Sys4:=seq(simplify(subs({x=pts4[j][1],y=pts4[j][2]},F=y),trig),j=1..4);
NiYvLCZJIkJHNiIiIiJJIkVHRiZGJ0YnLywmSSJBR0YmRidGKCEiIiIiIy8sJComRi0jRidGLSwoRitGJyomSSJDR0YmRidGLUYxRidGJUYnRidGMSEiIy8sJkYrRixGKEYsRic=
soln4:=solve({Sys4},{A,B,C,E});
PCYvSSJBRzYiIyIiIiIiIy9JIkJHRiUjIiImRigvSSJDR0YlLCYqJEYoRiYjISIkRighIiNGJy9JIkVHRiVGMQ==
The answer is yes!!! Will that always be the answer. Definitely not, but in this case it works.
FSoln4:=subs(soln4,F);
LCotSSRzaW5HNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHNiI2I0kieEdGKCMiIiIiIiMtSSRjb3NHRiVGKSMiIiZGLSomLCYqJEYtRisjISIkRi0hIiNGLEYsLUYkNiMsJEYqRi1GLEYsLUYvRjlGNQ==
plot7:=pointplot(pts4,symbol=circle,symbolsize=30):
plot8:=plot(FSoln4,x=-Pi/2..2*Pi,color=green):
display({plot7,plot8});
NictJSdQT0lOVFNHNic3JCQiIiEhIiIkIiM1ISIiNyQkIitGanpxOiEiKiQiIz8hIiI3JCQiK047KVImeSEjNSQhIz8hIiI3JCQiKyIpKilRN1ohIiokIiM1ISIiLSUnU1lNQk9MRzYkJSdDSVJDTEVHIiNJLSUnQ1VSVkVTRzYkN2ZyNyQkIS1sQmp6cTohIzYkIjJYXCV5TCsrKzUhIzs3JCQhMjg/cyYqNCo+JlsiISM7JCIxeFFKXTh1Jio9ISM6NyQkITJFUyVcdj1nKlIiISM7JCIxVUdDZGMjKUhGISM6NyQkITA9NSYpPUNeSyIhIzkkIjFkR3cyJnArUiQhIzo3JCQhMncmZl8sbGtdNyEjOyQiMXknZlc3eHMoUiEjOjckJCEyJ3pTOGAmKSlvOyIhIzskIjE3NjEhKlE+UVghIzo3JCQhMjk/VVpnSUozIiEjOyQiMSdmNzpdSEkpXCEjOjckJCExQzYqSGF1IikpKiohIzskIjAuW0tDcmBJJiEjOTckJCExTS1jUUkvWCIqISM7JCIxWzlqKlxacVwmISM6NyQkITFwKG9tKD5FTiopISM7JCIwNzhRYjhVXyYhIzk3JCQhMS90eDk0W0QoKSEjOyQiL15RMS5AVmIhIzg3JCQhMVJlKUcmKSpwOiYpISM7JCIxY1c/K20xYWIhIzo3JCQhMXRWKjR6PWZJKSEjOyQiMUdXdSUqWyNvYiYhIzo3JCQhMTNINUh4OCc0KSEjOyQiMSNScFU7VDpiJiEjOjckJCExVjlAbm1OJyl5ISM7JCIxcnhhc1pHUWIhIzo3JCQhMXgqPmBndmxuKCEjOyQiMWleTEMkUXJeJiEjOjckJCExNiZHTWElem11ISM7JCIweCtSXCg+KVsmISM5NyQkITF1V2JdKj0pKW8nISM7JCIyRGxpLiRRdzlgISM7NyQkITFPL29kTCUzImYhIzskIjJDJXpQK2wjRy8mISM7NyQkITI5SVshSGd1KkglISM8JCIxJlEjNEcleSUzVSEjOjckJCExXkwxW2trbU0hIzskIjJFMTtbOmE7bSQhIzs3JCQhLyV5cSdvYUxFISM5JCIxOiEqR0kkKXBlSSEjOjckJCExcGAsWiE+SiE9ISM7JCIxJjQ+I2VuXz9DISM6NyQkITBRQiZwQSJwcyohIzskIjFBKDM6O0VZdyIhIzo3JCQhMlxsTVlUVFw9IiEjPSQiMi0hMyopXEMpPjQiISM7NyQkIjAyYS1XSHFOKCEjOyQiMllEaTQnM3E+VyEjPDckJCIxckooPT8hMylbIiEjOyQhMU04ZmMpPitoKiEjPDckJCIyWyQ0c2Z1WFNBISM8JCExJG9BPXBhJikpZSEjOzckJCIyV1FjVWthVyRSISM8JCEyWChHJTRtbHBaIiEjOzckJCIxND43JTRKXHklISM7JCEyPT9SWyNmYSR5IiEjOzckJCIxTXUpUmEyYWomISM7JCExMGY1bHBsJik+ISM6NyQkIjFEclZvTD9YZyEjOyQhMlZzTy08SVYvIyEjOzckJCIxPG8pRz4qKlxYJyEjOyQhMkU6bW1TI1x4PyEjOzckJCIxajs2MHIqKWZtISM7JCEyRCZvMk4hKVwlMyMhIzs3JCQiMTNsTDxdemtvISM7JCEydmltNzpSXjMjISM7NyQkIjFgOGNISHBwcSEjOyQhMmp5S3gpXFd6PyEjOzckJCIxKj4neVQzZnVzISM7JCExIjQwVmdjdTEjISM6NyQkIjE9K0Z4JmUpPSEpISM7JCEyMGtGSHc0Pig+ISM7NyQkIjFQUXY3ajdqKCkhIzskITJFTChSbzsoenoiISM7NyQkIjEtMk9QRThbJyohIzskITJqJillZk1nclwiISM7NyQkIjJudic+JypRSmA1ISM7JCEydlYjKTQnPmYxNiEjOzckJCIxLSR6Wz1kSz8iISM6JCExPUJtaz5wakchIzs3JCQiMi49Jlx2b254OCEjOyQiMXM+Mi5jd3UhKSEjOzckJCIxVyVwZiZ5MUs6ISM6JCIyZG8yPCYpPUl4IiEjOzckJCIxbVA/SiNmOXEiISM6JCIyamdtRmJzUW8jISM7NyQkIjE3MmdhJjRAeSIhIzokIjJkTCFcKVxISC0kISM7NyQkIjFldyp6KClmRic9ISM6JCIxdyEqNFpnIkdHJCEjOjckJCIyImZGJVEmXCRbIT4hIzskIjFAQGlHeUgkUSQhIzo3JCQiMCd5b0grInAlPiEjOSQiMSh6aioqcCMqelgkISM6NyQkIjIxVDV3Y1p6Jz4hIzskIjJOTi5SbUBgWyQhIzs3JCQiMWhIYDBeKSopKT4hIzokIjImKW9HZiYpNGVdJCEjOzckJCIyOF5iTWtBKywjISM7JCIyJilbVStRXiQ+TiEjOzckJCIxaSF5OD1nNS4jISM6JCIxdUJCa0gmZV8kISM6NyQkIjInXCYqUjMjei4wIyEjOyQiMic9PylmWmBiXyQhIzs3JCQiMncuQGFCKXBwPyEjOyQiMiV6KGYwZSc+Pk4hIzs3JCQiMmNfVUNFPCEqMyMhIzskIjIxQkcmZipRbl0kISM7NyQkIjE4U1kqR08kM0AhIzokIjJ2KTRUPG45KVskISM7NyQkIjEqKXBdVlYocDkjISM6JCIxMkhrNC5aS00hIzo3JCQiMlcnKlx2UjdjPSMhIzskIjFKXGxocDNfTCEjOjckJCIyd0NnKD5eJypvQSEjOyQiMjFSXnhBJTQmNCQhIzs3JCQiMUowKD4leUpfQiEjOiQiMmE7Xl0mcCVvcyMhIzs3JCQiMTQ3XCUqeSopUUMhIzokIjBudHNocE1CIyEjOTckJCIxKCk9LFp6WkREISM6JCIyQjE1bSpld1A7ISM7NyQkIjAjKVtzOFkzZyMhIzkkIjJZWW4hSG55WTUhIzs3JCQiMWBkW0ZWQHdFISM6JCIxejdnOiZlKyslISM7NyQkIjJMbz94Ljl3diMhIzskITI4MFwheT0pcFkkISM8NyQkIjJOaGJ6dTghUkchIzskITJMKltTdiRwazciISM7NyQkIjJQbiZmOHc1QkghIzskITJuSCpvRD43Wj4hIzs3JCQiMU1kQno5PzJJISM6JCExVyhvVkplQHcjISM6NyQkIjIyQ3kwIzRaKjMkISM7JCExKEcpXGksZExOISM6NyQkIjJ1dT8+T1M8PCQhIzskITFXZE80ZmJmVSEjOjckJCIxIltObUZTOEQkISM6JCExWEgjZUspNCxcISM6NyQkIjE6LU4iPlM0TCQhIzokITFVK1QqUk5vWSYhIzo3JCQiMSopcEFEQks+TSEjOiQhMU47NndkMyopZiEjOjckJCIyTXcuImZXcTJOISM7JCExPDpHPChcU1EnISM6NyQkIjJzK0QjKTQ3dWEkISM7JCEyYmYhZVBBWjtsISM7NyQkIjI7RFl0dD5yZSQhIzskITEjeWAmRyE0Jz5tISM6NyQkIjEnXG5rUEZvaSQhIzokITFoU2AlbyxGcCchIzo3JCQiMihSKCllOl1gbU8hIzskITFfMTJsRjtObiEjOjckJCIxUUVvMEt0KG8kISM6JCExRlMiKTMnUl51JyEjOjckJCIybWB3ZFJKKjNQISM7JCEwYixEbUJpdSchIzk3JCQiMU4vKGVlSCx0JCEjOiQhMiYpUnY/TyhRUW4hIzs3JCQiMk5Ma2Z4Rjh2JCEjOyQhMXVgTk54aEBuISM6NyQkIjIvOF9oOkNQeiQhIzskITEyJXkzdypIaG0hIzo3JCQiMnQjKlJqYD9oJFEhIzskITFBLncwJmVhYychIzo3JCQiMXpCIyp6YCZIIlIhIzokITF5LSRHTWRGSSchIzo3JCQiMUpbXUIteiopUiEjOiQhMlh2OShcIVEkSGYhIzs3JCQiMDgjRy01enRTISM5JCEwJ2VQR0VdLWEhIzk3JCQiMFZmNXkiemRUISM5JCExMSEpR2d5WWpaISM6NyQkIjEsXmJuKkhvQiUhIzokITF5TzB1NyhIMiUhIzo3JCQiMXQyMGEib2VKJSEjOiQhMjIoKW8+azkxSiQhIzs3JCQiMS5xdktlWylSJSEjOiQhMmstK2VYblVYIyEjOzckJCIxTEtZNk41IlslISM6JCEyXnVKLCdwYGM6ISM7NyQkIjFzYGAkUng9YyUhIzokITEwX0lkKnpoZSchIzs3JCQiMTh2Z3Y3bFVZISM6JCIyTHU4WzddW1IjISM8NyQkIjI6RS5XLzlzcyUhIzskIjJgMlgpPm5qZTYhIzs3JCQiMiY0ISo+OG94NlshIzskIjFSZihvJWUjWy4jISM6NyQkIjIwenMrIjNBJCpbISM7JCIwVnklUm0pKj1HISM5NyQkIjFybCVwIVttdVwhIzokIjJVY1c7Iyp6d18kISM7NyQkIjFBcidwMGB6MCYhIzokIjEyKGYicGA3ZlQhIzo3JCQiMXR3KXBJVDc5JiEjOiQiMmIkKlI3QFRLbyUhIzs3JCQiMUNuQykpKVJRQSYhIzokIjFzKioqMyQ0KXAzJiEjOjckJCIxdmRdcCVRa0kmISM6JCIxaU0+PWhHb2AhIzo3JCQiMjAiZlRHIylRV2AhIzskIjBeVF1OaWBYJiEjOTckJCIyYi9FdCl6TCNRJiEjOyQiMUo1WCYpPm06YiEjOjckJCIxODZ5bUdKLGEhIzokIjFHcSYzbnRkYCYhIzo3JCQiMjA9T2l1KEc/YSEjOyQiMUMuYztMQFxiISM6NyQkIjFbN3BERUVSYSEjOiQiMWJoKyRSMmdiJiEjOjckJCIxO2o5MHZCZWEhIzokIjFjdyZmWCI+Y2IhIzo3JCQiMXNlSmVbKSp6YSEjOiQiMURZQVI1TFtiISM6NyQkIjFGYVs2QXQsYiEjOiQiMW9VL24jRz5gJiEjOjckJCIxIylcbGsmek1fJiEjOiQiMiYzTUEyKnpxXSYhIzs3JCQiMVJYI3kicEFYYiEjOiQiMTE0Jnp1KCpRWiYhIzo3JCQiMGxqVGlAKCllJiEjOSQiMUg3cjBHMCRRJiEjOjckJCIxaEZdSWpAS2MhIzokIjE0KGY3RnkwRSYhIzo3JCQiMFs1IylwPytyJiEjOSQiMW47byhSZ3InXCEjOjckJCIxKj49ZjFEeXkmISM6JCIxcz1Bbnclb2UlISM6NyQkIjFec19Ma3RgZiEjOiQiMil6ZWFHMDBUTiEjOzckJCIxSFwkeWdOSi4nISM6JCIxUHInNGlaJWVIISM6NyQkIjEzRTkjeU1ENichIzokIjE6LkY1Y3ZXQiEjOjckJCIxLzhRUSsneT4nISM6JCIxMms3PSR6KnA7ISM6NyQkIi1pJUgmPSRHJyEjNiQiMmFPV3M0KysrIiEjOy0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIzUhIiIkIiIhISIiLSUlVklFV0c2JDskIStGanpxOiEiKiQiKzNgPSRHJyEiKiUoREVGQVVMVEctJStBWEVTTEFCRUxTRzYnLUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NVEieDYiLyUnZmFtaWx5R1EhNiIvJSVzaXplR1EjMTA2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EldHJ1ZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ2l0YWxpYzYiUSE2Ii0lJUZPTlRHNiUlKERFRkFVTFRHJShERUZBVUxURyIjNSUrSE9SSVpPTlRBTEclK0hPUklaT05UQUxHLSUlUk9PVEc2Jy0lKUJPVU5EU19YRzYjJCIjcSEiIi0lKUJPVU5EU19ZRzYjJCIjXSEiIi0lLUJPVU5EU19XSURUSEc2IyQiJSEpUSEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiUhKVEhIiItJSlDSElMRFJFTkc2Ig==