MAT331, Feb 25th, 2010 If a curve fits, wear it (Chapter 2 of the online book) restart; pts:=[[9.461625032, -15.00000412], [-4.026835626, 5.229076781], [-8.417983819, 17.28331107], [7.335537549, -10.75702150], [.667622007, 1.352242696], [-3.642662167, 8.570158849], [4.108941669, -6.649973580], [-3.592756472, 7.661525316], [0.121773384e-1, -.5948721361], [3.588065808, -7.386147022], [9.371019578, -17.50451275], [6.303804064, -11.16461096], [9.396335873, -15.62444142], [-.807790200, 2.555503241], [.782176902, -4.501842488], [3.640727999, -6.449402329], [1.668897236, -3.014773741], [3.704029489, -3.221704164], [-11.29937335, 17.53901720], [8.475940602, -17.91640911], [-3.666884970, 8.440428354], [-5.617368657, 12.00375535], [-5.626834452, 7.337908700], [-2.138751424, 5.625807468], [-2.268585475, 4.613703098]]: with(plots): distancePointLine:=P->(m*P[1]+b-P[2])^2;

Zio2I0kiUEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlKiQsKComSSJtR0YlIiIiJjkkNiNGLUYtRi1JImJHRiVGLSZGLzYjIiIjISIiRjRGJUYlRiU=

distancePointsLine:=pts->sum(distancePointLine(pts[i]),i=1..nops(pts));

Zio2I0kkcHRzRzYiRiU2JEkpb3BlcmF0b3JHRiVJJmFycm93R0YlRiUtSSRzdW1HRiU2JC1JMmRpc3RhbmNlUG9pbnRMaW5lR0YlNiMmOSQ2I0kiaUdGJS9GMjsiIiItSSVub3BzRyUqcHJvdGVjdGVkRzYjRjBGJUYlRiU=

DD:=distancePointsLine(pts): simplify(DD);

LC5JIm1HNiIkIisuei54RyEiJ0kiYkdGJCQiK1JhbDlWISIpKiRGIyIiIyQiK11MZyY+KSEiKComRiMiIiJGKEYyJCIrMlxAI1skRisqJEYoRi0kIiNEIiIhJCIrS0VrOEVGJ0Yy

eq1:=diff(DD,m); eq2:=diff(DD,b);

LChJIm1HNiIkIitwMTdSOyEiJ0kiYkdGJCQiKzRcQCNbJCEiKSQiKy16LnhHRiciIiI=

LChJIm1HNiIkIis0XEAjWyQhIilJImJHRiQiI10kIitPYWw5VkYnIiIi

plot3d(DD, m=-2..-1,b=-0.8..2,axes=boxed,style=patchcontour);

