MAT331,
Feb 16th, 2010
If a curve fits, wear it (Chapter 2 of the online book)
EXERCISE: Is it possible to find a linear function passing through [1,-1/2], [2,1] and [3,1/3]? Why?
pts:=[ [1,-1/2], [2,1] , [3,1/3]];
NyU3JCIiIiMhIiIiIiM3JEYnRiQ3JCIiJCNGJEYq
with(plots):
graph1:=pointplot(pts,color=blue):
graph1;
LSUlUExPVEc2JS0lJ1BPSU5UU0c2JTckJCIjNSEiIiQhIiYhIiI3JCQiIz8hIiIkIiM1ISIiNyQkIiNJISIiJCIrTExMTEwhIzUtJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIjNSEiIi0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJCFSISIiLSUpQk9VTkRTX1lHNiMkIiQ/IiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJTVOISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJV1QISIiLSUpQ0hJTERSRU5HNiI=
It is not possible because these points are not colinear.
EXERCISE: Find the linear function passing closest to the points [1,-1/2], [2,1] and [3,1/3]
display(plot(x/3,x=0..4),graph1);
LSUlUExPVEc2Jy0lJ0NVUlZFU0c2JDdTNyQkIiIhISIiJCIiISEiIjckJCIxbW1tbTsnKT0oKSEjPCQiMUFBQUFzRzFIISM8NyQkIjJNTEwkZSc0MGoiISM8JCIyWVdXVz5LXVYmISM9NyQkIjJubW1tNmhPWyMhIzwkIjEpKSkpKSkpUXEpeSMpISM8NyQkIjJtbW1tInlZVUwhIzwkIjIqKSkpKSkpUWY6OTYhIzw3JCQiMUxMTGVGPig+JSEjOyQiMnZ4eEZEayEqUiIhIzw3JCQiMmttbTs+SycqKVwhIzwkIjJiYmIwdDVLbSIhIzw3JCQiMSoqKioqXEtkLCJlISM7JCIxTExMMyI+biQ+ISM7NyQkIjFtbW0iZlgoZW0hIzskIjIoKSkpKSlRJz1lPkEhIzw3JCQiMSoqKioqXFU3WV0oISM7JCIyKCoqKioqXFpQOl0jISM8NyQkIjFNTExML3B1JCkhIzskIjJ3eHh4Wmo6eiMhIzw3JCQiMW5tbTtjMFQiKiEjOyQiMSopKSkpKVEmPXEvJCEjOzckJCIqOCFRKzUhIikkIjFMTExML2dNTCEjOzckJCIqJiozcTMiISIpJCIybW1tbTtqTGkkISM8NyQkIiooPVxxNiEiKSQiMmttbW1jUjshUiEjPDckJCIybm1tImZCSVk3ISM7JCIxQUFBKD5UVjolISM7NyQkIjJMTExMT1trTCIhIzskIjF4eHh4eSNbWCUhIzs3JCQiMktMTEwmUSJHVCIhIzskIjIuNjY2XnokNFohIzw3JCQiLEQyWDtdIiEjNSQiMSoqKioqXGQkWzBdISM7NyQkIjJKTExMdnYtZSIhIzskIjFWV1dXZWVuXyEjOzckJCIsRDJZbG0iISM1JCIrdk46YmIhIzU3JCQiLHYiZXBbPCEjNSQiMUxMTGVnKSpHZSEjOzckJCIyTUxMJGUvVE09ISM7JCIxeHh4RjpxOWghIzs3JCQiMkxMTGVEQkoiPiEjOyQiMVdXVz52MnhqISM7NyQkIjJtbW07a0QhKSo+ISM7JCIxQEFBc2EzZ20hIzs3JCQiMmptbSJmYEAnMyMhIzskIjFAQUEoPl5TJnAhIzs3JCQiLHZ3JSlIOyMhIzUkIjEqKioqKlxBXCo0cyEjOzckJCIybG1tOyR5KmVDIyEjOyQiMWFiYjBoSydbKCEjOzckJCIqOWI6TCMhIikkIjFsbW1tLyY9eCghIzs3JCQiMicqKioqKlw1YWBUIyEjOyQiMSkqKioqKlwuPV4hKSEjOzckJCIyKCoqKipcN1JWJ1wjISM7JCIxKioqKipcUFk5SykhIzs3JCQiMicqKioqKlxAZmtlIyEjOyQiMWxtbTsyYEAnKSEjOzckJCIxTExMYDRObkUhIzokIjFWV1dXKXA2KikpISM7NyQkIio6P1B2IyEiKSQiMW1tbW1yMXoiKiEjOzckJCIybG1tInpNKT4kRyEjOyQiMWJiYklcJSpSJSohIzs3JCQiKihmYTxIISIpJCIxS0xMTEs6RCgqISM7NyQkIjFMTCRlZ2AhKSpIISM6JCIxVldXPj9eJCoqKiEjOzckJCIyKSoqKipcI0cyQTMkISM7JCIyKSoqKipcRkNTRjUhIzs3JCQiMUxMTCQpR1trSiEjOiQiMnh4eHhpRlswIiEjOzckJCIsRCJ5aF1LISM1JCIySkxMM0ZSTjMiISM7NyQkIjJvbW1tKWZkTEwhIzskIjJjYmJiKj4+NjYhIzs3JCQiMm1tbTtxNyU9TSEjOyQiMmFiYjAhNFpSNiEjOzckJCIyTExMZSNwYS1OISM7JCIyVldXPms6djsiISM7NyQkIiphZCl6TiEiKSQiMmxtbW0lZUckPiIhIzs3JCQiMkxMTCRHVVlvTyEjOyQiMlZXVyU0OSNHQSIhIzs3JCQiMmxtbW01Onh1JCEjOyQiMikpKSkpKSlvJFEjXDchIzs3JCQiLEQyOEEkUSEjNSQiLHZOL3VGIiEjNTckJCIsdlMpMzhSISM1JCIsRCFHTy84ISM1NyQkIiNTISIiJCIyTExMTExMTEwiISM7LSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJSdQT0lOVFNHNiY3JCQiIzUhIiIkISImISIiNyQkIiM/ISIiJCIjNSEiIjckJCIjSSEiIiQiK0xMTExMISM1LSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIzUhIiItJSVWSUVXRzYkOyQiIiEhIiIkIiNTISIiJShERUZBVUxURy0lK0FYRVNMQUJFTFNHNiQtSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1USJ4NiIvJSdmYW1pbHlHUSE2Ii8lJXNpemVHUSMxMDYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiJRITYiLSUlUk9PVEc2Jy0lKUJPVU5EU19YRzYjJCIkPyYhIiItJSlCT1VORFNfWUc2IyQiI10hIiItJS1CT1VORFNfV0lEVEhHNiMkIiUrTSEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiVTUSEiIi0lKUNISUxEUkVORzYi
pt1:=[(pts[1][1]+pts[2][1])/2,(pts[1][2]+pts[2][2])/2 ];
pt2:=[(pts[3][1]+pts[2][1])/2,(pts[3][2]+pts[2][2])/2 ];
line:=(x-pt1[1])*(pt2[2]-pt1[2])/(pt2[1]-pt1[1])+pt1[2];
NyQjIiIkIiIjIyIiIiIiJQ==
NyQjIiImIiIjI0YlIiIk
LCZJInhHNiIjIiImIiM3IyEiJCIiKSIiIg==
display(plot([x/3,line],x=0..4,thickness=3),graph1,plot([pt1,pt2],style=point));
LSUlUExPVEc2Ky0lJ0NVUlZFU0c2JTdTNyQkIiIhISIiJCIiISEiIjckJCIxbW1tbTsnKT0oKSEjPCQiMUFBQUFzRzFIISM8NyQkIjJNTEwkZSc0MGoiISM8JCIyWVdXVz5LXVYmISM9NyQkIjJubW1tNmhPWyMhIzwkIjEpKSkpKSkpUXEpeSMpISM8NyQkIjJtbW1tInlZVUwhIzwkIjIqKSkpKSkpUWY6OTYhIzw3JCQiMUxMTGVGPig+JSEjOyQiMnZ4eEZEayEqUiIhIzw3JCQiMmttbTs+SycqKVwhIzwkIjJiYmIwdDVLbSIhIzw3JCQiMSoqKioqXEtkLCJlISM7JCIxTExMMyI+biQ+ISM7NyQkIjFtbW0iZlgoZW0hIzskIjIoKSkpKSlRJz1lPkEhIzw3JCQiMSoqKioqXFU3WV0oISM7JCIyKCoqKioqXFpQOl0jISM8NyQkIjFNTExML3B1JCkhIzskIjJ3eHh4Wmo6eiMhIzw3JCQiMW5tbTtjMFQiKiEjOyQiMSopKSkpKVEmPXEvJCEjOzckJCIqOCFRKzUhIikkIjFMTExML2dNTCEjOzckJCIqJiozcTMiISIpJCIybW1tbTtqTGkkISM8NyQkIiooPVxxNiEiKSQiMmttbW1jUjshUiEjPDckJCIybm1tImZCSVk3ISM7JCIxQUFBKD5UVjolISM7NyQkIjJMTExMT1trTCIhIzskIjF4eHh4eSNbWCUhIzs3JCQiMktMTEwmUSJHVCIhIzskIjIuNjY2XnokNFohIzw3JCQiLEQyWDtdIiEjNSQiMSoqKioqXGQkWzBdISM7NyQkIjJKTExMdnYtZSIhIzskIjFWV1dXZWVuXyEjOzckJCIsRDJZbG0iISM1JCIrdk46YmIhIzU3JCQiLHYiZXBbPCEjNSQiMUxMTGVnKSpHZSEjOzckJCIyTUxMJGUvVE09ISM7JCIxeHh4RjpxOWghIzs3JCQiMkxMTGVEQkoiPiEjOyQiMVdXVz52MnhqISM7NyQkIjJtbW07a0QhKSo+ISM7JCIxQEFBc2EzZ20hIzs3JCQiMmptbSJmYEAnMyMhIzskIjFAQUEoPl5TJnAhIzs3JCQiLHZ3JSlIOyMhIzUkIjEqKioqKlxBXCo0cyEjOzckJCIybG1tOyR5KmVDIyEjOyQiMWFiYjBoSydbKCEjOzckJCIqOWI6TCMhIikkIjFsbW1tLyY9eCghIzs3JCQiMicqKioqKlw1YWBUIyEjOyQiMSkqKioqKlwuPV4hKSEjOzckJCIyKCoqKipcN1JWJ1wjISM7JCIxKioqKipcUFk5SykhIzs3JCQiMicqKioqKlxAZmtlIyEjOyQiMWxtbTsyYEAnKSEjOzckJCIxTExMYDRObkUhIzokIjFWV1dXKXA2KikpISM7NyQkIio6P1B2IyEiKSQiMW1tbW1yMXoiKiEjOzckJCIybG1tInpNKT4kRyEjOyQiMWJiYklcJSpSJSohIzs3JCQiKihmYTxIISIpJCIxS0xMTEs6RCgqISM7NyQkIjFMTCRlZ2AhKSpIISM6JCIxVldXPj9eJCoqKiEjOzckJCIyKSoqKipcI0cyQTMkISM7JCIyKSoqKipcRkNTRjUhIzs3JCQiMUxMTCQpR1trSiEjOiQiMnh4eHhpRlswIiEjOzckJCIsRCJ5aF1LISM1JCIySkxMM0ZSTjMiISM7NyQkIjJvbW1tKWZkTEwhIzskIjJjYmJiKj4+NjYhIzs3JCQiMm1tbTtxNyU9TSEjOyQiMmFiYjAhNFpSNiEjOzckJCIyTExMZSNwYS1OISM7JCIyVldXPms6djsiISM7NyQkIiphZCl6TiEiKSQiMmxtbW0lZUckPiIhIzs3JCQiMkxMTCRHVVlvTyEjOyQiMlZXVyU0OSNHQSIhIzs3JCQiMmxtbW01Onh1JCEjOyQiMikpKSkpKSlvJFEjXDchIzs3JCQiLEQyOEEkUSEjNSQiLHZOL3VGIiEjNTckJCIsdlMpMzhSISM1JCIsRCFHTy84ISM1NyQkIiNTISIiJCIyTExMTExMTEwiISM7LSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJSpUSElDS05FU1NHNiMiIiQtJSdDVVJWRVNHNiU3UzckJCIiISEiIiQhJHYkISIkNyQkIjFtbW1tOycpPSgpISM8JCExQUFBKDQ5blEkISM7NyQkIjJNTEwkZSc0MGoiISM8JCEyWlclcHY0aXFJISM8NyQkIjJubW1tNmhPWyMhIzwkITIoKSkpKVEsNzk6RiEjPDckJCIybW1tbSJ5WVVMISM8JCExKikpKVF3XUlkQiEjOzckJCIxTExMZUY+KD4lISM7JCExeXgtJW9wNisjISM7NyQkIjJrbW07PksnKilcISM8JCEyZGIwb2UnKTRuIiEjPDckJCIxKioqKipcS2QsImUhIzskITJNTCRlOTY1SDghIzw3JCQiMW1tbSJmWChlbSEjOyQhMSMqKSlRLG5BYigqISM8NyQkIjEqKioqKlxVN1ldKCEjOyQhMS8rXWk6eUlpISM8NyQkIjFNTExML3B1JCkhIzskITFzeHhGbFgwRSEjPDckJCIxbm1tO2MwVCIqISM7JCIxTTY2aHRKeGUhIz03JCQiKjghUSs1ISIpJCIxbG1tO2FdIz0lISM8NyQkIiomKjNxMyIhIikkIjFLTEwkZVI/eighIzw3JCQiKig9XHE2ISIpJCIyTkxMJGUlXHE3IiEjPDckJCIybm1tImZCSVk3ISM7JCIyeHhGbFxFSFciISM8NyQkIjJMTExMT1trTCIhIzskIjI9QUFzJVtgPT0hIzw3JCQiMktMTEwmUSJHVCIhIzskIjIoKSkpKSkpKVFDbjgjISM8NyQkIixEMlg7XSIhIzUkIjIlKioqXChvVyZvXSMhIzw3JCQiMkpMTEx2di1lIiEjOyQiMlliYmJJI1tNRyEjPDckJCIsRDJZbG0iISM1JCIyLytdKG8+JVI+JCEjPDckJCIsdiJlcFs8ISM1JCIybW07SGRLaWAkISM8NyQkIjJNTEwkZS9UTT0hIzskIjJFQUEoNHBQJCpRISM8NyQkIjJMTExlREJKIj4hIzskIjFjYkkqKm9NQFUhIzs3JCQiMm1tbTtrRCEpKj4hIzskIjIleXhGU281dlghIzw3JCQiMmptbSJmYEAnMyMhIzskIjF4eF8nKlJjVVwhIzs3JCQiLHZ3JSlIOyMhIzUkIi1ESmxWaV8hIzc3JCQiMmxtbTskeSplQyMhIzskIjFXVyU+ajJ6ZyYhIzs3JCQiKjliOkwjISIpJCIxTExMJDM4WydmISM7NyQkIjInKioqKipcNWFgVCMhIzskIjEqKioqXFBhKFJKJyEjOzckJCIyKCoqKipcN1JWJ1wjISM7JCIxKioqXChveiE9bCchIzs3JCQiMicqKioqKlxAZmtlIyEjOyQiMUtMJGVSOHAtKCEjOzckJCIxTExMYDRObkUhIzokIjFjYmIwQidSTyghIzs3JCQiKjo/UHYjISIpJCIxTExMZVIkUXMoISM7NyQkIjJsbW0iek0pPiRHISM7JCIxV1c+ajYkKlwhKSEjOzckJCIqKGZhPEghIikkIjFubW1UOlcxJSkhIzs3JCQiMUxMJGVnYCEpKkghIzokIjFhYklDKyo9dSkhIzs3JCQiMikqKioqXCNHMkEzJCEjOyQiLXZWLmAjNCohIzc3JCQiMUxMTCQpR1trSiEjOiQiMUFBQVpgTU4lKiEjOzckJCIsRCJ5aF1LISM1JCIxbm1UJikzQyV6KiEjOzckJCIyb21tbSlmZExMISM7JCIyWVdXVyoqKilSLCIhIzs3JCQiMm1tbTtxNyU9TSEjOyQiMldXJXBEJ1EkXDUhIzs3JCQiMkxMTGUjcGEtTiEjOyQiMmNiSUNiJVIlMyIhIzs3JCQiKmFkKXpOISIpJCIyTUxMJDN0Zzs2ISM7NyQkIjJMTEwkR1VZb08hIzskIjJjYjA9d0VOOiIhIzs3JCQiMmxtbW01Onh1JCEjOyQiMjc2NjYnemEnPSIhIzs3JCQiLEQyOEEkUSEjNSQiLnZvV2I8QSIhIzc3JCQiLHZTKTM4UiEjNSQiLkRKXWBhRCIhIzc3JCQiI1MhIiIkIjJubW1tbW07SCIhIzstJSZDT0xPUkc2JiUkUkdCRyQiMidcXydlTSV5Zz4hIzwkIjFpajkuZUBSISkhIzskIjInXF8nZU0leWc+ISM8LSUqVEhJQ0tORVNTRzYjIiIkLSUnUE9JTlRTRzYmNyQkIiM1ISIiJCEiJiEiIjckJCIjPyEiIiQiIzUhIiI3JCQiI0khIiIkIitMTExMTCEjNS0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiM1ISIiLSUnQ1VSVkVTRzYlNyQ3JCQiIzohIiIkIiNEISIjNyQkIiNEISIiJCIxbW1tbW1tbW0hIzstJSZDT0xPUkc2JiUkUkdCRyQiIzUhIiIkIiIhISIiJCIiISEiIi0lJlNUWUxFRzYjJSZQT0lOVEctJSVWSUVXRzYkOyQiIiEhIiIkIiNTISIiJShERUZBVUxURy0mJSZfQVhJU0c2IyIiIjYmLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiItJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSYlJl9BWElTRzYjIiIjNiYtJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIS0lKlRISUNLTkVTU0c2IyIiIi0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJStBWEVTTEFCRUxTRzYlLUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NVEieDYiLyUnZmFtaWx5R1EoREVGQVVMVDYiLyUlc2l6ZUdRIzEwNiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2IlEhNiItJSVGT05URzYlJShERUZBVUxURyUoREVGQVVMVEciIzUtJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiQ/JiEiIi0lKUJPVU5EU19ZRzYjJCIjISkhIiItJS1CT1VORFNfV0lEVEhHNiMkIiUrTSEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiU/USEiIi0lKUNISUxEUkVORzYi
EXERCISE: Find the quadratic function passing closest to the points [1,-1/2], [2,1] and [3,1/3]
EXERCISE: Find the quadratic function passing closest to the points [1,2], [2.4,1.5] , [3,1] and [-0.4,3]
EXERCISE: Find the linear function passing closest to the points [1,-1/2], [2,1], [2.4,1] and [3,1/3] . You can use the graph below to study the distance.
pts:=[ [1,-1/2], [2,1], [2.4,1] , [3,1/3] ];
NyY3JCIiIiMhIiIiIiM3JEYnRiQ3JCQiI0NGJkYkNyQiIiQjRiRGLQ==
Note: The graph below works only if all the vertical blue segments have length larger than 0.
L:=x->(x-1);
display(plot(pts, x = 0 .. 5, style = point, symbol = box, color = green), plot([L(x)], x = 0 .. 5, color = red), seq(plot([pts[i][1], t, t = pts[i][2] .. L(pts[i][1])], thickness = 3, color = blue), i = 1 .. nops(pts)));
Zio2I0kieEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCY5JCIiIiEiIkYrRiVGJUYl
LSUlUExPVEc2LS0lJ0NVUlZFU0c2JjcmNyQkIiM1ISIiJCEiJiEiIjckJCIjPyEiIiQiIzUhIiI3JCQiI0MhIiIkIiM1ISIiNyQkIiNJISIiJCIxTExMTExMTEwhIzstJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiM1ISIiJCIiISEiIi0lJ1NZTUJPTEc2IyUkQk9YRy0lJlNUWUxFRzYjJSZQT0lOVEctJSdDVVJWRVNHNiQ3UzckJCIiISEiIiQhIzUhIiI3JCQiMkxMTCQzeCYpKjMiISM8JCExbm1tIkhVLCIqKSEjOzckJCIybm07SDJQIlE/ISM8JCExTEwzRkgnPSd6ISM7NyQkIjJMTExlUndYNSQhIzwkITFtbTsvT1UmKm8hIzs3JCQiMU1MJDN4JTN5VCEjOyQhMWxtO0hfIj4jZSEjOzckJCIxbm0ieiU0XFlfISM7JCEyTUwkM18hNE52JSEjPDckJCIxTUxlUi0vUGkhIzskITJsbTsvd2ZIdyQhIzw3JCQiLURjbXBpcyEjNyQhMjArXVBNLnR0IyEjPDckJCIxTUxlKik+VkIkKSEjOyQhMmVtOy8sb2xuIiEjPDckJCItREpidyFRKiEjNyQhMmEqKipcKG9XQj4nISM9NyQkIjJvbW1USU9vLyIhIzskIjElb21tVElPbyUhIzw3JCQiMk5MJDNfPmpVNiEjOyQiMllMTDNfPmpVIiEjPDckJCIyLCsrRDt2L0QiISM7JCIxLCsrRDt2L0QhIzs3JCQiLHY9aChlOCEjNSQiMSwrK3Y9aChlJCEjOzckJCIsdiRbNmo5ISM1JCIyLysrXVBbNmolISM8NyQkIjJNTGUqW3ooeWIiISM7JCIxTUxlKlt6KHliISM7NyQkIjJvbW1UWGcwbiIhIzskIjFvbW1UWGcwbiEjOzckJCIyb21tbUo8Z3ciISM7JCIxb21tbUo8Z3chIzs3JCQiLkQxTWNxKD0hIzckIi1EMU1jcSgpISM3NyQkIjJtbW07cFdgKD4hIzskIjFtbW07cFdgKCohIzs3JCQiMi4rRDFmIz0kMyMhIzskIjIuK0QxZiM9JDMiISM7NyQkIi52PXhwZT0jISM3JCIyLSt2PXhwZT0iISM7NyQkIjFubSJIMjhJSCMhIzokIjJybTtIMjhJSCIhIzs3JCQiMm9tInpwU1MiUiMhIzskIjJvbSJ6cFNTIlIiISM7NyQkIjJMTCQzXz9gKFwjISM7JCIyTEwkM18/YChcIiEjOzckJCIyT0xlKik+cHhnIyEjOyQiMk9MZSopPnB4ZyIhIzs3JCQiLnYkZjR0LkYhIzckIi52JGY0dC48ISM3NyQkIjJQTCRlKkdzdCFHISM7JCIyUEwkZSpHc3QhPSEjOzckJCIrRFJXOUghIiokIitEUlc5PiEiKjckJCItREpFPj5JISM2JCItREpFPj4/ISM2NyQkIjIuK0QxUlUwNyQhIzskIjIuK0QxUlUwNyMhIzs3JCQiLXY9UzJMSyEjNiQiLXY9UzJMQSEjNjckJCIyam1tO3ApPU1MISM7JCIyam1tO3ApPU1CISM7NyQkIjItKyt2PV1AVyQhIzskIjItKyt2PV1AVyMhIzs3JCQiMUwkZSpbJHoqUk4hIzokIjFMJGUqWyR6KlJEISM6NyQkIjIxKytEWUtwayQhIzskIjIxKytEWUtwayMhIzs3JCQiMmxtIkgycWNaUCEjOyQiMmxtIkgycWNaRiEjOzckJCIyLytESjVmRiZRISM7JCIyLytESjVmRiZHISM7NyQkIjJsbW1UZy5jJlIhIzskIjJsbW1UZy5jJkghIzs3JCQiLkRjRXNLMSUhIzckIjIpKipcaWxBRmpJISM7NyQkIjFNTEwkKSpwcDslISM6JCIxTUxMJCkqcHA7JCEjOjckJCIxTUwzeGUsdFUhIzokIjFNTDN4ZSx0SyEjOjckJCIxbjtIZE89eVYhIzokIjJ0bSJIZE89eUwhIzs3JCQiMjArK10jPiNbWiUhIzskIjIwKytdIz4jW1okISM7NyQkIjFubVQmRyFlJmUlISM6JCIydW07YUchZSZlJCEjOzckJCIxTUxMJClRayVvJSEjOiQiMU1MTCQpUWslbyQhIzo3JCQiLkQxTW0teiUhIzckIi5EMU1tLXokISM3NyQkIjIxK3YkNDBPIipbISM7JCIyMSt2JDQwTyIqUSEjOzckJCIjXSEiIiQiI1MhIiItJSZDT0xPUkc2JiUkUkdCRyQiIzUhIiIkIiIhISIiJCIiISEiIi0lJ0NVUlZFU0c2JTdTNyQkIiM1ISIiJCEiJiEiIjckJCIjNSEiIiQhMW5tO0hVLCIqWyEjOzckJCIjNSEiIiQhMUwkM0ZIJz0neiUhIzs3JCQiIzUhIiIkITJrbTsvT1UmKm8lISM8NyQkIiM1ISIiJCEybG07SF8iPiNlJSEjPDckJCIjNSEiIiQhMk5MM18hNE52VyEjPDckJCIjNSEiIiQhMmttVGcoZkh3ViEjPDckJCIjNSEiIiQhLnZWTElQRiUhIzg3JCQiIzUhIiIkITFuOy8sb2xuVCEjOzckJCIjNSEiIiQhLnZvV0I+MSUhIzg3JCQiIzUhIiIkITFMTCRlcGpKJlIhIzs3JCQiIzUhIiIkITJubTt6L290JlEhIzw3JCQiIzUhIiIkISx2JFtfXFAhIzY3JCQiIzUhIiIkISxEIilRN2skISM2NyQkIiM1ISIiJCEsRDsmKW9gJCEjNjckJCIjNSEiIiQhMm5tVDUwQUBXJCEjPDckJCIjNSEiIiQhMlBMTGVhUiVITCEjPDckJCIjNSEiIiQhMk5MTExvIylSQiQhIzw3JCQiIzUhIiIkITIvK3YkZk8lSDckISM8NyQkIiM1ISIiJCExTUxMM2BsQ0khIzs3JCQiIzUhIiIkIS52JDR1Im8iSCEjODckJCIjNSEiIiQhLkQiRy04OUchIzg3JCQiIzUhIiIkITJOTCQzRnApcHEjISM8NyQkIiM1ISIiJCEyTUwzLSRmZjNFISM8NyQkIiM1ISIiJCEybm07eiV6WS1EISM8NyQkIiM1ISIiJCExbjsvLDNCI1IjISM7NyQkIiM1ISIiJCEyLitEMS9waUgjISM8NyQkIiM1ISIiJCExbm1UNXhpIz4jISM7NyQkIiM1ISIiJCEyLSsrXTJjYjMjISM8NyQkIiM1ISIiJCEyMCtdKG90ITMpPiEjPDckJCIjNSEiIiQhMi8rdiQ0d1h6PSEjPDckJCIjNSEiIiQhMi8rXTcpZiNwdyIhIzw3JCQiIzUhIiIkITJQTEwkMzgiZW0iISM8NyQkIiM1ISIiJCEsRCIpXHliIiEjNjckJCIjNSEiIiQhMW47L14xLWc5ISM7NyQkIiM1ISIiJCEsdmBuSU4iISM2NyQkIiM1ISIiJCEyUkwzRipIVl83ISM8NyQkIiM1ISIiJCEyLit2byozQ1o2ISM8NyQkIiM1ISIiJCEyTkxMZVInUlc1ISM8NyQkIiM1ISIiJCExLSt2VnRGbiQqISM8NyQkIiM1ISIiJCExbG1tbSxJSSQpISM8NyQkIiM1ISIiJCExb207SDclKXBzISM8NyQkIiM1ISIiJCExTUwzRk07PWkhIzw3JCQiIzUhIiIkITEsKytdMnleXyEjPDckJCIjNSEiIiQhMU1MJGU5KD5XVCEjPDckJCIjNSEiIiQhMiZvbW1tNmNgSiEjPTckJCIjNSEiIiQhMikpKipcUGZPdDQjISM9NyQkIiM1ISIiJCEyKikqKlxpIVxSJzMiISM9NyQkIiM1ISIiJCIiISEiIi0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiM1ISIiLSUqVEhJQ0tORVNTRzYjIiIkLSUnQ1VSVkVTRzYlN1M3JCQiIz8hIiIkIiM1ISIiNyQkIiM/ISIiJCIjNSEiIjckJCIjPyEiIiQiIzUhIiI3JCQiIz8hIiIkIiM1ISIiNyQkIiM/ISIiJCIjNSEiIjckJCIjPyEiIiQiIzUhIiI3JCQiIz8hIiIkIiM1ISIiNyQkIiM/ISIiJCIjNSEiIjckJCIjPyEiIiQiIzUhIiI3JCQiIz8hIiIkIiM1ISIiNyQkIiM/ISIiJCIjNSEiIjckJCIjPyEiIiQiIzUhIiI3JCQiIz8hIiIkIiM1ISIiNyQkIiM/ISIiJCIjNSEiIjckJCIjPyEiIiQiIzUhIiI3JCQiIz8hIiIkIiM1ISIiNyQkIiM/ISIiJCIjNSEiIjckJCIjPyEiIiQiIzUhIiI3JCQiIz8hIiIkIiM1ISIiNyQkIiM/ISIiJCIjNSEiIjckJCIjPyEiIiQiIzUhIiI3JCQiIz8hIiIkIiM1ISIiNyQkIiM/ISIiJCIjNSEiIjckJCIjPyEiIiQiIzUhIiI3JCQiIz8hIiIkIiM1ISIiNyQkIiM/ISIiJCIjNSEiIjckJCIjPyEiIiQiIzUhIiI3JCQiIz8hIiIkIiM1ISIiNyQkIiM/ISIiJCIjNSEiIjckJCIjPyEiIiQiIzUhIiI3JCQiIz8hIiIkIiM1ISIiNyQkIiM/ISIiJCIjNSEiIjckJCIjPyEiIiQiIzUhIiI3JCQiIz8hIiIkIiM1ISIiNyQkIiM/ISIiJCIjNSEiIjckJCIjPyEiIiQiIzUhIiI3JCQiIz8hIiIkIiM1ISIiNyQkIiM/ISIiJCIjNSEiIjckJCIjPyEiIiQiIzUhIiI3JCQiIz8hIiIkIiM1ISIiNyQkIiM/ISIiJCIjNSEiIjckJCIjPyEiIiQiIzUhIiI3JCQiIz8hIiIkIiM1ISIiNyQkIiM/ISIiJCIjNSEiIjckJCIjPyEiIiQiIzUhIiI3JCQiIz8hIiIkIiM1ISIiNyQkIiM/ISIiJCIjNSEiIjckJCIjPyEiIiQiIzUhIiI3JCQiIz8hIiIkIiM1ISIiLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIzUhIiItJSpUSElDS05FU1NHNiMiIiQtJSdDVVJWRVNHNiU3UzckJCIjQyEiIiQiIzUhIiI3JCQiI0MhIiIkIjJubW07Jyk9KDM1ISM7NyQkIiNDISIiJCIyTEwkZSc0MGosIiEjOzckJCIjQyEiIiQiMm5tbTZoT1stIiEjOzckJCIjQyEiIiQiMm5tbSJ5WVVMNSEjOzckJCIjQyEiIiQiMktMJGVGPig+LyIhIzs3JCQiI0MhIiIkIjJtbTs+SycqKVw1ISM7NyQkIiNDISIiJCItRHQ6NWU1ISM2NyQkIiNDISIiJCIybW07ZlgoZW01ISM7NyQkIiNDISIiJCItRENoL3Y1ISM2NyQkIiNDISIiJCIyS0xMViFwdSQzIiEjOzckJCIjQyEiIiQiMm5tbWhiNTk0IiEjOzckJCIjQyEiIiQiKzghUSs1IiEiKjckJCIjQyEiIiQiKyYqM3EzNiEiKjckJCIjQyEiIiQiMioqKioqKnA9XHE2IiEjOzckJCIjQyEiIiQiMm5tO2ZCSVk3IiEjOzckJCIjQyEiIiQiMkxMTGokW2tMNiEjOzckJCIjQyEiIiQiMktMTGBRIkdUNiEjOzckJCIjQyEiIiQiLUQyWDtdNiEjNjckJCIjQyEiIiQiMkxMTGBkRiFlNiEjOzckJCIjQyEiIiQiLUQyWWxtNiEjNjckJCIjQyEiIiQiLXYiZXBbPCIhIzY3JCQiI0MhIiIkIjJMTExlL1RNPSIhIzs3JCQiI0MhIiIkIjJMTCRlREJKIj4iISM7NyQkIiNDISIiJCIybG1tVGMtKSo+IiEjOzckJCIjQyEiIiQiMmxtO2ZgQCczNyEjOzckJCIjQyEiIiQiLXZ3JSlIOzchIzY3JCQiI0MhIiIkIjJtbW1KeSplQzchIzs3JCQiI0MhIiIkIis5YjpMNyEiKjckJCIjQyEiIiQiLDBUTjpDIiEjNTckJCIjQyEiIiQiLUQiUlYnXDchIzY3JCQiI0MhIiIkIjIqKioqKlxAZmtlNyEjOzckJCIjQyEiIiQiMktMTGA0Tm5FIiEjOzckJCIjQyEiIiQiMioqKioqKlwsc2BGIiEjOzckJCIjQyEiIiQiMm1tO3pNKT4kRyIhIzs3JCQiI0MhIiIkIisoZmE8SCIhIio3JCQiI0MhIiIkIjJMTCRlZ2AhKSpIIiEjOzckJCIjQyEiIiQiLURHMkEzOCEjNjckJCIjQyEiIiQiMkxMTCQpR1trSiIhIzs3JCQiI0MhIiIkIjIpKioqXDd5aF1LIiEjOzckJCIjQyEiIiQiMm1tbScpZmRMTCIhIzs3JCQiI0MhIiIkIjJubW0sRlQ9TSIhIzs3JCQiI0MhIiIkIjJNTCRlI3BhLU4iISM7NyQkIiNDISIiJCIrYWQpek4iISIqNyQkIiNDISIiJCIyS0xMR1VZb08iISM7NyQkIiNDISIiJCIybm1tMV5yWlAiISM7NyQkIiNDISIiJCItRDI4QSRRIiEjNjckJCIjQyEiIiQiMikqKipcMiUpMzhSIiEjOzckJCIjQyEiIiQiIzkhIiItJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIjNSEiIi0lKlRISUNLTkVTU0c2IyIiJC0lJ0NVUlZFU0c2JTdTNyQkIiNJISIiJCIrTExMTEwhIzU3JCQiI0khIiIkIjJPLyZ5TiM+bXAkISM8NyQkIiNJISIiJCIxViI+dE43RiwlISM7NyQkIiNJISIiJCIyJTM9akpAPm9WISM8NyQkIiNJISIiJCIxbCpRbURHZ3MlISM7NyQkIiNJISIiJCIxKj4yIVxPOyMzJiEjOzckJCIjSSEiIiQiMUNnQlluTTdhISM7NyQkIiNJISIiJCIxJTNsJT1BQmFkISM7NyQkIiNJISIiJCIxLG0iSG01eTUnISM7NyQkIiNJISIiJCIxUUQibzxiLVknISM7NyQkIiNJISIiJCIxNj9IIW8oeUFvISM7NyQkIiNJISIiJCIxKEcobzFsNVVyISM7NyQkIiNJISIiJCIxSyt2dVFlLHYhIzs3JCQiI0khIiIkIjE8UlUiSFBEJ3khIzs3JCQiI0khIiIkIjF2M1YieiNRNSMpISM7NyQkIiNJISIiJCIxUDtqSClmaV8pISM7NyQkIiNJISIiJCIxI2ZMLj1vPSEqKSEjOzckJCIjSSEiIiQiMUJtK0F4MD8jKiEjOzckJCIjSSEiIiQiMU9eKD0heT0hZiohIzs3JCQiI0khIiIkIjFXc29RYyJ5IioqISM7NyQkIiNJISIiJCIxKSkpKT1JdnNGNSEjOjckJCIjSSEiIiQiMiJSaWchZmM+MSIhIzs3JCQiI0khIiIkIjE0dkdDNW4oNCIhIzo3JCQiI0khIiIkIjIkKVxZSy1vLzgiISM7NyQkIiNJISIiJCIxR1dNPFMlZTsiISM6NyQkIiNJISIiJCIxaiwoSHQqZS03ISM6NyQkIiNJISIiJCIxRElXJylwZE03ISM6NyQkIiNJISIiJCIyLjs4bDRDInA3ISM7NyQkIiNJISIiJCIyR3dfO2s5W0kiISM7NyQkIiNJISIiJCIxRywyeDN0UjghIzo3JCQiI0khIiIkIjJPIWU+SVRedDghIzs3JCQiI0khIiIkIjJzKVshSG5DNVQiISM7NyQkIiNJISIiJCIyTHl4UWNIWlciISM7NyQkIiNJISIiJCIxIkdjIkhuciFbIiEjOjckJCIjSSEiIiQiMlYvVidca0s4OiEjOzckJCIjSSEiIiQiMid6NFwoW3gqWzohIzs3JCQiI0khIiIkIjIiUmJ2TkJfI2UiISM7NyQkIiNJISIiJCIxJm8reE8nZTw7ISM6NyQkIiNJISIiJCIyR2ZbIW95Jz1sIiEjOzckJCIjSSEiIiQiMWIobz1VZHhvIiEjOjckJCIjSSEiIiQiMlZBc3hLQkJzIiEjOzckJCIjSSEiIiQiMjgkSC1mPm5kPCEjOzckJCIjSSEiIiQiMk51ZmQpeXMjeiIhIzs3JCQiI0khIiIkIjJiO2o7a1NcIz0hIzs3JCQiI0khIiIkIjJpN09eNGc9Jz0hIzs3JCQiI0khIiIkIjFVQlclSCIpWyo9ISM6NyQkIiNJISIiJCIxTnA/IXkpM0k+ISM6NyQkIiNJISIiJCIyIzR3WE9veWo+ISM7NyQkIiNJISIiJCIjPyEiIi0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiM1ISIiLSUqVEhJQ0tORVNTRzYjIiIkLSUlVklFV0c2JDskIiIhISIiJCIjXSEiIiUoREVGQVVMVEctJiUmX0FYSVNHNiMiIiI2Ji0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqTElORVNUWUxFRzYjIiIhLSUqVEhJQ0tORVNTRzYjIiIiLSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0mJSZfQVhJU0c2IyIiIzYmLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiItJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUrQVhFU0xBQkVMU0c2JS1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRIng2Ii8lJ2ZhbWlseUdRKERFRkFVTFQ2Ii8lJXNpemVHUSMxMDYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiJRITYiLSUlRk9OVEc2JSUoREVGQVVMVEclKERFRkFVTFRHIiM1LSUlUk9PVEc2Jy0lKUJPVU5EU19YRzYjJCIkSSUhIiItJSlCT1VORFNfWUc2IyQiJD8iISIiLSUtQk9VTkRTX1dJRFRIRzYjJCIlIVwkISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJV1QISIiLSUpQ0hJTERSRU5HNiI=
distancePointLine:=P->(m*P[1]+b-P[2])^2;
Zio2I0kiUEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlKiQsKComSSJtR0YlIiIiJjkkNiNGLUYtRi1JImJHRiVGLSZGLzYjIiIjISIiRjRGJUYlRiU=
distancePointsLine:=pts->sum(distancePointLine(pts[i]),i=1..nops(pts));
Zio2I0kkcHRzRzYiRiU2JEkpb3BlcmF0b3JHRiVJJmFycm93R0YlRiUtSSRzdW1HRiU2JC1JMmRpc3RhbmNlUG9pbnRMaW5lR0YlNiMmOSQ2I0kiaUdGJS9GMjsiIiItSSVub3BzRyUqcHJvdGVjdGVkRzYjRjBGJUYlRiU=
DD:=distancePointsLine(pts);
LCoqJCwoSSJtRzYiIiIiSSJiR0YmRicjRiciIiNGJ0YqRicqJCwoRiVGKkYoRichIiJGJ0YqRicqJCwoRiUkIiNDRi1GKEYnRi1GJ0YqRicqJCwoRiUiIiRGKEYnI0YtRjRGJ0YqRic=
simplify(DD);
LC4qJEkibUc2IiIiIyQiKysrK3c+ISIpKiZGJCIiIkkiYkdGJUYrJCIrKysrIW8iRilGJCQhKysrKyspKiEiKiokRixGJiQiIiUiIiFGLCQhK25tbW1PRjEkIis2NjZoQkYxRis=
eq1:=diff(DD,m);
eq2:=diff(DD,b);
LChJIm1HNiIkIiVfUiEiI0kiYkdGJCQiJG8iISIiJCEjKSpGKyIiIg==
LChJIm1HNiIkIiRvIiEiIkkiYkdGJCIiKSMhIzYiIiQiIiI=
sol:=solve({eq1,eq2},{m,b});
PCQvSSJiRzYiJCEraitoPGUhIzUvSSJtR0YlJCIrKj1JRyZcRig=
foundedLine:=x->subs(sol, m*x+b);
Zio2I0kieEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUklc3Vic0clKnByb3RlY3RlZEc2JEkkc29sR0YlLCYqJkkibUdGJSIiIjkkRjFGMUkiYkdGJUYxRiVGJUYl
L:=foundedLine;
display(plot(pts, x = 0 .. 5, style = point, symbol = box, color = green), plot([L(x)], x = 0 .. 5, color = red), seq(plot([pts[i][1], t, t = pts[i][2] .. L(pts[i][1])], thickness = 3, color = blue), i = 1 .. nops(pts)));
SSxmb3VuZGVkTGluZUc2Ig==
LSUlUExPVEc2LS0lJ0NVUlZFU0c2JjcmNyQkIiM1ISIiJCEiJiEiIjckJCIjPyEiIiQiIzUhIiI3JCQiI0MhIiIkIiM1ISIiNyQkIiNJISIiJCIxTExMTExMTEwhIzstJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiM1ISIiJCIiISEiIi0lJ1NZTUJPTEc2IyUkQk9YRy0lJlNUWUxFRzYjJSZQT0lOVEctJSdDVVJWRVNHNiQ3UzckJCIiISEiIiQhK2oraDxlISM1NyQkIjJMTEwkM3gmKSozIiEjPCQhMUtfL1o/I3lGJiEjOzckJCIybm07SDJQIlE/ISM8JCEyViNRZiFRYiIzWyEjPDckJCIyTExMZVJ3WDUkISM8JCEyYilmRSQ0bSp6VSEjPDckJCIxTUwkM3glM3lUISM7JCExXzpbQ2NGW1AhIzs3JCQiMW5tInolNFxZXyEjOyQhMjpKJW8oPTciPkshIzw3JCQiMU1MZVItL1BpISM7JCExKVFRVCUqNCZHRiEjOzckJCItRGNtcGlzISM3JCEyX2NQeHQ+MEEjISM8NyQkIjFNTGUqKT5WQiQpISM7JCEyTUwhZlJiOiZwIiEjPDckJCItREpidyFRKiEjNyQhMiRRJSozIj53OTwiISM8NyQkIjJvbW1USU9vLyIhIzskITFwOSN5N3chR2ohIzw3JCQiMk5MJDNfPmpVNiEjOyQhMltsYyFlLlskZSIhIz03JCQiMiwrK0Q7di9EIiEjOyQiMS1vQ0MyInl2JCEjPDckJCIsdj1oKGU4ISM1JCIwKlsiNG1MNTcqISM7NyQkIix2JFs2ajkhIzUkIjIjSFE5dSNcKkc5ISM8NyQkIjJNTGUqW3ooeWIiISM7JCIyWEJ5TChbSCkqPSEjPDckJCIyb21tVFhnMG4iISM7JCIycEBZcD0jUmNDISM8NyQkIjJvbW1tSjxndyIhIzskIjFSTUc8UTxISCEjOzckJCIuRDFNY3EoPSEjNyQiMW44SVo3OHpNISM7NyQkIjJtbW07cFdgKD4hIzskIjFpdmNmbiRmJ1IhIzs3JCQiMi4rRDFmIz0kMyMhIzskIjFIbi9oJlIrXSUhIzs3JCQiLnY9eHBlPSMhIzckIjFWXU9LPGozXSEjOzckJCIxbm0iSDI4SUgjISM6JCIxRkxKM1ZIUmIhIzs3JCQiMm9tInpwU1MiUiMhIzskIjBqaHo/M20tJyEjOjckJCIyTEwkM18/YChcIyEjOyQiMV9NYiMzVUFiJyEjOzckJCIyT0xlKik+cHhnIyEjOyQiMSM9WDomekEpNCghIzs3JCQiLnYkZjR0LkYhIzckIjFLa19iLV50diEjOzckJCIyUEwkZSpHc3QhRyEjOyQiMSZRNVk7R28zKSEjOzckJCIrRFJXOUghIiokIjFedylHek5yaCkhIzs3JCQiLURKRT4+SSEjNiQiMVg4WVokUWY4KiEjOzckJCIyLitEMVJVMDckISM7JCIwdFRgXDF6aiohIzo3JCQiLXY9UzJMSyEjNiQiMWJPLmFjXz41ISM6NyQkIjJqbW07cCk9TUwhIzskIjJ4ITRcTWdncDUhIzs3JCQiMi0rK3Y9XUBXJCEjOyQiMjY/cytgeEk3IiEjOzckJCIxTCRlKlskeipSTiEjOiQiMidSZFxfMWByNiEjOzckJCIyMSsrRFlLcGskISM7JCIyOFReTXItWEEiISM7NyQkIjJsbSJIMnFjWlAhIzskIjJbWCxZSFhWRiIhIzs3JCQiMi8rREo1ZkYmUSEjOyQiMms/IVJgXldFOCEjOzckJCIybG1tVGcuYyZSISM7JCIxKFFWJClHI1F4OCEjOjckJCIuRGNFc0sxJSEjNyQiMkJSSiFbKDMyViIhIzs3JCQiMU1MTCQpKnBwOyUhIzokIjJpLDZuWW8/WyIhIzs3JCQiMU1MM3hlLHRVISM6JCIyW1k0cj4iZk06ISM7NyQkIjFuO0hkTz15ViEjOiQiMilmP3k3IXptZSIhIzs3JCQiMjArK10jPiNbWiUhIzskIjIoNFJ2MEJhTTshIzs3JCQiMW5tVCZHIWUmZSUhIzokIjIxXHAzLypSKm8iISM7NyQkIjFNTEwkKVFrJW8lISM6JCIyRnYzKGVOWVE8ISM7NyQkIi5EMU1tLXolISM3JCIyZSEpPiNvbHgheiIhIzs3JCQiMjErdiQ0ME8iKlshIzskIjI9XDZMelkzJT0hIzs3JCQiI10hIiIkIiwjKTNhWSo9ISM1LSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJSdDVVJWRVNHNiU3UzckJCIjNSEiIiQhIiYhIiI3JCQiIzUhIiIkITJgWislcFInKTRcISM8NyQkIiM1ISIiJCEyYWNJNjRQOSRbISM8NyQkIiM1ISIiJCExNFw1a3lCVlohIzs3JCQiIzUhIiIkITIkPS5OJipSWGFZISM8NyQkIiM1ISIiJCEyTydmIzM1IzRtWCEjPDckJCIjNSEiIiQhMWo/XDgkcFRbJSEjOzckJCIjNSEiIiQhMVUzTzdJTSpSJSEjOzckJCIjNSEiIiQhMmRsJlFJYmg2ViEjPDckJCIjNSEiIiQhMiRleGIiUnBUQSUhIzw3JCQiIzUhIiIkITIlM2xGKkc/VTglISM8NyQkIiM1ISIiJCEybidvLXJJKlwwJSEjPDckJCIjNSEiIiQhMXdMPikpPiFlJ1IhIzs3JCQiIzUhIiIkITFvVz95WUN3USEjOzckJCIjNSEiIiQhMWNUWTonUioqeSQhIzs3JCQiIzUhIiIkITA9Qi5ObTpyJCEjOjckJCIjNSEiIiQhMDFZdyZIUD1PISM6NyQkIiM1ISIiJCExKDReKUhmVVJOISM7NyQkIiM1ISIiJCExPUJfbzxmWk0hIzs3JCQiIzUhIiIkITJ0YVldKEhJbUwhIzw3JCQiIzUhIiIkITJqJEgtJkc7ckYkISM8NyQkIiM1ISIiJCEyV1wwYFkqPSM+JCEjPDckJCIjNSEiIiQhMXooUSRRc2QuSiEjOzckJCIjNSEiIiQhMVEpUkFiLkEtJCEjOzckJCIjNSEiIiQhMSUpKnkmUTVWTUghIzs3JCQiIzUhIiIkITFANCFcMWdLJUchIzs3JCQiIzUhIiIkITJEaytWWSYqUXcjISM8NyQkIiM1ISIiJCEwMSlRQCZ6InlFISM6NyQkIiM1ISIiJCEyT0pbXWNFJyplIyEjPDckJCIjNSEiIiQhMUBvUnRaKkhdIyEjOzckJCIjNSEiIiQhMiRwLVhnUzw+QyEjPDckJCIjNSEiIiQhMilmTXdcYTVFQiEjPDckJCIjNSEiIiQhMip5d21qInpDQyMhIzw3JCQiIzUhIiIkITJjMShlYC0+YEAhIzw3JCQiIzUhIiIkITJbaTcsQiJHcz8hIzw3JCQiIzUhIiIkITI9WCJlJ0hFUSk+ISM8NyQkIiM1ISIiJCExJHkjZSg0KGYrPiEjOzckJCIjNSEiIiQhMiZlYkEuJylmOD0hIzw3JCQiIzUhIiIkITJLUDlKblQmRzwhIzw3JCQiIzUhIiIkITJvKFImKlxcXFI7ISM8NyQkIiM1ISIiJCExQms0c0J0YDohIzs3JCQiIzUhIiIkITEjXHQjXHgtbTkhIzs3JCQiIzUhIiIkITJGazNdT1whejghIzw3JCQiIzUhIiIkITJ2eSU0al83KkgiISM8NyQkIiM1ISIiJCEyXkJyR0FCdj8iISM8NyQkIiM1ISIiJCEyNlcqNG1LZkQ2ISM8NyQkIiM1ISIiJCEyaVVwUSUpUSNRNSEjPDckJCIjNSEiIiQhMXo8QWNeSFkmKiEjPDckJCIjNSEiIiQhKnUpelonKSEjNS0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiM1ISIiLSUqVEhJQ0tORVNTRzYjIiIkLSUnQ1VSVkVTRzYlN1M3JCQiIz8hIiIkIis6LjApMyUhIzU3JCQiIz8hIiIkIjFBJjRQKlIicEAlISM7NyQkIiM/ISIiJCIyMkxDOmZQIUhWISM8NyQkIiM/ISIiJCIyazUzUklLXlglISM8NyQkIiM/ISIiJCIyJkdtKlEmRzEjZSUhIzw3JCQiIz8hIiIkIjFFeThELFIzWiEjOzckJCIjPyEiIiQiMi51JWZ3O15EWyEjPDckJCIjPyEiIiQiMWJMPWZVeVlcISM7NyQkIiM/ISIiJCIxQih5dmAvQTImISM7NyQkIiM/ISIiJCIxWTFiImZBcz4mISM7NyQkIiM/ISIiJCIxUkh0RSE+ZUsmISM7NyQkIiM/ISIiJCIuUW8kbzNSYSEjODckJCIjPyEiIiQiMFczT2MqZm1iISM6NyQkIiM/ISIiJCIwPCcpKXBlaiVwJiEjOjckJCIjPyEiIiQiMSJSTmJkQSE9ZSEjOzckJCIjPyEiIiQiMUxEI1xfcCskZiEjOzckJCIjPyEiIiQiMTxvOiY9L0wxJyEjOzckJCIjPyEiIiQiMFQlej05PHdoISM6NyQkIiM/ISIiJCIxV3FMVkdZMmohIzs3JCQiIz8hIiIkIjFNSWArIXlPVSchIzs3JCQiIz8hIiIkIjBmKTNaVz1ebCEjOjckJCIjPyEiIiQiMSdbZG90K0VuJyEjOzckJCIjPyEiIiQiMTloI3kqZUcqeichIzs3JCQiIz8hIiIkIjE4RFVBQ2k6cCEjOzckJCIjPyEiIiQiMSRIPTIzMjYvKCEjOzckJCIjPyEiIiQiMSQpbydRUF05PCghIzs3JCQiIz8hIiIkIjFiLnMkZjlcRyghIzs3JCQiIz8hIiIkIjElKWZoZiFmdVMoISM7NyQkIiM/ISIiJCIxJHlxUVtmU2AoISM7NyQkIiM/ISIiJCIxXlgwK0wiemwoISM7NyQkIiM/ISIiJCIxQnBsTiNbeHgoISM7NyQkIiM/ISIiJCIxXGhNK1ghMyJ6ISM7NyQkIiM/ISIiJCIxaWUzczlPSSEpISM7NyQkIiM/ISIiJCIxakVNZVMsZSIpISM7NyQkIiM/ISIiJCIxV2hSdWlvdCMpISM7NyQkIiM/ISIiJCIwankmZmw5KyUpISM6NyQkIiM/ISIiJCIxc1tsS2U4PiYpISM7NyQkIiM/ISIiJCIxWCxBM1JeVicpISM7NyQkIiM/ISIiJCIxKGZGOEA7XncpISM7NyQkIiM/ISIiJCIwPm1KYkFDKikpISM6NyQkIiM/ISIiJCIxIlJ1NXBMXSwqISM7NyQkIiM/ISIiJCIxaVhxKUdALzkqISM7NyQkIiM/ISIiJCIxQiY+UGhxWkUqISM7NyQkIiM/ISIiJCIxcHApKkdaLnokKiEjOzckJCIjPyEiIiQiMWFBeSpIJSoqNCYqISM7NyQkIiM/ISIiJCIxTj4zLGg3RicqISM7NyQkIiM/ISIiJCIxOVtxSkksXygqISM7NyQkIiM/ISIiJCIxJT5DJmZkYXIpKiEjOzckJCIjPyEiIiQiIzUhIiItJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIjNSEiIi0lKlRISUNLTkVTU0c2IyIiJC0lJ0NVUlZFU0c2JTdTNyQkIiNDISIiJCIrKFEjPXBnISM1NyQkIiNDISIiJCIxNDteQ0knWzonISM7NyQkIiNDISIiJCIxImV5ISpHOCVIaSEjOzckJCIjQyEiIiQiMTlIYmVHRDhqISM7NyQkIiNDISIiJCIxdWY2ciJbd1InISM7NyQkIiNDISIiJCIxPS5cbkJrIlsnISM7NyQkIiNDISIiJCIxeVtgUmReZmwhIzs3JCQiI0MhIiIkIjF1X28wIlwsaychIzs3JCQiI0MhIiIkIjFOX0lbLWFCbiEjOzckJCIjQyEiIiQiMTpLWmFSbTFvISM7NyQkIiNDISIiJCIxSiNwTSFwOyMqbyEjOzckJCIjQyEiIiQiMVteW1p6Wm5wISM7NyQkIiNDISIiJCIxMltuWC5FX3EhIzs3JCQiI0MhIiIkIjFmcmBvM1JQciEjOzckJCIjQyEiIiQiMUV0aCwqSCU+cyEjOzckJCIjQyEiIiQiMVlSOC4jSFJIKCEjOzckJCIjQyEiIiQiMUFST3lnXiNRKCEjOzckJCIjQyEiIiQiMW5LayJ5Z3ZYKCEjOzckJCIjQyEiIiQiMTpXJz1oYlthKCEjOzckJCIjQyEiIiQiMXgpXCJHajdBdyEjOzckJCIjQyEiIiQiMVdabV1YIXBxKCEjOzckJCIjQyEiIiQiMTRAIW9ZTHd5KCEjOzckJCIjQyEiIiQiMTw/eEBjJz0oeSEjOzckJCIjQyEiIiQiMSRSblcuPCNceiEjOzckJCIjQyEiIiQiMENzRyk0bEshKSEjOjckJCIjQyEiIiQiMWszZTFhSj4iKSEjOzckJCIjQyEiIiQiMUZcPVVxdiU+KSEjOzckJCIjQyEiIiQiMFlHYTJPaUYpISM6NyQkIiNDISIiJCIxcCIpXCopPVRnJCkhIzs3JCQiI0MhIiIkIjFqOVIsOndVJSkhIzs3JCQiI0MhIiIkIjF2akImKSpRQ18pISM7NyQkIiNDISIiJCIxUXFTWXMhNGgpISM7NyQkIiNDISIiJCIxXlAmUSIqKlIhcCkhIzs3JCQiI0MhIiIkIjFeS3MuYEZ2KCkhIzs3JCQiI0MhIiIkIjFRJDM+LSY9XykpISM7NyQkIiNDISIiJCIxI0guJWV3RU8qKSEjOzckJCIjQyEiIiQiMWM+am1HUTohKiEjOzckJCIjQyEiIiQiMW4+Q2I1Myk0KiEjOzckJCIjQyEiIiQiMUI0WV1PJCp5IiohIzs3JCQiI0MhIiIkIjFMbWlDKXlORSohIzs3JCQiI0MhIiIkIjEoemx3PS1eTSohIzs3JCQiI0MhIiIkIjFpTjAsO1pHJSohIzs3JCQiI0MhIiIkIjFvaGd0MTo2JiohIzs3JCQiI0MhIiIkIjE5ZGNhVDcoZSohIzs3JCQiI0MhIiIkIjEiR1t0TCk+dScqISM7NyQkIiNDISIiJCIxeHlXQyZ5P3YqISM7NyQkIiNDISIiJCIxaHpgWl02TikqISM7NyQkIiNDISIiJCIxInljYGYiZjkqKiEjOzckJCIjQyEiIiQiIzUhIiItJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIjNSEiIi0lKlRISUNLTkVTU0c2IyIiJC0lJ0NVUlZFU0c2JTdTNyQkIiNJISIiJCIrTExMTEwhIzU3JCQiI0khIiIkIjIwP2wnPTt1ZE0hIzw3JCQiI0khIiIkIjEuXSJHZykpZmMkISM7NyQkIiNJISIiJCIxOCxwZE9zKG8kISM7NyQkIiNJISIiJCIxVytgXmNFNVEhIzs3JCQiI0khIiIkIjJtJ1JtQ19BS1IhIzw3JCQiI0khIiIkIjFqNigzTyhIWFMhIzs3JCQiI0khIiIkIjImeSo+JjNxUGlUISM8NyQkIiNJISIiJCIxXkRnWjRZJEclISM7NyQkIiNJISIiJCIyMiZvJipvbDovVyEjPDckJCIjSSEiIiQiMTd1KjMmb0lHWCEjOzckJCIjSSEiIiQiMmsxMm9bZXdqJSEjPDckJCIjSSEiIiQiMThvKCpHRHdnWiEjOzckJCIjSSEiIiQiMjtJciJbP1AlKVshIzw3JCQiI0khIiIkIjE+QD9DRlwuXSEjOzckJCIjSSEiIiQiMTZHKzVkbTZeISM7NyQkIiNJISIiJCIxQl0kRyZRSFNfISM7NyQkIiNJISIiJCIxel4qPndlI1xgISM7NyQkIiNJISIiJCIvWVpiMyx3YSEjOTckJCIjSSEiIiQiMShmPmIkej8pZSYhIzs3JCQiI0khIiIkIjFXJjRZImZJNmQhIzs3JCQiI0khIiIkIjFPJVIzSUMmR2UhIzs3JCQiI0khIiIkIjE6VCYpKlFIMyZmISM7NyQkIiNJISIiJCIxMWozVE85amchIzs3JCQiI0khIiIkIjBgb1R0KkclPSchIzo3JCQiI0khIiIkIjF2PFJ4azc1aiEjOzckJCIjSSEiIiQiMTx1R0t4bT5rISM7NyQkIiNJISIiJCIxYGFiKUd2emAnISM7NyQkIiNJISIiJCIxV0llcSYpPmdtISM7NyQkIiNJISIiJCIxODAwMi14em4hIzs3JCQiI0khIiIkIjEmeSJmOD5ZJipvISM7NyQkIiNJISIiJCIxeCk0KUd5IlItKCEjOzckJCIjSSEiIiQiMW80JT49VCRSciEjOzckJCIjSSEiIiQiMCUpSCFbLWVpcyEjOjckJCIjSSEiIiQiL0okKWVKRHV0ISM5NyQkIiNJISIiJCIxWkAtWjZNJ1woISM7NyQkIiNJISIiJCIxUSdvJEhrQDZ3ISM7NyQkIiNJISIiJCIxdChHIyo+JUhKeCEjOzckJCIjSSEiIiQiMU1eW25AcFt5ISM7NyQkIiNJISIiJCIxNEBQZ3BmcnohIzs3JCQiI0khIiIkIi9NXSQpKW8qKjMpISM5NyQkIiNJISIiJCIxQjN1cjctNiMpISM7NyQkIiNJISIiJCIxMmMzJkdyNUwpISM7NyQkIiNJISIiJCIxKVJQJXpaUVQlKSEjOzckJCIjSSEiIiQiMS1nJEdcO3ljKSEjOzckJCIjSSEiIiQiMTZKXUMmKSozbykhIzs3JCQiI0khIiIkIjEzJEdnYm45ISkpISM7NyQkIiNJISIiJCIxKXBmTV1ubyIqKSEjOzckJCIjSSEiIiQiKzIwKTMvKiEjNS0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiM1ISIiLSUqVEhJQ0tORVNTRzYjIiIkLSUlVklFV0c2JDskIiIhISIiJCIjXSEiIiUoREVGQVVMVEctJiUmX0FYSVNHNiMiIiI2Ji0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqTElORVNUWUxFRzYjIiIhLSUqVEhJQ0tORVNTRzYjIiIiLSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0mJSZfQVhJU0c2IyIiIzYmLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiItJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUrQVhFU0xBQkVMU0c2JS1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRIng2Ii8lJ2ZhbWlseUdRKERFRkFVTFQ2Ii8lJXNpemVHUSMxMDYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiJRITYiLSUlRk9OVEc2JSUoREVGQVVMVEclKERFRkFVTFRHIiM1LSUlUk9PVEc2Jy0lKUJPVU5EU19YRzYjJCIkPyYhIiItJSlCT1VORFNfWUc2IyQiI10hIiItJS1CT1VORFNfV0lEVEhHNiMkIiUrTSEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiUhKVEhIiItJSlDSElMRFJFTkc2Ig==
plot3d(DD,m=-10..10,b=-10..10);
LSUnUExPVDNERzYlLSUlR1JJREc2JTskISQrIiEiIiQiJCsiISIiOyQhJCsiISIiJCIkKyIhIiI3Ozc7JCIxeXh4eHgtJD4lISM3JCIxTExMTEwzJylSISM3JCIxV1dXV1dwJXkkISM3JCIyODY2NjZoKSllJCEjOCQiMUxMTExMZSlSJCEjNyQiMjg2NjY2aFFAJCEjOCQiMlZXV1dXJXBNSSEjOCQiMk5MTExMJDNoRyEjOCQiMiR5eHh4eC0kcCMhIzgkIjF5eHh4eF9JRCEjNyQiMk5MTExMJGV0QiEjOCQiMlZXV1dXJT5BQSEjOCQiMjg2NjY2aGoyIyEjOCQiMk5MTExMJDNPPiEjOCQiMjQ2NjY2aDghPSEjOCQiMlZXV1dXJT5zOyEjOCQiMUxMTExMZVs6ISM3JCIyd3h4eHhGMFYiISM4JCIyeXh4eHhGIT04ISM4JCIyTkxMTEwkMzY3ISM4JCIyW1dXV1clcDQ2ISM4JCIyNjY2NjZoUSwiISM4JCIxTkxMTEwkZUIqISM4JCIxNzY2NjZoKVEpISM4JCIxVldXV1clcGYoISM4NzskIidESFAhIiMkIjFBQUFBQShSYCQhIzckIjEqKioqKioqKipcVU0kISM3JCIxTExMTEwzZ0ohIzckIjFBQUFBQVoiKUghIzckIjJsbW1tbTslM0chIzgkIjJsbW1tbTs0ayMhIzgkIjI8QUFBQXMqeUMhIzgkIjFMTExMTGVBQiEjNyQiJ3ZyQCEiIyQiMkFBQUFBc2stIyEjOCQiMiYqKioqKioqKipcbik9ISM4JCIxTExMTExlXzwhIzckIjJBQUFBQXNSaSIhIzgkIjJqbW1tbTs0XSIhIzgkIjJsbW1tbTtNUSIhIzgkIjI8QUFBQXM5RiIhIzgkIjFMTExMTDNsNiEjNyQiMikqKioqKioqKipcVTEiISM4JCIxQUFBQUFzKm8qISM4JCIxKioqKioqKioqKlwjeikhIzgkIjFKTExMTCQzJnohIzgkIjFBQUFBQXNrciEjOCQiMW1tbW1tO01rISM4JCIxbG1tbW07ZmQhIzg3OyQiMmxtbW1tO0hIJCEjOCQiMWJiYmJiSTRKISM3JCIxKioqKioqKioqXDckSCEjNyQiMSoqKioqKioqKlwoZUYhIzckIjJaYmJiYjA9ZiMhIzgkIjFtbW1tbVRJQyEjNyQiMUxMTExMZXVBISM3JCIyWmJiYmIwVjcjISM4JCIyRUxMTEwkZXo+ISM4JCIxbW1tbW1UUz0hIzckIjJfYmJiYjBvcSIhIzgkIjIkKioqKioqKioqXCh5OiEjOCQiMiYqKioqKioqKipcaVgiISM4JCIyX2JiYmIwJFI4ISM4JCIxbW1tbW0iekEiISM3JCIxTExMTEwzQTYhIzckIjJfYmJiYjA9LSIhIzgkIjFHTExMTCQzRiohIzgkIjFqbW1tbTt6JCkhIzgkIjFgYmJiYjBWdiEjOCQiMSgqKioqKioqKipcaW4hIzgkIjEoKioqKioqKioqXFBnISM4JCIxYGJiYmIwb2AhIzgkIjFsbW1tbTthWiEjOCQiMjlMTExMTGU+JSEjOTc7JCIydXh4eHhGUylHISM4JCIxTExMTEwzN0YhIzckIjJNV1dXVyVwWEQhIzgkIjIvNjY2NmhbUSMhIzgkIjFMTExMTGVIQSEjNyQiMi82NjY2aCl6PyEjOCQiMVdXV1dXcE4+ISM3JCIyR0xMTEwkMyh6IiEjOCQiMnJ4eHh4RlNtIiEjOCQiMnV4eHh4RmxgIiEjOCQiMUxMTExMZTk5ISM3JCIxV1dXV1c+KUgiISM3JCIyMTY2NjZodD0iISM4JCIxTExMTEwzIzMiISM3JCIxMTY2NjZoQikqISM4JCIxVFdXV1clPikpKSEjOCQiMExMTExMZSp6ISM3JCIxdXh4eHhGbHIhIzgkIjF1eHh4eEYhUichIzgkIjFKTExMTCQzbiYhIzgkIjJCV1dXV1dwKyYhIzkkIjA2NjY2NicpUiUhIzckIjJFTExMTExlJVEhIzkkIjA2NjY2NidbTCEjNyQiMVZXV1dXJXAhSCEjODc7JCIxTExMTExlLUQhIzckIjFiYmJiYklVQiEjNyQiMUxMTExMZSg9IyEjNyQiMmptbW1tOyVRPyEjOCQiMWJiYmJiIVsqPSEjNyQiMikqKioqKioqKipcbnYiISM4JCInREM7ISIjJCIyYWJiYmIwdFwiISM4JCIybG1tbW07ZlAiISM4JCIyTExMTEwkM2c3ISM4JCIyYWJiYmIwKVw2ISM4JCIxTExMTEwzWDUhIzckIjFsbW1tbTtmJSohIzgkIjFhYmJiYjBCJikhIzgkIjEpKioqKioqKioqXFV3ISM4JCIxKSoqKioqKioqKlw8byEjOCQiMWBiYmJiMFtnISM4JCIxbG1tbW07TWAhIzgkIjJFTExMTExlbiUhIzkkIjJhYmJiYmJJMiUhIzkkIjJQTExMTExlXyQhIzkkIjFubW1tbTtNSSEjOCQiMWNiYmJiMClmIyEjOCQiJnZAIyEiIyQiMigqKioqKioqKioqXCMqPSEjOTc7JCIxTExMTExlW0AhIzckIjFBQUFBQSgqKio+ISM3JCIyam1tbW07cCY9ISM4JCIxbW1tbW1UPjwhIzckIjI8QUFBQXN1ZSIhIzgkIjFMTExMTDNoOSEjNyQiJ0RTOCEiIyQiMkFBQUFBc1xBIiEjOCQiMikqKioqKioqKipcXzYiISM4JCIyTExMTEwkMzY1ISM4JCIxQUFBQUFzQyIqISM4JCIxbG1tbW07JT4pISM4JCIxbG1tbW07PnQhIzgkIjFBQUFBQXMqXCchIzgkIjFMTExMTCRldCYhIzgkIjIlKioqKioqKioqKlxGXSEjOSQiMUFBQUFBc3VWISM4JCImdngkISIjJCIxTExMTEwkZUIkISM4JCIyQkFBQUFBKFxGISM5JCIycG1tbW1tIj5CISM5JCIybW1tbW1tVCU+ISM5JCIyQkFBQUFBWmkiISM5JCIyTUxMTExMM08iISM5JCIyKCoqKioqKioqKipcXzYhIzk3OyQiMnl4eHh4Rj8jPSEjOCQiMUxMTExMMyZvIiEjNyQiMlZXV1dXJXBgOiEjOCQiMTY2NjY2J3lVIiEjNyQiMkdMTExMJGUyOCEjOCQiMjQ2NjY2aEc+IiEjOCQiMlZXV1dXJXAkMyIhIzgkIjFNTExMTCQzISkqISM4JCIxeXh4eHhGPykpISM4JCIxeHh4eHhGJip5ISM4JCIxTExMTEwkZS0oISM4JCIxWVdXV1clPkAnISM4JCIxNjY2NjZoYGEhIzgkIjFMTExMTCQzdiUhIzgkIjE2NjY2NmguVCEjOCQiMllXV1dXVz5eJCEjOSQiMUxMTExMJGUoSCEjOCQiMnV4eHh4eF9cIyEjOSQiMnh4eHh4eC0yIyEjOSQiMlBMTExMTDNxIiEjOSQiMllXV1dXV3BRIiEjOSQiMjk2NjY2NidHNiEjOSQiMU9MTExMTGUjKiEjOSQiMTY2NjY2Nid5KCEjOSQiMVdXV1dXV3BvISM5NzskIjJubW1tbTtIXyIhIzgkIjImKSkpKSkpKSkpUXdSIiEjOCQiMmptbW1tO3pGIiEjOCQiJ3ZqNiEiIyQiMigpKSkpKSkpKSlRXjAiISM4JCIxTExMTEwkM18qISM4JCIxTUxMTEwkZWEpISM4JCIxKikpKSkpKSkpKVFFdyEjOCQiJkR3JyEiIyQiMW5tbW1tO2FmISM4JCIxKikpKSkpKSkpKVEsXyEjOCQiMnVtbW1tbVRdJSEjOSQiJkQnUSEiIyQiMicpKSkpKSkpKSkpUXdLISM5JCIyRUxMTExMZXUjISM5JCIyTUxMTExMM0YjISM5JCIyKikpKSkpKSkpKSlRXj0hIzkkIiZ2WyIhIiMkIjJtbW1tbW0iejYhIzkkIjAqKSkpKSkpKSkpUUUqISM4JCIxcG1tbW1tIkgoISM5JCIyOSsrKysrXShlISM6JCIxKikpKSkpKSkpKSlRLCYhIzkkIjFMTExMTEwzWiEjOSQiMUxMTExMTGVcISM5NzskIidEXjchIiMkIjIoKSkpKSkpKSkpUXc4IiEjOCQiMkxMTExMJGVINSEjOCQiMUxMTExMJDNGKiEjOCQiMSgpKSkpKSkpKSlRLCQpISM4JCIxKSoqKioqKioqKlwoUSghIzgkIjFubW1tbTtIbCEjOCQiMSopKSkpKSkpKSlRRWQhIzgkIjJ1bW1tbW0ielwhIzkkIiZ2RyUhIiMkIjEqKSkpKSkpKSkpUV5PISM4JCIyUExMTExMMzIkISM5JCIySkxMTExMZWEjISM5JCIyJykpKSkpKSkpKSlRdz8hIzkkIiZEbSIhIiMkIjJtbW1tbW1USSIhIzkkIjIoKSkpKSkpKSkpKVEsNSEjOSQiMW1tbW1tbVR2ISM5JCIlRGMhIiMkIjEqKSkpKSkpKSkpKVFFJSEjOSQiMU1MTExMTGVNISM5JCIyT0xMTExMJDNLISM6JCIxKCkpKSkpKSkpKSlRXiQhIzkkIjInKSoqKioqKioqKipcUCUhIzokIjFtbW1tbW0ieiYhIzk3OyQiMnh4eHh4RnErIiEjOCQiMUpMTExMJDMwKiEjOCQiMVZXV1dXJXAzKSEjOCQiMTY2NjY2aHlyISM4JCIxTExMTEwkZUsnISM4JCIxNjY2NjZoR2IhIzgkIjJZV1dXV1dweSUhIzkkIjFNTExMTCQzNSUhIzgkIjF5eHh4eEZxTSEjOCQiMnh4eHh4eF8qRyEjOSQiMkpMTExMTGVQIyEjOSQiMllXV1dXVz4iPiEjOSQiMjY2NjY2Nk9dIiEjOSQiMUxMTExMJDM6IiEjOCQiMDY2NjY2aGApISM4JCIxV1dXV1dXPmghIzkkIjFLTExMTExlVSEjOSQiMnZ4eHh4eEYmSCEjOiQiMXl4eHh4eC1BISM5JCIyT0xMTExMJDM/ISM6JCIyT1dXV1dXJXBCISM6JCIwNjY2NjZoRyQhIzgkIjBMTExMTCRlWiEjOCQiMTQ2NjY2Nid5JyEjOSQiMVdXV1dXV3AkKiEjOTc7JCImRCF6ISIjJCIxQUFBQUFzKipwISM4JCIxKSoqKioqKioqKlxfaCEjOCQiMUxMTExMJDNPJiEjOCQiMUFBQUFBc0NZISM4JCIyam1tbW1tVCVSISM5JCIyam1tbW1tIj5MISM5JCIxQkFBQUFzXEYhIzgkIjJKTExMTExlQiMhIzkkIiZ2eCIhIiMkIjFBQUFBQXN1OCEjOCQiMigqKioqKioqKioqXEY1ISM5JCIxSkxMTExMZXQhIzkkIjFAQUFBQUEoKlwhIzkkIjJkbW1tbW07PiQhIzokIjJkbW1tbW07JT4hIzokIjI9QUFBQUFzQyIhIzokIjJNTExMTEwkMzYhIzokIjItKysrKytdXyIhIzokIjI5QUFBQUFzXCMhIzokIiVEUyEiIyQiMUpMTExMTDNoISM5JCIxPEFBQUFBWigpISM5JCIybm1tbW1tVD4iISM5JCIycG1tbW1tInA6ISM5NzskIjFubW1tbTs0ZyEjOCQiMWJiYmJiMEJfISM4JCIxKioqKioqKioqKlwjXCUhIzgkIjIlKioqKioqKioqKlw8USEjOSQiMmFiYmJiYiEpPiQhIzkkIjJqbW1tbW1UaiMhIzkkIjJNTExMTExlNyMhIzkkIjJhYmJiYmJJbiIhIzkkIjJKTExMTExlRiIhIzkkIjFtbW1tbW1UJCohIzkkIjFhYmJiYmIhWychIzkkIjInKSoqKioqKioqKipcPCUhIzokIjInKSoqKioqKioqKipcVSMhIzokIjJZYmJiYmIwQiIhIzokIjJEbW1tbW1tImYhIzskIjFOTExMTEwkMyYhIzokIjFtYmJiYmIwKSohIzokIjJWTExMTEwkMz8hIzokIjFubW1tbW0iZiQhIzkkIjFjYmJiYmJJZCEjOSQiMSoqKioqKioqKioqXFUpISM5JCIyKioqKioqKioqKipcbjYhIzkkIjJeYmJiYmIhWzohIzkkIjJtbW1tbW1UKT4hIzkkIjJNTExMTExlWiMhIzk3OyQiMnh4eHh4eC1SJSEjOSQiMUxMTExMJDNzJCEjOCQiMk1XV1dXV3A1JCEjOSQiMi42NjY2NidbRCEjOSQiMUxMTExMJGUvIyEjOCQiMjE2NjY2NicpZiIhIzkkIjJUV1dXV1dwPyIhIzkkIjFKTExMTEwzKCkhIzkkIjF4eHh4eHgtZiEjOSQiMmt4eHh4eEZsJCEjOiQiMUtMTExMTGU+ISM5JCIxUVdXV1dXJT4pISM6JCIxMzY2NjY2aEIhIzokIjJgTExMTExMMyMhIzskIjE8NjY2NjZodCEjOiQiMmRXV1dXVyU+PSEjOiQiMU9MTExMTGVNISM5JCIxInl4eHh4RmwmISM5JCIweXh4eHhGUykhIzgkIjJPTExMTEwzPCIhIzkkIjJWV1dXV1dwYiIhIzkkIjI2NjY2NjYnKSo+ISM5JCIyRUxMTExMZVwjISM5JCIyLjY2NjY2J1tJISM5JCIyWVdXV1dXcGwkISM5NzskIjJFTExMTExlLyQhIzkkIjFiYmJiYjAkXCMhIzgkIjJCTExMTExlKj4hIzkkIjFtbW1tbTthOiEjOCQiMWJiYmJiMG82ISM4JCIxKCoqKioqKioqKipcUCkhIzkkIjInKSoqKioqKioqKipcaSYhIzokIjFjYmJiYmJJTSEjOSQiMW1tbW1tbSJ6IiEjOSQiMUhMTExMTCQzKCEjOiQiMlViYmJiYmIhPSEjOyQiMnZMTExMTEwzIyEjOyQiMXVtbW1tbTt6ISM6JCIya2JiYmJiMCQ+ISM6JCIyOSsrKysrXWkkISM6JCIxLysrKysrdmUhIzkkIjBjYmJiYjBvKSEjOCQiMnJtbW1tbVQ/IiEjOSQiMlBMTExMTGVmIiEjOSQiMmRiYmJiYkkvIyEjOSQiMkpMTExMTGVhIyEjOSQiMmptbW1tbVQ1JCEjOSQiMWJiYmJiMD1QISM4JCImdlElISIjJCIyMSsrKysrRDYmISM5NzskIjFMTExMTCRlKD4hIzgkIjI8QUFBQUEoUjohIzkkIjFtbW1tbTtmNiEjOCQiMWhtbW1tbVQkKSEjOSQiMic9QUFBQUFaYyEjOiQiMUpMTExMTDNOISM5JCIyJCoqKioqKioqKioqXCM+ISM6JCIwQUFBQUFBKCopISM5JCIkRCUhIiMkIjFQTExMTEwkMyYhIzokIjJLQUFBQUFzOSIhIzokIjInb21tbW1tVEIhIzokIjInb21tbW1tIjQlISM6JCIya0FBQUFBc1InISM6JCIxUExMTExMZSMqISM5JCIyMisrKysrdkUiISM5JCIySEFBQUFBWm0iISM5JCIyMSsrKysrdjYjISM5JCIyUExMTExMZWkjISM5JCIyQkFBQUFBKCo9JCEjOSQiMW5tbW1tOzRRISM4JCIyam1tbW1tVFslISM5JCIxQUFBQUFzOV8hIzgkIjFNTExMTCQzKychIzgkIiZEJW8hIiM3OyQiMnV4eHh4eC09IiEjOSQiMUhMTExMTDMnKSEjOSQiMVJXV1dXV3BmISM5JCIxMjY2NjY2JylRISM5JCIyPUxMTExMJGVCISM6JCIyMjY2NjY2aFEiISM6JCIxW1dXV1dXJXAqISM6JCIyVExMTExMJDM2ISM6JCIxenh4eHh4LT0hIzkkIjIueXh4eHhGMCQhIzokIjFQTExMTExlWyEjOSQiMFhXV1dXJT5zISM4JCIyPDY2NjY2TywiISM5JCIxTUxMTEwkM08iISM4JCIyOTY2NjY2T3ciISM5JCIxWFdXV1clPkEjISM4JCIxTkxMTEwkZXQjISM4JCIxeXh4eHhGMEwhIzgkIjF5eHh4eEZJUiEjOCQiMU1MTExMJDNoJSEjOCQiMVhXV1dXJXBNJiEjOCQiMTY2NjY2aFFoISM4JCIxTExMTEwkZSlwISM4JCIxNzY2NjZoKSl5ISM4JCIxWFdXV1clcCUpKSEjODc7JCIxa21tbW1tImYnISM5JCIxJykpKSkpKSkpKSlRYyUhIzkkIjJWbW1tbW07NCQhIzokIjInKSoqKioqKioqKipcPCMhIzokIjIkKikpKSkpKSkpKSlRIj0hIzokIjFOTExMTEwzPyEjOSQiMmRMTExMTCRlRiEjOiQiMSMqKSkpKSkpKSkpUTElISM5JCIxLisrKysrRGYhIzkkIjFybW1tbW1UJCkhIzkkIjInKikpKSkpKSkpKVFKNiEjOSQiMW5tbW1tOyVbIiEjOCQiMjErKysrK0QqPSEjOSQiMi4qKSkpKSkpKSlRY0IhIzkkIjFNTExMTCRlKEchIzgkIjJQTExMTEwzWCQhIzkkIjAqKSkpKSkpKSlRIjMlISM3JCIxLCsrKytdblohIzgkIjFubW1tbTs0YiEjOCQiMCopKSkpKSkpKVExaiEjNyQiMW5tbW1tO2ZyISM4JCImdjEpISIjJCIwKikpKSkpKSkpUUohKiEjNyQiMk5MTExMJDMwNSEjOCQiMk5MTExMJGU3NiEjODc7JCIxKioqKioqKioqKipcNyUhIzkkIjInKSkpKSkpKSkpKSlRRSQhIzokIjFMTExMTExlSCEjOSQiMU1MTExMTDNLISM5JCIwKikpKSkpKSkpKVEsJSEjOCQiMjkrKysrK11QJiEjOiQiMW5tbW1tbSJIKCEjOSQiMSgpKSkpKSkpKSkpUXcqISM5JCIya21tbW1tIno3ISM5JCImdmoiISIjJCIxKikpKSkpKSkpKVFePyEjOCQiMkpMTExMTDNfIyEjOSQiMlBMTExMTGUvJCEjOSQiMSopKSkpKSkpKSlRRU8hIzgkIiZERSUhIiMkIjFubW1tbTthXCEjOCQiMSopKSkpKSkpKSlRLGQhIzgkIjFubW1tbTsvbCEjOCQiJkRPKCEiIyQiMSopKSkpKSkpKSlRdyMpISM4JCIxTExMTEwkZUMqISM4JCIyTkxMTEwkM0Y1ISM4JCIyKCkpKSkpKSkpKVFeOCIhIzgkIid2WzchIiMkIjJubW1tbTt6TyIhIzg3OyQiMid5eHh4eHgtVyEjOiQiMU1MTExMTDNaISM5JCIya1dXV1dXJXBiISM6JCIxOTY2NjY2JylwISM5JCIxTUxMTExMZSopISM5JCIyNjY2NjY2J1s2ISM5JCIyVldXV1dXcFgiISM5JCIxTExMTEwkMyM9ISM4JCIydXh4eHh4LUMjISM5JCIyeHh4eHh4X3IjISM5JCIyUExMTExMZUMkISM5JCIxV1dXV1clPiRRISM4JCIxNjY2NjZodFchIzgkIjFNTExMTCQzPCYhIzgkIjE3NjY2NmhCZiEjOCQiMVhXV1dXJT50JyEjOCQiMU5MTExMJGVmKCEjOCQiMXl4eHh4RjomKSEjOCQiMXl4eHh4RiFcKiEjOCQiMkxMTExMJDNfNSEjOCQiMlZXV1dXJXBnNiEjOCQiMTY2NjY2J1tGIiEjNyQiMkxMTExMJGUlUiIhIzgkIjI0NjY2NmgpPjohIzgkIjJWV1dXVyVwXTshIzg3OyQiJUR1ISIjJCIxQEFBQUFBKCopKSEjOSQiJkQ0IiEiIyQiMkpMTExMTDNOIiEjOSQiMUFBQUFBc2s7ISM4JCIyam1tbW1tVC4jISM5JCIyam1tbW1tImZDISM5JCIxQEFBQUFzUkghIzgkIjJFTExMTExlWiQhIzkkIjEqKioqKioqKioqXG5TISM4JCIxQUFBQUFzOVohIzgkIiZ2VCYhIiMkIjFMTExMTCRlPCchIzgkIjFBQUFBQXMqKXAhIzgkIjFtbW1tbTtmeSEjOCQiMW5tbW1tOyV5KSEjOCQiMUNBQUFBc2soKiEjOCQiMkxMTExMJDMhMyIhIzgkIidEKj0iISIjJCIyQUFBQUFzUkkiISM4JCIyKSoqKioqKioqKlxVVSIhIzgkIjFMTExMTDNdOiEjNyQiMUFBQUFBWiJvIiEjNyQiMmxtbW1tOyU9PSEjOCQiMmxtbW1tOzQnPiEjODc7JCIyam1tbW1tIj44ISM5JCIyYWJiYmJiSWUiISM5JCImRCE+ISIjJCImdkYjISIjJCIyYWJiYmJiITNGISM5JCIyam1tbW1tVD4kISM5JCIyRUxMTExMZXQkISM5JCIyVmJiYmJiSUwlISM5JCIxTExMTEwkZSlcISM4JCIxbG1tbW07JXAmISM4JCIxYmJiYmIwZWshIzgkIiZ2RighIiMkIiZEOikhIiMkIjFhYmJiYjAkMyohIzgkIjJtbW1tbTtwKyIhIzgkIjJMTExMTCQzNjYhIzgkIjJkYmJiYjAzQSIhIzgkIjJOTExMTCQzTzghIzgkIjJsbW1tbTtwWCIhIzgkIjJhYmJiYjBMZSIhIzgkIjIpKioqKioqKioqXF9yIiEjOCQiMikqKioqKioqKipcRiY9ISM4JCIyX2JiYmIwZSo+ISM4JCIybG1tbW07VzkjISM4JCIyTkxMTEwkZSlIIyEjODc7JCIyanh4eHh4LTwjISM5JCIxS0xMTEwkM2IjISM4JCIyTVdXV1dXcClIISM5JCIxNDY2NjZoeU0hIzgkIjI5TExMTExlLSUhIzkkIjA2NjY2NidHWSEjNyQiMVVXV1dXJXBHJiEjOCQiMExMTExMMysnISM3JCIxd3h4eHhGcW4hIzgkIjF3eHh4eEYmZighIzgkIjBMTExMTGVaKSEjNyQiMVZXV1dXJT5UKiEjOCQiMjQ2NjY2aC4vIiEjOCQiMUxMTExMM1g2ISM3JCIyNDY2NjZoYEQiISM4JCIyVldXV1clPnI4ISM4JCIyR0xMTEwkZSNcIiEjOCQiMnV4eHh4RiY+OyEjOCQiMnV4eHh4Rj92IiEjOCQiMkdMTExMJDMhKj0hIzgkIjJUV1dXVyVwTD8hIzgkIjIvNjY2NmhHPSMhIzgkIjFLTExMTGVQQiEjNyQiMTY2NjY2J3lcIyEjNyQiMlZXV1dXJXBqRSEjODc7JCIyOUxMTExMZUgkISM5JCIxYWJiYmIwJHokISM4JCIxS0xMTEwkZU0lISM4JCIxa21tbW07YVwhIzgkIjFgYmJiYjA9YyEjOCQiMSoqKioqKioqKipcUGohIzgkIjEpKioqKioqKioqXDdyISM4JCIxX2JiYmIwVnohIzgkIjFsbW1tbTtIKSkhIzgkIjBMTExMTDN4KiEjNyQiMl9iYmJiMG8yIiEjOCQiMUxMTExMMyM9IiEjNyQiMmptbW1tO0hIIiEjOCQiMWJiYmJiSTQ5ISM3JCIyKSoqKioqKioqKlw3YCIhIzgkIjImKioqKioqKioqXChlOyEjOCQiMl9iYmJiMD16IiEjOCQiMmxtbW1tOy8kPiEjOCQiMUxMTExMZXU/ISM3JCIyWmJiYmIwVkEjISM4JCIxTExMTExlekIhIzckIjJjbW1tbTsvYSMhIzgkIjJaYmJiYjBvcSMhIzgkIjEqKioqKioqKipcKHlHISM3JCIxKioqKioqKioqXGkwJCEjNzc7JCIyRUxMTExMZXAlISM5JCIxQkFBQUFzNGAhIzgkIjFubW1tbTt6ZiEjOCQiMW1tbW1tOy9uISM4JCIxQ0FBQUFzJVsoISM4JCIxTUxMTEwkM0spISM4JCIxKioqKioqKioqKlw3IyohIzgkIjJAQUFBQXNmLCIhIzgkIidEOzYhIiMkIjFMTExMTDNBNyEjNyQiMUFBQUFBWkw4ISM3JCIybm1tbW07L1giISM4JCIybG1tbW07SGQiISM4JCIxQUFBQUEoNHEiISM3JCIyTkxMTEwkZU09ISM4JCIndnQ+ISIjJCIxQUFBQUFaPUAhIzckIid2b0EhIiMkIjFMTExMTGVDQyEjNyQiMjxBQUFBc2ZlIyEjOCQiMW1tbW1tIkh2IyEjNyQiMW1tbW1tVERIISM3JCIyOEFBQUFzTTUkISM4JCIxTExMTEwzKEckISM3JCIyJioqKioqKioqKlxpWiQhIzg3OyQiMXl4eHh4RnFqISM4JCIxTkxMTEwkMzUoISM4JCIxWVdXV1clcCl5ISM4JCIxNjY2NjZoRygpISM4JCIxTkxMTEwkZWkqISM4JCIyODY2NjZoeTAiISM4JCIyVldXV1clcGU2ISM4JCIyTExMTEwkM2w3ISM4JCIxeXh4eHgteDghIzckIjJ5eHh4eEZYXCIhIzgkIjJMTExMTCRlPDshIzgkIjJbV1dXVyU+WTwhIzgkIjI4NjY2NmguKT0hIzgkIjJMTExMTCQzPz8hIzgkIjI4NjY2NmhgOyMhIzgkIjJWV1dXVyU+O0IhIzgkIjJOTExMTCRlc0MhIzgkIjJ1eHh4eEZYaiMhIzgkIjJ1eHh4eEY/IUchIzgkIjFMTExMTDN2SCEjNyQiMVdXV1dXcGBKISM3JCIxNjY2NjYneUwkISM3JCIyRUxMTEwkZUZOISM4JCIyLzY2NjZoR3MkISM4JCIyVldXV1clcEJSISM4LSUrUFJPSkVDVElPTkc2LCQiMXduKD5xKioqKipcISM7JCExLyRcM3BuNTIoISM7JCIxd24oPnEqKioqKlwhIzskIjF3big+cSoqKioqXCEjOyQiMS8kXDNwbjUyKCEjOyQiMXduKD5xKioqKipcISM7JCExLyRcM3BuNTIoISM7JCIiISEiIiQiMS8kXDNwbjUyKCEjOyQiIzUhIiItJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiQrIiEiIi0lKUJPVU5EU19ZRzYjJCIkKyIhIiItJS1CT1VORFNfV0lEVEhHNiMkIiUrUiEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiUrUiEiIi0lKUNISUxEUkVORzYi
EXERCISE: Consider the points
newpts := [seq([i, 1.5*i+1], i = -1 .. 5), [8, 1.5*8+.5], seq([i, 1.5*i+1], i = 9 .. 11)];
a. Find a polynomial passing through those points.
b. Fit the points with a line using the least squares method. Use the graph of the previous exercise to visualize the line you found, and see how "close" the line and the points are.
c. Fit the points with a quadratic function using the least squares method.
d. Find a piecewise linear function passing through all the points.
e. Plot the points and all the functions together.
f. Compare the results. In particular, compute the value of the function you found in each case when x=12.3
Problems 11 and 13
Fitting points to a circle
Here is a program that produces a list of points approximately in a circle. You do not have to worry about the details of this program (but you are wellcome to do it if you want), just run it to obtain the points.
Seed:=readlib(randomize)():
with(stats):
myrand := (a,b) -> evalf(a+(b-a)*(rand()/10^12)):
points_wanted := 50:
circlePts := proc()
local fuzz, a, b, r,s,x,i;
a:= myrand(-3,3);
b:= myrand(-3,3);
r:= myrand(1,2);
fuzz := [ random[normald[0,0.1]](points_wanted)];
x := [ random[uniform[a-r,a+r]](points_wanted)];
for i from 1 to points_wanted do if (myrand(-10,10) > 0) then s[i]:=1 else s[i]:=-1 fi:od:
return [seq([x[i],b+s[i]*sqrt(r^2-(x[i]-a)^2)+fuzz[i]],i=1..points_wanted)];
end:
myCirclePts:=circlePts():
JSFH