LSUnUExPVDNERzYoLSUlR1JJREc2JTskISM/ISIiJCEjNSEiIjskISIpISIiJCIyJCoqKioqKioqKioqKioqPiEjOzc7NzskIjEqUXYqNE8xXTwhIzgkIjIoKj4+bm8mKWVuIiEjOSQiMmtjPSE+TF4zOyEjOSQiMiMpW3RvXVl6YSIhIzkkIjFtUkddXz0lXCIhIzgkIjImKioqXCNcJkhzVyIhIzkkIjIpKWVyUFN6cVMiISM5JCIyS3RbUSJbdHQ4ISM5JCIyRlYiW3pkPlo4ISM5JCIyKSlvcDFJaXVLIiEjOSQiMiUqXDh1UE1YSiIhIzkkIjJlJ0dyND9UMzghIzkkIjImKXludj4mNDQ4ISM5JCIxbSN5NCVSZTs4ISM4JCIyJCpIVypSI3kzTCIhIzkkIjIjKSllWSU0eT5OIiEjOSQiMUxJYS9OKSl6OCEjOCQiMUxkPHFXZjk5ISM4JCIyKSkpUk8iKjQ2YzkhIzkkIi95NW9JVi86ISM2JCIyazsyL3FnJmY6ISM5JCIyJikzaSMpKVFcQDshIzkkIjJtY3M7akstcCIhIzkkIjImKmZRMSRweGw8ISM5JCIxKT1nXnlFIls9ISM4NzskIjF4MUExbV4mZSIhIzgkIjEqKj5LRTkuODohIzgkIjJrKCl5PiE9Tlo5ISM5JCIyKDQ4Pkx4WilRIiEjOSQiMicpSGYqPiM0a0wiISM5JCIxVkdHaWk5IkgiISM4JCIyTCU+O2cpKW9fNyEjOSQiMikpZidmOHEuQDchIzkkIjIoNG9lQTI+Jz4iISM5JCIybGRLcikqXCJ5NiEjOSQiMigpKlFCMlsicDsiISM5JCIyaXghKkc9JltpNiEjOSQiMic0SzU5NidbOyIhIzkkIjIpKT5yM2dVUzwiISM5JCIyS3UlPlYnSCs+IiEjOSQiMVZRMlRBI0dAIiEjOCQiMicpXDNYUj9DQyIhIzkkIjIpNCgpXC5UIyl5NyEjOSQiMm5aVyFvTC5BOCEjOSQiMikpelgiKT1bP1AiISM5JCIyam4tUWNvKUc5ISM5JCIwNjpdXCVcI1wiISM3JCIyKCk0JHkiKWYjSGMiISM5JCIyR2sxVC1qLGsiISM5JCIxVmQpPmkwVXMiISM4NzskIjFSPWFlbFVcOSEjOCQiMXQxK0FUank4ISM4JCIxaF0sVHNrOTghIzgkIjJkKyZlOmZZZDchIzkkIjJkXTVkOSE0MjchIzkkIjI6YyJSSio+TjsiISM5JCIyQjxHRUZibjciISM5JCIxUi5VcGh6JzQiISM4JCIyNzFvPGlVTzIiISM5JCIyIlI4bkhZSGQ1ISM5JCIyRDxJSj5feC8iISM5JCIyN2NXQEo6XS8iISM5JCIyY105bilSM1w1ISM5JCIyZCtTb0BlKmY1ISM5JCIyNzFARCtRdzIiISM5JCIyQG5kUE1CQDUiISM5JCIyKlEpXDBDOU04IiEjOSQiMjljKCpHcDU6PCIhIzkkIjIkUjMhM3E3a0AiISM5JCIyRG5mVUU/Im83ISM5JCIyPDF1S1FMbUsiISM5JCIyYStXeTBfPlIiISM5JCIxMSZwekd3U1kiISM4JCIyN2NdTzIxSWEiISM5JCIyRjwoKVtUVChHOyEjOTc7JCIyZCgpUXBZJHpUOCEjOSQiMi5BYlB4JHBzNyEjOSQiMio+cjdPJypSNTchIzkkIjJjZGFTMDZcOiIhIzkkIjJtZVB2LUdpNSIhIzkkIjJQOndsYl1WMSIhIzkkIjJkRnE2a3kjSDUhIzkkIjJNJio+OEc3NSsiISM5JCIxbT1EcVpeJnoqISM5JCIxZShmR0dpKlsnKiEjOSQiL0stXmBZcSYqISM3JCIxKT5VWihSLWcmKiEjOSQiMWRuLGEiUXdoKiEjOSQiMW5vJSkpKXlJVigqISM5JCIxS0RCekouUCoqISM5JCIyY1A8RFMiKSk+NSEjOSQiMkswbkUvbEcwIiEjOSQiMShHcyRRVWwjNCIhIzgkIjF3SWoqKSpbI1I2ISM4JCIyLFVca0hcRT4iISM5JCIyLUtAKWVeJkdEIiEjOSQiMXcoW25kbSk+OCEjOCQiMnV5Si1iJG8kUiIhIzkkIjJNTnEjemdJdTkhIzkkIjJkWmtRO008YyIhIzk3OyQiMnJ5NjlMPEVFIiEjOSQiMkZrJmUiUTVfPiIhIzkkIjJOMDp0KSozWTgiISM5JCIwLSsnW0oiMzMiISM3JCIyRGFTYSdHI1EuIiEjOSQiMTJpT3k4UU8qKiEjOSQiMU5EKXlsKmUtJyohIzkkIjFDVyZIXGBvTCohIzkkIjFzPWUkKUc8UiIqISM5JCIwKFt3SHlhNCEqISM4JCIxQ01dSiR5eiUqKSEjOSQiMVB2eilRa1cmKikhIzkkIjEuc2ssZytIISohIzkkIjFEQzBxSmdyIiohIzkkIjEwSywlKmVEI1EqISM5JCIxUSZITjxrNG0qISM5JCIyRjlnMyFHeCs1ISM5JCIycilHIypSWkRVNSEjOSQiMnQ9VFhCVTA0IiEjOSQiMkcvOlpHTmM5IiEjOSQiMlRYVy8qUWAyNyEjOSQiMjFVSDwwUWlGIiEjOSQiMkYlKnAmb3h1XjghIzkkIjIyLW00L2pTViIhIzkkIjJSbDwqb1E9QjohIzk3OyQiMkRkZz46KSo9QCIhIzkkIjIjUj5cWFI9WTYhIzkkIjI3Jyl5WEh2czMiISM5JCIxUjhBKj5zXi4iISM4JCIxQlA+JWZZKCkqKSohIzkkIjE4J0g8dkVRXiohIzkkIjFlNSNbWWlwPiohIzkkIjAxb010YCJbKikhIzgkIjE6MW5kMFNuKCkhIzkkIjFCKEd1JEhxYScpISM5JCIxJlJVRigzMTUnKSEjOSQiMTs7aGpWWkwnKSEjOSQiMSJSTytUVlxzKSEjOSQiMUNuLDchb1cpKSkhIzkkIjE6RWJwIls/NiohIzkkIjFjU2sjKVFvMiUqISM5JCIwMSJIXl5QcigqISM4JCIyO09cdj43Li0iISM5JCIySDxEYlYjSHE1ISM5JCIyIlJsMEgjeXE3IiEjOSQiMjlZViJ5JnAxPiIhIzkkIjImUmZ5I1ttNUUiISM5JCIxdFIpSCUqbyNROCEjOCQiMjpjUChlcEZBOSEjOSQiMTFuL0kwNDg6ISM4NzskIjJDQiZlR2ZqKj0iISM5JCIyLDZ1YVk5YzciISM5JCIyTGE9emIpUm81ISM5JCIyQmA9Zj8pKXosIiEjOSQiMW4ydSU0TVF1KiEjOSQiMW88JmVvVGVQKiEjOSQiMUEkPUQkWyFmMiohIzkkIjFOL3VNTi1XKSkhIzkkIjEtIj1Eeig+IW8pISM5JCIxQTgmZWdGV2UpISM5JCIxLCx1dUhyYyYpISM5JCIxTlc9KipRMChmKSEjOSQiMUJWPXouWDAoKSEjOSQiMW8oUlpULT4pKSkhIzkkIjB4XWUrNWs3KiEjOCQiMUJ0Xl9KKCpRJSohIzkkIjFNJVJaJj1mPikqISM5JCIyLXJeN2hFby0iISM5JCIyRy4mZSNmKlx5NSEjOSQiMi4iUlpIInlwOCIhIzkkIjJOTT0+QWlBPyIhIzkkIjJBTD0qcD1OdTchIzkkIjJzKFFadHFDYDghIzkkIjJvKFxlS3klKlE5ISM5JCIyQWpeczlhOWAiISM5NzskIjFuZEdoMSRlPiIhIzgkIjJiOks6JT5dTDYhIzkkIjIqKjRNdHh5ejIiISM5JCIyLmcicG82RUg1ISM5JCIxY2wvYzZcdCkqISM5JCIxbkV0IT1FQ18qISM5JCIxTFYoNHc7JVIjKiEjOSQiMWU6eCcqR1lDISohIzkkIjFOVjcpZWt2KCkpISM5JCIwbktdJD1zKXopISM4JCIxZmxcUFkkenkpISM5JCIuO2IqSD9YKSkhIzYkIjEsNTQ0cF9xKikhIzkkIjFkOkF5aiFSOyohIzkkIjFudyFIU1RgVSohIzkkIjFPJFxKKT4kW3YqISM5JCIxY1kqPSJ5Qjo1ISM4JCIyTyQqSDUpenpoNSEjOSQiMnJ3P2Rxal42IiEjOSQiMWNyJ2YpXEx2NiEjOCQiMjA1cDwjPUpVNyEjOSQiMjFnRUpAJTQ7OCEjOSQiMmpsUis7I28nUiIhIzkkIjFuIzNEbXZTWyIhIzgkIjJRVkswc3UjeTohIzk3OyQiMl48aStOI1tJNyEjOSQiMmAybU9QWSlwNiEjOSQiMjRgREcmZiw7NiEjOSQiMVUwYSgzIioqbzUhIzgkIjIqMzYieXhyKEc1ISM5JCIxM0JQTy1lYCoqISM5JCIxKDMqPV0jKVwobyohIzkkIjFCOWM+PVoqWyohIzkkIjE2JCpbVzRdZiQqISM5JCIxYUYoXGkmZShIKiEjOSQiMWA8LGhlcy4kKiEjOSQiMTNqZ187I3pQKiEjOSQiMFVjKCoqSDw/JiohIzgkIjEpM2lDISp6L3QqISM5JCIyM0xzZ0IlKTMrIiEjOSQiMiUzU1hQZ19ONSEjOSQiMj5DIlIlUnRwMiIhIzkkIjI4LiUpb0lFXzciISM5JCIyZFBLXHglRyE9IiEjOSQiMmJGTyYpeltAQyIhIzkkIjI4dCZweCQ9M0oiISM5JCIyQnU1Q14kSCdRIiEjOSQiMiU0OG8tVWRvOSEjOSQiMUp1XVsvbWQ6ISM4JCIyJzMiKikpXEFiYDshIzk3OyQiMiZlVyJcKjRmJEgiISM5JCIyJnBlKD13WllCIiEjOSQiMmUkR1IlMzVEPSIhIzkkIjIlZWBZaXo8UDYhIzkkIjJoViQ0J1JeJyk0IiEjOSQiMidwcUYmUUlwMSIhIzkkIjIkZWksSVwsVTUhIzkkIjEuNUpJXSFSLSIhIzgkIjEuODsnbytFLCIhIzgkIjInZXJjKCo9NTM1ISM5JCIyKnAmR1htMy8sIiEjOSQiMmlgWHEpNF8+NSEjOSQiMihlIT1eJylRYS4iISM5JCIyazhZKClIaSJlNSEjOSQiMidwKEh6RyJwKDMiISM5JCIyJGUqb0UkZS1DNiEjOSQiMkdxakgkZjtuNiEjOSQiMkorOSkpZTZyQCIhIzkkIjIpZSk+LSFHJ1FGIiEjOSQiMiVwNz1uJj51TCIhIzkkIjJsQilwKik9eTI5ISM5JCIyI2YyeG4oXFxbIiEjOSQiMVApKVIsSyMqbzohIzgkIjImcENlIT4tKGY7ISM5JCIyKmU7S05uR2Q8ISM5NzskIjJlaFVlZmNeUSIhIzkkIjFROjsxaCF6SyIhIzgkIjJmLE8/PGh1RiIhIzkkIjIiXGdZJHpAUUIiISM5JCIxUTtYcXopcD4iISM4JCIySHkjKkhxZnA7IiEjOSQiMkZbKjMiKnB0VjYhIzkkIjIjUTx1TSk+dDciISM5JCIyI1wmXFJCM3g2IiEjOSQiMmMiSHIpPS1cNiIhIzkkIjIkUT0uKnAsKj02ISM5JCIyZUoxXHcxKEg2ISM5JCIyJFxqTCdRPHQ5IiEjOSQiMVE+S2pOdHI2ISM4JCIyRzNqZUhiSD8iISM5JCIyRXlmUmUjKTRDIiEjOSQiMiNRP2hGYSJlRyIhIzkkIjIpXCk+byNRWFA4ISM5JCIybEAkZSJ5KCplUiIhIzkkIjJ5OC0+SFo2WSIhIzkkIjE7bXhkQj9MOiEjOCQiMiZcbT96SDE3OyEjOSQiMilRQT5jIkh4cCIhIzkkIjJGUUwoKSkzPyF6IiEjOSQiMSQzSW88eSUqKT0hIzg3OyQiMnlrWUc6el5dIiEjOSQiMjkzQmxTQCdcOSEjOSQiMDJiZEBwM1MiISM3JCIyWGhVMGVBKmU4ISM5JCIxOmQpM10ieUI4ISM4JCIyLFAleXdmVyZIIiEjOSQiMjhlUSMzZyJSRiIhIzkkIjIiWyRbX2YiPmY3ISM5JCIyLG44eXRzN0QiISM5JCIxW1gkZlZmLEQiISM4JCIyOik0aCpvXmVEIiEjOSQiMi0oSCUpKVxcJG83ISM5JCIyWl5JTydHbChHIiEjOSQiMllodFJ5aFBKIiEjOSQiMi9GcylmaW5ZOCEjOSQiMjlbRThIJ1InUSIhIzkkIjIoW2lMeT0jSFYiISM5JCIyM2QsNC1gaVsiISM5JCIyJ1tDLT4oKlFZOiEjOSQiMjgpKSlwcz5MODshIzkkIjAoMyQ+eXpxbyIhIzckIjE6JT1uOUx3dyIhIzgkIjJbXmhxMSMqXCY9ISM5JCIyLzxnSGFjIlw+ISM5JCIyOVE5V2RFLDAjISM5NzskIjJOYkVmbWVPbCIhIzkkIjImKVxnSG0keipmIiEjOSQiMicpKipcYkBNRmIiISM5JCIyUTAmcEIuWzc6ISM5JCIyYG0mUigpPi56OSEjOSQiMjwkPWwxIypRXzkhIzkkIjFhTlkiKT5iSzkhIzgkIjI7JDMkPUo/Jj45ISM5JCIyYm1geD4lSDg5ISM5JCIxYT9CUk8oUVQiISM4JCIyKSkqZkVPJ2U3VSIhIzkkIjInKVxiKSk9XGFWIiEjOSQiMWEwKyhIWGtYIiEjOCQiMmE7LDIncEMlWyIhIzkkIjJCTGQqelQmKT06ISM5JCIyVTBwWiZwRWc6ISM5JCIxS2o4JkcmWzM7ISM4JCIxbCJmNTw0Tm0iISM4JCIyWGJQRGhRYHMiISM5JCIyKilcciY0TyhSeiIhIzkkIjEqKjQ7aVRUcD0hIzgkIjJbMDEuRmc7Jj4hIzkkIjJkbTFTJD5yUz8hIzkkIjJDJEdFYCJwbDgjISM5JCIxYVgyRz5CUkEhIzg3OyQiMldMIzNOXmZJPSEjOSQiMHp0YShHVXk8ISM3JCIyOSEzVXJoMEw8ISM5JCIxb0wjSC0mXCVwIiEjOCQiMikqWyIpKkglUkZtIiEjOSQiMiNvXmYjUip5UDshIzkkIjEsV3c1XGs+OyEjOCQiMi4+KlslKWZJMzshIzkkIjJXYHBQaHNQZyIhIzkkIjJbVjAnKXpXZ2ciISM5JCIyLipvKipRRDc6OyEjOSQiMjohUiVcJGUrSjshIzkkIjIob2tXJ28lcGA7ISM5JCIwZi9ONCo9JG8iISM3JCIxcCM9aDAqWz48ISM4JCIyN10oR3VYZmk8ISM5JCIyLUg3IVtjXTc9ISM5JCIxTkVIeEFBcD0hIzgkIjJiYEdAWVdGJD4hIzkkIjIyKio+RD9zSSsjISM5JCIyOCtuJSlcMC0zIyEjOSQiMidvJnAqXFY5a0AhIzkkIjEicEZxdikpW0QjISM4JCIyKG84az4oUUNOIyEjOSQiMS0xInlCJXpjQyEjODc7JCIyKSopUkpnJikpZi4jISM5JCIyaCZIMVchNGIpPiEjOSQiMiN5dU8kM049JT4hIzkkIjJsYkYjeW0nXCE+ISM5JCIyKCo9VidHUSFcKD0hIzkkIjIoeVZoTWxrXj0hIzkkIjJMN1RoeiU+Tj0hIzkkIjJHVUJLaFtiIz0hIzkkIjImeTcnZSl6cUE9ISM5JCIwcGFTIkhuRT0hIzckIjJtbC55UlZ1JD0hIzkkIjIleSIzclY+XSY9ISM5JCIybERvPi4sJXo9ISM5JCIyKCopUVEjPSllNT4hIzkkIjIlel1OKSkzZVs+ISM5JCIyTiM9KSlcInpMKj4hIzkkIjFCVCdwJ0gpXC8jISM4JCIyJHo+Z1JCUi5AISM5JCIyMFImem5zZ29AISM5JCIxZFZhXnhpU0EhIzgkIjF6KVszemAlPkIhIzgkIjJ5JiozZFEmMzBDISM5JCIyMmZDaGBBdlwjISM5JCIweSY0VV93J2YjISM3JCIxQ0RpLk4iR3EjISM4NzskIjJ5XkA7JSpRKXBBISM5JCIxJyp6c29AMEBBISM4JCIyJ0grUl40MnpAISM5JCIyJj13ZypHJipROSMhIzkkIjJHdyFRJD1EYjYjISM5JCIySFk0RmpnUjQjISM5JCIyJz1QZlA7P3o/ISM5JCIyKUhOLik+WzcyIyEjOSQiMmsqKUdTSissMiMhIzkkIjE+KXpiKXp2dj8hIzgkIjEoSCdvNzdBKTMjISM4JCIyJ0gkW2AqKipbMkAhIzkkIjIpPWZjTFZjTEAhIzkkIjJNMVJ0QVdrOyMhIzkkIjJQd25tbkhoPyMhIzkkIjE+P2Iibz9FRCMhIzgkIjIuJD0qPkM8ZkkjISM5JCIyaz4oKXpOPmdPIyEjOSQiMik+IlEmSHEjSFYjISM5JCIyd2ZXbURTbV0jISM5JCIwajEkUiFmcmUjISM3JCIwQUN2UCRbdUUhIzckIjFrdEhyS2hvRiEjOCQiMWtnaT8oWyZwRyEjOCQiMT4uXkQoKkd4SCEjODc7JCIyRCNcK3ppOUtEISM5JCIyOCIqbyVcQTAmWyMhIzkkIjJjWClbdlB3V0MhIzkkIjFjTjFkM0c2QyEjOCQiMTZVPiVcLllRIyEjOCQiMkFVISlvb0paTyMhIzkkIjEqPUFeVmw7TiMhIzgkIjE2Jj4qUVpTWEIhIzgkIjEqUXMjKWZcZk0jISM4JCIxQjM9OCtJYEIhIzgkIjI5IltrJClmWG5CISM5JCIxY1ZtNHZUKVEjISM4JCIyZFhSN2YlPTtDISM5JCIyOTZxJEdzdl1DISM5JCIySEtjNVVOQFwjISM5JCIyJSozKUhwIj4uYSMhIzkkIjE3YTR0JTNgZiMhIzgkIjEiSFtDTC5ybCMhIzgkIjFDbk5aUHFERiEjOCQiMjpyP3lyNDYhRyEjOSQiMWMtJVFDQEwpRyEjOCQiMm1OOmFLUUIoSCEjOSQiMkIsWUQnNDtvSSEjOSQiMUJBQmIiKnlxSiEjOCQiMCpSWi5IQSFHJCEjNzc7JCIyMD9rQ2Q1SCNHISM5JCIvZEcnRzR2eCMhIzYkIjJsdmljYjgqUUYhIzkkIjJsTyZmIVFCcnEjISM5JCIyT2AkM2goUUBvIyEjOSQiMWNzNyhwZlJtIyEjOCQiMk9gRSgpPSdlX0UhIzkkIjFuOCllQj0hW0UhIzgkIjJrdiJmUWVEXUUhIzkkIjIycWRvKiopSGZFISM5JCIxLCN6MXJaXm4jISM4JCIxY2kwISk+IXlwIyEjOCQiMW4pKSlcIT1FRkYhIzgkIjJWLnhhPUZOdyMhIzkkIjFkMl9AIilmMUchIzgkIjFOKzc4WVpjRyEjOCQiMW9bRmdtOjhIISM4JCIxY18pSEVXbShIISM4JCIyOT9eN1VQcC8kISM5JCIxLEYyTmguQ0ohIzgkIjFjKFxXU1N6PyQhIzgkIjImb0JRSC1sKUgkISM5JCIyY2BxKTRjOydSJCEjOSQiMWRVImZhJ1srTiEjOCQiMU1OXlBJaDZPISM4NzskIjJFTioqPiM9OFVKISM5JCIxaiR5InpLVSk0JCEjOCQiMUpIIj5IPzoxJCEjOCQiMWBJP2dHVUpJISM4JCIyOHRbUyk0ODNJISM5JCIxaypcTW1XOypIISM4JCIxYG5TKSpRJz4pSCEjOCQiMSk0PiopbykzekghIzgkIjJ1KnApXC4+SSlIISM5JCIyUFg1bSRcdiQqSCEjOSQiMlBZKnkkUidINkkhIzkkIjI4LkNsU1ZjLiQhIzkkIjFgVCJbKGZ6bUkhIzgkIjJDJClmJyk0YVo1JCEjOSQiMW01MXl4XlxKISM4JCIyTiZ5LDhxMyxLISM5JCIxKT5JTiE9WWZLISM4JCIxKTQpZlxAa0NMISM4JCIyYGJANzBHbVIkISM5JCIxbDBTMyY+YVokISM4JCIwOE42Xzs1YyQhIzckIjFiX1UqMz5NbCQhIzgkIjJCJDRGOHNpX1AhIzkkIjFuQG4jKjNrZVEhIzgkIjFhKkd3N2c5KFIhIzg3OyQiMih5LmhGKyIpKlskISM5JCIyQyFwOUdVelpNISM5JCIwKSpRVSlSZTdNISM3JCIyUGgnKWVIelRRJCEjOSQiMkQhKSozaixlaUwhIzkkIjFZJltlZSd5WkwhIzgkIjJqJUc7ayYpelJMISM5JCIyRXFLISk0OydRTCEjOSQiMTgiZXU9UlVNJCEjOCQiMDNSQyR5bWNMISM3JCIxLWMoSC4tZlAkISM4JCIwb24hKnlUPlMkISM3JCIxOGByK3J5TU0hIzgkIjEuJj16J3pWdU0hIzgkIjJrQ3gxUiUqM18kISM5JCIxWzoqKm9qOnVOISM4JCIxLTknRyFSQU1PISM4JCIxOG9HIypwNCxQISM4JCIyL3luc2p2WngkISM5JCIyQ0kveSQpZl8mUSEjOSQiMSJRJypRZlxEJVIhIzgkIjE5U2EwXGtPUyEjOCQiMTBzdXNkYVBUISM4JCIxW2ZdJj5fX0MlISM4JCIydkM/UTxrKGZWISM5NzskIjIwRyhIKj1YZidRISM5JCIyV0oiPkxAaURRISM5JCIyVyE0a0tZNSN6JCEjOSQiMVtnayhvI1JsUCEjOCQiMVtuPylIJ1tYUCEjOCQiMS9JS2thUUtQISM4JCIyViJbKmY9ITRFUCEjOSQiMSI9QUtZK21zJCEjOCQiMS5eKydIO1J0JCEjOCQiMSJlVlZvUCFbUCEjOCQiMTp3QkdZJypvUCEjOCQiMkw/KG9GcnAneiQhIzkkIjJ2TSNwIz1ONyRRISM5JCIxXElEJHp5RChRISM4JCIxMCRwJGZ6cz9SISM4JCIxOjYvIm8jb3ZSISM4JCIyOVtvI2VIV1BTISM5JCIxMDkwInozZzUlISM4JCIxIykpKlF6LFEiPSUhIzgkIjE7UkdCcmJqVSEjOCQiMk1dTEZpUkROJSEjOSQiMGxReG5GJFtXISM3JCIyJlwkKkgpR0A0YiUhIzkkIjJNZzpXWD8ubSUhIzkkIjE6dTN3Xl93WiEjODc7JCIyazBncUlQMEYlISM5JCIyPGc2VipwIT5CJSEjOSQiMjhxPXJCIzMrVSEjOSQiMmpOIltOSTF2VCEjOSQiMnVjKlIqUVxvOiUhIzkkIjFNTCgpKUhUYTklISM4JCIxZEUhUndRMzklISM4JCIyYGAoWyV5VEk5JSEjOSQiMW96aWcuMF9UISM4JCIydSZSSyNca3k7JSEjOSQiMSxiZHpUWyE+JSEjOCQiMS1FUUElNCo+VSEjOCQiMWRfdT8tOWNVISM4JCIyJW9NbXVsPCpIJSEjOSQiMU5zOCVbPSFcViEjOCQiMWVsO1xmbTBXISM4JCIxTzl2cCo9InBXISM4JCIyKHA9KmVheCRSWCEjOSQiMiVleWV4O1c7WSEjOSQiMSwlUltPNi5xJSEjOCQiMjtdWXdnJyk0eiUhIzkkIjFkIjRnU24lKSlbISM4JCIyJXB0IypmUHYjKlwhIzkkIjFONlNwYyVRNSYhIzgkIjFlL1ZNSnVAXyEjODc7JCIxMSgpKjNRJ2UuWiEjOCQiMWp4XTYpW21tJSEjOCQiMk1Qc3d6O2xqJSEjOSQiMVREUlIuPjhZISM4JCIxaiNvbVZwbWYlISM4JCIwYSpcKjNhcGUlISM3JCIyWVAnKXlIV1NlJSEjOSQiMWooRz0xU3plJSEjOCQiMTNuSyJRVCcpZiUhIzgkIjEyLVFjI1toaCUhIzgkIjJWRSopcG9nL2slISM5JCIyTyhROnQneTpuJSEjOSQiMjEvdVtALSY0WiEjOSQiMkR3XEBKSlV2JSEjOSQiMVQ1KVwnZncwWyEjOCQiMlAoeU90aDVrWyEjOSQiMk5FNXQkPkRIXCEjOSQiMTQjM29ELjcrJiEjOCQiMiUzPCc9OGcqel0hIzkkIjFsMlppRF9sXiEjOCQiMXRgalswKnlEJiEjOCQiMVRiTiE0a3FOJiEjOCQiMWo3aig9VklZJiEjOCQiMVREWVN5I2VkJiEjOCQiMXYkXClbIT1hcCYhIzg3OyQiMEI4M1QjNGxeISM3JCIxKHp6WmVaKUheISM4JCIyJj4+STkkMzk1JiEjOSQiMnVmeiQqZnUoel0hIzkkIjFLRyxTayVcMSYhIzgkIjIkPjs/T1EjcDAmISM5JCIxbWYlenkxZDAmISM4JCIxbGVDJkgmSGhdISM4JCIyLkssImUkKm90XSEjOSQiMEw3bCgqKSlHNCYhIzckIjEoKil5LzolKik9XiEjOCQiMUA1KyEpW3FeXiEjOCQiMShweV07QDg+JiEjOCQiMUs+cjBJdVBfISM4JCIxQjIhPlNxNEgmISM4JCIxbV1rYEwrXmAhIzgkIjFtXCU0Jz0leVQmISM4JCIxQC8hUSNmWyJcJiEjOCQiMU05QFViJD5kJiEjOCQiMSkqejw7Mj5mYyEjOCQiMUAscVg5RGBkISM4JCIxKnp4MnQ8VCZlISM4JCIxSzVUciYqeWhmISM4JCIxQSkqZm5wRXdnISM4JCIxbFRNPipcdj4nISM4NzskIjFGTyFvUmJdbCYhIzgkIjExeDc5TF1AYyEjOCQiMVR0Kyh5Y1pmJiEjOCQiMUhEVzplIltkJiEjOCQiMXNLVipSIW9oYiEjOCQiMXYmeipRME5iYiEjOCQiMFYiM01pI2ViJiEjNyQiMVQpUVpbMkpjJiEjOCQiMTM9JjRIJT54YiEjOCQiMEw/RmwnMylmJiEjNyQiMTNXL3FYeURjISM4JCIxVFMjSC8pR2djISM4JCIxSCNmODIoZixkISM4JCIxdipcYGw2KFxkISM4JCIxdmkqW3pKWSFlISM4JCIxSCIpKioqW2RqJ2UhIzgkIjFUYmxTKCkpWyRmISM4JCIxMyZvb2FELSwnISM4JCIwLlAnM3pPIzQnISM3JCIxMjYnZiNlSiI9JyEjOCQiMHVTKSlIcHFGJyEjNyQiMCRmRkYkRyd6aiEjNyQiMXVtRTZIKiopWychIzgkIjF3SCIzMGpeZychIzgkIjFKWyJmdVEiR24hIzgtJStBWEVTTEFCRUxTRzYlLUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NVEibTYiLyUnZmFtaWx5R1EhNiIvJSVzaXplR1EjMTA2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EldHJ1ZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ2l0YWxpYzYiLUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NVEiYjYiLyUnZmFtaWx5R1EhNiIvJSVzaXplR1EjMTA2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EldHJ1ZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ2l0YWxpYzYiUSE2Ii0lK1BST0pFQ1RJT05HNiwkIjF3big+cSoqKioqXCEjOyQhMS8kXDNwbjUyKCEjOyQiMXduKD5xKioqKipcISM7JCIxd24oPnEqKioqKlwhIzskIjEvJFwzcG41MighIzskIjF3big+cSoqKioqXCEjOyQhMS8kXDNwbjUyKCEjOyQiIiEhIiIkIjEvJFwzcG41MighIzskIiM1ISIiLSUmU1RZTEVHNiMlL1NVUkZBQ0VDT05UT1VSRy0lKkFYRVNTVFlMRUc2IyUkQk9YRy0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJCsiISIiLSUpQk9VTkRTX1lHNiMkIiQrIiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJStSISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJStSISIiLSUpQ0hJTERSRU5HNiI=

sol:=solve({eq1,eq2},{m,b});

PCQvSSJiRzYiJCIrWyhwKVtPISM1L0kibUdGJSQhK0dWKUh3IiEiKg==

sol[1];

L0kiYkc2IiQiK1socClbTyEjNQ==

foundedLine:=x->subs(sol, m*x+b);

Zio2I0kieEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUklc3Vic0clKnByb3RlY3RlZEc2JEkkc29sR0YlLCYqJkkibUdGJSIiIjkkRjFGMUkiYkdGJUYxRiVGJUYl

foundedLine(x);

LCZJInhHNiIkIStHVilIdyIhIiokIitbKHApW08hIzUiIiI=

Exercise: For the problems above, define a distance function by taking the absolute value of the vertical distances. Then try to minimize that distance. What happens?JSFHJSFH absValueDist:=(p,m,b)->abs(m*p[1]+b-p[2] );

Zio2JUkicEc2IkkibUdGJUkiYkdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUkkYWJzRyUqcHJvdGVjdGVkRzYjLCgqJjklIiIiJjkkNiNGMkYyRjI5JkYyJkY0NiMiIiMhIiJGJUYlRiU=

E:=(data,m,b)->sum(absValueDist(data[i],m,b) , i=1..nops(data) );

Zio2JUklZGF0YUc2IkkibUdGJUkiYkdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUkkc3VtRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliR0YlNiQtSS1hYnNWYWx1ZURpc3RHRiU2JSY5JDYjSSJpR0YlOSU5Ji9GNzsiIiItSSVub3BzR0YuNiNGNUYlRiVGJQ==

E(pts,m,b);

LFQtSSRhYnNHJSpwcm90ZWN0ZWRHNiMsKEkibUc2IiQhK0VjJG8tJSEiKkkiYkdGKSIiIiQhKyJ5dyFIX0YsRi5GLi1GJDYjLChGKCQiK0tdaWglKkYsRi1GLiQiKzcvKys6ISIpRi5GLi1GJDYjLChGKCQhKz5RKXpUKUYsRi1GLiQhKzI2TEc8RjhGLkYuLUYkNiMsKEYoJCIrXHZgTnRGLEYtRi4kIitdQHF2NUY4Ri5GLi1GJDYjLChGKCQiKjI/aW4nRixGLUYuJCErJ3BVQU4iRixGLkYuLUYkNiMsKEYoJCErbkBtVU9GLEYtRi4kIStcKWUsZClGLEYuRi4tRiQ2IywoRigkIitwOyUqM1RGLEYtRi4kIishZXQqXG1GLEYuRi4tRiQ2IywoRigkIStza3YjZiRGLEYtRi4kISs7YF9od0YsRi5GLi1GJDYjLChGKCQiKiVRdDw3ISM1Ri1GLiQiK2g4c1tmRmJvRi5GLi1GJDYjLChGKCQiKzNlMSllJEYsRi1GLiQiK0FxOSdRKEYsRi5GLi1GJDYjLChGKCQiK3kmPjVQKkYsRi1GLiQiK3Y3WF08RjhGLkYuLUYkNiMsKEYoJCIra1MhUUknRixGLUYuJCIrJzRoazYiRjhGLkYuLUYkNiMsKEYoJCIrdGVMJ1IqRixGLUYuJCIrVVRXaTpGOEYuRi4tRiQ2IywoRigkISorLXoyKUYsRi1GLiQhK1RLXWJERixGLkYuLUYkNiMsKEYoJCIqLXA8I3lGLEYtRi4kIispW1U9XSVGLEYuRi4tRiQ2IywoRigkIisqKnpzU09GLEYtRi4kIitIQlNca0YsRi5GLi1GJDYjLChGKCQiK09zKilvO0YsRi1GLiQiK1RQeDlJRixGLkYuLUYkNiMsKEYoJCIrKltIU3EkRixGLUYuJCIra1RxQEtGLEYuRi4tRiQ2IywoRigkIStOdCQqSDZGOEYtRi4kISs/PCFSdiJGOEYuRi4tRiQ2IywoRigkIistMSVmWilGLEYtRi4kIis2NGsieiJGOEYuRi4tRiQ2IywoRigkIStxXClvbSRGLEYtRi4kISthJEcvVylGLEYuRi4tRiQ2IywoRigkIStkJ290aCZGLEYtRi4kIStOYlArN0Y4Ri5GLi1GJDYjLChGKCQhK19XJG9pJkYsRi1GLiQhKysoM3pMKEYsRi5GLi1GJDYjLChGKCQhK0M5dlFARixGLUYuJCErb3UhZWkmRixGLkYuLUYkNiMsKEYoJCErdmFlb0FGLEYtRi4kISspNC5QaCVGLEYuRi4=

plot3d(E(pts,m,b),m=-3..2,b=0..2,style=patchcontour);

LSUnUExPVDNERzYmLSUlR1JJREc2JTskISNJISIiJCIjPyEiIjskIiIhISIiJCIjPyEiIjc7NzskIjItKz5KJTNtI1wiISM5JCIybm0meTR2Iz1cIiEjOSQiMmxtbHAmW28iXCIhIzkkIjItKypIIT49RFwiISM5JCIxTEJqQjpOJFwiISM4JCIybG1scCZbPSVcIiEjOSQiLipIIT49XVwiISM1JCIxTEJqQjomZVwiISM4JCIybG1scCZbbydcIiEjOSQiLipIIT49dlwiISM1JCIxTEJqQjpOKVwiISM4JCIybG1scCZbPSpcIiEjOSQiLipIIT49K10iISM1JCIxTEJqQjomM10iISM4JCIybG1scCZbbyw6ISM5JCIuKkghPj1EXSIhIzUkIjFMQmpCOk4uOiEjOCQiMmxtbHAmWz0vOiEjOSQiLipIIT49XV0iISM1JCIyTExLT19eZV0iISM5JCIyb21scCZbbzE6ISM5JCIuKkghPj12XSIhIzUkIjJMTEtPX14kMzohIzkkIjJvbWxwJls9NDohIzkkIjIpKioqZTQkKXosXiIhIzk3OyQiMmxtIXpHQzddNyEjOSQiMktMZGE0KkdcNyEjOSQiMioqKlI3aWRYWzchIzkkIjJsbSF6R0NpWjchIzkkIjIpKipSIW9KWnVDIiEjOSQiMk1MUCxsIUdbNyEjOSQiMmxtcU0pUjZcNyEjOSQiMikqKlIhb0paKlw3ISM5JCIyTUxQLGwheV03ISM5JCIybG1xTSlSaF43ISM5JCIyKSoqUiFvSlpDRCIhIzkkIjJNTFAsbCFHYDchIzkkIjJsbXFNKVI2YTchIzkkIjIpKipSIW9KWlxEIiEjOSQiMk1MUCxsIXliNyEjOSQiMmltcU0pUmhjNyEjOSQiMikqKlIhb0padUQiISM5JCIyTUxQLGwhR2U3ISM5JCIybG1xTSlSNmY3ISM5JCIyKCoqUiFvSloqZjchIzkkIi45RU9IN0UiISM1JCIyKCoqUmhpJEhQRSIhIzkkIjIoKipSaGkkSGlFIiEjOSQiLjlFT0gobzchIzUkIjIoKipSaGkkSDdGIiEjOTc7JCIybW06cCspZjQ1ISM5JCIybW06cCspNDI1ISM5JCIyaW0melp0ImYrIiEjOSQiMktMaVcsJTMwNSEjOSQiMikqKipHNm9dVSsiISM5JCIyaW0melp0VC41ISM5JCIyTExVbXg0SysiISM5JCIya212KjRKLy81ISM5JCIyJioqKjNMV3dbKyIhIzkkIjJMTFVteDRkKyIhIzkkIjJrbXYqNEphMTUhIzkkIjImKioqM0xXd3QrIiEjOSQiMkxMVW14NCMzNSEjOSQiMmttdio0Si80NSEjOSQiMiYqKiozTFd3KTQ1ISM5JCIyTExVbXg0MiwiISM5JCIyKSoqKm9VKil5QSwiISM5JCIyJioqKm9VKil5WiwiISM5JCIyJioqKm9VKil5cywiISM5JCIyKSoqKm9VKil5KD41ISM5JCIyJioqKm9VKil5QS0iISM5JCIyKSoqKm9VKil5Wi0iISM5JCIyJioqKm9VKil5cy0iISM5JCIxTEInendeLy4iISM4JCIyKSoqKkdZVj1ZLiIhIzk3OyQiMSkqKlImbzByZnghIzkkIjEoKipSJm8wck14ISM5JCIxIyoqUiZvMHI0eCEjOSQiMSkqKlImbzByJW8oISM5JCIxKCoqUiZvMHJmdyEjOSQiMSMqKlImbzByTXchIzkkIjEmKipSJm8wcjR3ISM5JCIxaW0rb0U3I2YoISM5JCIwbTFbT0EoKmUoISM4JCIxJyoqUiIpcGIhKWYoISM5JCIxSExaSiEqUTF3ISM5JCIwbTFbT0FaaCghIzgkIjEoKipSI2ZVR0x3ISM5JCIxJyoqUiNmVUdldyEjOSQiMUZMPFFkXipvKCEjOSQiMScqKlJbUyM9SnghIzkkIjFobV1yIVxHeCghIzkkIjFGTDxRZF45eSEjOSQiMScqKlJbUyM9Y3khIzkkIjFqbV1yIVx5KnkhIzkkIjFHTDxRZF5SeiEjOSQiMSkqKlJbUyM9Iil6ISM5JCIxaG1dciFcRy0pISM5JCIxR0w8UWReayEpISM5JCIxKSoqUltTIz0xIikhIzk3OyQiMSYqKipbIzMpM2VjISM5JCIyOW1jIlxadipmJiEjOiQiMmtLQmVUQDlhJiEjOiQiMSYqKipbIzMpMyRbJiEjOSQiMjltYyJcWnZDYSEjOiQiMjFtY25qYkBQJiEjOiQiMSYqKiozOTV3VGAhIzkkIjFqbXYhb0ZNTCYhIzkkIjJjKioqb0VPPlFgISM6JCIxaG1ESEApNE8mISM5JCIxSExzbjBmIlImISM5JCIyYioqKlFNc0RMYSEjOiQiMlkqKio9InBhbFomISM6JCIxR0xfQyEpKVtgJiEjOSQiMWltJnlOQEtmJiEjOSQiMlkqKio9InBhOmwmISM6JCIxR0xfQyEpKSk0ZCEjOSQiMWltJnlOQCNvZCEjOSQiMUdMVWRqeVFlISM5JCIyJ0dMVWRqeThmISM6JCIxR0xVZGp5KSlmISM5JCIxR0xVZGp5amchIzkkIjFHTFVkanlRaCEjOSQiMmxLQnVOJ3k4aSEjOiQiMUdMVWRqeSlHJyEjOTc7JCIxS0xfej9UIlslISM5JCIxbW0mR1RYKCpRJSEjOSQiMicpKioqPVkoeSEpSCUhIzokIjInKSpcJzR4K0FBJSEjOiQiMUwkKVJvMjxqVCEjOSQiMicpKlwxTnUkWzUlISM6JCIxbDt0LFRdWVMhIzkkIjFMJClSbzI8KSlSISM5JCIxazt0YydHIVFSISM5JCIxbW1EYzFINVIhIzkkIjEqKlxPcihRQipRISM5JCIxbDspKj5gI3oqUSEjOSQiMSoqXDYrRGE2UiEjOSQiMWxtNWFUZl9SISM5JCIxSkx4PzNFJSpSISM5JCIxSkxQYShmKVxTISM5JCIwTEJvcHFINyUhIzgkIjFrOzhmSyc+PyUhIzkkIjFqbT9fTyMpPVYhIzkkIjFrbT9fTyNRVyUhIzkkIjFrbT9fTyMpb1ghIzkkIjFqbT9fTyNRcCUhIzkkIjFqbT9fTyMpPVshIzkkIjJDbTFBbEJRJVwhIzokIjFqbT9fTyMpb10hIzk3OyQiMicqKipSK3cwTjElISM6JCIyJyoqKlIrdzAmUVMhIzokIi0vZ2RdOFMhIzUkIjInKioqUit3MCYpKVIhIzokIjInKioqUit3ME4nUiEjOiQiLS9nZF1RUiEjNSQiMicqKipSK3cwTiJSISM6JCIya21jJmV2LigqUSEjOiQiMW5tP2I9SSQqUSEjOSQiMjErUyYpPU47IVIhIzokIjIxKyEqcCJmIUciUiEjOiQiLSpwImYheSRSISM1JCIyMSshKnAiZiFHJ1IhIzokIjIxKyEqcCJmIXkpUiEjOiQiMUxMN3NedzxTISM5JCIxLCt6UT1WZlMhIzkkIjJ2bWNhXSk0LFQhIzokIjJtbWN1P2U3OiUhIzokIi16UzpmNFUhIzUkIjEqKioqZS5PLXRVISM5JCIxKioqKmUuTy1bViEjOSQiLWYuTy1CVyEjNSQiMSoqKiplLk8tKVwlISM5JCIyQ0xCPFdlKHpYISM6JCIyJioqKipRM15VclkhIzo3OyQiMXA7WzozNU5ZISM5JCIyRCs6KVtUVlZZISM6JCIyOytsXSMzRmhZISM6JCIxLF0xRDNGJ28lISM5JCIyOytsXSMzRjZaISM6JCIyOytsXSMzRk9aISM6JCIxLF0xRDNGaFohIzkkIjFwbVhaNWIjeiUhIzkkIjFOTDc5eEBNWyEjOSQiMjkrIXohUSUpZShbISM6JCIxcG1YWjViPFwhIzkkIjJrTEJUcjwjZlwhIzokIjI5KyF6IVElKTMrJiEjOiQiMXBtWFo1YlVdISM5JCIxLF1KXDRIJTQmISM5JCIxTiRbRUdDRTomISM5JCIxTyRbJXpBZjpfISM5JCIxTyRbJXpBZiFIJiEjOSQiMU0kWyV6QWZsYCEjOSQiMU8kWyV6QWZTYSEjOSQiMlhMWyV6QWY6YiEjOiQiMU0kWyV6QWYhZiYhIzkkIjJYTFslekFmbGMhIzokIjJYTFslekFmU2QhIzokIjFNJFslekFmOmUhIzk3OyQiMU9MI1I4bUNWJyEjOSQiMU9MI1I4bXVYJyEjOSQiMU9MI1I4bUNbJyEjOSQiMU9MI1I4bXVdJyEjOSQiMU9MI1I4bUNgJyEjOSQiMU9MI1I4bXViJyEjOSQiMU9MI1I4bUNlJyEjOSQiMU9MI1I4bXVnJyEjOSQiMU9MI1I4bUNqJyEjOSQiMU9MI1I4bXVsJyEjOSQiMVBMI1I4bUNvJyEjOSQiMU9MI1I4bXVxJyEjOSQiMU9MI1I4bUN0JyEjOSQiMU9MI1I4bXV2JyEjOSQiMU1MI1I4bUN5JyEjOSQiMU9MI1I4bXUhbyEjOSQiMU1MI1I4bUMkbyEjOSQiMU1MI1I4bXUmbyEjOSQiMU9MczBpeCcpbyEjOSQiLVJzR1dHcCEjNSQiMXBtWHokZWMpcCEjOSQiMU9MX1s4Zl1xISM5JCIxTUxfWzhmRHIhIzkkIjFRTF9bOGYrcyEjOSQiMVFMX1s4ZnZzISM5NzskIi1IRCNSOm8pISM1JCIwbmM+KmU/dCcpISM4JCIxTkxpZUQoW20pISM5JCItSEQjUmxsKSEjNSQiMG5jPiplP1snKSEjOCQiMU5MaWVEKClSJykhIzkkIi1IRCNSOmopISM1JCIwbmM+KmU/QicpISM4JCIxTkxpZUQoW2gpISM5JCIxLCslXCk9dzsnKSEjOSQiMU9MRj1fNEQnKSEjOSQiMG4xO2JHTWopISM4JCIxLV1yIzMncF8nKSEjOSQiMW87M2goPUdwKSEjOSQiMU4kW3hVJltNKCkhIzkkIjEtXVQlNF9oeCkhIzkkIjFvOzNoKD15IikpISM5JCIxTiRbeFUmW2YpKSEjOSQiMS1dVCU0XzYhKikhIzkkIjFvOzNoKD1HJSopISM5JCIxbjspeSE0WikpKikhIzkkIjEuXUBUVSFvLyohIzkkIjFNJFtYZFBeNSohIzkkIjFuOyl5ITRaaiIqISM5JCIxLl1AVFUhPUEqISM5NzskIjFuYzdDRjoqNCIhIzgkIjJOTCN6IVI+JCk0IiEjOSQiMjArZnUwJ1soNCIhIzkkIjJPTDIsZ3oiKTQiISM5JCIybm1TTSRILCo0IiEjOSQiMjErdW5FWSkqNCIhIzkkIjJPTDIsZ3oxNSIhIzkkIjJubVNNJEheLDYhIzkkIjIxK3VuRVlCNSIhIzkkIjJPTDIsZ3pKNSIhIzkkIjJubVNNJEgsLzYhIzkkIjIxK3VuRVlbNSIhIzkkIjJPTDIsZ3pjNSIhIzkkIjFuMVdMSF4xNiEjOCQiMjErdW5FWXQ1IiEjOSQiMk9MMixneiIzNiEjOSQiMnBtU00kSCw0NiEjOSQiMjErdW5FWSk0NiEjOSQiMk9MMixnejE2IiEjOSQiMk1MKEglUUJDNiIhIzkkIjJUTChIJVFCXDYiISM5JCIyTkwoSCVRQnU2IiEjOSQiLldKWipHPzYhIzUkIjJ0bTUpUmhYQzYhIzkkIjJPTGlkQFEkSDYhIzk3OyQiMW4iM3omPTZVOCEjOCQiMi9dVDc+WEhNIiEjOSQiMk8kW2RDJnlQTSIhIzkkIjFuIjN6Jj1oVzghIzgkIjIvXVQ3PlhhTSIhIzkkIjJPJFtkQyZ5aU0iISM5JCIxbiIzeiY9Nlo4ISM4JCIyL11UNz5Yek0iISM5JCIyTyRbZEMmeShbOCEjOSQiMW4iM3omPWhcOCEjOCQiMi9dVDc+WC9OIiEjOSQiMk8kW2RDJnk3TiIhIzkkIjFuIjN6Jj02XzghIzgkIjIvXVQ3PlhITiIhIzkkIjJSJFtkQyZ5UE4iISM5JCIxbiIzeiY9aGE4ISM4JCIyL11UNz5YYU4iISM5JCIyUiRbZEMmeWlOIiEjOSQiMW4iM3omPTZkOCEjOCQiMi9dVDc+WHpOIiEjOSQiMlIkW2RDJnkoZTghIzkkIjFuIjN6Jj1oZjghIzgkIjJvbWJiLTQxTyIhIzkkIjFuY2JEITRKTyIhIzgkIjFuY0JiRSVmTyIhIzg3OyQiMjErNGQ2VyYpZSIhIzkkIjJQTFUhXHVQKmUiISM5JCIxbmNQI3k1LWYiISM4JCIyMSs0ZDZXNWYiISM5JCIyUExVIVx1KD1mIiEjOSQiMW5jUCN5NUZmIiEjOCQiMjErNGQ2V05mIiEjOSQiMlBMVSFcdVAlZiIhIzkkIjFuY1AjeTVfZiIhIzgkIjIxKzRkNldnZiIhIzkkIjJQTFUhXHUob2YiISM5JCIxbmNQI3k1eGYiISM4JCIyMSs0ZDZXJilmIiEjOSQiMlBMVSFcdVAqZiIhIzkkIjFuY1AjeTUtZyIhIzgkIjIxKzRkNlc1ZyIhIzkkIjJQTFUhXHUoPWciISM5JCIxbmNQI3k1RmciISM4JCIyLis0ZDZXTmciISM5JCIyUExVIVx1UC87ISM5JCIxbmNQI3k1X2ciISM4JCIyLys+PC9CdGciISM5JCIyLys+PC9CKTQ7ISM5JCIyLys+PC9CQmgiISM5JCIyLys+PC9CW2giISM5NzskIjJOJCk0TlB3XCQ9ISM5JCIxbkolb3E0ZSQ9ISM4JCIyMl13LC9WbSQ9ISM5JCIyTiQpNE5Qd3UkPSEjOSQiMW5KJW9xNCRRPSEjOCQiMjJddywvViJSPSEjOSQiMk4kKTROUHcqUj0hIzkkIjFuSiVvcTQzJT0hIzgkIjIyXXcsL1Y7JT0hIzkkIjJOJCk0TlB3QyU9ISM5JCIxbkolb3E0TCU9ISM4JCIyMl13LC9WVCU9ISM5JCIyTiQpNE5Qd1wlPSEjOSQiMW5KJW9xNGUlPSEjOCQiMjJddywvVm0lPSEjOSQiMk4kKTROUHd1JT0hIzkkIjFuSiVvcTQkWz0hIzgkIjIyXXcsL1YiXD0hIzkkIjJOJCk0TlB3Klw9ISM5JCIyUEwjKXkwUDomPSEjOSQiMlBMIyl5MFBTJj0hIzkkIjJNTCMpeTBQbCY9ISM5JCIxTUIpeTBQIWY9ISM4JCIyVUwjKXkwUDonPSEjOSQiMk1MIyl5MFBTJz0hIzk3OyQiMW4xSkonMzkzIyEjOCQiMi0rV1knPkMjMyMhIzkkIjJVTHh6SHZJMyMhIzkkIjFuMUpKJzNSMyMhIzgkIjItK1dZJz51JTMjISM5JCIyVUx4ekh2YjMjISM5JCIxbjFKSiczazMjISM4JCIyLStXWSc+QygzIyEjOSQiMlVMeHpIdiEpMyMhIzkkIjFuMUpKJzMqKTMjISM4JCIyLStXWSc+dSozIyEjOSQiMlVMeHpIdjA0IyEjOSQiMW4xSkonMzk0IyEjOCQiMi0rV1knPkMjNCMhIzkkIjFNdCh6SHZJNCMhIzgkIjFuMUpKJzNSNCMhIzgkIjItK1dZJz51JTQjISM5JCIydm1YUzJeZDQjISM5JCIxbmMvdTVEKTQjISM4JCIxbmMvdTV2K0AhIzgkIjJ2bVhTMl5LNSMhIzkkIjFuYy91NXYwQCEjOCQiMW5jL3U1RDNAISM4JCIydm1YUzJeMjYjISM5JCIxbmMvdTVEOEAhIzg3OyQiMi9dNiIqKTMleUsjISM5JCIyTiRbV0FVbkdCISM5JCIydjt5ZGIyJkhCISM5JCIyL102IiopM01JQiEjOSQiMk4kW1dBVTxKQiEjOSQiMnY7eWRiMj9MIyEjOSQiMi9dNiIqKTMlR0wjISM5JCIyTiRbV0FVbkxCISM5JCIydjt5ZGIyWEwjISM5JCIyL102IiopM01OQiEjOSQiMk4kW1dBVTxPQiEjOSQiMnY7eWRiMnFMIyEjOSQiMi9dNiIqKTMleUwjISM5JCIyTiRbV0FVblFCISM5JCIydjt5ZGIyJlJCISM5JCIyL102IiopM01TQiEjOSQiMi0rNC00bENNIyEjOSQiMjMrNC00bFxNIyEjOSQiMi0rNC00bHVNIyEjOSQiMi0rNC00bCpcQiEjOSQiMjMrNC00bENOIyEjOSQiMi0rNC00bFxOIyEjOSQiLjQtNGx1TiMhIzUkIjIzKzQtNGwqZkIhIzkkIjItKzQtNGxDTyMhIzk3OyQiMkxMN3A5dFVkIyEjOSQiMndtWC1bMV5kIyEjOSQiMSwhek4iKVJmZCMhIzgkIjJMTDdwOXRuZCMhIzkkIjJ3bVgtWzF3ZCMhIzkkIjEsIXpOIilSJXlEISM4JCIyTEw3cDl0I3pEISM5JCIyd21YLVsxLGUjISM5JCIxLCF6TiIpUjRlIyEjOCQiMkxMN3A5dDxlIyEjOSQiMW5jQyFbMUVlIyEjOCQiMSwhek4iKVJNZSMhIzgkIjJMTDdwOXRVZSMhIzkkIjJ3bVgtWzFeZSMhIzkkIjFNQlAxInptZSMhIzgkIjFNQlAxInoiKmUjISM4JCIyVkxzajV6O2YjISM5JCIxTUJQMSJ6VGYjISM4JCIxTUJQMSJ6bWYjISM4JCIxTUJQMSJ6IipmIyEjOCQiMU1CUDEiejtnIyEjOCQiMU1CUDEielRnIyEjOCQiMU1CUDEiem1nIyEjOCQiMU1CUDEieiI0RSEjOCQiMU1CUDEiejtoIyEjODc7JCIxbkpyL2FxP0chIzgkIi9sL1EoUTojRyEjNiQiMlYkKXo4MnNCI0chIzkkIjFuSnIvYT9CRyEjOCQiL2wvUShRUyNHISM2JCIyViQpejgyc1sjRyEjOSQiMW5Kci9hcURHISM4JCIvbC9RKFFsI0chIzYkIjJWJCl6ODJzdCNHISM5JCIxbkpyL2E/R0chIzgkIi9sL1EoUSFIRyEjNiQiMlYkKXo4MnMpSEchIzkkIjJ0bU5ENyQqMyRHISM5JCIxbmNgQUpSTEchIzgkIjJ0bU5ENyQqZSRHISM5JCIydG1ORDckUlFHISM5JCIydG1ORDckKjMlRyEjOSQiMnRtTkQ3JFJWRyEjOSQiMnRtTkQ3JCplJUchIzkkIjJ0bU5ENyRSW0chIzkkIjJ0bU5ENyQqMyZHISM5JCIydG1ORDckUmBHISM5JCIydG1ORDckKmUmRyEjOSQiMnRtTkQ3JFJlRyEjOSQiMW5jYEFKKjMnRyEjODc7JCIyMCs5RG1QcjEkISM5JCIyT0xaZSo0KHoxJCEjOSQiMnRtIT1IViEpb0khIzkkIjIwKzlEbVAncEkhIzkkIjJPTFplKjRacUkhIzkkIjJ0bSE9SFZJckkhIzkkIjIwKzlEbVBAMiQhIzkkIjJPTFplKjQoSDIkISM5JCIydG0hPUhWIVEyJCEjOSQiMjArOURtUFkyJCEjOSQiMk9MWmUqNFp2SSEjOSQiMSwhKnBRcmd4SSEjOCQiMSwhKnBRcjUhMyQhIzgkIjIuKypwUXJnIzMkISM5JCIxLCEqcFFyNSYzJCEjOCQiMi4rKnBRcmcoMyQhIzkkIjIuKypwUXI1ITQkISM5JCIxLCEqcFFyZyM0JCEjOCQiMi4rKnBRcjUmNCQhIzkkIjIuKypwUXJnKDQkISM5JCIxLCEqcFFyNStKISM4JCIyLisqcFFyZy1KISM5JCIxLCEqcFFyNTBKISM4JCIxLCEqcFFyZzJKISM4JCIyLisqcFFyNTVKISM5NzskIjFNW0o/KnBOSiQhIzgkIjIkbyJbT0QuV0okISM5JCIyLl0iKXBlT19KJCEjOSQiMU1bSj8qcGdKJCEjOCQiMiRvIltPRC5wSiQhIzkkIjIuXSIpcGVPeEokISM5JCIxTVtKPypwJj1MISM4JCIyJG8iW09ELiU+TCEjOSQiMi5dIilwZU8tSyQhIzkkIjJXTGlbOkA9SyQhIzkkIjFOQidbOkBWSyQhIzgkIjJMTGlbOkBvSyQhIzkkIjJXTGlbOkAkSEwhIzkkIjFOQidbOkA9TCQhIzgkIjJMTGlbOkBWTCQhIzkkIjFNQidbOkBvTCQhIzgkIjFOQidbOkAkUkwhIzgkIjJMTGlbOkA9TSQhIzkkIjFNQidbOkBWTSQhIzgkIjFOQidbOkBvTSQhIzgkIjJMTGlbOkAkXEwhIzkkIjFNQidbOkA9TiQhIzgkIjFOQidbOkBWTiQhIzgkIjJMTGlbOkBvTiQhIzkkIjFNQidbOkAkZkwhIzg3OyQiMW5jNnlAK2dOISM4JCIyMitcOV5OM2MkISM5JCIxTUJ5VylvO2MkISM4JCIxbmM2eUBdaU4hIzgkIjIyK1w5Xk5MYyQhIzkkIjFNQnlXKW9UYyQhIzgkIjFuYzZ5QCtsTiEjOCQiMnVtRDU8TmdjJCEjOSQiMnVtRDU8TiZvTiEjOSQiMW9jLXJeLnJOISM4JCIydW1ENTxOTmQkISM5JCIydW1ENTxOZ2QkISM5JCIxb2Mtcl5geU4hIzgkIjJ1bUQ1PE41ZSQhIzkkIjJ1bUQ1PE5OZSQhIzkkIjFvYy1yXi4nZSQhIzgkIjJ1bUQ1PE4mKWUkISM5JCIydW1ENTxONWYkISM5JCIydW1ENTxOTmYkISM5JCIydW1ENTxOZ2YkISM5JCIydW1ENTxOJilmJCEjOSQiMnVtRDU8TjVnJCEjOSQiMnVtRDU8Tk5nJCEjOSQiMnVtRDU8TmdnJCEjOSQiMW5jLXJeYDNPISM4NzskIi9sImZWTWshUSEjNiQiMU0pXCNweEUyUSEjOCQiMW1KZS02NTNRISM4JCIvbCJmVk0qM1EhIzYkIjFNKVwjcHh3NFEhIzgkIjFtSmUtNmc1USEjOCQiMi8rKj0oPVxGIlEhIzkkIjIvKyo9KD1cXyJRISM5JCIuKj0oPVx4IlEhIzUkIjIvKyo9KD1cLSNRISM5JCIyLysqPSg9XEYjUSEjOSQiLio9KD1cXyNRISM1JCIyLysqPSg9XHgjUSEjOSQiMi8rKj0oPVwtJFEhIzkkIi4qPSg9XEYkUSEjNSQiMi8rKj0oPVxfJFEhIzkkIjI6Kyo9KD1ceCRRISM5JCIuKj0oPVwtJVEhIzUkIjIvKyo9KD1cRiVRISM5JCIyOisqPSg9XF8lUSEjOSQiLio9KD1ceCVRISM1JCIyLysqPSg9XC0mUSEjOSQiMSwhKj0oPVxGJlEhIzgkIi4qPSg9XF8mUSEjNSQiMi8rKj0oPVx4JlEhIzk3OyQiMk5MPFBwbUcwJSEjOSQiMW8xMEYrcWBTISM4JCIuJVFnTGBhUyEjNSQiMk5MPFBwbWAwJSEjOSQiMU1CTi5LJ3AwJSEjOCQiMk1MX0w/aiVmUyEjOSQiMlhMX0w/aj4xJSEjOSQiMU1CTi5LWWtTISM4JCIyTUxfTD9qcDElISM5JCIyWExfTD9qJXBTISM5JCIxTUJOLksnPjIlISM4JCIyTUxfTD9qVzIlISM5JCIyWExfTD9qcDIlISM5JCIxTUJOLktZelMhIzgkIjFMQk4uSyc+MyUhIzgkIjJYTF9MP2pXMyUhIzkkIjJYTF9MP2pwMyUhIzkkIjFMQk4uS1kqMyUhIzgkIjJNTF9MP2o+NCUhIzkkIjJYTF9MP2pXNCUhIzkkIjFMQk4uSydwNCUhIzgkIjJNTF9MP2olKjQlISM5JCIyWExfTD9qPjUlISM5JCIxTEJOLktZL1QhIzgkIjJNTF9MP2pwNSUhIzk3OyQiMW4iPTomKilIKkglISM4JCIxLDomW0dLLEklISM4JCIxbGNePnM8LFYhIzgkIjJ2bTomPnNuLlYhIzkkIjJ2bTomPnM8MVYhIzkkIjFsY14+c24zViEjOCQiMnZtOiY+czw2ViEjOSQiMnZtOiY+c244ViEjOSQiMWxjXj5zPDtWISM4JCIydm06Jj5zbj1WISM5JCIydm06Jj5zPEBWISM5JCIxbGNePnNuQlYhIzgkIjJ2bTomPnM8RVYhIzkkIjJ2bTomPnNuR1YhIzkkIjFsY14+czxKViEjOCQiMnZtOiY+c25MViEjOSQiMW9jXj5zPE9WISM4JCIxbGNePnNuUVYhIzgkIjJ2bTomPnM8VFYhIzkkIjFvY14+c25WViEjOCQiMWxjXj5zPFlWISM4JCIydm06Jj5zbltWISM5JCIxb2NePnM8XlYhIzgkIjFsY14+c25gViEjOCQiMmttOiY+czxjViEjOTc7JCIuPiQ0N3RYWCEjNSQiMSwhemNCIip5YSUhIzgkIjEsIXpjQiJSXVghIzgkIi56Y0IiKkdiJSEjNSQiMSwhemNCIlJiWCEjOCQiMSwhemNCIip5YiUhIzgkIi56Y0IiUmdYISM1JCIxLCF6Y0IiKkdjJSEjOCQiMSwhemNCIlJsWCEjOCQiLnpjQiIqeWMlISM1JCIxLCF6Y0IiUnFYISM4JCIxLCF6Y0IiKkdkJSEjOCQiLnpjQiJSdlghIzUkIjEsIXpjQiIqeWQlISM4JCIxLCF6Y0IiUiFlJSEjOCQiLnpjQiIqR2UlISM1JCIxLCF6Y0IiUiZlJSEjOCQiMSwhemNCIip5ZSUhIzgkIi56Y0IiUiFmJSEjNSQiMSwhemNCIipHZiUhIzgkIjIwK3pjQiJSJmYlISM5JCIuemNCIip5ZiUhIzUkIjEsIXpjQiJSK1khIzgkIjIwK3pjQiIqR2clISM5JCIuemNCIlIwWSEjNS0lK1BST0pFQ1RJT05HNiwkIjFhYUAhZUQoXEYhIzskITFpI291QD1FZyohIzskIjFOYiVHKj41JXklISM8JCIxWkdhNE83T3ohIzskIjJQIVFGJz4jelpEISM8JCIwYmpLVSxeXyYhIzokITFwTyNlcEV1VSYhIzskITIpPU8+d3hkUjYhIzwkIjEnR29vZUM4SykhIzskIiM1ISIiLSUmU1RZTEVHNiMlL1NVUkZBQ0VDT05UT1VSRy0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJCsiISIiLSUpQk9VTkRTX1lHNiMkIiQrIiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJStSISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJStSISIiLSUpQ0hJTERSRU5HNiI=

First try "blindly" to repeat what we did before. eqm:=diff(E(pts,m,b),m); eqb:=diff(E(pts,m,b),b);

LFQtSSRhYnNHJSpwcm90ZWN0ZWRHNiQiIiIsKEkibUc2IiQiK0tdaWglKiEiKkkiYkdGKkYnJCIrNy8rKzohIilGJ0YrLUYkNiRGJywoRikkIStFYyRvLSVGLUYuRickISsieXchSF9GLUYnRjUtRiQ2JEYnLChGKSQhKz5RKXpUKUYtRi5GJyQhKzI2TEc8RjFGJ0Y8LUYkNiRGJywoRikkIitcdmBOdEYtRi5GJyQiK11AcXY1RjFGJ0ZDLUYkNiRGJywoRikkIioyP2luJ0YtRi5GJyQhKydwVUFOIkYtRidGSi1GJDYkRicsKEYpJCErbkBtVU9GLUYuRickIStcKWUsZClGLUYnRlEtRiQ2JEYnLChGKSQiK3A7JSozVEYtRi5GJyQiKyFldCpcbUYtRidGWC1GJDYkRicsKEYpJCErc2t2I2YkRi1GLkYnJCErO2BfaHdGLUYnRmluLUYkNiRGJywoRikkIiolUXQ8NyEjNUYuRickIitoOHNbZkZib0YnRmBvLUYkNiRGJywoRikkIiszZTEpZSRGLUYuRickIitBcTknUShGLUYnRmhvLUYkNiRGJywoRikkIit5Jj41UCpGLUYuRickIit2N1hdPEYxRidGX3AtRiQ2JEYnLChGKSQiK2tTIVFJJ0YtRi5GJyQiKyc0aGs2IkYxRidGZnAtRiQ2JEYnLChGKSQiK3RlTCdSKkYtRi5GJyQiK1VUV2k6RjFGJ0ZdcS1GJDYkRicsKEYpJCEqKy16MilGLUYuRickIStUS11iREYtRidGZHEtRiQ2JEYnLChGKSQiKi1wPCN5Ri1GLkYnJCIrKVtVPV0lRi1GJ0Zbci1GJDYkRicsKEYpJCIrKip6c1NPRi1GLkYnJCIrSEJTXGtGLUYnRmJyLUYkNiRGJywoRikkIitPcyopbztGLUYuRickIitUUHg5SUYtRidGaXItRiQ2JEYnLChGKSQiKypbSFNxJEYtRi5GJyQiK2tUcUBLRi1GJ0Zgcy1GJDYkRicsKEYpJCErTnQkKkg2RjFGLkYnJCErPzwhUnYiRjFGJ0Zncy1GJDYkRicsKEYpJCIrLTElZlopRi1GLkYnJCIrNjRrInoiRjFGJ0ZedC1GJDYkRicsKEYpJCErcVwpb20kRi1GLkYnJCErYSRHL1cpRi1GJ0ZldC1GJDYkRicsKEYpJCErZCdvdGgmRi1GLkYnJCErTmJQKzdGMUYnRlx1LUYkNiRGJywoRikkIStfVyRvaSZGLUYuRickISsrKDN6TChGLUYnRmN1LUYkNiRGJywoRikkIStDOXZRQEYtRi5GJyQhK291IWVpJkYtRidGanUtRiQ2JEYnLChGKSQhK3ZhZW9BRi1GLkYnJCErKTQuUGglRi1GJ0Zhdg==

LFQtSSRhYnNHJSpwcm90ZWN0ZWRHNiQiIiIsKEkibUc2IiQiK0tdaWglKiEiKkkiYkdGKkYnJCIrNy8rKzohIilGJ0YnLUYkNiRGJywoRikkIStFYyRvLSVGLUYuRickISsieXchSF9GLUYnRictRiQ2JEYnLChGKSQhKz5RKXpUKUYtRi5GJyQhKzI2TEc8RjFGJ0YnLUYkNiRGJywoRikkIitcdmBOdEYtRi5GJyQiK11AcXY1RjFGJ0YnLUYkNiRGJywoRikkIioyP2luJ0YtRi5GJyQhKydwVUFOIkYtRidGJy1GJDYkRicsKEYpJCErbkBtVU9GLUYuRickIStcKWUsZClGLUYnRictRiQ2JEYnLChGKSQiK3A7JSozVEYtRi5GJyQiKyFldCpcbUYtRidGJy1GJDYkRicsKEYpJCErc2t2I2YkRi1GLkYnJCErO2BfaHdGLUYnRictRiQ2JEYnLChGKSQiKiVRdDw3ISM1Ri5GJyQiK2g4c1tmRmJvRidGJy1GJDYkRicsKEYpJCIrM2UxKWUkRi1GLkYnJCIrQXE5J1EoRi1GJ0YnLUYkNiRGJywoRikkIit5Jj41UCpGLUYuRickIit2N1hdPEYxRidGJy1GJDYkRicsKEYpJCIra1MhUUknRi1GLkYnJCIrJzRoazYiRjFGJ0YnLUYkNiRGJywoRikkIit0ZUwnUipGLUYuRickIitVVFdpOkYxRidGJy1GJDYkRicsKEYpJCEqKy16MilGLUYuRickIStUS11iREYtRidGJy1GJDYkRicsKEYpJCIqLXA8I3lGLUYuRickIispW1U9XSVGLUYnRictRiQ2JEYnLChGKSQiKyoqenNTT0YtRi5GJyQiK0hCU1xrRi1GJ0YnLUYkNiRGJywoRikkIitPcyopbztGLUYuRickIitUUHg5SUYtRidGJy1GJDYkRicsKEYpJCIrKltIU3EkRi1GLkYnJCIra1RxQEtGLUYnRictRiQ2JEYnLChGKSQhK050JCpINkYxRi5GJyQhKz88IVJ2IkYxRidGJy1GJDYkRicsKEYpJCIrLTElZlopRi1GLkYnJCIrNjRrInoiRjFGJ0YnLUYkNiRGJywoRikkIStxXClvbSRGLUYuRickISthJEcvVylGLUYnRictRiQ2JEYnLChGKSQhK2Qnb3RoJkYtRi5GJyQhK05iUCs3RjFGJ0YnLUYkNiRGJywoRikkIStfVyRvaSZGLUYuRickISsrKDN6TChGLUYnRictRiQ2JEYnLChGKSQhK0M5dlFARi1GLkYnJCErb3UhZWkmRi1GJ0YnLUYkNiRGJywoRikkISt2YWVvQUYtRi5GJyQhKyk0LlBoJUYtRidGJw==

plot3d(E(pts,m,b),m=-2..-1.6,b=0..2,axes=boxed, style=patchcontour); plot3d(E(pts,m,b),m=-2..-1.6,b=0..2,axes=boxed, style=patchcontour); plot3d(E(pts,m,b),m=-1.95..-1.8,b=0.5..1.5,axes=boxed, style=patchcontour); plot3d(E(pts,m,b),m=-1.95..-1.8,b=0.7..1.2,axes=boxed, style=patchcontour); plot3d(E(pts,m,b),m=-1.905..-1.85,b=0.8..1.05,axes=boxed, style=patchcontour); plot3d(E(pts,m,b),m=-1.905..-1.86,b=0.8..0.95,axes=boxed, style=patchcontour); plot3d(E(pts,m,b),m=-1.904..-1.886,b=0.8..0.9,axes=boxed, style=patchcontour); From the above picture we deduce that the minimum of E is when m is between -1.904 and -1.886 and b is between 0.8 and 0.9. 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 newpts := [seq([i, 1.5*i+1], i = -1 .. 5), [8, 1.5*8+.5], seq([i, 1.5*i+1], i = 9 .. 11)];

Ny03JCEiIiQhIiZGJDckIiIhJCIiIkYoNyRGKiQiI0RGJDckIiIjJCIjU0YkNyQiIiQkIiNiRiQ3JCIiJSQiI3FGJDckIiImJCIjJilGJDckIiIpJCIkRCJGJDckIiIqJCIkWCJGJDckIiM1JCIkZyJGJDckIiM2JCIkdiJGJA==

plot(newpts,style=point,symbol=solidbox,symbolsize=16);

LSUlUExPVEc2Jy0lJ0NVUlZFU0c2JDctNyQkISM1ISIiJCEiJiEiIjckJCIiISEiIiQiIzUhIiI3JCQiIzUhIiIkIiNEISIiNyQkIiM/ISIiJCIjUyEiIjckJCIjSSEiIiQiI2IhIiI3JCQiI1MhIiIkIiNxISIiNyQkIiNdISIiJCIjJikhIiI3JCQiIyEpISIiJCIkRCIhIiI3JCQiIyEqISIiJCIkWCIhIiI3JCQiJCsiISIiJCIkZyIhIiI3JCQiJDUiISIiJCIkdiIhIiItJSZDT0xPUkc2JiUkUkdCRyQiIzUhIiIkIiIhISIiJCIiISEiIi0lK0FYRVNMQUJFTFNHNiRRITYiUSE2Ii0lJ1NZTUJPTEc2JCUqX1NPTElEQk9YRyIjOy0lJlNUWUxFRzYjJSZQT0lOVEctJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiRTIiEiIi0lKUJPVU5EU19ZRzYjJCIkSSIhIiItJS1CT1VORFNfV0lEVEhHNiMkIiVJUCEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiVdTiEiIi0lKUNISUxEUkVORzYi

with(CurveFitting): pol:=PolynomialInterpolation(newpts,x);

LDgqJEkieEc2IiIjNSQiK2FsKUdmJSEjOyokRiQiIiokISslM3EzLSMhIzkqJEYkIiIpJCIrI1woeUJPISM4KiRGJCIiKCQhK3BZMzFNISM3KiRGJCIiJyQiKyMpcCwkeiIhIzYqJEYkIiImJCErIjMmKmU7JkY9KiRGJCIiJSQiK2M8JTQjb0Y9KiRGJCIiJCQiKD92QCYhIzUqJEYkIiIjJCErTURDXScpRj1GJCQiKyNcalhiIiEiKiRGJiEiIiIiIg==

line:=LeastSquares(newpts,x);

LCYkIitDLF8lKSoqISM1IiIiSSJ4RzYiJCIrdT9yIVwiISIq

display(plot({line,pol},x=-1..20,thickness=3,color=[blue,green]),plot(newpts,style=point,symbol=solidbox,symbolsize=16));

LSUlUExPVEc2Ki0lJ0NVUlZFU0c2JTdTNyQkISM1ISIiJCExKioqKipmaCtFI1whIzs3JCQhMSwrK0R3ZkFhISM7JCIyYHhvPlApKTQhPiEjPDckJCEyJSoqKlxQSEMpUjkhIzwkIjFgJG9bbWMiUXkhIzs3JCQiMjErK0QnM0FSSSEjPCQiMlg5R2ZMNzpYIiEjOzckJCIxLCtdUGcmemEoISM7JCIyW0M3QiRcakJAISM7NyQkIi5EIik+RU4/IiEjNyQiMWEtS1pJYyN6IyEjOjckJCIyKCoqXGkrcGI+OyEjOyQiMTA8cFFddTdNISM6NyQkIi5EY2ZLLjAjISM3JCIxWlJUb3YhXDAlISM6NyQkIi5EY1ZUZVwjISM3JCIxIm92XyhILj5aISM6NyQkIjIpKipcN0I6IypSSCEjOyQiMTNcYEImRzVRJiEjOjckJCIyLysrdlo3blIkISM7JCIxL2ReOUIoPjEnISM6NyQkIjIvK10oKT5hISp6JCEjOyQiMUBxRHd6dWhtISM6NyQkIjIwKytEbyYqPkQlISM7JCIxQCcpKltKYHBMKCEjOjckJCIsdilwejFaISM1JCIxeElgNjYkXCwpISM6NyQkIix2SiMzWF4hIzUkIjFfNiNbTSlHbycpISM6NyQkIi5EY1EoM1ZiISM3JCIxI0cpUTAkKmZoIyohIzo3JCQiLHYhUk47ZyEjNSQiMUh0d0JNNW4qKiEjOjckJCIyMCsrK3Rzc1QnISM7JCIyVT47YHp2azAiISM6NyQkIi5EMWpPTylvISM3JCIybjhxXUEoKmY3IiEjOjckJCIyJikqKioqXHFaa0goISM7JCIyKTNgTnFhYCg9IiEjOjckJCIxLF1pISlvT1x4ISM6JCIyL00/KHlFMGI3ISM6NyQkIi52PS9gMT0pISM3JCIxKnk3di1YJD44ISM5NyQkIjEsK0QxXGxJJykhIzokIjI6VW4oW3RVJ1EiISM6NyQkIi5ESjQoKlEvKiEjNyQiMSVcODAnKUghWzkhIzk3JCQiLXY9WWoqWyohIzYkIjEsPmsiXHdXXiIhIzk3JCQiMSoqXGlOMWpfKiohIzokIjFbKSkqbydlXCRlIiEjOTckJCIvdiRIK25iLiIhIzckIjJZI1wmZlZ4TmsiISM6NyQkIi5EO08nNHo1ISM2JCIxZVFfd3JZMzwhIzk3JCQiLCZbazFDNiEiKiQiMVxzKnppL2J4IiEjOTckJCIyKSoqXDcwNDFvNiEjOiQiMmQ/eGNwKDNUPSEjOjckJCIvRDEveWk1NyEjNyQiMXZnKGVvVVghPiEjOTckJCIudnkzIip5RCIhIzYkIjJ2RDdQaikqXCg+ISM6NyQkIjIqKioqKlxdI2YuSSIhIzokIjIvJXpcT2tJUT8hIzo3JCQiLXZ5SXFYOCEjNSQiMnlmXihHNSFmNSMhIzo3JCQiMioqXGlsSyJ6J1EiISM6JCIyYj5sZXhecjsjISM6NyQkIi1ETTtySjkhIzUkIjJvdS00LTpUQiMhIzo3JCQiMikqXGlJOXlSWiIhIzokIjJ3TzpmQEFySCMhIzo3JCQiMi1dN0wjKWUiPTohIzokIjJtZCwlKilHKUhPIyEjOjckJCIyKSoqKlxQXmA4YyIhIzokIjIoXCVwZnR0dFUjISM6NyQkIjIpKlxpOk51bGciISM6JCIxUkgqWylceSVcIyEjOTckJCIrJFJGLGwiISIpJCIxVktpSCtyZkQhIzk3JCQiLnYkb21tJXAiISM2JCIyaic+PmxfNUVFISM6NyQkIi9EMU9yJClRPCEjNyQiMidcem85MiY+cCMhIzo3JCQiLCYzX1V6PCEiKiQiMnYiNExXZVhfRiEjOjckJCIudik+UCVmIz0hIzYkIjFPQCE0Yiw9I0chIzk3JCQiK0ovYm49ISIpJCIxL2NgeF4jUSlHISM5NyQkIi9EMWo9Ij4iPiEjNyQiMS42JCoqSGIqXEghIzk3JCQiMi1dUFJUclYmPiEjOiQiMSQzeF1FXUssJCEjOTckJCIkKyMhIiIkIjIpKioqUiNcJHA3MyQhIzotJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIjNSEiIi0lKlRISUNLTkVTU0c2IyIiJC0lJ0NVUlZFU0c2JTdmbjckJCEjNSEiIiQhMVlJVSgqKioqKioqXCEjOzckJCExLCsrRHdmQWEhIzskIjJHISoqMylbIkguOSEjPDckJCEyJSoqKlxQSEMpUjkhIzwkIjEiUWM6OCEzV3ghIzs3JCQiMjErK0QnM0FSSSEjPCQiMiV5VVlpPyZcWSIhIzs3JCQiMSwrXVBnJnphKCEjOyQiMVtvayhIIW9PQCEjOjckJCIuRCIpPkVOPyIhIzckIjFGQSdRPWZGIUchIzo3JCQiMigqKlxpK3BiPjshIzskIjJ4PiJvcl9QRU0hIzs3JCQiLkRjZksuMCMhIzckIjByNl0nKjNmMiUhIzk3JCQiLkRjVlRlXCMhIzckIjE5YVZGXDVZWiEjOjckJCIyKSoqXDdCOiMqUkghIzskIjFcbChvPEAuVCYhIzo3JCQiMi8rK3ZaN25SJCEjOyQiMWF6ImY0IWYjNCchIzo3JCQiMi8rXSgpPmEhKnokISM7JCIxZl44cUhwJ3AnISM6NyQkIjIwKytEbyYqPkQlISM7JCIxeVF0LCJlN1EoISM6NyQkIix2KXB6MVohIzUkIjElKTRra3BYbSEpISM6NyQkIix2SiMzWF4hIzUkIjEiSDdhOkotcikhIzo3JCQiLkRjUSgzVmIhIzckIjE8d3kpZidlbSMqISM6NyQkIix2IVJOO2chIzUkIjE8bWRETF0iKSkqISM6NyQkIjIwKysrdHNzVCchIzskIjJRTkQ2Iyk+bS4iISM6NyQkIi5EMWpPTylvISM3JCIxIkgwYC9JOjQiISM5NyQkIjImKSoqKipccVprSCghIzskIjJ5SyQ9MyVHSTkiISM6NyQkIjEsXWkhKW9PXHghIzokIjJqLyFba1pGMzchIzo3JCQiLnY9L2AxPSkhIzckIjJeKWZCOmR4I0ciISM6NyQkIjEsK0QxXGxJJykhIzokIjJsTHUjeXNvczghIzo3JCQiLkRKNCgqUS8qISM3JCIyYiYqKVFMcigqZTkhIzo3JCQiLXY9WWoqWyohIzYkIjFhdDw2WCo0YSIhIzk3JCQiMSoqXGlOMWpfKiohIzokIjJGSWtnMVtqZiIhIzo3JCQiL3YkSCtuYi4iISM3JCIyTydwOTsiMzBpIiEjOjckJCIuRDtPJzR6NSEjNiQiMThwJFFKNUluIiEjOTckJCIsJltrMUM2ISIqJCIyPTtfZWshM04+ISM6NyQkIjIpKipcNzA0MW82ISM6JCIybF9RXE4zRiJHISM6NyQkIi9EMS95aTU3ISM3JCIxRjpxbCdcOy0mISM5NyQkIi52eTMiKnlEIiEjNiQiMkg6Z0U4YzMyIiEjOTckJCIyKioqKipcXSNmLkkiISM6JCIxPDUicGglb0JAISM4NyQkIi12eUlxWDghIzUkIjFIOiJ6KEdpXFUhIzg3JCQiMioqXGlsSyJ6J1EiISM6JCIxQFAwQU4jXGsoISM4NyQkIi1ETTtySjkhIzUkIjJLdCJ5Z0MoNFIiISM4NyQkIjIpKlxpSTl5UloiISM6JCIxRCkqbyZHak1OIyEjNzckJCIyLV03TCMpZSI9OiEjOiQiMkMqXCEpZm9lVFIhIzg3JCQiMikqKipcUF5gOGMiISM6JCIvOid5P05nTCchIzU3JCQiMikqXGk6TnVsZyIhIzokIjIpKXApeWVAPzg1ISM3NyQkIiskUkYsbCIhIikkIjJucCQqNGAlUWI6ISM3NyQkIi52JG9tbSVwIiEjNiQiMVQoKj0oPSdHZ0IhIzY3JCQiL0QxT3IkKVE8ISM3JCIxalRiOjQwLU4hIzY3JCQiLCYzX1V6PCEiKiQiMVUjZSkqKVIlcCZcISM2NyQkIi92PWtXby09ISM3JCIxLjpRdShSRSwnISM2NyQkIi52KT5QJWYjPSEjNiQiMWFyUWMpNEpFKCEjNjckJCIvdlZ2cXVZPSEjNyQiMTRsUW9ia3ImKSEjNjckJCIrSi9ibj0hIikkIjJMPUolUSRSJjM1ISM2NyQkIjItREpxOUooKik9ISM6JCIydDglW1FlbSY+IiEjNjckJCIvRDFqPSI+Ij4hIzckIjJXJylISTVLSVQiISM2NyQkIjBEInldbl9BPiEjOCQiMT1ELU1rJypHOiEjNTckJCIsJlE7OUw+ISIqJCIyVXZhRilRR2A7ISM2NyQkIjB2PWlfY1AlPiEjOCQiMkZ5NGxXN2x5IiEjNjckJCIyLV1QUlRyViY+ISM6JCIyKVxeP2cnMyNIPiEjNjckJCIxREpYZyZ5ZCc+ISM5JCIxXGc2cih6UDQjISM1NyQkIjB2b3BxJj14PiEjOCQiMjdXLVgzSTJGIyEjNjckJCIxdlZbYEdmKSk+ISM5JCIwWFhoNXIzWSMhIio3JCQiJCsjISIiJCIwLiRIOSsxbEUhIiotJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiM1ISIiJCIiISEiIi0lKlRISUNLTkVTU0c2IyIiJC0lJ0NVUlZFU0c2JjctNyQkISM1ISIiJCEiJiEiIjckJCIiISEiIiQiIzUhIiI3JCQiIzUhIiIkIiNEISIiNyQkIiM/ISIiJCIjUyEiIjckJCIjSSEiIiQiI2IhIiI3JCQiI1MhIiIkIiNxISIiNyQkIiNdISIiJCIjJikhIiI3JCQiIyEpISIiJCIkRCIhIiI3JCQiIyEqISIiJCIkWCIhIiI3JCQiJCsiISIiJCIkZyIhIiI3JCQiJDUiISIiJCIkdiIhIiItJSZDT0xPUkc2JiUkUkdCRyQiIzUhIiIkIiIhISIiJCIiISEiIi0lJ1NZTUJPTEc2JCUqX1NPTElEQk9YRyIjOy0lJlNUWUxFRzYjJSZQT0lOVEctJSVWSUVXRzYkOyQhIzUhIiIkIiQrIyEiIiUoREVGQVVMVEctJiUmX0FYSVNHNiMiIiI2Ji0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqTElORVNUWUxFRzYjIiIhLSUqVEhJQ0tORVNTRzYjIiIiLSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0mJSZfQVhJU0c2IyIiIzYmLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiItJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUrQVhFU0xBQkVMU0c2JS1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRIng2Ii8lJ2ZhbWlseUdRKERFRkFVTFQ2Ii8lJXNpemVHUSMxMDYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiJRITYiLSUlRk9OVEc2JSUoREVGQVVMVEclKERFRkFVTFRHIiM1LSUlUk9PVEc2Jy0lKUJPVU5EU19YRzYjJCIkcSUhIiItJSlCT1VORFNfWUc2IyQiIyEpISIiLSUtQk9VTkRTX1dJRFRIRzYjJCIlP00hIiItJS5CT1VORFNfSEVJR0hURzYjJCIlK04hIiItJSlDSElMRFJFTkc2Ig==

quad:=LeastSquares(newpts,x,curve=a*x^2+b*x+c);

LCgkIitQUXBANSEiKiIiIkkieEc2IiQiK2tOJ29ZIkYlKiRGJyIiIyQiKywlXEdMIyEjNw==

spl:=Spline(newpts,x,degree=1);

LUkqcGllY2V3aXNlRyUqcHJvdGVjdGVkRzY1MkkieEc2IiIiISwmJCIiIkYpRixGJyQiIzohIiIyRidGLEYqMkYnIiIjLCYkIiM1Ri9GLEYnRi0yRiciIiRGMzJGJyIiJUYzMkYnIiImRjMyRiciIiksJiQiK0xMTEw9ISIqRixGJyQiK0xMTEw4RkEyRiciIiosJiQhI05GL0YsRickIiM/Ri8yRidGNUYzRjM=

display(plot({line,pol,quad,spl},x=-1..11,thickness=3,color=[blue,green,cyan,brown]),plot(newpts,style=point,symbol=solidbox,symbolsize=16));

LSUlUExPVEc2LC0lJ0NVUlZFU0c2JTdTNyQkISM1ISIiJCExKioqKipmaCtFI1whIzs3JCQhKTpNJVEoISIpJCEyL0YqZTNwU0I1ISM8NyQkIStENVozXiEjNSQiMVhMTihvZyNwQiEjOzckJCEqbDshXEQhIiokIjFweixGLW4lPSchIzs3JCQiMiUpMysrK1gudSMhIz4kIjJkVXl4MlBEKyIhIzs3JCQiMiUqKioqKlxGeTpmIyEjPCQiMiNmZSNcciJ5JVEiISM7NyQkIjEqKioqKlxkJyopb1whIzskIjJZKDRHQzk8UjwhIzs3JCQiMSoqKioqXCg+Wkl1ISM7JCIyQzZsN1dAaDUjISM7NyQkIit2bkJ3KiohIzUkIjFLdnZ0O2kmWyMhIzo3JCQiLHZzJFFeNyEjNSQiMmRQIj5XWyFSJ0chIzs3JCQiMi0rKys4MkNeIiEjOyQiMSlvJSo9ZTpJRCQhIzo3JCQiKyZvO0J1IiEiKiQiMlg6aDlgV2RmJCEjOzckJCIqUlM2KyMhIikkIjIkcD9hYGhkIilSISM7NyQkIiomby1oQSEiKSQiMSg9THQuIyoqb1YhIzo3JCQiKmh2OV4jISIpJCIxOk55RiFSQnUlISM6NyQkIix2MjIqUUYhIzUkIjEvTFJpUVAiMyYhIzo3JCQiKjRYJDRJISIpJCIxO2MqRzI+WFsmISM6NyQkIipjVCVRSyEiKSQiMSZSTm8lZi5FZSEjOjckJCIsdkBOXF0kISM1JCIxbV1yZjdJQmkhIzo3JCQiMigqKioqKipmcyMzdSQhIzskIjFQLmp2JFtcZCchIzo3JCQiMi4rK3ZAUScqKlIhIzskIjEpM2ZbKjR3Z3AhIzo3JCQiLERYKDNZVSEjNSQiMV87NVhlOUd0ISM6NyQkIit2OEIuWCEiKiQiMUEiR15SdDlyKCEjOjckJCIsdndwJFJaISM1JCIxQkcnenUoW2ohKSEjOjckJCIrRHAyJSpcISIqJCIxJEgjKUh5Iz1WJSkhIzo3JCQiLHYyWSdlXyEjNSQiMCo+LTgjenYkKSkhIzk3JCQiLERJYSopWyYhIzUkIjEpKjR4a0MhND0qISM6NyQkIismXCRwUGQhIiokIjFmeEkjKTRxXiYqISM6NyQkIipVbFkqZiEiKSQiMEclZUxAeE0qKiEjOTckJCIrOkIxWWkhIiokIjEqb0s/RGA0LiIhIzk3JCQiLHZ0LCQqWychIzUkIjJIL0trQzhzMSIhIzo3JCQiK1h3UGZuISIqJCIyJT1GMCl5dHU1IiEjOjckJCIxKioqKioqZkcwLXEhIzokIjFCKWVgblxPOSIhIzk3JCQiKlhnNkUoISIpJCIydGlZbTp2QT0iISM6NyQkIjImKioqKlxQLyZmXCghIzskIjEkb29OZXZzQCIhIzk3JCQiMSwrKzV6al94ISM6JCIyRXI9XUhTYkQiISM6NyQkIjEqKioqXDwzOyUqeiEjOiQiMmZcYWtTVzpIIiEjOjckJCIsdiU9aVkjKSEjNSQiMk8vPXEyeiJIOCEjOjckJCIqbFtNXCkhIikkIi9hW1lRKGZPIiEjNzckJCIsdlZgPXYpISM1JCIyYSdmOmdYXC85ISM6NyQkIionenMrISohIikkIjFSaHImZSVmVDkhIzk3JCQiKzAiUV9EKiEiKiQiMVFvSz9aYHo5ISM5NyQkIix2eFN3XSohIzUkIjIpKlJuOXBncl4iISM6NyQkIippcyZSKCohIikkIjI6UWo3ME48YiIhIzo3JCQiLCZvI1IwKyIhIiokIjJ5Mi8kcDZPImYiISM6NyQkIitLWEpDNSEiKSQiMiNmZ21USyFvaSIhIzo3JCQiLXZAUm1cNSEjNSQiMmsiXC4kKT1mazshIzo3JCQiLURBbCNSMiIhIzUkIjJbIXBhTS93KzwhIzo3JCQiJDUiISIiJCItayNbRydSPCEjNS0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiM1ISIiLSUqVEhJQ0tORVNTRzYjIiIkLSUnQ1VSVkVTRzYlN1M3JCQhIzUhIiIkITI5KyEqZnhvJEdXISM8NyQkISk6TSVRKCEiKSQhMm0qW3BGV2lAZyEjPTckJCErRDVaM14hIzUkIjEqKnpjcWlmSEYhIzs3JCQhKmw7IVxEISIqJCIxTy9wKFslUnprISM7NyQkIjIlKTMrKytYLnUjISM+JCIyMCE9dm1OckQ1ISM7NyQkIjIlKioqKipcRnk6ZiMhIzwkIjIuP05AcCoqPlMiISM7NyQkIjEqKioqKlxkJyopb1whIzskIjFDTCgpXCJSNnYiISM6NyQkIjEqKioqKlwoPlpJdSEjOyQiMmB2VyIzMiRINiMhIzs3JCQiK3ZuQncqKiEjNSQiMkZ5eGFWJFIoWyMhIzs3JCQiLHZzJFFeNyEjNSQiMS14PCQ9YzQnRyEjOjckJCIyLSsrKzgyQ14iISM7JCIyTig+NF9bX1hLISM7NyQkIismbztCdSIhIiokIjE2TTdha14lZSQhIzo3JCQiKlJTNisjISIpJCIyKDQsaC5lVm1SISM7NyQkIiomby1oQSEiKSQiMTMnXEAkeUJdViEjOjckJCIqaHY5XiMhIikkIjBjLDpQKy9zJSEjOTckJCIsdjIyKlFGISM1JCIxUSd5XidwemNdISM6NyQkIio0WCQ0SSEiKSQiMVV0W2AiPnJYJiEjOjckJCIqY1QlUUshIikkIjE5JT5CYDZseiYhIzo3JCQiLHZATlxdJCEjNSQiMSRbWmJsODs+JyEjOjckJCIyKCoqKioqKmZzIzN1JCEjOyQiMSYqKVJfV0E7YSchIzo3JCQiMi4rK3ZAUScqKlIhIzskIjEqRyJHWWkkZiNwISM6NyQkIixEWCgzWVUhIzUkIjFRPillTiU9I0goISM6NyQkIit2OEIuWCEiKiQiMGV6MiR6aXV3ISM5NyQkIix2d3AkUlohIzUkIjEwZWFiQTVFISkhIzo3JCQiK0RwMiUqXCEiKiQiMUtSYWFrXTAlKSEjOjckJCIsdjJZJ2VfISM1JCIxOyNHXUdAKip6KSEjOjckJCIsRElhKilbJiEjNSQiMHVtJFJrX1YiKiEjOTckJCIrJlwkcFBkISIqJCIwbiR5J1IyXF4qISM5NyQkIipVbFkqZiEiKSQiMSZISnZVJCkpKSopKiEjOjckJCIrOkIxWWkhIiokIjJpIT5mP0ZbRjUhIzo3JCQiLHZ0LCQqWychIzUkIjJWYzI+RyYpUTEiISM6NyQkIitYd1BmbiEiKiQiMkZOSTVZT1Y1IiEjOjckJCIxKioqKioqZkcwLXEhIzokIjJtNz9GbjcyOSIhIzo3JCQiKlhnNkUoISIpJCIxUSVlUFAjZXo2ISM5NyQkIjImKioqKlxQLyZmXCghIzskIjJwJCozQjpKW0AiISM6NyQkIjEsKys1empfeCEjOiQiMiwpcCszb1JgNyEjOjckJCIxKioqKlw8MzslKnohIzokIjB5QGI6NygqRyIhIzg3JCQiLHYlPWlZIykhIzUkIjEnW0UzRCx4SyIhIzk3JCQiKmxbTVwpISIpJCIydSk9dSVHcltPIiEjOjckJCIsdlZgPXYpISM1JCIybydcJj1LOlFTIiEjOjckJCIqJ3pzKyEqISIpJCIyKD0ncCxeXzhXIiEjOjckJCIrMCJRX0QqISIqJCIxOygpWyVmcCh6OSEjOTckJCIsdnhTd10qISM1JCIyLT5qUVQpKnleIiEjOjckJCIqaXMmUigqISIpJCIyVFBsVy9oSGIiISM6NyQkIiwmbyNSMCsiISIqJCIyJ3k+WFV3PCRmIiEjOjckJCIrS1hKQzUhIikkIjJGKHpLeGQ8SDshIzo3JCQiLXZAUm1cNSEjNSQiMiZSRCN5TSdlbjshIzo3JCQiLURBbCNSMiIhIzUkIjIwTlZkQ3lWcSIhIzo3JCQiJDUiISIiJCIvQCY9eVlSdSIhIzctJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiM1ISIiJCIiISEiIi0lKlRISUNLTkVTU0c2IyIiJC0lJ0NVUlZFU0c2JTdWNyQkISM1ISIiJCExWUlVKCoqKioqKipcISM7NyQkISk6TSVRKCEiKSQhMSxZKj43J1wuOyEjOzckJCErRDVaM14hIzUkIjI9dip6SChlKyE+ISM8NyQkISpsOyFcRCEiKiQiMT1CWSlwKGUlKWYhIzs3JCQiMiUpMysrK1gudSMhIz4kIi96ciFSZlUrIiEjODckJCIyJSoqKioqXEZ5OmYjISM8JCIyZj9HKSk0SHRSIiEjOzckJCIxKioqKipcZCcqKW9cISM7JCIyJil5XmIzIypSdiIhIzs3JCQiMSoqKioqXCg+Wkl1ISM7JCIyeEZqKm9TRz5AISM7NyQkIit2bkJ3KiohIzUkIjJbLG14SnRrXCMhIzs3JCQiLHZzJFFeNyEjNSQiMmxmbS0+d1QoRyEjOzckJCIyLSsrKzgyQ14iISM7JCIyVztpYmw8X0UkISM7NyQkIismbztCdSIhIiokIjJgSCNlPzJLNk8hIzs3JCQiKlJTNisjISIpJCIxaGU/JjM/PCslISM6NyQkIiomby1oQSEiKSQiMjpQKUc+JUdMUiUhIzs3JCQiKmh2OV4jISIpJCIxLF5QIWVdJnBaISM6NyQkIix2MjIqUUYhIzUkIjAmKnBvUF0rNiYhIzk3JCQiKjRYJDRJISIpJCIxY0wnKkgkW1JeJiEjOjckJCIqY1QlUUshIikkIjEyTGBNPyhmJmUhIzo3JCQiLHZATlxdJCEjNSQiMiZcKFs6c3BZRCchIzs3JCQiMigqKioqKipmcyMzdSQhIzskIjF3VXkpbyYqKjNtISM6NyQkIjIuKyt2QFEnKipSISM7JCIxc3g2ZUpYKipwISM6NyQkIixEWCgzWVUhIzUkIjFrWiZlMDxCUCghIzo3JCQiK3Y4Qi5YISIqJCIxV1glZmhxM3coISM6NyQkIix2d3AkUlohIzUkIjFvaFFMUDI6IikhIzo3JCQiK0RwMiUqXCEiKiQiMCdcVEl6TSJcKSEjOTckJCIsdjJZJ2VfISM1JCIwJz1yWytFcykpISM5NyQkIixESWEqKVsmISM1JCIxRWpKQyMzSD4qISM6NyQkIismXCRwUGQhIiokIjFaLUg/NiFlXyohIzo3JCQiKlVsWSpmISIpJCIxaSFcYClwV2EpKiEjOjckJCIrOkIxWWkhIiokIjJDPjpPQ05pLCIhIzo3JCQiLHZ0LCQqWychIzUkIjJQIypbKXBTNVg1ISM6NyQkIitYd1BmbiEiKiQiMlRTJ3AjXCp6dzUhIzo3JCQiMSoqKioqKmZHMC1xISM6JCIyTHIpKSpvYyFlNSIhIzo3JCQiKlhnNkUoISIpJCIyUV9edl0oUlE2ISM6NyQkIjImKioqKlxQLyZmXCghIzskIjJkXHVHakcuPCIhIzo3JCQiMSwrKzV6al94ISM6JCIyZEIvX1oiejM3ISM6NyQkIjEqKioqXDwzOyUqeiEjOiQiMjtgI3A+InkqWzchIzo3JCQiLHYlPWlZIykhIzUkIjJSY0lvI1FHJkgiISM6NyQkIipsW01cKSEiKSQiMnMoPV06SEhXOCEjOjckJCIsdlZgPXYpISM1JCIyRSc9KWYkKUciKVIiISM6NyQkIionenMrISohIikkIjJHeihmaShcLFgiISM6NyQkIiswIlFfRCohIiokIjJhYzRTNSpSKzohIzo3JCQiLHZ4U3ddKiEjNSQiMiNwXXMpeTJRYSIhIzo3JCQiKmlzJlIoKiEiKSQiMSNIQV8xKlJ2OiEjOTckJCIsJm8jUjArIiEiKiQiMkdIOCUpPSZSKzshIzo3JCQiK0tYSkM1ISIpJCIyLkU9QWE3WWgiISM6NyQkIi52b0EqKXAuIiEjNiQiMmlsNWdiLDhpIiEjOjckJCItdkBSbVw1ISM1JCIxciwwLkhxSDshIzk3JCQiK0Ffemg1ISIpJCIyWDp5WGcnKT1rIiEjOjckJCItREFsI1IyIiEjNSQiMSk9ISk+T0s5bSIhIzk3JCQiMiwrRDZFanAzIiEjOiQiMXNURHJeYiZwIiEjOTckJCIkNSIhIiIkIjJPNT4nNDErXTwhIzotJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiM1ISIiJCIjNSEiIi0lKlRISUNLTkVTU0c2IyIiJC0lJ0NVUlZFU0c2JTdTNyQkISM1ISIiJCEiJiEiIjckJCEpOk0lUSghIikkITItKysrRDdsMiIhIzw3JCQhK0Q1WjNeISM1JCIyLisrRFkkSFBCISM8NyQkISpsOyFcRCEiKiQiK0RdWndoISM1NyQkIjIlKTMrKytYLnUjISM+JCIsdl41VCsiISM1NyQkIjIlKioqKipcRnk6ZiMhIzwkIjIpKioqXDd1TygpUSIhIzs3JCQiMSoqKioqXGQnKilvXCEjOyQiLUQnW01gdSIhIzY3JCQiMSoqKioqXCg+Wkl1ISM7JCIyKSoqKlxpenFYNiMhIzs3JCQiK3ZuQncqKiEjNSQiLUQ7YlYnXCMhIzY3JCQiLHZzJFFeNyEjNSQiMigqKipcN2Z2cShHISM7NyQkIjItKysrODJDXiIhIzskIjIwKytdcDUnb0shIzs3JCQiKyZvO0J1IiEiKiQiLHYtdk1oJCEjNTckJCIqUlM2KyMhIikkIjIwKytdZTU8KyUhIzs3JCQiKiZvLWhBISIpJCIrdi1hIlIlISIqNyQkIipodjleIyEiKSQiMjArK11UOHN3JSEjOzckJCIsdjIyKlFGISM1JCItRDsxTzNeISM2NyQkIio0WCQ0SSEiKSQiK053LDliISIqNyQkIipjVCVRSyEiKSQiKk1pdyZlISIpNyQkIix2QE5cXSQhIzUkIi1ERUdTZGkhIzY3JCQiMigqKioqKipmcyMzdSQhIzskIjImKioqKioqKikzQzZtISM7NyQkIjIuKyt2QFEnKipSISM7JCIxLCtERXRYKipwISM6NyQkIixEWCgzWVUhIzUkIjEsK3Z5NjhwdCEjOjckJCIrdjhCLlghIiokIjIxKytEMVpbdighIzs3JCQiLHZ3cCRSWiEjNSQiLUReWTA0IikhIzY3JCQiK0RwMiUqXCEiKiQiMSoqKipcKFE6NlwpISM6NyQkIix2MlknZV8hIzUkIjBYIWVNOSdbJSkpISM5NyQkIixESWEqKVsmISM1JCIxOnEkeTFSPjoqISM6NyQkIismXCRwUGQhIiokIjFvMlVDOGYkWyohIzo3JCQiKlVsWSpmISIpJCIxVyVvd2I/aSMpKiEjOjckJCIrOkIxWWkhIiokIjJYWWU8a1RoLCIhIzo3JCQiLHZ0LCQqWychIzUkIjJsTl1aY3QmWzUhIzo3JCQiK1h3UGZuISIqJCIxYTh1Jm8kZSUzIiEjOTckJCIxKioqKioqZkcwLXEhIzokIjEpZm01UVNwNiIhIzk3JCQiKlhnNkUoISIpJCIyemlDKGYhKVteNiEjOjckJCIyJioqKipcUC8mZlwoISM7JCIyI29NUSJSJHojPSIhIzo3JCQiMSwrKzV6al94ISM6JCIyYUMzeFE9cUAiISM6NyQkIjEqKioqXDwzOyUqeiEjOiQiMnlfTD9XQCNcNyEjOjckJCIsdiU9aVkjKSEjNSQiLCZwVksqSCIhIio3JCQiKmxbTVwpISIpJCIyKioqKioqKkgoKm9bOCEjOjckJCIsdlZgPXYpISM1JCIyKioqKipcKG9xLlMiISM6NyQkIionenMrISohIikkIjItKytTPjQsWCIhIzo3JCQiKzAiUV9EKiEiKiQiLXY6ZEcpWyIhIzU3JCQiLHZ4U3ddKiEjNSQiMiwrRG02WWhfIiEjOjckJCIqaXMmUigqISIpJCIrJCplJDRjIiEiKTckJCIsJm8jUjArIiEiKiQiLXYtKjMzZyIhIzU3JCQiK0tYSkM1ISIpJCIrKXpya2oiISIpNyQkIi12QFJtXDUhIzUkIi5ERSllXHU7ISM2NyQkIi1EQWwjUjIiISM1JCIudkx5KikzciIhIzY3JCQiJDUiISIiJCIkdiIhIiItJSZDT0xPUkc2JiUkUkdCRyQiMV9NbVglKWVxayEjOyQiMndtb1YoKWVxayIhIzwkIjJ3bW9WKCllcWsiISM8LSUqVEhJQ0tORVNTRzYjIiIkLSUnQ1VSVkVTRzYmNy03JCQhIzUhIiIkISImISIiNyQkIiIhISIiJCIjNSEiIjckJCIjNSEiIiQiI0QhIiI3JCQiIz8hIiIkIiNTISIiNyQkIiNJISIiJCIjYiEiIjckJCIjUyEiIiQiI3EhIiI3JCQiI10hIiIkIiMmKSEiIjckJCIjISkhIiIkIiREIiEiIjckJCIjISohIiIkIiRYIiEiIjckJCIkKyIhIiIkIiRnIiEiIjckJCIkNSIhIiIkIiR2IiEiIi0lJkNPTE9SRzYmJSRSR0JHJCIjNSEiIiQiIiEhIiIkIiIhISIiLSUnU1lNQk9MRzYkJSpfU09MSURCT1hHIiM7LSUmU1RZTEVHNiMlJlBPSU5URy0lJVZJRVdHNiQ7JCEjNSEiIiQiJDUiISIiJShERUZBVUxURy0mJSZfQVhJU0c2IyIiIjYmLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiItJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSYlJl9BWElTRzYjIiIjNiYtJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIS0lKlRISUNLTkVTU0c2IyIiIi0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJStBWEVTTEFCRUxTRzYlLUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NVEieDYiLyUnZmFtaWx5R1EoREVGQVVMVDYiLyUlc2l6ZUdRIzEwNiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2IlEhNiItJSVGT05URzYlJShERUZBVUxURyUoREVGQVVMVEciIzUtJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiRTIiEiIi0lKUJPVU5EU19ZRzYjJCIkSSIhIiItJS1CT1VORFNfV0lEVEhHNiMkIiVJUCEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiVdTiEiIi0lKUNISUxEUkVORzYi

subs(x=12.3,{line,pol,quad}); subs(x=12.3,1.0+1.5*x);

PCUkIitfNVVMPiEiKSQiK1pfcVQ+RiUkIis1KFJyIW9GJQ==

JCIlWD4hIiM=

Problem 11 pts:=[[1.02,-4.30], [1.00,-2.12], [0.99, 0.52], [1.03, 2.51], [1.00, 3.34], [1.02, 5.30]];

Nyg3JCQiJC0iISIjJCEkSSVGJjckJCIkKyJGJiQhJDcjRiY3JCQiIyoqRiYkIiNfRiY3JCQiJC4iRiYkIiReI0YmNyRGKiQiJE0kRiY3JEYkJCIkSSZGJg==

plot(pts,style=point,symbol=solidbox,symbolsize=16);

LSUlUExPVEc2Jy0lJ0NVUlZFU0c2JDcoNyQkIiQtIiEiIyQhI1YhIiI3JCQiIzUhIiIkISQ3IyEiIzckJCIjKiohIiMkIiNfISIjNyQkIiQuIiEiIyQiJF4jISIjNyQkIiM1ISIiJCIkTSQhIiM3JCQiJC0iISIjJCIjYCEiIi0lJkNPTE9SRzYmJSRSR0JHJCIjNSEiIiQiIiEhIiIkIiIhISIiLSUrQVhFU0xBQkVMU0c2JFEhNiJRITYiLSUnU1lNQk9MRzYkJSpfU09MSURCT1hHIiM7LSUmU1RZTEVHNiMlJlBPSU5URy0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJEklISIiLSUpQk9VTkRTX1lHNiMkIiRJIiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJT9NISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJVNQISIiLSUpQ0hJTERSRU5HNiI=

E:=(m,b,data)->sum((m*data[i][1]+b-data[i][2])^2 , i=1..nops(data) );

Zio2JUkibUc2IkkiYkdGJUklZGF0YUdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUkkc3VtRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliR0YlNiQqJCwoKiY5JCIiIiYmOSY2I0kiaUdGJTYjRjVGNUY1OSVGNSZGNzYjIiIjISIiRj8vRjo7RjUtSSVub3BzR0YuNiNGOEYlRiVGJQ==

E(m,b,pts);

LC4qJCwoSSJtRzYiJCIkLSIhIiNJImJHRiYiIiIkIiRJJUYpRisiIiNGKyokLChGJSQiJCsiRilGKkYrJCIkNyNGKUYrRi5GKyokLChGJSQiIyoqRilGKkYrJCEjX0YpRitGLkYrKiQsKEYlJCIkLiJGKUYqRiskISReI0YpRitGLkYrKiQsKEYlRjFGKkYrJCEkTSRGKUYrRi5GKyokLChGJUYnRipGKyQhJEkmRilGK0YuRis=

solution:=solve({diff(E(m,b,pts),m),diff(E(m,b,pts),b)},{m,b});

PCQvSSJiRzYiJCErbm07eEkhIikvSSJtR0YlJCIrTExMTEpGKA==

line:=subs(solution,m*x+b);

LCZJInhHNiIkIitMTExMSiEiKSQhK25tO3hJRiciIiI=

with(plots): display(plot(pts,style=point,symbol=solidbox,symbolsize=16),plot(line, x=0..2,color=blue,thickness=3));

LSUlUExPVEc2KS0lJ0NVUlZFU0c2JjcoNyQkIiQtIiEiIyQhI1YhIiI3JCQiIzUhIiIkISQ3IyEiIzckJCIjKiohIiMkIiNfISIjNyQkIiQuIiEiIyQiJF4jISIjNyQkIiM1ISIiJCIkTSQhIiM3JCQiJC0iISIjJCIjYCEiIi0lJkNPTE9SRzYmJSRSR0JHJCIjNSEiIiQiIiEhIiIkIiIhISIiLSUnU1lNQk9MRzYkJSpfU09MSURCT1hHIiM7LSUmU1RZTEVHNiMlJlBPSU5URy0lJ0NVUlZFU0c2JTdTNyQkIiIhISIiJCErbm07eEkhIik3JCQiMUxMTEwzVmZWISM8JCEybzNxdjtyMCVIISM6NyQkIjFubW0iSFtEOikhIzwkITEnRykpUTo/PCNHISM5NyQkIjJMTEwkZTAkPUMiISM8JCEya2hqKGUoZiEpbyMhIzo3JCQiMkxMTCQzUkJyOyEjPCQhMCR6eHZMXmBEISM4NyQkIjJsbW0iempmKTQjISM8JCEydlJrI29rZz5DISM6NyQkIjJLTExlNDtbXCMhIzwkITIkXD9zQndYJkgjISM6NyQkIjIoKioqKlxpJ3ldIUghIzwkITImcCw9J28zcDsjISM6NyQkIjFMTCRlenMkSEwhIzskITEsSzMiPmpSLiMhIzk3JCQiMiYqKioqXDdpSV9QISM8JCExeF0oUTBXOSE+ISM5NyQkIjFubW07X00oPSUhIzskITJBLSVlSz04bDwhIzo3JCQiMk1MTCQzeV9xWCEjPCQhMXRYKFEmejFYOyEjOTckJCIqbCs+KyYhIiokITJNMStOWS8qNDohIzo3JCQiKnZXXVYmISIqJCEyWiR5JSl5Zj11OCEjOjckJCIqTmZDJmUhIiokITIoWzwnMzInUlY3ISM6NyQkIjFMTCRlejY6QichIzskITJERjhYSUVZNyIhIzo3JCQiMW1tbTs9QyNvJyEjOyQhMCc9blljKFIkKSohIzk3JCQiMW1tbW0jcFMxKCEjOyQhMFpLLChIZVAnKSEjOTckJCIsRE9EIzN2ISM2JCExVUYrcVEqZUMoISM6NyQkIjFtbW1tKHk4IXohIzskITExXHUuSyw5ZyEjOjckJCIsRE9JRkwpISM2JCExdnZGcVtXaVkhIzo3JCQiLHYzek11KSEjNiQhMUZ5dSk+S2FQJCEjOjckJCIxbm1tIkhfPzwqISM7JCEyR3osRFxwRC4jISM7NyQkIjFubTt6aWhsJiohIzskITE6bCp6SkJTKnohIzs3JCQiMUxMTDMjRywqKiohIzskIjE2UzloL04yYCEjOzckJCIyS0wkZXp3NVY1ISM7JCIyQVRuTSZSUDc+ISM7NyQkIi12JFEjXCIzIiEjNiQiMik+XSpRbyU0OkohIzs3JCQiMktMTGUiKltINyIhIzskIjFXek9BZzE5VyEjOjckJCIqZHhkOyIhIikkIjEudVc+MC5jZCEjOjckJCIyKSoqKipcX3FuMjchIzskIjEsVyg0KDQpKW9xISM6NyQkIjIpKioqXGkmcEBbNyEjOyQiMEZSJClHSiJSJCkhIzk3JCQiMikqKioqXDInSEtIIiEjOyQiMWtELSVwRiZcKCohIzo3JCQiMmxtbW1adk9MIiEjOyQiMExXJWZKbyw2ISM4NyQkIit2KydvUCIhIiokIjIqelYoMyFcKnBCIiEjOjckJCIyS0wkZVI8KmZUIiEjOyQiMlE2KikpXHVnZjghIzo3JCQiKyYpSHhlOSEiKiQiMTQvPV9gbCRcIiEjOTckJCIybG07SCFvLSpcIiEjOyQiMTgjUjspUnk+OyEjOTckJCIyKioqKlw3ay42YSIhIzskIjIoKUhPPHVDO3YiISM6NyQkIjJtbW07V1RBZSIhIzskIjFUIkdJPUIwKT0hIzk3JCQiLUQxKjNgaSIhIzYkIjJvKltQMHpZOj8hIzo3JCQiMk1MTEwqenltOyEjOyQiMTpiYjZkVlhAISM5NyQkIjJMTEwzTjEjNDwhIzskIjJjUFRdY1kkeUEhIzo3JCQiMm1tO0hZdDd2IiEjOyQiMUpeYlxvOjVDISM5NyQkIip4RyoqeSIhIikkIjIvcE88Im9GSkQhIzo3JCQiMm1tbVQ2S1UkPSEjOyQiMXl1eGNSNHFFISM5NyQkIjJMTExMYmRRKD0hIzskIjFlSjptT0QleiMhIzk3JCQiLURPbDU7PiEjNiQiMSlINkUiUWpFSCEjOTckJCItdi5VYWM+ISM2JCIyJj15NC8mUUwwJCEjOjckJCIjPyEiIiQiKyoqKipcKj0kISIpLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIzUhIiItJSpUSElDS05FU1NHNiMiIiQtJSVWSUVXRzYkOyQiIiEhIiIkIiM/ISIiJShERUZBVUxURy0mJSZfQVhJU0c2IyIiIjYmLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiItJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSYlJl9BWElTRzYjIiIjNiYtJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIS0lKlRISUNLTkVTU0c2IyIiIi0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJStBWEVTTEFCRUxTRzYlLUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NVEieDYiLyUnZmFtaWx5R1EoREVGQVVMVDYiLyUlc2l6ZUdRIzEwNiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2IlEhNiItJSVGT05URzYlJShERUZBVUxURyUoREVGQVVMVEciIzUtJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiQhXCEiIi0lKUJPVU5EU19ZRzYjJCIjISkhIiItJS1CT1VORFNfV0lEVEhHNiMkIiVJTSEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiVTUSEiIi0lKUNISUxEUkVORzYi