<Text-field style="Heading 1" layout="Heading 1">Initializations</Text-field>

<Text-field style="Heading 1" layout="Heading 1">Example 1: Parameter Estimation with Noisy Data</Text-field>

<Text-field style="Heading 2" layout="Heading 2">The Model</Text-field>

Consider the nonlinear damped oscillator in which the damping force is assumed proportional to the velocity, a system governed by the nonlinear differential equation

DE := diff(y(t),t,t) + a*diff(y(t),t)*abs(diff(y(t),t)) + b*y(t) = 0;

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSNERUdGKC8sKC1JJWRpZmZHRiY2JC1JInlHRig2I0kidEdGKC1JIiRHRiY2JEY2IiIjIiIiKihJImFHRihGOy1GMTYkRjNGNkY7LUkkYWJzR0YmNiNGPkY7RjsqJkkiYkdGKEY7RjNGO0Y7IiIhNyMvLCgtRjE2JEY+RjZGO0Y8RjtGQ0Y7RkU=

The unknown parameters are the damping coefficient LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JlEiYUYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnRjIvRjZRJ25vcm1hbEYn and the spring constant LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JlEiYkYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnRjIvRjZRJ25vcm1hbEYn.

To simulate observational data, we start with the model

LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzY0LUkjbWlHRiQ2JlEieUYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2LlEiJ0YnRjIvRjZRJ25vcm1hbEYnLyUmZmVuY2VHRjQvJSpzZXBhcmF0b3JHRjQvJSlzdHJldGNoeUdGNC8lKnN5bW1ldHJpY0dGNC8lKGxhcmdlb3BHRjQvJS5tb3ZhYmxlbGltaXRzR0Y0LyUnYWNjZW50R0Y0LyUnbHNwYWNlR1EsMC4xMTExMTExZW1GJy8lJ3JzcGFjZUdRJjAuMGVtRidGOC1GOTYuUScmcGx1cztGJ0YyRjxGPkZARkJGREZGRkhGSi9GTVEsMC4yMjIyMjIyZW1GJy9GUEZWLUkjbW5HRiQ2JVEiMkYnRjJGPC1GOTYtUTEmSW52aXNpYmxlVGltZXM7RidGPEY+RkBGQkZERkZGSEZKL0ZNRlFGT0YrRjhGZm4tSShtZmVuY2VkR0YkNiktRiM2JkYrRjhGMkY8RjJGPC9JK21zZW1hbnRpY3NHRiRRJGFic0YnLyUlb3BlbkdRKSZ2ZXJiYXI7RicvJSZjbG9zZUdGZG9GX29GUi1GWTYlUSI0RidGMkY8RmZuRistRjk2LlEiPUYnRjJGPEY+RkBGQkZERkZGSEZKL0ZNUSwwLjI3Nzc3NzhlbUYnL0ZQRl5wLUZZNiVRIjBGJ0YyRjxGMkY8

LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYpLUkjbWlHRiQ2I1EhRictRiM2Jy1GLDYmUSJ5RicvJSdpdGFsaWNHUSV0cnVlRicvJStleGVjdXRhYmxlR1EmZmFsc2VGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUTAmQXBwbHlGdW5jdGlvbjtGJy9GO1Enbm9ybWFsRicvJSZmZW5jZUdGOS8lKnNlcGFyYXRvckdGOS8lKXN0cmV0Y2h5R0Y5LyUqc3ltbWV0cmljR0Y5LyUobGFyZ2VvcEdGOS8lLm1vdmFibGVsaW1pdHNHRjkvJSdhY2NlbnRHRjkvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZTLUkobWZlbmNlZEdGJDYlLUYjNiUtSSNtbkdGJDYlUSIwRidGN0ZBRjdGQUY3RkFGN0ZBLUY+Ni5RIj1GJ0Y3RkFGQ0ZFRkdGSUZLRk1GTy9GUlEsMC4yNzc3Nzc4ZW1GJy9GVUZdby1JJm1mcmFjR0YkNigtRmZuNiVRIjFGJ0Y3RkEtRmZuNiVRIjJGJ0Y3RkEvJS5saW5ldGhpY2tuZXNzR0Zkby8lK2Rlbm9tYWxpZ25HUSdjZW50ZXJGJy8lKW51bWFsaWduR0ZccC8lKWJldmVsbGVkR0Y5RitGN0ZB

LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYpLUkjbWlHRiQ2JlEieUYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2LlEiJ0YnRjIvRjZRJ25vcm1hbEYnLyUmZmVuY2VHRjQvJSpzZXBhcmF0b3JHRjQvJSlzdHJldGNoeUdGNC8lKnN5bW1ldHJpY0dGNC8lKGxhcmdlb3BHRjQvJS5tb3ZhYmxlbGltaXRzR0Y0LyUnYWNjZW50R0Y0LyUnbHNwYWNlR1EsMC4xMTExMTExZW1GJy8lJ3JzcGFjZUdRJjAuMGVtRictSShtZmVuY2VkR0YkNiUtRiM2JS1JI21uR0YkNiVRIjBGJ0YyRjxGMkY8RjJGPC1GOTYuUSI9RidGMkY8Rj5GQEZCRkRGRkZIRkovRk1RLDAuMjc3Nzc3OGVtRicvRlBGaW4tRlg2JVEiMkYnRjJGPEYyRjw=

compute equi-spaced displacement LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictRiM2Jy1GLDYmUSJ5RicvJSdpdGFsaWNHUSV0cnVlRicvJStleGVjdXRhYmxlR1EmZmFsc2VGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUTAmQXBwbHlGdW5jdGlvbjtGJy9GO1Enbm9ybWFsRicvJSZmZW5jZUdGOS8lKnNlcGFyYXRvckdGOS8lKXN0cmV0Y2h5R0Y5LyUqc3ltbWV0cmljR0Y5LyUobGFyZ2VvcEdGOS8lLm1vdmFibGVsaW1pdHNHRjkvJSdhY2NlbnRHRjkvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZTLUkobWZlbmNlZEdGJDYlLUYjNiUtRiw2JlEidEYnRjRGN0Y6RjdGQUY3RkFGN0ZBRitGN0ZB, then add random noise to the observations.

<Text-field style="Heading 1" layout="Heading 1">Generating the Observations</Text-field>

Starting with the differential equation

q := eval(DE, [a=2, b=4]);

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSJxR0YoLywoLUklZGlmZkdGJjYkLUkieUdGKDYjSSJ0R0YoLUkiJEdGJjYkRjYiIiMiIiIqJi1GMTYkRjNGNkY7LUkkYWJzR0YmNiNGPUY7RjpGMyIiJSIiITcjLywoLUYxNiRGPUY2RjtGPEY6RjNGQkZD

and the initial conditions

ic := y(0)=1/2, D(y)(0)=2;

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSNpY0dGKDYkLy1JInlHRig2IyIiISMiIiIiIiMvLS1JIkRHRiU2I0YxRjJGNjckRi9GNw==

we obtain a numeric solution via

Q := dsolve({q,ic}, y(t), numeric, output=listprocedure):

The function

Yt := rhs(Q[2]):

is a Maple-generated procedure that provides LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictRiM2Jy1GLDYmUSJ5RicvJSdpdGFsaWNHUSV0cnVlRicvJStleGVjdXRhYmxlR1EmZmFsc2VGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUTAmQXBwbHlGdW5jdGlvbjtGJy9GO1Enbm9ybWFsRicvJSZmZW5jZUdGOS8lKnNlcGFyYXRvckdGOS8lKXN0cmV0Y2h5R0Y5LyUqc3ltbWV0cmljR0Y5LyUobGFyZ2VvcEdGOS8lLm1vdmFibGVsaW1pdHNHRjkvJSdhY2NlbnRHRjkvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZTLUkobWZlbmNlZEdGJDYlLUYjNiUtRiw2JlEidEYnRjRGN0Y6RjdGQUY3RkFGN0ZBRitGN0ZB numerically. Choosing

N := 50; d := .2;

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSJOR0YoIiNdNyNGLg==

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSJkR0YoJCIiIyEiIjcjRi4=

as the number of observations and the time interval between observations, we write

Y := [seq(Yt(d*k),k=0..N)]:

to obtain a list of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2JlEiTkYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2LlEnJnBsdXM7RidGMi9GNlEnbm9ybWFsRicvJSZmZW5jZUdGNC8lKnNlcGFyYXRvckdGNC8lKXN0cmV0Y2h5R0Y0LyUqc3ltbWV0cmljR0Y0LyUobGFyZ2VvcEdGNC8lLm1vdmFibGVsaW1pdHNHRjQvJSdhY2NlbnRHRjQvJSdsc3BhY2VHUSwwLjIyMjIyMjJlbUYnLyUncnNwYWNlR0ZOLUkjbW5HRiQ2JVEiMUYnRjJGPEYyRjw= observations, and

P := [seq([d*(k-1),Y[k]], k=1..N+1)]:

to obtain a list of corresponding points.

Figure 1 contains a graph of the numeric solution of the governing initial value problem, with the computed points superimposed.

p1 := odeplot(Q,[t,y(t)],0..10, numpoints=100, color=red): p2 := plot(P, style=point, color=black): display([p1,p2]);

NistJSdDVVJWRVNHNiQ3YHE3JCQiIiEhIiIkIiImISIiNyQkIjAsLCwsLCwsIiEjOiQiMUBCXmUhekxnJyEjOzckJCIwLS0tLS0tLSMhIzokIjFCJSk+P3JIK3chIzs3JCQiMC4uLi4uLi4kISM6JCIwc18tI2ZFaiIpISM6NyQkIjAvLy8vLy8vJSEjOiQiMXgpNGluRzlPKSEjOzckJCIwMDAwMDAwMCYhIzokIjEoW0czS2ojPSMpISM7NyQkIjAxMTExMTExJyEjOiQiMTw1biVRekN3KCEjOzckJCIwMjIyMjIyMighIzokIjF5QGpdMkRocSEjOzckJCIwMzMzMzMzMykhIzokIjFWR3Q/cSV5PichIzs3JCQiMDQ0NDQ0NDQqISM6JCIxeGNeKG9EJVtfISM7NyQkIjAsLCwsLCwsIiEjOSQiMToiKSo9IUcoPkYlISM7NyQkIjA2NjY2NjY2IiEjOSQiMCIqUlpTJzM1TCEjOjckJCIwQEBAQEBAQCIhIzkkIjEweUIlPipvIVIjISM7NyQkIjBKSkpKSkpKIiEjOSQiMjg/cUZ2SUBgIiEjPDckJCIwVFRUVFRUVCIhIzkkIjEkXEoqPkc8a3UhIzw3JCQiMF9eXl5eXl4iISM5JCIxJyoqKW8sTSRbOSUhIz03JCQiMGloaGhoaGgiISM5JCEwdiN5ZlBCdmQhIzs3JCQiMHNycnJycnIiISM5JCEyYVlAOWdrcDUiISM8NyQkIjAjPT09PT09PSEjOSQhMiZcZFF2dlpXOiEjPDckJCIwIz4+Pj4+Pj4hIzkkITIwQiI0JXBKJSkpPSEjPDckJCIwLS0tLS0tLSMhIzkkITIxVEI4Ij50UEAhIzw3JCQiMDc3Nzc3NzcjISM5JCExPkpMJypSbiJIIyEjOzckJCIwQUFBQUFBQSMhIzkkITIoekhpIjNkKVxCISM8NyQkIjBLS0tLS0tLIyEjOSQhMDF5akNJQkojISM6NyQkIjBVVVVVVVVVIyEjOSQhMkEqPXN6PlgjPSMhIzw3JCQiMF9fX19fX18jISM5JCEyZVJGI1x0JSpwPiEjPDckJCIwamlpaWlpaSMhIzkkITJCN3A1ZjIoKm8iISM8NyQkIjB0c3Nzc3NzIyEjOSQhMkclSCM9K3kmZjghIzw3JCQiMCRHR0dHR0dHISM5JCExJSl6RXBhQCMpKiohIzw3JCQiMCRISEhISEhIISM5JCExOVMhPVEqeU1pISM8NyQkIjAuLi4uLi4uJCEjOSQhMk14OEpPI2U4RCEjPTckJCIwODg4ODg4OCQhIzkkIjIzPilRMFhCVzUhIz03JCQiMEJCQkJCQkIkISM5JCIxcHVmZnU2RFYhIzw3JCQiMExMTExMTEwkISM5JCIxPVJuL1psUHMhIzw3JCQiMFZWVlZWVlYkISM5JCIwd0FCJWU1NSgqISM7NyQkIjBhYGBgYGBgJCEjOSQiMmwqKWY2I294bzYhIzw3JCQiMGtqampqamokISM5JCIyaUZvP1lnSUoiISM8NyQkIjB1dHR0dHR0JCEjOSQiMiRbRyp5Vkw2UyIhIzw3JCQiMCVRUVFRUVFRISM5JCIyeUljbTIxOVYiISM8NyQkIjAlUlJSUlJSUiEjOSQiMV5nJFIsXE1TIiEjOzckJCIwLy8vLy8vLyUhIzkkIjJHIikpKipII1wiPjghIzw3JCQiMDk5OTk5OTklISM5JCIyW285TSFIciQ9IiEjPDckJCIwQ0NDQ0NDQyUhIzkkIjJzKUhPJ2VgXisiISM8NyQkIjBNTU1NTU1NJSEjOSQiMSlHYDlIeFMkeiEjPDckJCIwV1dXV1dXVyUhIzkkIjJDPnQ+I3AiUmYmISM9NyQkIjBhYWFhYWFhJSEjOSQiMVhBaGRGZlRKISM8NyQkIjBsa2tra2trJSEjOSQiMU5YLj9WO0RvISM9NyQkIjB2dXV1dXV1JSEjOSQhMHZYdi9xdG8iISM7NyQkIjAmW1tbW1tbWyEjOSQhMGZzbiVHLSUpUSEjOzckJCIwJlxcXFxcXFwhIzkkITIlPU9vJG95aSRlISM9NyQkIjAwMDAwMDAwJiEjOSQhMDwnKVJNKillWyghIzs3JCQiMDo6Ojo6OjomISM5JCEwPlVMVE5xeSkhIzs3JCQiMEREREREREQmISM5JCExJClwKnBJT2NxKiEjPDckJCIwTk5OTk5OTiYhIzkkITIkZTdkR3InPS0iISM8NyQkIjBYWFhYWFhYJiEjOSQhMiYpXHpaeGs4LiIhIzw3JCQiMGNiYmJiYmImISM5JCExSyU9PzM7KSopKiohIzw3JCQiMG1sbGxsbGwmISM5JCExL29dXTxzbCMqISM8NyQkIjB3dnZ2dnZ2JiEjOSQhMEhkSyNHaiI9KSEjOzckJCIwJ2VlZWVlZWUhIzkkITE3YlRzN3UmeichIzw3JCQiMCdmZmZmZmZmISM5JCExTi5sISpmIyp5XiEjPDckJCIwMTExMTExMSchIzkkITBaJT1ETzk0TSEjOzckJCIwOzs7Ozs7OychIzkkITEnKlxgIj1JaWMiISM8NyQkIjBFRUVFRUVFJyEjOSQiMidlOlxjLipRcyMhIz43JCQiME9PT09PT08nISM5JCIxJVFqJ1I3KVsuIyEjPDckJCIwWVlZWVlZWSchIzkkIjE6Y1ZKIT1ybCQhIzw3JCQiMGNjY2NjY2MnISM5JCIxLXRBSCcpKlEzJiEjPDckJCIwbm1tbW1tbSchIzkkIjFKW0FELl9waSEjPDckJCIweHd3d3d3dychIzkkIjEqXCk9bUsjejwoISM8NyQkIjAob29vb29vbyEjOSQiMSpcYVY/L0V5KCEjPDckJCIwKHBwcHBwcHAhIzkkIjE7USU9bFZrMSkhIzw3JCQiMDIyMjIyMjIoISM5JCIxOG8pR04wOy0pISM8NyQkIjA8PDw8PDw8KCEjOSQiMTY/JWZ1NTxsKCEjPDckJCIwRkZGRkZGRighIzkkIjFAeSN5KXBSd3AhIzw3JCQiMFBQUFBQUFAoISM5JCIxJz5vUnV4Mi4nISM8NyQkIjBaWlpaWlpaKCEjOSQiMllBM1oiUVtpWyEjPTckJCIwZWRkZGRkZCghIzkkIjBuZVgtcXpfJCEjOzckJCIwb25ubm5ubighIzkkIjJBdiRRImZGJykzIyEjPTckJCIweXh4eHh4eCghIzkkIjFoeFk8J0c6MichIz03JCQiMCl5eXl5eXl5ISM5JCExd0g/LmtLYiYpISM9NyQkIjApenp6enp6eiEjOSQhMTE5Wm5jblVBISM8NyQkIjAzMzMzMzMzKSEjOSQhMllmZyIzOlMuTiEjPTckJCIwPT09PT09PSkhIzkkITFhbl81QiRRZiUhIzw3JCQiMEdHR0dHR0cpISM5JCExXGEvPVB1eGEhIzw3JCQiMFFRUVFRUVEpISM5JCEydkhlWywqKm83JyEjPTckJCIwW1tbW1tbWykhIzkkITBXMmsmPj1AbCEjOzckJCIwZWVlZWVlZSkhIzkkITFMaUxKWWlbbSEjPDckJCIwcG9vb29vbykhIzkkITFtIXk8NXNjXSchIzw3JCQiMHp5eXl5eXkpISM5JCExaVpmP1omKSo0JyEjPDckJCIwKikpKSkpKSkpKSkpKSkpISM5JCExSip5aWI3PFgmISM8NyQkIjAqKikqKSopKikqKSopKikhIzkkITEnSCRbdyR6TGYlISM8NyQkIjA0NDQ0NDQ0KiEjOSQhMkRSJ2ZfJ3BoYyQhIz03JCQiMD4+Pj4+Pj4qISM5JCExYk52U3RzPEMhIzw3JCQiMEhISEhISEgqISM5JCExb3o+NyIpMyo+IiEjPDckJCIwUlJSUlJSUiohIzkkIjFhS3lkIykpZiFRISM+NyQkIjBcXFxcXFxcKiEjOSQiMl8nPTB6VHNWNyEjPTckJCIvJ2ZmZmZmZiohIzgkIjJEclkzank5UCMhIz03JCQiLyhwcHBwcHAqISM4JCIydiUpUnNDYSl6TCEjPTckJCIvKXp6enp6eiohIzgkIjJFeik9MmRJTFUhIz03JCQiLyoqKSopKikqKSopKikqISM4JCIxSmI4USFRRyFcISM8NyQkIiQrIiEiIiQiMSdlOjB1IVFtYCEjPC0lJkNPTE9SRzYmJSRSR0JHJCIjNSEiIiQiIiEhIiIkIiIhISIiLSUnQ1VSVkVTRzYlN1U3JCQiIiEhIiIkIiImISIiNyQkIiIjISIiJCIvJioqKSlcYlxlKCEjOTckJCIiJSEiIiQiMSlbUigpKSlHLE8pISM7NyQkIiInISIiJCIxZik+aUNgdnooISM7NyQkIiIpISIiJCIxKDQiNCYpWyo0RichIzs3JCQiIzUhIiIkIjBKPi5VRiZwViEjOjckJCIjNyEiIiQiMm0xcjQqbzQpXCMhIzw3JCQiIzkhIiIkIjFlNEJQYi48JikhIzw3JCQiIzshIiIkITF2SW1yV0dXWyEjPDckJCIjPSEiIiQhMmw2Kik+dSFmczkhIzw3JCQiIz8hIiIkITBNIlFpRlsmNCMhIzo3JCQiI0EhIiIkITJuX19eMCVHWEIhIzw3JCQiI0MhIiIkITE6dWUmKWVpQEEhIzs3JCQiI0UhIiIkITI5OVEoeS03bzwhIzw3JCQiI0chIiIkITFMRUluLFwsNiEjOzckJCIjSSEiIiQhMT4nZlslb2M9TyEjPDckJCIjSyEiIiQiMTBzIlEpZU42TCEjPDckJCIjTSEiIiQiMTohKXAxRmRBKikhIzw3JCQiI08hIiIkIjEwIjQ0JW9abjchIzs3JCQiI1EhIiIkIjE3ajU5endFOSEjOzckJCIjUyEiIiQiMlUzUillJG8kZjghIzw3JCQiI1UhIiIkIjJNI1wqRy8neiUzIiEjPDckJCIjVyEiIiQiMXImeTI2KClRaychIzw3JCQiI1khIiIkIjIpeTlyTGY6Mz0hIz03JCQiI1shIiIkITJWZjBQKlEwY0chIz03JCQiI10hIiIkITFQUzY5OzEtbiEjPDckJCIjXyEiIiQhMSVvKzExdXRGKiEjPDckJCIjYSEiIiQhMmFiYkhmbDkuIiEjPDckJCIjYyEiIiQhMWRYIypvVmg9KCohIzw3JCQiI2UhIiIkITFmXCo+KillQ2ooISM8NyQkIiNnISIiJCEyYTksLT0hKltbJSEjPTckJCIjaSEiIiQhMS1AaC4qeitqKSEjPTckJCIjayEiIiQiMSIqUilROUQiUUUhIzw3JCQiI20hIiIkIjJ2VVRQJlwyO2IhIz03JCQiI28hIiIkIjF3PV1IQlowdSEjPDckJCIjcSEiIiQiMSxxSSVlNnczKSEjPDckJCIjcyEiIiQiMTEhekI+a0BcKCEjPDckJCIjdSEiIiQiMDkmZjlnR1lkISM7NyQkIiN3ISIiJCIxVz50OjRzKj0kISM8NyQkIiN5ISIiJCIyJ3lpWSgqZiFwIkchIz43JCQiIyEpISIiJCEyRW9WNydcVzFEISM9NyQkIiMjKSEiIiQhMV42Um9UNXBaISM8NyQkIiMlKSEiIiQhMTdJTjMtXzJpISM8NyQkIiMnKSEiIiQhMSczKSk+W2JbaychIzw3JCQiIykpISIiJCEyJm9yb0whUVYuJyEjPTckJCIjISohIiIkITFWRCw5N14oXCUhIzw3JCQiIyMqISIiJCEyeEM1cDBLQEsjISM9NyQkIiMlKiEiIiQiMnRnQSJSYHE8NiEjPjckJCIjJyohIiIkIjJaJylwSmNpVlQjISM9NyQkIiMpKiEiIiQiMk4tRlYoUWVbVSEjPTckJCIkKyIhIiIkIjEnZTowdSFRbWAhIzwtJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lJlNUWUxFRzYjJSZQT0lOVEctJSVWSUVXRzYkOyQiIiEhIiIkIiQrIiEiIjskITIoekhpIjNkKVxCISM8JCIxeCk0aW5HOU8pISM7LSYlJl9BWElTRzYjIiIiNiYtJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIS0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJiUmX0FYSVNHNiMiIiM2Ji0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqTElORVNUWUxFRzYjIiIhLSUqVEhJQ0tORVNTRzYjIiIhLSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0lK0FYRVNMQUJFTFNHNiQtSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1USJ0NiIvJSdmYW1pbHlHUSE2Ii8lJXNpemVHUSMxMDYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiItSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1USJ5NiIvJSdmYW1pbHlHUSE2Ii8lJXNpemVHUSMxMDYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiItJSlfVklTSUJMRUc2IyIiIi0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiIiEhIiItJSlCT1VORFNfWUc2IyQiIiEhIiItJS1CT1VORFNfV0lEVEhHNiMkIiUrXSEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiUrXSEiIi0lKUNISUxEUkVORzYiLSUrQU5OT1RBVElPTkc2Jy0lKUJPVU5EU19YRzYjJCIiISEiIi0lKUJPVU5EU19ZRzYjJCIiISEiIi0lLUJPVU5EU19XSURUSEc2IyQiJStdISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJStdISIiLSUpQ0hJTERSRU5HNiI=NiI=

Figure 1 Numeric solution of governing initial value problem

<Text-field style="Heading 1" layout="Heading 1">Adding Random Noise</Text-field>

We distinguish between two types of random noise appearing in experimental data. True experimental noise is relative to the data-collection process itself. External interference is additive.

A sample of size LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2JlEiTkYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2LlEnJnBsdXM7RidGMi9GNlEnbm9ybWFsRicvJSZmZW5jZUdGNC8lKnNlcGFyYXRvckdGNC8lKXN0cmV0Y2h5R0Y0LyUqc3ltbWV0cmljR0Y0LyUobGFyZ2VvcEdGNC8lLm1vdmFibGVsaW1pdHNHRjQvJSdhY2NlbnRHRjQvJSdsc3BhY2VHUSwwLjIyMjIyMjJlbUYnLyUncnNwYWNlR0ZOLUkjbW5HRiQ2JVEiMUYnRjJGPEYyRjw= taken from a normal distribution with mean 1 and standard deviation LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYmLUkjbWlHRiQ2I1EhRictRiw2JlEoJnNpZ21hO0YnLyUnaXRhbGljR1EmZmFsc2VGJy8lK2V4ZWN1dGFibGVHRjQvJSxtYXRodmFyaWFudEdRJ25vcm1hbEYnRjVGNw== provides factors normally distributed about 1. Multiplying each computed value of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictRiM2Jy1GLDYmUSJ5RicvJSdpdGFsaWNHUSV0cnVlRicvJStleGVjdXRhYmxlR1EmZmFsc2VGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUTAmQXBwbHlGdW5jdGlvbjtGJy9GO1Enbm9ybWFsRicvJSZmZW5jZUdGOS8lKnNlcGFyYXRvckdGOS8lKXN0cmV0Y2h5R0Y5LyUqc3ltbWV0cmljR0Y5LyUobGFyZ2VvcEdGOS8lLm1vdmFibGVsaW1pdHNHRjkvJSdhY2NlbnRHRjkvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZTLUkobWZlbmNlZEdGJDYlLUYjNiUtRiw2JlEidEYnRjRGN0Y6RjdGQUY3RkFGN0ZBRitGN0ZB by one such factor perturbes the computed values with noise of magnitude relative to that of the computed value.

Choosing LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2I1EhRictRiw2JlEoJnNpZ21hO0YnLyUnaXRhbGljR1EmZmFsc2VGJy8lK2V4ZWN1dGFibGVHRjQvJSxtYXRodmFyaWFudEdRJ25vcm1hbEYnLUkjbW9HRiQ2LlEiPUYnRjVGNy8lJmZlbmNlR0Y0LyUqc2VwYXJhdG9yR0Y0LyUpc3RyZXRjaHlHRjQvJSpzeW1tZXRyaWNHRjQvJShsYXJnZW9wR0Y0LyUubW92YWJsZWxpbWl0c0dGNC8lJ2FjY2VudEdGNC8lJ2xzcGFjZUdRLDAuMjc3Nzc3OGVtRicvJSdyc3BhY2VHRk4tSSNtbkdGJDYlUSQwLjNGJ0Y1RjdGNUY3, we have

S := Sample(Normal(1,.3),N+1): Yn := zip(`*`, S, Vector(Y))^%T:

for our noisy data. Figure 2 shows the noisy observations as red dots, and the original computed displacements as black dots.

p3 := plot([seq([d*(k-1),Yn[k]],k=1..N+1)],style=point, color=red): display([p2,p3]);

NiotJSdDVVJWRVNHNiU3VTckJCIiISEiIiQiIiYhIiI3JCQiIiMhIiIkIi8mKiopKVxiXGUoISM5NyQkIiIlISIiJCIxKVtSKCkpKUcsTykhIzs3JCQiIichIiIkIjFmKT5pQ2B2eighIzs3JCQiIikhIiIkIjEoNCI0JilbKjRGJyEjOzckJCIjNSEiIiQiMEo+LlVGJnBWISM6NyQkIiM3ISIiJCIybTFyNCpvNClcIyEjPDckJCIjOSEiIiQiMWU0QlBiLjwmKSEjPDckJCIjOyEiIiQhMXZJbXJXR1dbISM8NyQkIiM9ISIiJCEybDYqKT51IWZzOSEjPDckJCIjPyEiIiQhME0iUWlGWyY0IyEjOjckJCIjQSEiIiQhMm5fX14wJUdYQiEjPDckJCIjQyEiIiQhMTp1ZSYpZWlAQSEjOzckJCIjRSEiIiQhMjk5USh5LTdvPCEjPDckJCIjRyEiIiQhMUxFSW4sXCw2ISM7NyQkIiNJISIiJCExPidmWyVvYz1PISM8NyQkIiNLISIiJCIxMHMiUSllTjZMISM8NyQkIiNNISIiJCIxOiEpcDFGZEEqKSEjPDckJCIjTyEiIiQiMTAiNDQlb1puNyEjOzckJCIjUSEiIiQiMTdqNTl6d0U5ISM7NyQkIiNTISIiJCIyVTNSKWUkbyRmOCEjPDckJCIjVSEiIiQiMk0jXCpHLyd6JTMiISM8NyQkIiNXISIiJCIxciZ5MjYoKVFrJyEjPDckJCIjWSEiIiQiMil5OXJMZjozPSEjPTckJCIjWyEiIiQhMlZmMFAqUTBjRyEjPTckJCIjXSEiIiQhMVBTNjk7MS1uISM8NyQkIiNfISIiJCExJW8rMTF1dEYqISM8NyQkIiNhISIiJCEyYWJiSGZsOS4iISM8NyQkIiNjISIiJCExZFgjKm9WaD0oKiEjPDckJCIjZSEiIiQhMWZcKj4qKWVDaighIzw3JCQiI2chIiIkITJhOSwtPSEqW1slISM9NyQkIiNpISIiJCExLUBoLip6K2opISM9NyQkIiNrISIiJCIxIipSKVE5RCJRRSEjPDckJCIjbSEiIiQiMnZVVFAmXDI7YiEjPTckJCIjbyEiIiQiMXc9XUhCWjB1ISM8NyQkIiNxISIiJCIxLHFJJWU2dzMpISM8NyQkIiNzISIiJCIxMSF6Qj5rQFwoISM8NyQkIiN1ISIiJCIwOSZmOWdHWWQhIzs3JCQiI3chIiIkIjFXPnQ6NHMqPSQhIzw3JCQiI3khIiIkIjIneWlZKCpmIXAiRyEjPjckJCIjISkhIiIkITJFb1Y3J1xXMUQhIz03JCQiIyMpISIiJCExXjZSb1Q1cFohIzw3JCQiIyUpISIiJCExN0lOMy1fMmkhIzw3JCQiIycpISIiJCExJzMpKT5bYltrJyEjPDckJCIjKSkhIiIkITImb3JvTCFRVi4nISM9NyQkIiMhKiEiIiQhMVZELDk3XihcJSEjPDckJCIjIyohIiIkITJ4QzVwMEtASyMhIz03JCQiIyUqISIiJCIydGdBIlJgcTw2ISM+NyQkIiMnKiEiIiQiMlonKXBKY2lWVCMhIz03JCQiIykqISIiJCIyTi1GVihRZVtVISM9NyQkIiQrIiEiIiQiMSdlOjB1IVFtYCEjPC0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUmU1RZTEVHNiMlJlBPSU5URy0lJ0NVUlZFU0c2JTdVNyQkIiIhISIiJCIybGYtITNRTyJSJCEjPDckJCIiIyEiIiQiMXJ5dyg+VmgkbyEjOzckJCIiJSEiIiQiMW0iNCxYUUMibyEjOzckJCIiJyEiIiQiMVl6KmUrWyMqSCkhIzs3JCQiIikhIiIkIjFFUCNvKXA6QmkhIzs3JCQiIzUhIiIkIjE5TDhBcXdObSEjOzckJCIjNyEiIiQiMicqPjMlPXgpR0YiISM8NyQkIiM5ISIiJCIyc0ciUUlhOmA3ISM8NyQkIiM7ISIiJCExPCgqSG9RJz00JiEjPDckJCIjPSEiIiQhMXM1LHkqSHMiPiEjOzckJCIjPyEiIiQhMS9VYSlmMldFIyEjOzckJCIjQSEiIiQhMiYpPUZcKT0lUUoiISM8NyQkIiNDISIiJCEyPWJ3TXIkeig+IyEjPDckJCIjRSEiIiQhMlBsIzRnMW5rPCEjPDckJCIjRyEiIiQhMmhnJCpHQCVmIzMiISM8NyQkIiNJISIiJCExLyRwLV5xIj1NISM8NyQkIiNLISIiJCIxZFNrPnMnUV4kISM8NyQkIiNNISIiJCIxZC05JlI6emwpISM8NyQkIiNPISIiJCIyX3JzU1lKZ0IiISM8NyQkIiNRISIiJCIxbTR5P3khW04iISM7NyQkIiNTISIiJCIyJj0nKVFaQldEOSEjPDckJCIjVSEiIiQiMVVaQi8nZSNbJykhIzw3JCQiI1chIiIkIjAqPTB2VElsZCEjOzckJCIjWSEiIiQiMTxHaiozKj4lSCMhIzw3JCQiI1shIiIkITEmPktTImUxXkshIzw3JCQiI10hIiIkITEjKUgleSplZXIjKiEjPDckJCIjXyEiIiQhMV9sVjY0QUkqKiEjPDckJCIjYSEiIiQhMmE4YWcjM0M/NSEjPDckJCIjYyEiIiQhMWxiOEEpZio0JSkhIzw3JCQiI2UhIiIkITEtalFtMltmXCEjPDckJCIjZyEiIiQhMCc0OGMiKT1ZRyEjOzckJCIjaSEiIiQhMFtjdnRhMGsmISM8NyQkIiNrISIiJCIwcylHXE8iXHMjISM7NyQkIiNtISIiJCIxPm1hYC9TIVwlISM8NyQkIiNvISIiJCIxPSl6JHk8NUspKSEjPDckJCIjcSEiIiQiMXNabUZIRnR3ISM8NyQkIiNzISIiJCIxL2dhIlF1KltjISM8NyQkIiN1ISIiJCIxQlsrW1RAbiEpISM8NyQkIiN3ISIiJCIyJHBiLHFYO1lBISM9NyQkIiN5ISIiJCIyT29uLyhHcSs/ISM+NyQkIiMhKSEiIiQhMjBcRU09UC5bIiEjPTckJCIjIykhIiIkITFDXkdlTCMqSHIhIzw3JCQiIyUpISIiJCEwZjolKm89enIpISM7NyQkIiMnKSEiIiQhMS1iajpxIjNnJSEjPDckJCIjKSkhIiIkITFpWUglb0IuWCchIzw3JCQiIyEqISIiJCExJCoqbzc+azVdJyEjPDckJCIjIyohIiIkITJfSnVgRCZHOzkhIz03JCQiIyUqISIiJCIxMVN1P0AlZSIqKiEjPjckJCIjJyohIiIkIjE4JilmKVI4ImZDISM8NyQkIiMpKiEiIiQiMSdlLWVudjFGJiEjPDckJCIkKyIhIiIkIjE3WCU+UUpyOiYhIzwtJSZDT0xPUkc2JiUkUkdCRyQiIzUhIiIkIiIhISIiJCIiISEiIi0lJlNUWUxFRzYjJSZQT0lOVEctJiUmX0FYSVNHNiMiIiI2Ji0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqTElORVNUWUxFRzYjIiIhLSUqVEhJQ0tORVNTRzYjIiIhLSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0mJSZfQVhJU0c2IyIiIzYmLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUrQVhFU0xBQkVMU0c2JFEhNiJRITYiLSUpX1ZJU0lCTEVHNiMiIiItJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiIhISIiLSUpQk9VTkRTX1lHNiMkIiIhISIiLSUtQk9VTkRTX1dJRFRIRzYjJCIlK10hIiItJS5CT1VORFNfSEVJR0hURzYjJCIlK10hIiItJSlDSElMRFJFTkc2Ii0lK0FOTk9UQVRJT05HNictJSlCT1VORFNfWEc2IyQiIiEhIiItJSlCT1VORFNfWUc2IyQiIiEhIiItJS1CT1VORFNfV0lEVEhHNiMkIiUrXSEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiUrXSEiIi0lKUNISUxEUkVORzYiRzYi

Figure 2 Computed displacements (black) and noisy observations (red)

<Text-field style="Heading 1" layout="Heading 1">The Sum of Squares of Deviations</Text-field>

The following Maple procedure is a function of the two parameters LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JlEiYUYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnRjIvRjZRJ25vcm1hbEYn and LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JlEiYkYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnRjIvRjZRJ25vcm1hbEYn appearing in the differential equation of our nonlinear model. For each pairLUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictSShtZmVuY2VkR0YkNiUtRiM2Jy1GLDYmUSJhRicvJSdpdGFsaWNHUSV0cnVlRicvJStleGVjdXRhYmxlR1EmZmFsc2VGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYuUSIsRidGOi9GPlEnbm9ybWFsRicvJSZmZW5jZUdGPC8lKnNlcGFyYXRvckdGOS8lKXN0cmV0Y2h5R0Y8LyUqc3ltbWV0cmljR0Y8LyUobGFyZ2VvcEdGPC8lLm1vdmFibGVsaW1pdHNHRjwvJSdhY2NlbnRHRjwvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR1EsMC4zMzMzMzMzZW1GJy1GLDYmUSJiRidGN0Y6Rj1GOkZERjpGREYrRjpGRA==, the differential equation is integrated numerically so that LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2JlEiTkYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2LlEnJnBsdXM7RidGMi9GNlEnbm9ybWFsRicvJSZmZW5jZUdGNC8lKnNlcGFyYXRvckdGNC8lKXN0cmV0Y2h5R0Y0LyUqc3ltbWV0cmljR0Y0LyUobGFyZ2VvcEdGNC8lLm1vdmFibGVsaW1pdHNHRjQvJSdhY2NlbnRHRjQvJSdsc3BhY2VHUSwwLjIyMjIyMjJlbUYnLyUncnNwYWNlR0ZOLUkjbW5HRiQ2JVEiMUYnRjJGPEYyRjw= equispaced nodes are computed. These nodes are compared to the noisy observations via the Euclidean norm of the difference between the vector of observations and the vector of computed solution values.

SS := proc(a,b) local F, V; if not type([a,b],[numeric,numeric]) then return 'SS'(a,b); elif a<0 or b<0 then return 1e100; end if; F := dsolve(eval({DE,ic},{:-a=a,:-b=b}), y(t), numeric, output=Array([seq(d*k,k=0..N)])); V := convert(Column(F[2,1],2),Vector); Norm(V-Yn,2); end proc:

(The penalty-function construct that keeps the parameters in the first quadrant was suggested by Dr. Allan Wittkopf of Maplesoft.)

Figure 3 contains a graph of the surface LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictRiM2Jy1GLDYmUSNTU0YnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2LVEwJkFwcGx5RnVuY3Rpb247RicvRjtRJ25vcm1hbEYnLyUmZmVuY2VHRjkvJSpzZXBhcmF0b3JHRjkvJSlzdHJldGNoeUdGOS8lKnN5bW1ldHJpY0dGOS8lKGxhcmdlb3BHRjkvJS5tb3ZhYmxlbGltaXRzR0Y5LyUnYWNjZW50R0Y5LyUnbHNwYWNlR1EmMC4wZW1GJy8lJ3JzcGFjZUdGUy1JKG1mZW5jZWRHRiQ2JS1GIzYnLUYsNiZRImFGJ0Y0RjdGOi1GPjYuUSIsRidGN0ZBRkMvRkZGNkZHRklGS0ZNRk9GUS9GVVEsMC4zMzMzMzMzZW1GJy1GLDYmUSJiRidGNEY3RjpGN0ZBRjdGQUY3RkFGK0Y3RkE=, which indicates that there is indeed at least one minimum for this function.

plot3d(SS,0..5, 2..6, axes=box, view=.3..3/2, labels=[a,b,""], orientation=[-65,70]);

LSUnUExPVDNERzYuLSUlR1JJREc2JTskIiIhISIiJCIjXSEiIjskIiM/ISIiJCIjZyEiIjdTN1MkIjEqKioqKlI4WDVQKCEjOiQiKzksPF9zISIqJCIrNGhDJDQoISIqJCIrTipwZSFwISIqJCIrM1ZYNG4hIiokIitPW09DbCEiKiQiKy08J2VPJyEiKiQiKkdhLUMnISIpJCIrTitAV2ghIiokIit2dDNuZyEiKiQiMiYqKioqKmYwP10qZiEjOyQiKykqUjM6ZiEiKiQiKyNHWCQ9ZSEiKiQiK1lOUCxkISIqJCIrRUVPbWIhIiokIisnKj0/P2EhIiokIismZU9HRiYhIiokIitYKDRfOCYhIiokIit6ISo0PF0hIiokIisnek1fI1whIiokIitEeDdpWyEiKiQiKyc9OGYjWyEiKiQiK0olMzciWyEiKiQiK0xDazVbISIqJCIramthO1shIiokIipZcEIjWyEiKSQiKypwKWVCWyEiKiQiK2hCMz1bISIqJCIqUVZnIVshIikkIisoel0mKnklISIqJCIqVDs+eCUhIikkIiojKj1wdiUhIikkIitOJlEhW1ohIiokIioqZiF5dSUhIikkIisjSG10diUhIiokIitlK0t3WiEiKiQiKzx5JUchWyEiKiQiKzEsMU1bISIqJCIrbSZ5aydbISIqJCIrPUxhJypbISIqJCIrIltpNSNcISIqJCIrakFiUFwhIiokIitMdFdXXCEiKiQiK3BQPVRcISIqJCIrPCRwIkdcISIqJCIrXzZtMVwhIiokIitDRWR5WyEiKiQiK1A6QllbISIqJCIrV002N1shIio3UyQiK3Q4JFE6JiEiKiQiK20kNCZcXSEiKiQiK1IxJDMkXCEiKiQiK2JCQyxbISIqJCIrOCZ5d20lISIqJCIrLiRHeGAlISIqJCIrWCllc1QlISIqJCIrLU4qKjNWISIqJCIrV3hBN1UhIiokIistXmBCVCEiKiQiKmghNFFTISIpJCIqUk82JlIhIikkIitCI1sjZlEhIiokIisuJFI2dyQhIiokIislKlsoeWwkISIqJCIrVmJiX04hIiokIissOmtcTSEiKiQiKjltU04kISIpJCIqJmZKcUshIikkIituT2ssSyEiKiQiKipbY1xKISIpJCIrMW50OEohIiokIitcJ0dBNCQhIiokIitOTzEjMyQhIiokIioyIyopekkhIikkIit5W2cjMyQhIiokIitpUHooMyQhIiokIitLdilSNCQhIiokIitRKWUxNSQhIiokIiskb2MhM0ohIiokIitIOiRwNiQhIiokIis3cUBHSiEiKiQiK1hCdFVKISIqJCIrWCZbNDskISIqJCIrZTsmRz0kISIqJCIrMC4iej8kISIqJCIrJlxNXkIkISIqJCIrJT05S0UkISIqJCIrWSZ6MUgkISIqJCIrI3puZ0okISIqJCIrS28zUUwhIiokIitsWnRiTCEiKiQiKzFDUW9MISIqJCIrcyIzZVAkISIqJCIrekE9eUwhIiokIis5dSx3TCEiKiQiK1hpM3FMISIqJCIrQ1ZLaEwhIiokIiomcHNdTCEiKTdTJCIrI0hOYzglISIqJCIrKmZNKVJTISIqJCIrVGNwT1IhIiokIitYJG92I1EhIiokIipRN2dyJCEiKSQiK1xBOzFPISIqJCIrWEVOLE4hIiokIit6R0cuTSEiKiQiKyRlNj1KJCEiKiQiK3d6RkRLISIqJCIrWyR6NjkkISIqJCIrN2sncDAkISIqJCIrKipveXFIISIqJCIreGIsIylHISIqJCIqKj5XInojISIpJCIqNkI2cSMhIikkIitwYys5RSEiKiQiK09uUUxEISIqJCIrRGdNaUMhIiokIitudEIuQyEiKiQiK2YjUXROIyEiKiQiKEFaSyMhIickIis7SVAvQiEiKiQiK3lqWyVIIyEiKiQiK1gwKUdIIyEiKiQiKyUqb1IoSCMhIiokIitjIz1pSSMhIiokIisqR2EhPUIhIiokIitzRzxLQiEiKiQiKyhcOCRbQiEiKiQiK3QkR2xPIyEiKiQiKzRzKnBRIyEiKiQiK2BqJik0QyEiKiQiKzIrMk5DISIqJCIrPidlQlkjISIqJCIrJ3AoPSJcIyEiKiQiKnotM18jISIpJCIraDJKXUQhIiokIituW3d5RCEiKiQiK24qXF9nIyEiKiQiK3haKCpHRSEiKiQiK3ZVTFxFISIqJCIrQE8mZm0jISIqJCIreG5xeUUhIiokIit5V3EobyMhIiokIitUXkUkcCMhIiokIithUSdlcCMhIiokIislUXdncCMhIiokIis4JT1YcCMhIio3UyQiK0EtZEFOISIqJCIrJipwcktNISIqJCIrblxkUUwhIiokIitIYWlTSyEiKiQiK0UkbzQ5JCEiKiQiKmVFQS8kISIpJCIqMDhtJUghIikkIitfdE1iRyEiKiQiKzpVW29GISIqJCIrPi4zJm8jISIqJCIrQ0BoLkUhIiokIitBZFtBRCEiKiQiK05WYVNDISIqJCIrUS1ZZEIhIiokIit3TiFSRiMhIiokIisselkiPiMhIiokIis4LVY3QCEiKiQiK0ImeSRSPyEiKiQiK04pKnp1PiEiKiQiK1dMcD8+ISIqJCIrYCV6I3k9ISIqJCIrUXgmeSU9ISIqJCIrTiJbKUc9ISIqJCIraCp6Kj49ISIqJCIrRCUqZT49ISIqJCIrTCpHZiM9ISIqJCIrLFtVUD0hIiokIistRiVHJj0hIiokIit6WExyPSEiKiQiKyI9I1IjKj0hIiokIitsU3U6PiEiKiQiK20xQlQ+ISIqJCIrMykqb28+ISIqJCIrNjsoeSo+ISIqJCIrPTpSRz8hIiokIitAZHNmPyEiKiQiKzhOQiI0IyEiKiQiKzI2QEFAISIqJCIrVWklPjojISIqJCIraCJwKHpAISIqJCIrJz02Xj8jISIqJCIoWHZBIyEiJyQiK3UxIm9DIyEiKiQiKztvI0dFIyEiKiQiK2sub3ZBISIqJCIqaTBjRyMhIikkIisiKilcSEgjISIqJCIrLEw4KUgjISIqJCIrYXVnLEIhIio3UyQiK2VwbS9KISIqJCIrKG8qUj5JISIqJCIrVExbSkghIiokIituRSY0JUchIiokIisuIT4iXEYhIiokIitUMnlkRSEiKiQiK2p0X29EISIqJCIrI1wmSCNbIyEiKiQiK24mPSMqUiMhIiokIitfbHI9QiEiKiQiKiR6eFJBISIpJCIrcyc9ODsjISIqJCIrWHpiIzMjISIqJCIrRnlLLj8hIiokIismb0ZVIz4hIiokIitRP2ZZPSEiKiQiK0gtTHM8ISIqJCIrKnBdT3EiISIqJCIrWChRRmsiISIqJCIoUjlmIiEiJyQiKyMpKSopNGIiISIqJCIrUmsnPV8iISIqJCIrSyJwUF0iISIqJCIra2R1JlwiISIqJCIrJGZAa1wiISIqJCIrdEZFLzohIiokIitsZCF5XiIhIiokIitBeSJlYCIhIiokIisjR2R0YiIhIiokIit4WXYiZSIhIiokIitPPGEzOyEiKiQiK2g5T1A7ISIqJCIrJylRKXltIiEiKiQiK09KdSpwIiEiKiQiKyIpR11LPCEiKiQiKyRHX2N3IiEiKiQiKyoqKj4nKXoiISIqJCIrKGU4MyQ9ISIqJCIrc1dsaD0hIiokIiswRGghKj0hIiokIit4YkM8PiEiKiQiK09xQVQ+ISIqJCIrbUZPaT4hIiokIitkemYhKT4hIiokIitkJDNnKj4hIiokIisjZiR5Mz8hIiokIit1Oj8+PyEiKiQiK0Z4ZkY/ISIqJCIrYFZMTT8hIio3UyQiK2ViXil6IyEiKiQiK1VvLDxGISIqJCIrI1tvUWojISIqJCIrQkIlKVtEISIqJCIrdCpvRlkjISIqJCIrJjRMcFAjISIqJCIrOSpHREgjISIqJCIrLT1JNUEhIiokIitqWVRJQCEiKiQiKmYlXF8/ISIpJCIrQzQkZSg+ISIqJCIrO3BrKio9ISIqJCIrc3BRQj0hIiokIipPd3B1IiEiKSQiK0FDJzRuIiEiKiQiKms2bGYiISIpJCIrKUcrYF8iISIqJCIraFpJZjkhIiokIitEOWErOSEiKiQiK0NGel04ISIqJCIqRWM4SiIhIikkIitpLihHRyIhIiokIitGbERsNyEiKiQiKyg9InpkNyEiKiQiKy53R2Y3ISIqJCIqJUhMbzchIikkIis9cF0kRyIhIiokIistaWEuOCEiKiQiKzMtVUY4ISIqJCIqb0tWTiIhIikkIitLUm4kUSIhIiokIitVQSZcVCIhIiokIitQJUd4VyIhIiokIis8WWMiWyIhIiokIisiPicqZl4iISIqJCIrOmhfXTohIiokIipyTFllIiEiKSQiKyY9K3loIiEiKiQiKW5gXDshIigkIit5XVN6OyEiKiQiKypRWHFxIiEiKiQiKipvPks8ISIpJCIrdCkzWnYiISIqJCIrM1dhdTwhIiokIis5UHgieiIhIiokIis3ImZsIT0hIiokIit1Zzg+PSEiKiQiK3p1eUg9ISIqJCIrZiwjKVE9ISIqN1MkIisqPSlSakQhIiokIisxLzgmWyMhIiokIis0ZyhlUyMhIiokIit6IlJfSyMhIiokIitMNnRWQSEiKiQiK3BzR2lAISIqJCIrMiNIPTMjISIqJCIrajwnSCsjISIqJCIrVFImZSM+ISIqJCIrJ2VkLSY9ISIqJCIrXlVsdjwhIiokIisnNGo5cSIhIiokIitfWEZGOyEiKiQiK256MWA6ISIqJCIrbG9NejkhIiokIitFVToyOSEiKiQiK3oyK1E4ISIqJCIrLGBxdDchIiokIismPiI9OzchIiokIitOaz5uNiEiKiQiK19GOEc2ISIqJCIpQ3oqNCIhIigkIitDdUgjMyIhIiokIitzVDZ2NSEiKiQiKnolPng1ISIpJCIqSig+KDMiISIpJCIreVtxLjYhIiokIitGNlNENiEiKiQiK2tRO142ISIqJCIrXGQzIT0iISIqJCIrWUpXNjchIiokIilKa1c3ISIoJCIrdCJvIno3ISIqJCIrbEpgOTghIiokIitGdkRdOCEiKiQiKmlmZVEiISIpJCIrdSVlM1UiISIqJCIrYTp6YTkhIiokIiteWEIoWyIhIiokIisxaCJ5XiIhIiokIis8YkJZOiEiKiQiKzImeUFkIiEiKiQiKyYqZSNlZiIhIiokIitDVSZvaCIhIiokIis8MVZOOyEiKiQiKiIpKnBeOyEiKSQiK1Aib2VtIiEiKiQiK2ZJPXk7ISIqJCIrI2Y4KilvIiEiKjdTJCIrRWRzd0IhIiokIitYQUosQiEiKiQiKydbeWBBIyEiKiQiKk49JVtAISIpJCIqXDsyMiMhIikkIitXJWVIKj4hIiokIis7VSdlIj4hIiokIitqKyQqUj0hIiokIisvekpsPCEiKiQiKyZ5ZT1wIiEiKiQiK1k3Oz47ISIqJCIrUXR4WTohIiokIit5UlJ1OSEiKiQiK2E4LS05ISIqJCIrXkM3STghIiokIitQKVInZjchIiokIitWVyY+PiIhIiokIitrW3dHNiEiKiQiKyRvOT4yIiEiKiQiK1kzPUI1ISIqJCIrOUxXUykqISM1JCIrUmF2YSYqISM1JCIrW1AmelAqISM1JCIrVy9ZMSQqISM1JCIqbWU0TCohIiokIis7P1hRJSohIzUkIitHJ1FYaCohIzUkIit3JVxgJSkqISM1JCIrRilwPSwiISIqJCIrI3BZQy8iISIqJCIrVFZXdjUhIiokIisiZSk9NTYhIiokIislRy5oOSIhIiokIitBMm4jPSIhIiokIituVFM+NyEiKiQiKzR3JGVEIiEiKiQiKzlSXyJIIiEiKiQiK2hkL0U4ISIqJCIrcysuZjghIiokIispPmMsUiIhIiokIitwJ28iPjkhIiokIitIIyopZVciISIqJCIqKD1BcTkhIikkIit3UzojXCIhIiokIissMnY2OiEiKiQiK1BVOUg6ISIqJCIrai5fVzohIiokIitpMjVlOiEiKiQiK3duN3E6ISIqN1MkIitAZCJbQSMhIiokIitcbyg+OiMhIiokIitgRycqeT8hIiokIitmPT4wPyEiKiQiK1E7eUk+ISIqJCIrdSpIaSY9ISIqJCIrTkE1I3kiISIqJCIrNG4hKTM8ISIqJCIrIj0jXE87ISIqJCIrQSwvbDohIiokIis3ZzklXCIhIiokIipoYU1VIiEiKSQiK3UqPUZOIiEiKiQiKkNsPkciISIpJCIrIylcaDY3ISIqJCIrVVdfVTYhIiokIitQUCdmMiIhIiokIitWMWA4NSEiKiQiK0FbNnEmKiEjNSQiKydcMz8zKiEjNSQiKyU+RW9vKSEjNSQiKzxNLSdSKSEjNSQiKlJMW0ApISIqJCIoZTk5KSEiKCQiK3llbW4iKSEjNSQiK2J2cyFHKSEjNSQiK1BTI2VZKSEjNSQiK2V4QjMoKSEjNSQiK0ZAdyUqKikhIzUkIisiZndVSiohIzUkIisjPShwZCcqISM1JCIrXjJ3LDUhIiokIis7KnooUTUhIiokIip0KEd3NSEiKSQiKzR0ejg2ISIqJCIrYyJmMzoiISIqJCIrRXAwKD0iISIqJCIrXjosQTchIiokIitkM1JiNyEiKiQiK3ljInBHIiEiKiQiK09sTzs4ISIqJCIrR0pmVjghIiokIismKnlebzghIiokIitkcDgiUiIhIiokIissVF42OSEiKiQiKik9eEg5ISIpJCIrZDAzWTkhIiokIiskNFUxWSIhIiokIissbG50OSEiKjdTJCIqInAjKSk0IyEiKSQiK1drTUc/ISIqJCIrTHAmeiY+ISIqJCIrKjNCcSk9ISIqJCIrXm1gOj0hIiokIitlKltRdSIhIiokIitFTFNzOyEiKiQiK1tjYCw7ISIqJCIrOHFQSjohIiokIioqPiU9WSIhIikkIitzeW8jUiIhIiokIitLXGlCOCEiKiQiK3ZCWGE3ISIqJCIqVi1fPSIhIikkIitbYEQ7NiEiKiQiKywqKlFbNSEiKiQiKyk9ZXgjKSohIzUkIitDLDw0IyohIzUkIisnZmNhaykhIzUkIitVcHBhIikhIzUkIit2MXNgeCEjNSQiK3JxdGJ1ISM1JCIraVozb3MhIzUkIitMOyQzPighIzUkIis5KHBvQCghIzUkIitAP25MdCEjNSQiK0BEKWVfKCEjNSQiKy0kSHl4KCEjNSQiMSwrK0tGQ3YhKSEjOyQiK2FqLTElKSEjNSQiKXlJZygpISIpJCIrM2o1SSIqISM1JCIrOV8oKTMmKiEjNSQiK1dHMyIqKSohIzUkIiskbyo9RjUhIiokIitNMHFrNSEiKiQiKzxDRCw2ISIqJCIrWSYpXE82ISIqJCIrdEI5cTYhIiokIitWeiQ+PyIhIiokIithXHBKNyEiKiQiK01oR2Y3ISIqJCIqaltZRyIhIikkIiotJnkyOCEiKSQiKz0kZShHOCEiKiQiKy9Rb1o4ISIqJCIrQiM9Wk8iISIqJCIreSNbK1EiISIqJCIqRHhRUiIhIik3UyQiKyI0ZEYqPiEiKiQiKj1oVyM+ISIpJCIrUCZvayY9ISIqJCIrPG81KXkiISIqJCIqZG8jPjwhIikkIitkKz5dOyEiKiQiKzRbQCJlIiEiKiQiK3J0aDc6ISIqJCIrUjheVzkhIiokIitGOiRvUCIhIiokIitkLFE0OCEiKiQiK01jIj5DIiEiKiQiKzAlcFU8IiEiKiQiKyV6dGs1IiEiKiQiKyVmcSlRNSEiKiQiK1ktbEAoKiEjNSQiK0dZTHUhKiEjNSQiKzw9KjNZKSEjNSQiK3JELykqeSEjNSQiK18uKlFTKCEjNSQiKipINycqcCEiKiQiKylHSicqbychIzUkIipjLlJcJyEiKiQiK2ozKzZrISM1JCIrVWcvTmshIzUkIitBXnRgbCEjNSQiKk8/NnYnISIqJCIxLCsrI0g5MCwoISM7JCIrKSpbYzt0ISM1JCIxLCsrOU9CY3chIzskIjEsKyspKSlSKj0hKSEjOyQiK2JIQidSKSEjNSQiK2YpMzh5KSEjNSQiK3QiZiZvIiohIzUkIisjW0pLYiohIzUkIitQaT1KKiohIzUkIismZnopSDUhIiokIitgKClHbDUhIiokIis6ejIqNCIhIiokIitcLy5KNiEiKiQiKjB6NDsiISIpJCIrVV8iKSk9IiEiKiQiKyVRKVs5NyEiKiQiKyNmMCFRNyEiKiQiK1oiSCVmNyEiKiQiKlluKXk3ISIpJCIrXHRZJ0giISIqJCIrLD9TNzghIiokIitfMCdvSyIhIio3UyQiK25mTy0+ISIqJCIqNzZoJD0hIikkIitKVUxxPCEiKiQiK0NtTC88ISIqJCIrLUokemoiISIqJCIraWpFcjohIiokIitqdWYvOiEiKiQiK2wiXCJROSEiKiQiK3RvLHM4ISIqJCIrQSlbaEkiISIqJCIrTiJ6LkMiISIqJCIrLFVddTYhIiokIisqUnkkMzYhIiokIit2QC5VNSEiKiQiK2JHb2QoKiEjNSQiKyoqKmVBNSohIzUkIis5IkdSWSkhIzUkIisqPUFmJnkhIzUkIitqTUglSCghIzUkIispeSgpcHonISM1JCIrYjRNI1EnISM1JCIqMmluMSchIiokIisoZmQ9J2UhIzUkIjEqKioqKjRLKHlyZCEjOyQiKzgzMSN6JiEjNSQiK2lFcDVmISM1JCIxKioqKipSeG42NichIzskIitMOXl2aiEjNSQiMSoqKioqKlw7NilvJyEjOyQiK25AQ01xISM1JCIqJXomSFMoISIqJCIpIz1heSghIikkIipPR1k8KSEiKiQiK2JpJlxjKSEjNSQiK2JOdF4qKSEjNSQiKzNZK0okKiEjNSQiKyNmayQqcCohIzUkIitkJCpRMDUhIiokIitcZ0BSNSEiKiQiKzlBQXI1ISIqJCIrV0BFLDYhIiokIitAMUNINiEiKiQiK0tzNmI2ISIqJCIqIkghKnk2ISIpJCIrciFlMT8iISIqJCIrNGBbPzchIiokIitkRV9RNyEiKiQiK3k5JFxEIiEiKiQiK2whKikpcDchIio3UyQiK1U+Y0M9ISIqJCIpakFnPCEiKCQiKzNmXidwIiEiKiQiKypHOUZqIiEiKiQiKztXZG86ISIqJCIrby47LzohIiokIitmbG1SOSEiKiQiK2wqcF9QIiEiKiQiK0YoWzVKIiEiKiQiKldmcEMiISIpJCIrbjknRz0iISIqJCIrVSp5Jj02ISIqJCIrKClmKVIwIiEiKiQiKztfMiIqKSohIzUkIitgODlVIyohIzUkIip0PikpZikhIiokIitYVTxxeiEjNSQiKyIqeWBvdCEjNSQiMSwrK00/QDRvISM7JCIrKHk/KTRqISM1JCIxKioqKipcUWIhKillISM7JCIrZCl6WGMmISM1JCIqcHMqXGAhIiokIispPUw6RCYhIzUkIis+Kj5qRSYhIzUkIitbaCRHUSYhIzUkIit5TjYlZSYhIzUkIisicCw6JmUhIzUkIitcISplbmghIzUkIistIlt3XichIzUkIitkYiIqKilvISM1JCIxLCsrLVFCdnMhIzskIit4JW9rbSghIzUkIitzTil6MCkhIzUkIisvPEBYJSkhIzUkIit0O05DKSkhIzUkIiszITNBPiohIzUkIjEsKytDQzJZJiohIzskIjEsKys9Ij5QKSkqISM7JCIrJCkpUi4tIiEiKiQiK1I3UV01ISIqJCIrT01UeTUhIiokIisxQ1MvNiEiKiQiK2E9T0c2ISIqJCIrRi9OXTYhIiokIitYX1lxNiEiKiQiK1FgJCkpPSIhIiokIiskNDhjPyIhIiokIis7Vyc0QSIhIio3UyQiK2AqSHF2IiEiKiQiKygzM1hwIiEiKiQiKyx6dEs7ISIqJCIrMT8qNGQiISIqJCIrKGVzKjM6ISIqJCIrMzpvWTkhIiokIio+YFVRIiEiKSQiK1JqIz1LIiEiKiQiKydvaiVmNyEiKiQiKzBoNyg+IiEiKiQiK3JKcE02ISIqJCIrYT8sczUhIiokIis9PSgqMzUhIiokIipRWWZYKiEiKiQiK1BTP0ApKSEjNSQiKzF1cSE+KSEjNSQiKyMqKWZGZCghIzUkIituKFEpeXAhIzUkIis8WVZCayEjNSQiK3lhbEJmISM1JCIrVW9GKVwmISM1JCIrJ2ZSZTsmISM1JCIrMVBqVFwhIzUkIitNYyxNWyEjNSQiKylcQTwlWyEjNSQiK3dgKFImXCEjNSQiKiUpZk46JiEiKiQiK3JTMEBhISM1JCIrWzo9UWQhIzUkIitRJT4mKjMnISM1JCIxKioqKip6Y1BHWSchIzskIjEsKys4KDQoW28hIzskIiskZSwqUnMhIzUkIitNLXpJdyEjNSQiMSwrKzpDKm8sKSEjOyQiMSwrK1V3YSVSKSEjOyQiK14nSDJ3KSEjNSQiK3MkNEg2KiEjNSQiK1oxLlwlKiEjNSQiKidIXG4oKiEiKiQiKy1JcjE1ISIqJCIrOzFzTTUhIiokIitLTXVnNSEiKiQiK01seiUzIiEiKiQiKiNmJHA1IiEiKSQiK0JOREY2ISIqJCIraDAoZTkiISIqJCIrXywkSDsiISIqJCIreigpZXk2ISIqN1MkIio9KCp6cCIhIikkIitaYT5QOyEiKiQiK0BqRHg6ISIqJCIrXG1XPDohIiokIitXdlVkOSEiKiQiKyJlWXJSIiEiKiQiKyJ6Km9POCEiKiQiKkZpaEYiISIpJCIrLC1oOjchIiokIisiMytdOiIhIiokIisqKWZBJTQiISIqJCIrOVs6TDUhIiokIitPZipvciohIzUkIispXClcKTQqISM1JCIrX21deSUpISM1JCIxKioqKipccmE7J3khIzskIityQ2RicyEjNSQiK1soMzRuJyEjNSQiK2NcTEBoISM1JCIqIlJKQmMhIiokIiowKVEmPiYhIiokIiowJ2ZjWyEiKiQiKyZHQ05pJSEjNSQiKyx6UDFYISM1JCIyKCoqKioqcGhCY10lISM8JCIrXiFcOWglISM1JCIrR0txMVshIzUkIitcQFhyXSEjNSQiK2MudydRJiEjNSQiK0hBa09kISM1JCIrOE1bM2ghIzUkIitibmojXCchIzUkIiszdHkiKW8hIzUkIistTUpxcyEjNSQiK0g1IVFsKCEjNSQiKyozLihHISkhIzUkIis5ejcjUikhIzUkIislKVFvVCgpISM1JCIqKCpcYTIqISIqJCIrM2YlPlIqISM1JCIrKW8qNCFwKiEjNSQiMSwrKzN6Q3AqKiEjOyQiK2ZMIkgtIiEiKiQiK0I2KnAvIiEiKiQiKlE2I3A1ISIpJCIrNC1tKjMiISIqJCIreUBYMzYhIiokIislUUBkNyIhIiokIis7YmhUNiEiKjdTJCIrdyh6Z2siISIqJCIqRzlwZSIhIikkIit5MXJHOiEiKiQiKzs+dHE5ISIqJCIrQXVnNzkhIiokIismeVFVTiIhIiokIislZXFjSCIhIiokIiszdihwQiIhIiokIisvRD55NiEiKiQiK0hcRz42ISIqJCIrdjY7ZzUhIiokIitwZHErNSEiKiQiKkZWJDMlKiEiKiQiK3p3bDApKSEjNSQiKzheJjQ/KSEjNSQiKyNvZCYpZighIzUkIjEsKysqPSFcMHEhIzskIisjNGU7VichIzUkIis8NSQqKillISM1JCIrJWUmKmZSJiEjNSQiK3lcJXonXCEjNSQiK2hKKltpJSEjNSQiKzZHRCVRJSEjNSQiKzAkPnhEJSEjNSQiK0ghSHVDJSEjNSQiK08weFdWISM1JCIrJ3pISWAlISM1JCIrb3E2I3olISM1JCIrQEtwLV4hIzUkIitPTVBbYSEjNSQiK0YqUmkiZSEjNSQiKnRzaz4nISIqJCIrdShwO2UnISM1JCIqIj4+bXAhIiokIiskRyNvWHQhIzUkIis+PWw7eCEjNSQiMSoqKioqXEFCajIpISM7JCIrJj4tQ1UpISM1JCIrJj1tSXYpISM1JCIrOSwjcDEqISM1JCIrMlonSE8qISM1JCIrd3dkUycqISM1JCIrbypIJioqKSohIzUkIitoaypSLCIhIiokIitSdkJPNSEiKiQiKXZ2YzUhIigkIitZSW12NSEiKiQiK1czMyQ0IiEiKiQiK2osOjQ2ISIqN1MkIitDKnosZyIhIiokIisoenREYSIhIiokIisjSD5nWyIhIiokIitbKHkoSDkhIiokIitYXlh0OCEiKiQiKig9InBKIiEiKSQiKzdoOmc3ISIqJCIrLHVCLjchIiokIisjM3hoOSIhIiokIilcJSopMyIhIigkIityJGY5LiIhIiokIitOJSo+TygqISM1JCIrQztiYCIqISM1JCIrSCwlb2MpISM1JCIrQmgpeSh6ISM1JCIrMyhcMVIoISM1JCIrKCkqXDsibyEjNSQiK2AhPixEJyEjNSQiKztNOD1kISM1JCIrSD5nSV8hIzUkIituKltcIVshIzUkIis3YSoqZlchIzUkIitQSlw4VSEjNSQiKy1LPHlTISM1JCIrJmVndzAlISM1JCIrcCI9WjklISM1JCIrYSVlTUslISM1JCIrciVcU2QlISM1JCIrWGQweFshIzUkIissRShlQCYhIzUkIitMPFF4YiEjNSQiKlAyOyZmISIqJCIrJjRiNUwnISM1JCIqMGwrcichIiokIjEsKysuP0YlMyghIzskIit1IVstWCghIzUkIjEqKioqKiozakQweSEjOyQiK3M/M1oiKSEjNSQiK0E6KFJaKSEjNSQiK0otZSV5KSEjNSQiKywzJ3oyKiEjNSQiK2NIXmAkKiEjNSQiKzBcLDYnKiEjNSQiK3AsZl0pKiEjNSQiK3gncHMrIiEiKiQiKyl6M3ktIiEiKiQiKyRId24vIiEiKiQiK3Y0SGs1ISIqJCIrWFRbITMiISIqN1MkIis8UlRmOiEiKiQiKjkoSC46ISIpJCIrd0tKWzkhIiokIit2bXMkUiIhIiokIitHKT0iUjghIiokIittMkslRyIhIiokIisoejIkSDchIiokIitOajV1NiEiKiQiKzp3cz02ISIqJCIrQzo5ajUhIiokIis+c0YyNSEiKiQiK3MieS9eKiEjNSQiKyVRU1IlKikhIzUkIis+OlF0JCkhIzUkIjEsKysuLF0reSEjOyQiK0NFKipHcyEjNSQiKzpzI1xtJyEjNSQiKyRlUXA2JyEjNSQiK0xsUCdmJiEjNSQiK25TTTxeISM1JCIrOHlaJ3AlISM1JCIrKHkyP04lISM1JCIrPT1iLFQhIzUkIiszMVRlUiEjNSQiK0JRVEZSISM1JCIrIilleC1TISM1JCIrPWt6cFQhIzUkIitzQUo0VyEjNSQiKiMqR0BxJSEiKiQiK207aUpdISM1JCIqZW9YUSYhIiokIispKilvM3YmISM1JCIxKioqKipIIylRSDcnISM7JCIrZTwwJlwnISM1JCIxLCsrLGUjRydvISM7JCIrLD4lR0EoISM1JCIrRVVQc3YhIzUkIispPmcjNHohIzUkIisqUnM8QikhIzUkIjEsKys8JzMnUSYpISM7JCIrKjNGKUcpKSEjNSQiK0hvJj01KiEjNSQiK20yWWQkKiEjNSQiKzJKdCZmKiEjNSQiKyJwIjQ8KSohIzUkIisneUFBKyIhIiokIitUeD9ANSEiKiQiKygzdihRNSEiKiQiKyNvW10wIiEiKjdTJCIpPTFCOiEiKCQiKykpKm8kbzkhIiokIiphJClbVCIhIikkIis8USg9TyIhIiokIitqUiEqMzghIiokIitqcXhiNyEiKiQiKygpM1ctNyEiKiQiK3A2ISpbNiEiKiQiK2ErOyY0IiEiKiQiK3N6PVQ1ISIqJCIrI1JJI3ApKiEjNSQiK1ktIkhLKiEjNSQiKztAV3MoKSEjNSQiK2BaND0jKSEjNSQiK05uW2h3ISM1JCIrPEg0MXIhIzUkIilXZ2RsISIpJCIreihHVC0nISM1JCIxLCsrWj1OO2IhIzskIitFM2VaXSEjNSQiMigqKioqKlxCc05qJSEjPCQiK1pHdyJIJSEjNSQiK0pQRVJTISM1JCIrai5cKilRISM1JCIrSEU6W1EhIzUkIitQLiozIlIhIzUkIisxK1hrUyEjNSQiK1RIcyFIJSEjNSQiK1FGMnJYISM1JCIrOiR5ISopWyEjNSQiK1tEXkpfISM1JCIrTyUqPSllJiEjNSQiMSoqKioqKioqeVVeZiEjOyQiKzhbVTpqISM1JCIrTmx1dm0hIzUkIiolWyUqR3EhIiokIislUSpHc3QhIzUkIisqNEVPcSghIzUkIisiKW9CQCEpISM1JCIrVlEjUUspISM1JCIrbTlZNScpISM1JCIrZE1jISkpKSEjNSQiK2pQKFE4KiEjNSQiKzgmZi9QKiEjNSQiKzR4byFmKiEjNSQiK0RqPSZ6KiEjNSQiK3puI1spKiohIzUkIit4WTE7NSEiKiQiKyJbIlFLNSEiKjdTJCIrNC5gIVwiISIqJCIrJz0rc1YiISIqJCIrQ2s5JlEiISIqJCIrdkFrTDghIiokIitDI1FBRyIhIiokIis2RnJJNyEiKiQiKzAwKip5NiEiKiQiK2MlZXE3IiEiKiQiK0YvIlwyIiEiKiQiK2xvXkE1ISIqJCIreS1DKXAqISM1JCIrbTZxbiIqISM1JCIrIypbPEwnKSEjNSQiK1gxKlw0KSEjNSQiKyYzO1piKCEjNSQiKyQpM2w6cSEjNSQiK1VwOSRbJyEjNSQiKzc9J1snZiEjNSQiK2FTJzNaJiEjNSQiKlsicDhdISIqJCIrMHY/M1khIzUkIis4KSkpM0YlISM1JCIqbyUzPVMhIiokIissbCVIJ1EhIzUkIis4QXU2USEjNSQiK1JTVWhRISM1JCIrbE9MK1MhIzUkIitVcHE2VSEjNSQiKz5cd3hXISM1JCIrPCQqWyN5JSEjNSQiKyJHZkY2JiEjNSQiKik+T2VhISIqJCIrMWNfNmUhIzUkIitmak5taCEjNSQiKnBTJD1sISIqJCIxLCsrW1IpUidvISM7JCIreURfK3MhIzUkIituMHlEdiEjNSQiK1U1LlF5ISM1JCIrSl8nZjgpISM1JCIrWmNrPSUpISM1JCIrMyRvYW8pISM1JCIrZCZbaCQqKSEjNSQiKzNKcnEiKiEjNSQiMSwrKzFaWipRKiEjOyQiK2FyKyRmKiEjNSQiK0csNiN5KiEjNSQiMSwrK3E8dmQqKiEjOyQiK19FNTc1ISIqN1MkIiteZ0toOSEiKiQiKzY2STQ5ISIqJCIrPGloZTghIiokIis5MmIzOCEiKiQiKyNvVydlNyEiKiQiK2RjbDM3ISIqJCIrNG5bZTYhIiokIitYITQiMzYhIiokIitKIzN2MCIhIiokIitzX2wxNSEiKiQiKyVvTF1iKiEjNSQiK0Q3LFMhKiEjNSQiKyY9QTdfKSEjNSQiK1RoMCoqeiEjNSQiK1QuL3Z1ISM1JCIrLyZHQiZwISM1JCIrKlwqKmZWJyEjNSQiK1JiSExmISM1JCIrZiI+UFgmISM1JCIqQnchNF0hIiokIitXRUs4WSEjNSQiKjE4PkclISIqJCIrKEhgLS4lISM1JCIrTE4oNChRISM1JCIrQ0BoNVEhIzUkIitKJmZzJVEhIzUkIis6UykzKFIhIzUkIitbKGZpOyUhIzUkIitCIipvO1chIzUkIiptWG5xJSEiKiQiKzI9YUJdISM1JCIqKCopKW9OJiEiKiQiKyc9WCopcCYhIzUkIitmM3RWZyEjNSQiKylbVm1RJyEjNSQiKzNxNUNuISM1JCIraTtKYHEhIzUkIitPZy5zdCEjNSQiK0xoYHl3ISM1JCIxLCsrXipwOSh6ISM7JCIxLCsrKkhzKVwjKSEjOyQiK3c8NjgmKSEjNSQiK1kocDN3KSEjNSQiKzpKNyQqKikhIzUkIit2Vzg1IyohIzUkIiphOENUKiEiKiQiK0tAcCsnKiEjNSQiK0c6KGV4KiEjNSQiKzwmeSpRKiohIzU3UyQiK1hXLk45ISIqJCIrazBFJVEiISIqJCIrR2MpW0wiISIqJCIrPm4+J0ciISIqJCIrPFdzUDchIiokIitOIzQjKj0iISIqJCIrKlJOMDkiISIqJCIrJSoqZTs0IiEiKiQiK20iZkQvIiEiKiQiKiczMUsqKiEiKiQiKncmZk4lKiEiKiQiK0QxeU4qKSEjNSQiKyM0YUNWKSEjNSQiK0MkZWcjeiEjNSQiK3NONj11ISM1JCIraSRROyJwISM1JCIxKioqKioqUVFaNmshIzskIituSFtDZiEjNSQiK01ea2ZhISM1JCIrJykzMkddISM1JCIrJ0hDR2slISM1JCIrOE9LPVYhIzUkIisjPlIqb1MhIzUkIistZGgxUiEjNSQiKzNOJHokUSEjNSQiKy41JD4nUSEjNSQiK1AjUSwoUiEjNSQiK2RlJSpbVCEjNSQiKzIqKipHUSUhIzUkIitsXktkWSEjNSQiK3VabmZcISM1JCIrT1smKXpfISM1JCIrWSEqKio0YyEjNSQiK3RRLldmISM1JCIrLGpFeGkhIzUkIitrbS4xbSEjNSQiMSwrK1chcHUjcCEjOyQiKy4wSFJzISM1JCIrJEcyKFJ2ISM1JCIrbEdNRnkhIzUkIisibyU+LCIpISM1JCIxLCsraEpmZyQpISM7JCIrblw8MCcpISM1JCIrIjRvWyQpKSEjNSQiKy4/KSlcISohIzUkIiouaDFEKiEiKiQiKy44KHlWKiEjNSQiK0VGTTcnKiEjNSQiK14iUl14KiEjNTdTJCIrQWJJNjkhIiokIio4SjxPIiEiKSQiKzEuaDg4ISIqJCIrcyhRaUUiISIqJCIrVCdSIj43ISIqJCIrVSpRPzwiISIqJCIrak4hWzciISIqJCIqaXd0MiIhIikkIitAKkgoSDUhIiokIitpUk09KSohIzUkIit4JlxsTCohIzUkIisuaWVeKSkhIzUkIjEsKysjUSFRaiQpISM7JCIrVnVWc3khIzUkIitFIm8tUSghIzUkIiopb3gqKW8hIiokIitfO2YwayEjNSQiK1RnQU1mISM1JCIrPD08JVsmISM1JCIrQCd5ZTEmISM1JCIrJXBfOnAlISM1JCIrKGYkZXVWISM1JCIrYXFIR1QhIzUkIittKGZRJ1IhIzUkIislM1F4KVEhIzUkIis7enMqKlEhIzUkIipHInojKlIhIiokIilzI1w6JSEiKSQiK3okMz9QJSEjNSQiK2pFQUlZISM1JCIrJik0WzxcISM1JCIrLFcmUUEmISM1JCIrLixbVGIhIzUkIitGYkFrZSEjNSQiKyY9M3Q9JyEjNSQiKyN6IilwXSchIzUkIis4JSpHP28hIzUkIjEsKytYaSpbNyghIzskIiszKmYqPXUhIzUkIissNzAseCEjNSQiKzBLN3F6ISM1JCIqTGZhQSkhIiokIis2OWxtJSkhIzUkIiskNHZOcCkhIzUkIisjKW9QMSopISM1JCIqcFNhNSohIiokIitYZ08iSCohIzUkIitVKT5cWSohIzUkIitUXypwaSohIzU3UyQiK002JSkqUSIhIiokIissc1RUOCEiKiQiKypwKFwlSCIhIiokIisoPihRWzchIiokIispPi9FPyIhIiokIitHLidvOiIhIiokIistJDM1NiIhIiokIis4JnpcMSIhIiokIitrKFEoPTUhIiokIio4b0RzKiEiKiQiK214LWIjKiEjNSQiK0cqPVh5KSEjNSQiMSwrK0AwMDYkKSEjOyQiKyMpKXleJHkhIzUkIiojKjQlZXQhIiokIispPlpOKW8hIzUkIipYQF1UJyEiKiQiKz0oPiFmZiEjNSQiK1ZzZEJiISM1JCIrKXk2Jj1eISM1JCIrR2pAYlohIzUkIileNFlXISIpJCIrW1hYLlUhIzUkIiszM21QUyEjNSQiK186J1wmUiEjNSQiK1s0dWJSISM1JCIrQ0NCTVMhIzUkIisiZmoqelQhIzUkIitZMDohUSUhIzUkIit6JEg+aSUhIzUkIitVSXUkKlshIzUkIit5ITNmPSYhIzUkIjEsKystMmghXCYhIzskIis2Tm8sZSEjNSQiKzkkcFU2JyEjNSQiMSwrKzVPYENrISM7JCIrQ2hXSG4hIzUkIit5d2ZFcSEjNSQiMSwrKzMnKjM5dCEjOyQiMSwrKz9tViFmKCEjOyQiK2p0YGF5ISM1JCIxLCsrWnlpMCIpISM7JCIrTiFcS00pISM1JCIrYnBAbiYpISM1JCIraUNpeCgpISM1JCIrOmt5dSopISM1JCIrMVBDZiIqISM1JCIrVVVwSiQqISM1JCIrOCpwSFwqISM1N1MkIitwX1FxOCEiKiQiK3peMUI4ISIqJCIrZm5IeDchIiokIitCVFJLNyEiKiQiKyUpPih5PSIhIiokIittLVZWNiEiKiQiK1slMyopNCIhIiokIitkckFhNSEiKiQiK25YTTQ1ISIqJCIrI1s0QmsqISM1JCIrb1ZmKT0qISM1JCIrOm42SygpISM1JCIrJyoqZkhGKSEjNSQiK0c4dDZ5ISM1JCIrTlQjKlx0ISM1JCIrelFEISpvISM1JCIrNylmcFYnISM1JCIrPG0jZipmISM1JCIrVVd2dWIhIzUkIitQc20jPSYhIzUkIittIipHSVshIzUkIis5UjVIWCEjNSQiKkZQL0glISIqJCIrJlsiKVE3JSEjNSQiKicpM2EuJSEiKiQiKSZbZS0lISIpJCIoTDA0JSEiKCQiK2smKVI/VSEjNSQiK0t5KFJTJSEjNSQiK3g9UUhZISM1JCIrMWdsJilbISM1JCIrJTQ/TTsmISM1JCIrd0EoXFgmISM1JCIrajA4YWQhIzUkIisnPSQpZjAnISM1JCIrRWhpY2ohIzUkIippUkhsJyEiKiQiK2xfWFVwISM1JCIrWWA/QnMhIzUkIjEqKioqKjRWV09cKCEjOyQiKyRwQkV2KCEjNSQiK2NhSioqeiEjNSQiK0EuQEwjKSEjNSQiMSwrKyEpRzJhJSkhIzskIipcQj5tKSEiKiQiK0opSHEmKSkhIzUkIitHLicpUiEqISM1JCIrJG9ZNUAqISM1JCIrOzFPciQqISM1N1MkIisjKnlyXzghIiokIitYc1gxOCEiKiQiK3I+emg3ISIqJCIrKFxZIT03ISIqJCIrZUN0dTYhIiokIiszLGFKNiEiKiQiK09vSCkzIiEiKiQiKzdJIlwvIiEiKiQiKyo0VDgrIiEiKiQiKylHMGJkKiEjNSQiKyx5O04iKiEjNSQiK0EjeUFwKSEjNSQiK0VKKXBDKSEjNSQiKzstJCoqeighIzUkIis7OmdfdCEjNSQiK2B6aTJwISM1JCIrJzRsIXBrISM1JCIxKioqKio+RTdELychIzskIitNQDpOYyEjNSQiMSoqKioqNCo9bGJfISM7JCIrWmYiUiJcISM1JCIrMVxkP1khIzUkIilVLydRJSEiKSQiKm0oPT5VISIqJCIrXy5uRFQhIzUkIitWVGwxVCEjNSQiKzQnMyVlVCEjNSQiKnhMSkYlISIqJCIrckNoU1chIzUkIisuKUgqXFkhIzUkIjIuKys1WCJ5ISpbISM8JCIrXi85YV4hIzUkIitQL1lLYSEjNSQiK2ImKmU+ZCEjNSQiKyU+djAsJyEjNSQiK29KWCxqISM1JCIqJClRISplJyEiKiQiKyV5KnlxbyEjNSQiKygpM25XciEjNSQiKyk+JjM0dSEjNSQiKyR6O0dtKCEjNSQiKzksKlwheiEjNSQiK3ZILk4iKSEjNSQiKmViRU4pISIqJCIrJD5CeWIpISM1JCIrL0p0XSgpISM1JCIrTFl6SiopISM1JCIxKioqKio0OiFlLCIqISM7JCIrdWx6ZyMqISM1N1MkIiszemtPOCEiKiQiKzxVUyJIIiEiKiQiKzpoelo3ISIqJCIqVGZePyIhIikkIistTCtqNiEiKiQiKjA1NTciISIpJCIrX1oqKnk1ISIqJCIrYilmby4iISIqJCIrcHheWCoqISM1JCIrXjJRPyYqISM1JCIrMysoSDQqISM1JCIrekdAaicpISM1JCIxKioqKioqenBJSiMpISM7JCIrY1olenooISM1JCIrVSZ5WE8oISM1JCIqKEh3THAhIiokIitsKXkkNGwhIzUkIispUVtuNCchIzUkIis5XW0tZCEjNSQiK1Q5Rk5gISM1JCIrO1MhUSsmISM1JCIqdic0PVohIiokIitfb3UoWyUhIzUkIitHViQ0SyUhIzUkIitRKj5JQSUhIzUkIitwWlMmPiUhIzUkIisielpeQiUhIzUkIitRL2ROViEjNSQiK2hcZyhbJSEjNSQiKzImKUgibyUhIzUkIjIuKyshRzEsMlwhIzwkIit2YTZjXiEjNSQiK3cpXDdVJiEjNSQiK3RQTCdwJiEjNSQiK0Q5VHdmISM1JCIrbUdZZGkhIzUkIitZW0NPbCEjNSQiK2ZHOTVvISM1JCIramIyeHEhIzUkIitjWFBOdCEjNSQiK19odyRlKCEjNSQiK0JCS0B5ISM1JCIqLTh1LykhIiokIislKioqb2gjKSEjNSQiK3QoZVNZKSEjNSQiKzp1bGEnKSEjNSQiKyVvSVEkKSkhIzUkIitrRTQtISohIzUkIipDIzRnIiohIio3UyQiKm80P0siISIpJCIrLT51eDchIiokIismKnk5TjchIiokIitbS2QkPiIhIiokIioocF9fNiEiKSQiKy1LbzY2ISIqJCIqKnklMzIiISIpJCIrUlciKkg1ISIqJCIrUG5CKSkpKiEjNSQiK3QqNGFaKiEjNSQiKz43WmchKiEjNSQiK2JAUVYnKSEjNSQiKk0kUkMjKSEiKiQiMSoqKioqZkA3VSF5ISM7JCIqaHpVUSghIiokIjEsKysmPmpxJ3AhIzskIiptd2liJyEiKiQiK0wkeXA6JyEjNSQiK2ZHZnZkISM1JCIrJFtwKD5hISM1JCIrUDY3KTQmISM1JCIrdkd3PlshIzUkIis9T2IkZiUhIzUkIjIoKioqKio0cVhxVSUhIzwkIit2IipIRFYhIzUkIit3RyopKkclISM1JCIrQUBhPVYhIzUkIislM15iUyUhIzUkIisqNCoqR2ElISM1JCIrIyk9YkBaISM1JCIyLisrZ1lLRCRcISM8JCIrai5rbl4hIzUkIitDQHY+YSEjNSQiK2B4KEdvJiEjNSQiKzI4Ml9mISM1JCIrdU9IQmkhIzUkIiswcUMkXCchIzUkIitNS0VmbiEjNSQiK1U3Pj5xISM1JCIrUW1KcnMhIzUkIit5QUk5diEjNSQiKzRrO1p4ISM1JCIxLCsrRzdBcHohIzskIjEsKysmKUgxIT0pISM7JCIxLCsraHpgeiQpISM7JCIrRyFHeGMpISM1JCIrcysiXHUpISM1JCIrOVdhNiopISM1JCIrN0FBbyEqISM1N1MkIipoZCczOCEiKSQiK2BvS2w3ISIqJCIrMV9xQjchIiokIis0IltKPSIhIiokIitOYztWNiEiKiQiKyNbQ001IiEiKiQiKy1Bc2o1ISIqJCIrJ3pWUi0iISIqJCIrTWVDUykqISM1JCIrVEVHUiUqISM1JCIreidmai4qISM1JCIrb0NbSicpISM1JCIxLCsrYDgjXEEpISM7JCIrUEdUPHkhIzUkIjEsKyt0UFA1dSEjOyQiK0txPTFxISM1JCIrJFIsJTNtISM1JCIrLS0kPUEnISM1JCIrMm5gX2UhIzUkIitUdHIyYiEjNSQiK3RjVCY+JiEjNSQiKyg+KDNDXCEjNSQiK0xvIj5xJSEjNSQiK3VCIWZgJSEjNSQiKy8oPTNWJSEjNSQiK21iTylRJSEjNSQiKydvM29TJSEjNSQiK1o+SSJbJSEjNSQiK2NEeC9ZISM1JCIrUiNSIXBaISM1JCIqSyF6bFwhIiokIitlMUUoPSYhIzUkIitfIyllRWEhIzUkIis4OyJ6biYhIzUkIisiUTtqJGYhIzUkIitRNncoPichIzUkIitzZiEqZWshIzUkIit6ay48biEjNSQiMSwrKzNgJSpwcCEjOyQiMSwrKzpOJ2VAKCEjOyQiK0cnKlJgdSEjNSQiKyo9PDpvKCEjNSQiK0tjWioqeSEjNSQiK3YsIm81KSEjNSQiKyRlOUxJKSEjNSQiKzpvKypbKSEjNSQiKztHNmsnKSEjNSQiK0YnRyFIKSkhIzUkIitwVEglKSopISM1N1MkIitNU1knSCIhIiokIitNSy5hNyEiKiQiKzFSTTg3ISIqJCIreEJ3dDYhIiokIitoISp6TTYhIiokIitvZzYnNCIhIiokIiskSCx2MCIhIiokIiteRCQpPTUhIiokIitLTlMrKSohIzUkIispM2YzVCohIzUkIitiKDMmPiEqISM1JCIra1ZQRScpISM1JCIrVUp3SiMpISM1JCIxLCsrRiU+ayR5ISM7JCIrSl90VHUhIzUkIitlSCxdcSEjNSQiK21UamttISM1JCIrdlg8IUgnISM1JCIrLjVQS2YhIzUkIis8Yyl6ZiYhIzUkIitwPWElSCYhIzUkIit0cyopSF0hIzUkIittIlE7IlshIzUkIitbbERZWSEjNSQiK0c5RVFYISM1JCIrQjxXKlslISM1JCIrOW9eKVwlISM1JCIqYHY4YyUhIiokIithSXpyWSEjNSQiK08leUIjWyEjNSQiK3B1WDBdISM1JCIrck9yOF8hIzUkIisnZmswVyYhIzUkIjEqKioqKmZfJ0ghbyYhIzskIitjcjFHZiEjNSQiK1xcIyl6aCEjNSQiK1QpPkFWJyEjNSQiKyY+JlwjbychIzUkIioyJlJHcCEiKiQiK1d1NG9yISM1JCIrW3U7K3UhIzUkIituPF1CdyEjNSQiK1BGSlB5ISM1JCIxLCsrcFwzVCEpISM7JCIrQEhiTSMpISM1JCIreSV6d1QpISM1JCIrJnlLMWYpISM1JCIrTEd2YCgpISM1JCIrekVfMiopISM1N1MkIitGd0omRyIhIiokIitzLnZWNyEiKiQiKz5hJlI/IiEiKiQiK0wnM2A7IiEiKiQiK3lFS0Y2ISIqJCIrJXlgJyozIiEiKiQiK0VGM181ISIqJCIrQi1bOTUhIiokIituTnJuKCohIzUkIitbWzoqUSohIzUkIisiR1AqMyEqISM1JCIrOD00RicpISM1JCIrOSRlUkMpISM1JCIrYCZ5LSd5ISM1JCIxKioqKipSakF1WighIzskIisvIilmKDQoISM1JCIrPSJRU3MnISM1JCIrJUckNGhqISM1JCIrbSs9OWchIzUkIisuPG4qbyYhIzUkIiskZSxZUiYhIzUkIitZc0hPXiEjNSQiKngnekBcISIqJCIrPWM5ZFohIzUkIis2QmhZWSEjNSQiKyMqKlE/ZiUhIzUkIisiW0VEZiUhIzUkIitgdGZXWSEjNSQiKy1xeFVaISM1JCIra0JTISlbISM1JCIrdlVTXV0hIzUkIisqcDpmQyYhIzUkIisxYmtnYSEjNSQiK1ZCMCpvJiEjNSQiK0Ikb2ojZiEjNSQiKzg/ZW9oISM1JCIrSCQ+QlQnISM1JCIrJUcnemFtISM1JCIrI1I8UCpvISM1JCIxLCsrRHlARnIhIzskIis5OCJRTighIzUkIitPXk1zdiEjNSQiK1ouKT55KCEjNSQiMSwrK0JyOSMpeiEjOyQiK0duX3MiKSEjNSQiK2MlR0lOKSEjNSQiKzdXd0ImKSEjNSQiKm1AXW8pISIqJCIrJT5JcyQpKSEjNTdTJCIrLyQ9XkYiISIqJCIrRTBRTTchIiokIisiZlZhPiIhIiokIis+SHBkNiEiKiQiKyJRVjE3IiEiKiQiK0ZxJVIzIiEiKiQiKy94UFo1ISIqJCIrTiIqejU1ISIqJCIrOSQ+OHUqISM1JCIrN3lKdCQqISM1JCIrNG16LiEqISM1JCIrIypcIUdqKSEjNSQiK0EqeTFFKSEjNSQiK18ieSIpKXkhIzUkIjEsKyswdGk7diEjOyQiKzdaOltyISM1JCIrKCp5JWV5JyEjNSQiK2I+JFFWJyEjNSQiK0x1Qig0JyEjNSQiKih5MCN5JiEiKiQiKnQlKltcJiEiKiQiK3JjZVVfISM1JCIrck1vSl0hIzUkIiswcyR5J1shIzUkIitZeTNiWiEjNSQiKyRROWBwJSEjNSQiKyU9TXpvJSEjNSQiKzBWLElaISM1JCIrSnl1O1shIzUkIitRVThVXCEjNSQiK2FAbio0JiEjNSQiK1dtJEhHJiEjNSQiK3drJGZbJiEjNSQiK2lXSi5kISM1JCIrWGZTSWYhIzUkIitrU0JqaCEjNSQiMSoqKioqejJJJSlSJyEjOyQiMSoqKioqUl4jPkxtISM7JCIrMy4+bG8hIzUkIittIz5ENCghIzUkIjEsKysiZVdPSighIzskIis6LVBGdiEjNSQiK0MsIkd0KCEjNSQiMSoqKioqKm8rTUh6ISM7JCIxLCsrKjMjZjsiKSEjOyQiKyczQFdIKSEjNSQiKz4iKSlHWSkhIzUkIitmSUFBJykhIzUkIitVJjNHeCkhIzU3UyQiK24leWRFIiEiKiQiK3J0JGVBIiEiKiQiK1lJcyg9IiEiKiQiKz06JDM6IiEiKiQiK3kmelk2IiEiKiQiK05lIip5NSEiKiQiK0B4SVY1ISIqJCIrKlI3eCsiISIqJCIrIltiL3MqISM1JCIqLytFTyohIiokIituMU8uISohIzUkIit5U3lVJykhIzUkIisoR0A3RykhIzUkIjEsKytWKj0lPnohIzskIjEqKioqKmYmR25ldiEjOyQiMSoqKioqNClcLSxzISM7JCIrLnZUXG8hIzUkIitPM3gybCEjNSQiKlZVND0nISIqJCIrYGhkdWUhIzUkIit0cihbZiYhIzUkIitbREJbYCEjNSQiK2NIeFNeISM1JCIyKCoqKioqXGh4eChcISM8JCIpRDRqWyEiKSQiKyY+OCcpeiUhIzUkIit1VS0leSUhIzUkIitBPCdvIlshIzUkIisxYyFIKlshIzUkIisoPmFuKyYhIzUkIitaKlxDOiYhIzUkIjEsKytTJ3pSSyYhIzskIipIaGNeJiEiKiQiKztgTEFkISM1JCIrciVbJVJmISM1JCIrXHgyamghIzUkIismUXkpKlEnISM1JCIrOlctPG0hIzUkIjEsKytzeDtVbyEjOyQiMSwrKyV5bkwxKCEjOyQiK3dNMHpzISM1JCIrTSp6elsoISM1JCIrXFRAKm8oISM1JCIrKlwqMyMpeSEjNSQiKjEjPW0hKSEiKiQiMSwrK0ImKkhUIykhIzskIitYMlgyJSkhIzUkIitpcyJbYykhIzUkIip4R1ByKSEiKjdTJCIrXCc+c0QiISIqJCIrKXlWIT03ISIqJCIrJzM+Mj0iISIqJCIrZjNsVzYhIiokIis+ImYkNDYhIiokIitrKCpbdTUhIiokIiteUSEpUjUhIiokIitNQzowNSEiKiQiK0FoWS8oKiEjNSQiK1knZmpOKiEjNSQiK3lXKnArKiEjNSQiK3cjPWtsKSEjNSQiMSwrK3NvKFxJKSEjOyQiK2BoVGB6ISM1JCIrdWopSGcoISM1JCIrY1RsYnMhIzUkIisqcEFVInAhIzUkIitFalQjZSchIzUkIit0bSNbRSchIzUkIisqeiV5bWYhIzUkIisjND1UcCYhIzUkIiokSCRHWCYhIiokIipAYSdbXyEiKiQiKyIqKVxsMyYhIzUkIitaXDxxXCEjNSQiKzheVixcISM1JCIqVVAtKVshIiokIitgVl8vXCEjNSQiKyozJWZxXCEjNSQiKywkKWV0XSEjNSQiKjpXITNfISIqJCIrYyNmJG9gISM1JCIrdC46XGIhIzUkIittIWZhdSYhIzUkIisiUWpHJmYhIzUkIispXC92OychIzUkIjEqKioqKnBPZ2dRJyEjOyQiK1g9cjBtISM1JCIxLCsrIVImM0NvISM7JCIrPmpAUnEhIzUkIis9ZlxccyEjNSQiKyhvU09YKCEjNSQiK2RCbl13ISM1JCIrPiJ6KVJ5ISM1JCIrJSpIej8hKSEjNSQiK1QocEo+KSEjNSQiK0ItKHBOKSEjNSQiKy1TSzcmKSEjNSQiKmw8JmYnKSEiKjdTJCIrY0FQXDchIiokIipZSjRAIiEiKSQiKz5UT3U2ISIqJCIrWWIzUjYhIiokIistJj1ZNSIhIiokIitAbGdxNSEiKiQiKWEhby4iISIoJCIrJj1nSSsiISIqJCIqQXJGcCohIiokIitAIUhTTiohIzUkIitpVjk5ISohIzUkIisvQzt0JykhIzUkIjEsKytZXVVKJCkhIzskIis2KWUnKil6ISM1JCIqNyczXHchIiokIitZVWU2dCEjNSQiK1AuJCl6cCEjNSQiK3BXTmRtISM1JCIrciteW2ohIzUkIiteQExlZyEjNSQiK1s0SiN6JiEjNSQiK0BaM2NiISM1JCIrXCFSXU4mISM1JCIrTmYlUT4mISM1JCIrL24qZjImISM1JCIrWFpSLl0hIzUkIjIuKytdTk1oKFwhIzwkIismZTVEKlwhIzUkIilwRlxdISIpJCIrVCdvPzkmISM1JCIqNndlRSYhIiokIitqRVw6YSEjNSQiKndBZWUmISIqJCIrQzc3c2QhIzUkIitgIio0cWYhIzUkIitbNChmPCchIzUkIitkW1gnUSchIzUkIip2UigpZichIiokIis4U1c1byEjNSQiKykpM2Q+cSEjNSQiMSoqKioqPjIjXENzISM7JCIrSzIpUVUoISM1JCIrVj5zO3chIzUkIisqZWRBIXkhIzUkIithaygqenohIzUkIitLIyllXCIpISM1JCIrL2wrNiQpISM1JCIrdGxKayUpISM1JCIrSlJ2NCcpISM1N1MkIisqKlE8VTchIiokIiswJFJXPyIhIiokIitkJylmbzYhIiokIisoPXhTOCIhIiokIis8MlMrNiEiKiQiKzkzQG41ISIqJCIrVyhlVS4iISIqJCIrRktRLDUhIiokIitGTSdbbyohIzUkIispNDheTiohIzUkIitfSExDISohIzUkIiokZWIjcCkhIiokIit4dTZnJCkhIzUkIiljc0YhKSEiKSQiMSwrK0RJYydwKCEjOyQiKnZHJW90ISIqJCIrVF4oZS8oISM1JCIrWl5ES24hIzUkIitaMHBKayEjNSQiKnhaKltoISIqJCIrKVEsIyopZSEjNSQiK1Bhd2RjISM1JCIrdGBxZmEhIzUkIitXNVcqSCYhIzUkIispWzAuPSYhIzUkIitwMT8vXiEjNSQiK3BYUHJdISM1JCIrJ2VGLzMmISM1JCIrJj0+JkdeISM1JCIrJlxDPEAmISM1JCIrLjhZRGAhIzUkIitKMylbWSYhIzUkIitPczxEYyEjNSQiMSoqKioqKikzQT0hZSEjOyQiKy5HbSEqZiEjNSQiKz1qKno9JyEjNSQiKkcmZiFSJyEiKiQiKzJFbCZmJyEjNSQiKzQnKnkrbyEjNSQiK0QkKSpSKyghIzUkIjEsKytpM2gucyEjOyQiMSoqKioqKj4qeiMpUighIzskIitUciVwZSghIzUkIitEbSIpb3ghIzUkIitsYkxWeiEjNSQiMSoqKioqUkltLDYpISM7JCIrajE9cCMpISM1JCIrLGJUPyUpISM1JCIreFExayYpISM1N1MkIisqZW9iQiIhIiokIipmNyYpPiIhIikkIistJXBMOyIhIiokIitmUmRINiEiKiQiK01gbCc0IiEiKiQiKyYpUkRrNSEiKiQiK0JpNks1ISIqJCIrbmIyKzUhIiokIis+Y0ghbyohIzUkIisnKj49ZiQqISM1JCIrZj85UCEqISM1JCIrdlg+OSgpISM1JCIrWXptIVIpISM1JCIrd1lDbiEpISM1JCIxLCsrNSVwXXUoISM7JCIxLCsrJ29uZVUoISM7JCIxLCsrTWUxN3IhIzskIismcF5vIW8hIzUkIipKTlReJyEiKiQiKkBBJVFpISIqJCIrPjVoJSlmISM1JCIrTSMzeHYmISM1JCIrOCtdaWIhIzUkIisqcHVKUyYhIzUkIitvTSNIRyYhIzUkIiooZmouXyEiKiQiK2ZgcGxeISM1JCIrZG8nejsmISM1JCIrQWInej8mISM1JCIrIypwOyNHJiEjNSQiKicqeWpRJiEiKiQiKzg+NDtiISM1JCIxKioqKipcPiZ5bWMhIzskIis6MDhNZSEjNSQiMSoqKioqcChcNzlnISM7JCIrJW9rSj8nISM1JCIrJD5pISlSJyEjNSQiMSoqKioqUktYZ2YnISM7JCIrPjV0JXonISM1JCIrNE01IypwISM1JCIrJ29zaz0oISM1JCIxLCsrbUJZd3QhIzskIitaIXo0YyghIzUkIitGZz5SeCEjNSQiK0UuXjV6ISM1JCIrYkZidSEpISM1JCIrTkw5SiMpISM1JCIxKioqKipSbSFHIVEpISM7JCIxLCsrS202QSYpISM7N1MkIitba11INyEiKiQiK1RBNSQ+IiEiKiQiK2kkRydlNiEiKiQiK2QjSGI3IiEiKiQiKldQTDQiISIpJCIrRkBwaDUhIiokIitGZExJNSEiKiQiKyJ6bTQqKiohIzUkIitIWW55JyohIzUkIisqemBlTyohIzUkIitfXT9fISohIzUkIitLQnRQKCkhIzUkIjEsKytRXnVBJSkhIzskIjEsKys0VSF6NSkhIzskIit1JD1WeighIzUkIisyNmokWyghIzUkIiprZiJ5ciEiKiQiMSwrK1RGJDQpbyEjOyQiK1lgbCZmJyEjNSQiK15FZ0VqISM1JCIrOm5TeWchIzUkIiswaCFlJmUhIzUkIitFZktqYyEjNSQiKyJHU1xdJiEjNSQiKyRlQVBRJiEjNSQiKig+YSxgISIqJCIrKlInKillXyEjNSQiK1c2KVtEJiEjNSQiK28nSHRHJiEjNSQiK1hcMmBgISM1JCIrOjpHW2EhIzUkIispUmwob2IhIzUkIitzeUU1ZCEjNSQiK0sncCdvZSEjNSQiK2RyNlNnISM1JCIqayg0QGkhIiokIitDRlwzayEjNSQiKzBfYipmJyEjNSQiKyxJIj56JyEjNSQiMSoqKioqNCMpUU4pcCEjOyQiMSoqKioqZkNORjwoISM7JCIrMz40ZXQhIzUkIjEqKioqKlJAJFxRdiEjOyQiK0t1MTh4ISM1JCIrJkglPSIpeSEjNSQiKzosVlUhKSEjNSQiMSwrK09eZSc+KSEjOyQiKjsyT00pISIqJCIrUFpnJFspISM1N1MkIitlPSVSQSIhIiokIit0VzspPSIhIiokIitoR0xhNiEiKiQiK0M7IT43IiEiKiQiK1xtUyE0IiEiKiQiKmgnW2Y1ISIpJCIrMyl6KUc1ISIqJCIrVz41JSkqKiEjNSQiK1BnbHonKiEjNSQiK1V3enUkKiEjNSQiKiplP3AhKiEiKiQiK0I6J0d3KSEjNSQiK2Y0MWMlKSEjNSQiK1QuV1wiKSEjNSQiMSwrK0slZVMleSEjOyQiKy5jXFR2ISM1JCIxLCsrIjNkUkMoISM7JCIrQDRLYXAhIzUkIisxWjV3bSEjNSQiK2ZvTzhrISM1JCIraUVccWghIzUkIisvYSg+JmYhIzUkIis5YjVpZCEjNSQiKzx5bS9jISM1JCIqYDpFWyYhIiokIisieS55UiYhIzUkIisqKmUjM04mISM1JCIrSiZ6NE0mISM1JCIrOil5ak8mISM1JCIxKioqKipcZiI9Q2EhIzskIioueTNeJiEiKiQiK0xoZUFjISM1JCIpVEliZCEiKSQiK2MvNjBmISM1JCIreDVJb2chIzUkIisjb3M5QychIzUkIitiMGNAayEjNSQiK2s7J2VnJyEjNSQiK05gLCN6JyEjNSQiK2ZfKnoocCEjNSQiKzt5M2lyISM1JCIxLCsrXnYnR00oISM7JCIrJmUkPT52ISM1JCIrXCpSLHAoISM1JCIrTyRwXSZ5ISM1JCIreiU9TiwpISM1JCIrPS1CbCIpISM1JCIxLCsrXzI3NSQpISM7JCIrRVtFWyUpISM1N1MkIislUk0pPTchIiokIitrI2ZPPSIhIiokIipRVy86IiEiKSQiKzxQbD02ISIqJCIreHIjeTMiISIqJCIrT0RnZDUhIiokIit6YHJGNSEiKiQiMSwrKy5AJSl6KiohIzskIitiQCRIbyohIzUkIitHNnMmUSohIzUkIitZaSd5MyohIzUkIispcEMkKnkpISM1JCIxLCsrWjhQIVwpISM7JCIrOiE+Oz4pISM1JCIqeSMzJSp5ISIqJCIraW9GKmYoISM1JCIxLCsrWnhINHQhIzskIisyNilvLSghIzUkIitrNlFibiEjNSQiKjpLJylcJyEiKiQiMSoqKioqejQ3M0UnISM7JCIrKFJ4aC8nISM1JCIrTnkhKWVlISM1JCIrLCQ9QnEmISM1JCIrTjBiemIhIzUkIitSTk0jXCYhIzUkIitESVBUYSEjNSQiKyMpSDZFYSEjNSQiKzRyI1xXJiEjNSQiK09oRSZcJiEjNSQiMSoqKioqXHBAUmQmISM7JCIrLSkqR3hjISM1JCIxKioqKipIJD5oLGUhIzskIityVztWZiEjNSQiMSoqKioqND4jUik0JyEjOyQiK05FKlJFJyEjNSQiMSoqKioqSEt4cFYnISM7JCIrIjN6WWgnISM1JCIrUiplWnonISM1JCIqeSM+dnAhIiokIitgJWZVOighIzUkIitoYV5JdCEjNSQiK1EhKnktdiEjNSQiK3pAOnF3ISM1JCIrKiplIT4keSEjNSQiMSoqKioqPlJndil6ISM7JCIraWgjbzgpISM1JCIxLCsrJHp3JnojKSEjOyQiKyQpbyVlVCkhIzU3UyQiKyopcDk5NyEiKiQiK080Yno2ISIqJCIrKDRHcDkiISIqJCIrbj92OjYhIiokIitkbGMmMyIhIiokIismSDRnMCIhIiokIis9QSJvLSIhIiokIis8ZSF6KCoqISM1JCIrajBCKW8qISM1JCIrSCNvJClSKiEjNSQiK183JXo1KiEjNSQiKigpKilvIikpISIqJCIrc25ZRCYpISM1JCIqdV1VQikhIiokIjEsKytnMUFXeiEjOyQiMSoqKioqZixAb2woISM7JCIrdyZcU1AoISM1JCIrXj1eKTQoISM1JCIrQztSTG8hIzUkIisoUVRCZSchIzUkIitROUtcaiEjNSQiK2xgUVFoISM1JCIrMSg9TSZmISM1JCIra1UoeXomISM1JCIqUClcdWMhIiokIitSQzYmZSYhIzUkIitcKmUvYCYhIzUkIjEsKytGdjs1YiEjOyQiK2lsI0dfJiEjNSQiK3IkXGhjJiEjNSQiMSoqKioqUi0uc2omISM7JCIrYHFrS2QhIzUkIitlPyYqW2UhIzUkIitPNWUjKWYhIzUkIistKUcsOCchIzUkIiswSlMpRychIzUkIit0RlthayEjNSQiK3R4dURtISM1JCIrT2kpKSp6JyEjNSQiKzIielsocCEjNSQiK0AtK1xyISM1JCIrOVx6P3QhIzUkIit6KG8hKlsoISM1JCIrW2MnR2woISM1JCIxLCsrLCxZNnkhIzskIitHIkdWJ3ohIzUkIitASjk2IikhIzUkIitYdHVeIykhIzUkIis7XTgnUSkhIzU3UyQiKylcWSk0NyEiKiQiKi8yZTwiISIpJCIqel9QOSIhIikkIisuazs4NiEiKiQiK3hlZiQzIiEiKiQiK1UpeVkwIiEiKiQiK2RXOUU1ISIqJCIxKioqKioqKkhPIXkqKiEjOyQiK2NeSiZwKiEjNSQiKyc0NURUKiEjNSQiK2guQUgiKiEjNSQiKiMzT1gpKSEiKiQiKyo+ZDZjKSEjNSQiKz5mO3gjKSEjNSQiMSoqKioqKj0xSyUqeiEjOyQiK3MjKSpScighIzUkIitzJTMiUXUhIzUkIislPUcicHIhIzUkIit3PzM1cCEjNSQiKyI9XldtJyEjNSQiKyhRLGdWJyEjNSQiK2RqZkdpISM1JCIxKioqKipmJCpRZi8nISM7JCIrcyVSOCplISM1JCIrYEpYbmQhIzUkIipNIjN3YyEiKiQiK1hzLT1jISM1JCIxKioqKipITmRJZiYhIzskIjEqKioqKj4/ZSoqZiYhIzskIis9TG9PYyEjNSQiK1YvYitkISM1JCIrJSlSWSl5JiEjNSQiK145NigqZSEjNSQiMSoqKioqSFpTSi0nISM7JCIrKkcpR2poISM1JCIxKioqKioqMzpaOWohIzskIjEqKioqKipHcSVRWichIzskIiswSyUpUW0hIzUkIitRNDwybyEjNSQiKyhRR28ocCEjNSQiK0gjKTNZciEjNSQiK0xWWzh0ISM1JCIxLCsrI3AreVooISM7JCIxLCsrdGoxUXchIzskIitYKj5OeighIzUkIiskKUdoVnohIzUkIislZiIpejMpISM1JCIrVkdWRSMpISM1JCIrKWVJKmUkKSEjNTdTJCIqYC1mPyIhIikkIitBIylSczYhIiokIisqMyEqMzkiISIqJCIrVCdwMzYiISIqJCIreSEqKT0zIiEiKiQiK1JqZWA1ISIqJCIrIjQpb0Q1ISIqJCIrKT05KykqKiEjNSQiK28icFJxKiEjNSQiK1QoW3pVKiEjNSQiMSoqKioqcCZvXl4iKiEjOyQiK0ZYY3UpKSEjNSQiKz5OSChmKSEjNSQiKypcRi1LKSEjNSQiKlRlVS8pISIqJCIrMiEzMngoISM1JCIxLCsrViMqUSx2ISM7JCIrQ09tUXMhIzUkIispKSozYSlwISM1JCIrKW9VXHUnISM1JCIrVkomM18nISM1JCIrMSg9b0onISM1JCIrV3hRT2ghIzUkIisvJ0dGKWYhIzUkIissLlVlZSEjNSQiKkhTX3cmISIqJCIrLFQvL2QhIzUkIipMPlpuJiEiKiQiK3QwQndjISM1JCIraGV1MWQhIzUkIitWdiJRdyYhIzUkIjEqKioqKipvWmRXZSEjOyQiK2Y7ImYlZiEjNSQiMSoqKioqenRdWTEnISM7JCIrMDRuKD4nISM1JCIraCgpKj5NJyEjNSQiKyp6bFtcJyEjNSQiK0tCdmBtISM1JCIrWSozayJvISM1JCIrYSJSMylwISM1JCIrIm8/YDkoISM1JCIraERRM3QhIzUkIisoKmV6b3UhIzUkIitwJmZiaSghIzUkIis0cip5eCghIzUkIitraEFEeiEjNSQiK0Q7Om4hKSEjNSQiMSwrK0o+WC4jKSEjOyQiK1ZEME0kKSEjNTdTJCIrNnlHLTchIiokIituISlIcDYhIiokIitFWEpRNiEiKiQiK1xxJCkzNiEiKiQiK0RGVSEzIiEiKiQiK2MmNEYwIiEiKiQiK00/VUQ1ISIqJCIrJiplaiQpKiohIzUkIitNLyw5KCohIzUkIis2aV1XJSohIzUkIittVW11IiohIzUkIiovW1YhKikhIiokIitpSXRMJykhIzUkIisnKT1KaiQpISM1JCIrKkdOUjQpISM1JCIrc1knbyN5ISM1JCIrNyJHUWMoISM1JCIrOVkyMnQhIzUkIitCIVwkZnEhIzUkIisqKlwhUSNvISM1JCIqQXhRZychIiokIit6YTIuayEjNSQiMSoqKioqekspeUNpISM7JCIrWSNwPzInISM1JCIrWSI9dSVmISM1JCIqTiNmX2UhIiokIisiNCJcKXkmISM1JCIrbSczXnYmISM1JCIqNHE6diYhIiokIisoPVNpeCYhIzUkIitNSClvI2UhIzUkIisjZlEzIWYhIzUkIitiVD4mKmYhIzUkIitISCVwNSchIzUkIip4KjRMaSEiKiQiKzspKXpxaiEjNSQiK2pITjxsISM1JCIrIlspSHFtISM1JCIxKioqKipSOjh1I28hIzskIjEsKytJWXMnKXAhIzskIjEsKyspbz1sOSghIzskIitaOUowdCEjNSQiKyVbbz1ZKCEjNSQiMSoqKioqUkBwXmgoISM7JCIxLCsrQ19Ua3ghIzskIitiYioqM3ohIzUkIjEsKys3ZltbISkhIzskIitcUGojPSkhIzUkIitXYEw2JCkhIzU3UyQiK0VyKCopPiIhIiokIitPQlttNiEiKiQiKyd6LWc4IiEiKiQiK3dsLzI2ISIqJCIrKlt2Ino1ISIqJCIrJj1HPzAiISIqJCIrKj1GYC0iISIqJCIxKioqKioqZUZzKSkqKiEjOyQiK2VpRUQoKiEjNSQiK2IqSD9ZKiEjNSQiK1VUXyk+KiEjNSQiKzxeZU0qKSEjNSQiK09FT3EnKSEjNSQiKyYpKT1qUykhIzUkIis3REVWIikhIzUkIioncFIjKXkhIiokIitRO1BEdyEjNSQiKlVHVlAoISIqJCIrUk4pPTgoISM1JCIrbSNcNSFwISM1JCIrbmk0Jm8nISM1JCIrUW9RKFsnISM1JCIxKioqKipcI2U8NmohIzskIitfX1JmaCEjNSQiK2tNWk1nISM1JCIrWydcIlFmISM1JCIraTZPcmUhIzUkIisleiY+TWUhIzUkIilQI2YjZSEiKSQiKyIzKDNYZSEjNSQiKyRlWScqKWUhIzUkIistTDhkZiEjNSQiK3VWI1svJyEjNSQiK1AnbylcaCEjNSQiKypwPSVwaSEjNSQiK1JNcitrISM1JCIrdXU5VGwhIzUkIitPJTMkKW8nISM1JCIrZmMsU28hIzUkIikpPlYqcCEiKSQiMSwrK28mNCZcciEjOyQiK0FmNS90ISM1JCIxLCsrI1Fmb1goISM7JCIxKioqKipwJCpRbmcoISM7JCIrQGgiSHYoISM1JCIrb153JSp5ISM1JCIrZiwkPS4pISM1JCIrWypHUTspISM1JCIrbSlHMUgpISM1N1MkIitFeCVmPiIhIiokIitrIUhSOyIhIiokIisqZkxSOCIhIiokIisuIXlhNSIhIiokIitKI0cieTUhIiokIitoVl9eNSEiKiQiK0xrUUQ1ISIqJCIrWls2JioqKiEjNSQiKzEmKmVQKCohIzUkIio3IVEhWyohIiokIitjbSdIQSohIzUkIis0MjtsKikhIzUkIitfKHlxcSkhIzUkIitGJWYiXCUpISM1JCIram47Iz4pISM1JCIqWFxzJHohIiokIisneiEpZm8oISM1JCIrWzVTU3UhIzUkIitwdysucyEjNSQiMSwrK15SbndwISM7JCIrIm9JWHcnISM1JCIxKioqKipIVyZ6cGwhIzskIipAKmUmUichIiokIitIXHVXaSEjNSQiK1goPic+aCEjNSQiK3hVJD4tJyEjNSQiKz1nbF9mISM1JCIxKioqKipmR2o+ImYhIzskIio4XSMqKmUhIiokIittR0E4ZiEjNSQiKzorLV9mISM1JCIrWSNmTCwnISM1JCIrRVdvJTQnISM1JCIrcztJJD4nISM1JCIrOTtcMWohIzUkIitBdmZKayEjNSQiK3lrNG1sISM1JCIrc2ZqMm4hIzUkIitMKW9TJm8hIzUkIjEsKys6SloucSEjOyQiK1clW1Q6KCEjNSQiK3IiPVlJKCEjNSQiMSwrKyd6QE9YKCEjOyQiK3JzNit3ISM1JCIrPiNmS3UoISM1JCIrbk1SIyl5ISM1JCIrclMvPCEpISM1JCIrPSoqKm85KSEjNSQiK0V3enIjKSEjNTdTJCIrQyF6Sj4iISIqJCIrSyQ9OzsiISIqJCIrIj0pM0s2ISIqJCIrIilINi82ISIqJCIrS1JFeDUhIiokIiteOT1eNSEiKiQiK19aZUQ1ISIqJCIrKilwRSs1ISIqJCIrWSNbM3YqISM1JCIrZWJWKlwqISM1JCIrWE4peUMqISM1JCIrKj52ZioqKSEjNSQiKz0lKnpWKCkhIzUkIjEsKysjZmk8XCkhIzskIis2S2ZTIykhIzUkIis7MVEiKnohIzUkIisqW0hjdSghIzUkIigiRzB2ISIoJCIrZkhzc3MhIzUkIit0NnFdcSEjNSQiMSwrKzwzQFVvISM7JCIrNFJMXW0hIzUkIitgOjJ5ayEjNSQiKiNIO0dqISIqJCIrYl0qRz8nISM1JCIrRUAoUjUnISM1JCIrI1wiUktnISM1JCIrJWUtJSkpZiEjNSQiKzkvX3JmISM1JCIrdSIpZiEpZiEjNSQiMSoqKioqXDxQUiwnISM7JCIrWHVVcGchIzUkIit2R25XaCEjNSQiMSoqKioqPltFckInISM7JCIrJWYkPldqISM1JCIrR2RLamshIzUkIitSNzIjZichIzUkIis7PDlHbiEjNSQiKzFhVnBvISM1JCIreiVbUywoISM1JCIrNyIqSGdyISM1JCIrayIzbkkoISM1JCIrNHYsX3UhIzUkIis2UTwmZighIzUkIjEsKys1dUlOeCEjOyQiK0ckWzwoeSEjNSQiKyVvKSpSKykhIzUkIjEsKytnJj08OCkhIzskIis3c3JhIykhIzU3UyQiKyR5XjE+IiEiKiQiK1I9YGY2ISIqJCIrJXlbLzgiISIqJCIrTF8kSDUiISIqJCIrSm1jdzUhIiokIityXSk0MCIhIiokIitcITNmLSIhIiokIitfaTcsNSEiKiQiKylHQlx3KiEjNSQiK1gmKTM+JiohIzUkIisqXHdKRiohIzUkIisuSSVwLSohIzUkIis0JVsveSkhIzUkIisoM29TYCkhIzUkIitmJCpbKUcpISM1JCIxLCsrU1F2VyEpISM7JCIxLCsrSERIL3khIzskIit0Wycqb3YhIzUkIitELi9UdCEjNSQiMSwrK1Q7OkJyISM7JCIrYFo7PXAhIzUkIitYWS9IbiEjNSQiK01DbmVsISM1JCIrLj9wNGshIzUkIittLU0lRychIzUkIitEXkglPSchIzUkIitGLWU1aCEjNSQiKiVRXmpnISIqJCIrX15yVWchIzUkIitDPzxaZyEjNSQiK1IsTXZnISM1JCIxKioqKipcYWpfNychIzskIisqbyxaPichIzUkIitvY0MiRychIzUkIismW0BDUSchIzUkIitoLHkmXCchIzUkIipEXyo9bSEiKiQiMSoqKioqZlo1KFxuISM7JCIrLzMqZilvISM1JCIrJClRI2YtKCEjNSQiK3pXJHk7KCEjNSQiMSwrK04iZS1KKCEjOyQiKmpHPlgoISIqJCIrdFF5ImYoISM1JCIrRkslKkd4ISM1JCIrcmFyaXkhIzUkIisvdmQjKnohIzUkIisvdzs9IikhIzUkIitFN0ZSIykhIzU3UyQiKz0pWyQpPSIhIiokIitHSGxkNiEiKiQiK2MnKioqRzYhIiokIis7JUg+NSIhIiokIitAQS13NSEiKiQiK2k7IzQwIiEiKiQiK3BVTUU1ISIqJCIrKlt4PysiISIqJCIrZjxyeigqISM1JCIrcE5DUiYqISM1JCIrJCozdylIKiEjNSQiKzB6KXowKiEjNSQiKzs1J3AiKSkhIzUkIitXKy13JikhIzUkIitqaCJlTCkhIzUkIjEsKytLQ00oNCkhIzskIjEsKytuSyc+J3khIzskIioiKVs5aighIiokIit1Xyh6UyghIzUkIitbWzAlPighIzUkIitqJ1FDKnAhIzUkIis4OShmIW8hIzUkIitBWlZQbSEjNSQiK0RLUSpbJyEjNSQiK0omKSpSTychIzUkIismKnkkSEUnISM1JCIrVnBDKD0nISM1JCIrVEdJUGghIzUkIitLRSNHNichIzUkIis6cCJINichIzUkIitQOD1PaCEjNSQiK0ppITM9JyEjNSQiKydRJnBXaSEjNSQiK0svZERqISM1JCIrTl0yQGshIzUkIitEKmUpR2whIzUkIjEqKioqKj5GS21rJyEjOyQiKmNFQXgnISIqJCIxLCsrdy9pLnAhIzskIjEsKytHLCkqUXEhIzskIjEqKioqKkhmVG08KCEjOyQiMSwrKz09Ojp0ISM7JCIxLCsrcEtDYHUhIzskIjEsKys8I1IpKmUoISM7JCIxLCsrME8wQ3ghIzskIitUKHpeJnkhIzUkIitwR24jKXohIzUkIitKVzkxIikhIzUkIitRZE5EIykhIzUtJSVWSUVXRzYlJShERUZBVUxURyUoREVGQVVMVEc7JCIiJCEiIiQiIzohIiItJiUmX0FYSVNHNiMiIiI2Ji0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqTElORVNUWUxFRzYjIiIhLSUqVEhJQ0tORVNTRzYjIiIhLSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0mJSZfQVhJU0c2IyIiIzYmLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSYlJl9BWElTRzYjIiIkNiYtJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIS0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJStBWEVTTEFCRUxTRzYlLUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NVEiYTYiLyUnZmFtaWx5R1EhNiIvJSVzaXplR1EjMTA2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EldHJ1ZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ2l0YWxpYzYiLUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NVEiYjYiLyUnZmFtaWx5R1EhNiIvJSVzaXplR1EjMTA2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EldHJ1ZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ2l0YWxpYzYiUSE2Ii0lK0xJR0hUTU9ERUxHNiNRKExJR0hUXzM2Ii0lK1BST0pFQ1RJT05HNiwkIiknUmFXIiEiKSQiKHlJMSohIigkIilFSnJSISIpJCEpYHYqNCQhIikkIilGPUVVISIpJCEoM2xeKSEiKCQhKEVwUiohIigkIiIhISIiJCIpOj8/TSEiKSQiIzUhIiItJSpBWEVTU1RZTEVHNiMlJEJPWEctJSlfVklTSUJMRUc2IyIiIi0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiIiEhIiItJSlCT1VORFNfWUc2IyQiIiEhIiItJS1CT1VORFNfV0lEVEhHNiMkIiUrcSEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiUrcSEiIi0lKUNISUxEUkVORzYiLSUrQU5OT1RBVElPTkc2Jy0lKUJPVU5EU19YRzYjJCIiISEiIi0lKUJPVU5EU19ZRzYjJCIiISEiIi0lLUJPVU5EU19XSURUSEc2IyQiJStxISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJStxISIiLSUpQ0hJTERSRU5HNiI=NiI=

Figure 3 The surface LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2I1EhRictRiw2JlEjU1NGJy8lJ2l0YWxpY0dRJXRydWVGJy8lK2V4ZWN1dGFibGVHUSZmYWxzZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYnLUkjbW9HRiQ2LVEwJkFwcGx5RnVuY3Rpb247RicvRjlRJ25vcm1hbEYnLyUmZmVuY2VHRjcvJSpzZXBhcmF0b3JHRjcvJSlzdHJldGNoeUdGNy8lKnN5bW1ldHJpY0dGNy8lKGxhcmdlb3BHRjcvJS5tb3ZhYmxlbGltaXRzR0Y3LyUnYWNjZW50R0Y3LyUnbHNwYWNlR1EmMC4wZW1GJy8lJ3JzcGFjZUdGUy1JKG1mZW5jZWRHRiQ2JS1GIzYnLUYsNiZRImFGJ0YyRjVGOC1GPjYuUSIsRidGNUZBRkMvRkZGNEZHRklGS0ZNRk9GUS9GVVEsMC4zMzMzMzMzZW1GJy1GLDYmUSJiRidGMkY1RjhGNUZBRjVGQUYrRjVGQUYrRjVGQQ== showing existence of a minimum

<Text-field style="Heading 1" layout="Heading 1">Minimization </Text-field>

What remains is to compute a minimum, a task Maple performs with its NLPSolve command from the Optimization package. By option, this command can implement the Nelder-Meade (simplex) algorithm that is Maple's only direct method - it does not differentiate the objective function. In two dimensions, the simplex is a triangle. Function values are obtained at the vertices of the triangle, and for a minimization, the triangle is reflected away from the vertex with the highest value. Various scaling and anti-cycling devices are included in the algorithm.

Using this command, we determine

R := NLPSolve(SS(a,b), method=nonlinearsimplex, initialpoint=[a=5,b=5], evaluationlimit=250);

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSJSR0YoNyQkIjMmKjRrN2pNLzJRISM9NyQvSSJhR0YoJCIwVTxqIilIai4jISM5L0kiYkdGKCQiMFxtMUI0WyslRjc3IzckRi83JC9GNF00MDAwNEE2NzQ5Njc5RkNCL0Y5XTQwMTAwNEVDQjYwMUI5RDY=

where the first number in the output list is the minimum function value, and the sublist contains the minimizing values of the parameters LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JlEiYUYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnRjIvRjZRJ25vcm1hbEYn and LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JlEiYkYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnRjIvRjZRJ25vcm1hbEYn. Surprisingly, they are very close toLUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictSShtZmVuY2VkR0YkNiUtRiM2Jy1JI21uR0YkNiVRIjJGJy8lK2V4ZWN1dGFibGVHUSZmYWxzZUYnLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGJy1JI21vR0YkNi5RIixGJ0Y4RjsvJSZmZW5jZUdGOi8lKnNlcGFyYXRvckdRJXRydWVGJy8lKXN0cmV0Y2h5R0Y6LyUqc3ltbWV0cmljR0Y6LyUobGFyZ2VvcEdGOi8lLm1vdmFibGVsaW1pdHNHRjovJSdhY2NlbnRHRjovJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR1EsMC4zMzMzMzMzZW1GJy1GNTYlUSI0RidGOEY7RjhGO0Y4RjtGK0Y4Rjs= used to generate the noisy obserations.

The optimum values of the parameters, namely

A := eval(a,R[2]); B := eval(b,R[2]);

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSJBR0YoJCIwVTxqIilIai4jISM5NyNdNDAwMDRBNjc0OTY3OUZDQg==

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSJCR0YoJCIwXG0xQjRbKyUhIzk3I100MDEwMDRFQ0I2MDFCOUQ2

cause the differential equation governing our model to become

qq := eval(DE, [a=A, b=B]);

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSNxcUdGKC8sKC1JJWRpZmZHRiY2JC1JInlHRig2I0kidEdGKC1JIiRHRiY2JEY2IiIjIiIiKiYtRjE2JEYzRjZGOy1JJGFic0dGJjYjRj1GOyQiMFU8aiIpSGouIyEjOUYzJCIwXG0xQjRbKyVGRCIiITcjLywoLUYxNiRGPUY2RjtGPF00MDAwNEE2NzQ5Njc5RkNCRjNdNDAxMDA0RUNCNjAxQjlENkZH

A numeric integration of this equation results in Figure 4 where the noisy observations, and two numeric solutions of the model can be seen. In green, we have the solution for the derived parameters, and in red, the solution for the original parameters,LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictSShtZmVuY2VkR0YkNiUtRiM2Jy1JI21uR0YkNiVRIjJGJy8lK2V4ZWN1dGFibGVHUSZmYWxzZUYnLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGJy1JI21vR0YkNi5RIixGJ0Y4RjsvJSZmZW5jZUdGOi8lKnNlcGFyYXRvckdRJXRydWVGJy8lKXN0cmV0Y2h5R0Y6LyUqc3ltbWV0cmljR0Y6LyUobGFyZ2VvcEdGOi8lLm1vdmFibGVsaW1pdHNHRjovJSdhY2NlbnRHRjovJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR1EsMC4zMzMzMzMzZW1GJy1GNTYlUSI0RidGOEY7RjhGO0Y4RjtGK0Y4Rjs=.

Q := dsolve({qq, ic},y(t),numeric): p4 := odeplot(Q,[t,y(t)],0..10, numpoints=100, color=green, thickness=3): display([p4,p1,p3]);

NiwtJSdDVVJWRVNHNiU3YHE3JCQiIiEhIiIkIiImISIiNyQkIjAsLCwsLCwsIiEjOiQiMSZvJW9BUykpKWYnISM7NyQkIjAtLS0tLS0tIyEjOiQiMXg/Xnd6RillKCEjOzckJCIwLi4uLi4uLiQhIzokIjFKWGJbJEdSOSkhIzs3JCQiMC8vLy8vLy8lISM6JCIwKHlKM11qTiQpISM6NyQkIjAwMDAwMDAwJiEjOiQiMEhqeSh6InA9KSEjOjckJCIwMTExMTExMSchIzokIjAqbzIlby11cyghIzo3JCQiMDIyMjIyMjIoISM6JCIxJ3AmeiRbJkdEcSEjOzckJCIwMzMzMzMzMykhIzokIjFKaHZyKzhraCEjOzckJCIwNDQ0NDQ0NCohIzokIjAjWztcVW4+XyEjOjckJCIwLCwsLCwsLCIhIzkkIjIkKmYxcChSOF1VISM8NyQkIjA2NjY2NjY2IiEjOSQiMSVRVVR4J0gnSCQhIzs3JCQiMEBAQEBAQEAiISM5JCIxQkIvVXRSJlEjISM7NyQkIjBKSkpKSkpKIiEjOSQiMU8pUSxzXl9gIiEjOzckJCIwVFRUVFRUVCIhIzkkIjF5ZTRcQnZ1diEjPDckJCIwX15eXl5eXiIhIzkkIjEocCZvJ3BxKG9mISM9NyQkIjBpaGhoaGhoIiEjOSQhMVpaZ1s3ZEliISM8NyQkIjBzcnJycnJyIiEjOSQhMmopKT04JXBPeDUhIzw3JCQiMCM9PT09PT09ISM5JCEyQEMhSE5GJTReIiEjPDckJCIwIz4+Pj4+Pj4hIzkkITEjXHYhekRBXz0hIzs3JCQiMC0tLS0tLS0jISM5JCExXCNlRiV6OitAISM7NyQkIjA3Nzc3Nzc3IyEjOSQhMW9CJ0grdVNEIyEjOzckJCIwQUFBQUFBQSMhIzkkITElcHREJFxlOEIhIzs3JCQiMEtLS0tLS0sjISM5JCEyLFUlNGg4c3lBISM8NyQkIjBVVVVVVVVVIyEjOSQhMlpaMmxFd0U6IyEjPDckJCIwX19fX19fXyMhIzkkITEoKSo0P2A+XCU+ISM7NyQkIjBqaWlpaWlpIyEjOSQhMiUqKlxaKEhtK24iISM8NyQkIjB0c3Nzc3NzIyEjOSQhMWsib3BxInBYOCEjOzckJCIwJEdHR0dHR0chIzkkITBzPz9jM0AhKiohIzs3JCQiMCRISEhISEhIISM5JCEyOl5WTkZqQ0AnISM9NyQkIjAuLi4uLi4uJCEjOSQhMjo4JXA9Yj5ZRCEjPTckJCIwODg4ODg4OCQhIzkkIjA9VylvJW8zaCohIzw3JCQiMEJCQkJCQkIkISM5JCIxV14qKW9bNyg+JSEjPDckJCIwTExMTExMTCQhIzkkIjFVOXgkW1M6MighIzw3JCQiMFZWVlZWVlYkISM5JCIxSmdFS0FTOCYqISM8NyQkIjBhYGBgYGBgJCEjOSQiMllISSIqW2hvOSIhIzw3JCQiMGtqampqamokISM5JCIyKT1bNENtdipHIiEjPDckJCIwdXR0dHR0dCQhIzkkIjJPR0kkW0pLeDghIzw3JCQiMCVRUVFRUVFRISM5JCIyIz4oUWsiMyl6UyIhIzw3JCQiMCVSUlJSUlJSISM5JCIyMkcyPyhvRyJRIiEjPDckJCIwLy8vLy8vLyUhIzkkIjFoJCk+YG8wKkgiISM7NyQkIjA5OTk5OTk5JSEjOSQiMidvKm9ZLCFRbTYhIzw3JCQiMENDQ0NDQ0MlISM5JCIxQSF6QSIpWzciKiohIzw3JCQiME1NTU1NTU0lISM5JCIxPGgiKTRubkl5ISM8NyQkIjBXV1dXV1dXJSEjOSQiMndZQCh5Z2dIYiEjPTckJCIwYWFhYWFhYSUhIzkkIjJfbSwzVyYpcDYkISM9NyQkIjBsa2tra2trJSEjOSQiMVt4U0xaXm5wISM9NyQkIjB2dXV1dXV1JSEjOSQhMnZcI1wiPXZraiIhIz03JCQiMCZbW1tbW1tbISM5JCEyY1JxK0Mpeip6JCEjPTckJCIwJlxcXFxcXFwhIzkkITEjSD4+KjQtQmQhIzw3JCQiMDAwMDAwMDAmISM5JCEwKil5JHlCb1t0ISM7NyQkIjA6Ojo6Ojo6JiEjOSQhMDd6KFFpZEonKSEjOzckJCIwRERERERERCYhIzkkITFuMmtJJm8hUSYqISM8NyQkIjBOTk5OTk5OJiEjOSQhMXZwNSRcVVgrIiEjOzckJCIwWFhYWFhYWCYhIzkkITIielteSCRHVCwiISM8NyQkIjBjYmJiYmJiJiEjOSQhMGhaUClbI1sjKSohIzs3JCQiMG1sbGxsbGwmISM5JCExPCN5c2gyVTYqISM8NyQkIjB3dnZ2dnZ2JiEjOSQhMG1HUGQ4IVwhKSEjOzckJCIwJ2VlZWVlZWUhIzkkITEjelYwcTdsbychIzw3JCQiMCdmZmZmZmZmISM5JCExQCxyUHNcJzQmISM8NyQkIjAxMTExMTExJyEjOSQhMShcXShIcHViTCEjPDckJCIwOzs7Ozs7OychIzkkITJHIilbM2EqKkdhIiEjPTckJCIwRUVFRUVFRSchIzkkIjJNelNRPEd5bCMhIz43JCQiME9PT09PT08nISM5JCIyKD1sMkNBYioqPiEjPTckJCIwWVlZWVlZWSchIzkkIjJsIz1LRW1EJmYkISM9NyQkIjBjY2NjY2NjJyEjOSQiMT0tUiR6VyYpKlwhIzw3JCQiMG5tbW1tbW0nISM5JCIxRElYLmlVa2ghIzw3JCQiMHh3d3d3d3cnISM5JCIxcCZHbik+VGRxISM8NyQkIjAob29vb29vbyEjOSQiMXVwLl4jRzlsKCEjPDckJCIwKHBwcHBwcHAhIzkkIjFpOmdYWGtIeiEjPDckJCIwMjIyMjIyMighIzkkIjE4RipRanRWKXkhIzw3JCQiMDw8PDw8PDwoISM5JCIxZC1lQWpCPnYhIzw3JCQiMEZGRkZGRkYoISM5JCIwWT9GQVZPJm8hIzs3JCQiMFBQUFBQUFAoISM5JCIxeit1LWJDQWYhIzw3JCQiMFpaWlpaWlooISM5JCIydyxuZ3I1P3glISM9NyQkIjBlZGRkZGRkKCEjOSQiMDgleW8oeSZlTSEjOzckJCIwb25ubm5ubighIzkkIjElKlFPWXhVVT8hIzw3JCQiMHl4eHh4eHgoISM5JCIxKkdfJHBWJUgmZSEjPTckJCIwKXl5eXl5eXkhIzkkITFZPyUzUW8mRyYpISM9NyQkIjApenp6enp6eiEjOSQhMiUzS0ZGIyk9O0EhIz03JCQiMDMzMzMzMzMpISM5JCEyV3ZQMyV6bWFNISM9NyQkIjA9PT09PT09KSEjOSQhMjBJWy91bV5fJSEjPTckJCIwR0dHR0dHRykhIzkkITFLeTBEPDUjUiYhIzw3JCQiMFFRUVFRUVEpISM5JCExQyUzNkhoeC0nISM8NyQkIjBbW1tbW1tbKSEjOSQhMWA7JEchUlY3ayEjPDckJCIwZWVlZWVlZSkhIzkkITFgJT41KFJWTWwhIzw3JCQiMHBvb29vb28pISM5JCEwKClvPlZ3LlInISM7NyQkIjB6eXl5eXl5KSEjOSQhMWkoW2ghUiV5KWYhIzw3JCQiMCopKSkpKSkpKSkpKSkpKSEjOSQhMkNtSWNqS3NNJiEjPTckJCIwKiopKikqKSopKikqKSopISM5JCExWT9gLUhMK1ghIzw3JCQiMDQ0NDQ0NDQqISM5JCExM00/cTYlelskISM8NyQkIjA+Pj4+Pj4+KiEjOSQhMWkiUVs2WHFOIyEjPDckJCIwSEhISEhISCohIzkkITIqKXpUcSVmI3o6IiEjPTckJCIwUlJSUlJSUiohIzkkIjEqZVA3KmZgYGUhIz43JCQiMFxcXFxcXFwqISM5JCIyazwpeihwbEpDIiEjPTckJCIvJ2ZmZmZmZiohIzgkIjJCMCgqeilcTF1CISM9NyQkIi8ocHBwcHBwKiEjOCQiMShvcWYhKj0kUkwhIzw3JCQiLyl6enp6enoqISM4JCIxaikpPUZYQ3ZUISM8NyQkIi8qKikqKSopKikqKSopKiEjOCQiMW9rKGVKKyhIWyEjPDckJCIkKyIhIiIkIjFFejM0OTUiRyYhIzwtJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiM1ISIiJCIiISEiIi0lKlRISUNLTkVTU0c2IyIiJC0lJ0NVUlZFU0c2JDdgcTckJCIiISEiIiQiIiYhIiI3JCQiMCwsLCwsLCwiISM6JCIxQEJeZSF6TGcnISM7NyQkIjAtLS0tLS0tIyEjOiQiMUIlKT4/ckgrdyEjOzckJCIwLi4uLi4uLiQhIzokIjBzXy0jZkVqIikhIzo3JCQiMC8vLy8vLy8lISM6JCIxeCk0aW5HOU8pISM7NyQkIjAwMDAwMDAwJiEjOiQiMShbRzNLaiM9IykhIzs3JCQiMDExMTExMTEnISM6JCIxPDVuJVF6Q3coISM7NyQkIjAyMjIyMjIyKCEjOiQiMXlAal0yRGhxISM7NyQkIjAzMzMzMzMzKSEjOiQiMVZHdD9xJXk+JyEjOzckJCIwNDQ0NDQ0NCohIzokIjF4Y14ob0QlW18hIzs3JCQiMCwsLCwsLCwiISM5JCIxOiIpKj0hRyg+RiUhIzs3JCQiMDY2NjY2NjYiISM5JCIwIipSWlMnMzVMISM6NyQkIjBAQEBAQEBAIiEjOSQiMTB5QiU+Km8hUiMhIzs3JCQiMEpKSkpKSkoiISM5JCIyOD9xRnZJQGAiISM8NyQkIjBUVFRUVFRUIiEjOSQiMSRcSio+RzxrdSEjPDckJCIwX15eXl5eXiIhIzkkIjEnKiopbyxNJFs5JSEjPTckJCIwaWhoaGhoaCIhIzkkITB2I3lmUEJ2ZCEjOzckJCIwc3JycnJyciIhIzkkITJhWUA5Z2twNSIhIzw3JCQiMCM9PT09PT09ISM5JCEyJlxkUXZ2Wlc6ISM8NyQkIjAjPj4+Pj4+PiEjOSQhMjBCIjQlcEolKSk9ISM8NyQkIjAtLS0tLS0tIyEjOSQhMjFUQjgiPnRQQCEjPDckJCIwNzc3Nzc3NyMhIzkkITE+SkwnKlJuIkgjISM7NyQkIjBBQUFBQUFBIyEjOSQhMih6SGkiM2QpXEIhIzw3JCQiMEtLS0tLS0sjISM5JCEwMXlqQ0lCSiMhIzo3JCQiMFVVVVVVVVUjISM5JCEyQSo9c3o+WCM9IyEjPDckJCIwX19fX19fXyMhIzkkITJlUkYjXHQlKnA+ISM8NyQkIjBqaWlpaWlpIyEjOSQhMkI3cDVmMigqbyIhIzw3JCQiMHRzc3Nzc3MjISM5JCEyRyVIIz0reSZmOCEjPDckJCIwJEdHR0dHR0chIzkkITElKXpFcGFAIykqKiEjPDckJCIwJEhISEhISEghIzkkITE5UyE9USp5TWkhIzw3JCQiMC4uLi4uLi4kISM5JCEyTXg4Sk8jZThEISM9NyQkIjA4ODg4ODg4JCEjOSQiMjM+KVEwWEJXNSEjPTckJCIwQkJCQkJCQiQhIzkkIjFwdWZmdTZEViEjPDckJCIwTExMTExMTCQhIzkkIjE9Um4vWmxQcyEjPDckJCIwVlZWVlZWViQhIzkkIjB3QUIlZTU1KCohIzs3JCQiMGFgYGBgYGAkISM5JCIybCopZjYjb3hvNiEjPDckJCIwa2pqampqaiQhIzkkIjJpRm8/WWdJSiIhIzw3JCQiMHV0dHR0dHQkISM5JCIyJFtHKnlWTDZTIiEjPDckJCIwJVFRUVFRUVEhIzkkIjJ5SWNtMjE5ViIhIzw3JCQiMCVSUlJSUlJSISM5JCIxXmckUixcTVMiISM7NyQkIjAvLy8vLy8vJSEjOSQiMkciKSkqKkgjXCI+OCEjPDckJCIwOTk5OTk5OSUhIzkkIjJbbzlNIUhyJD0iISM8NyQkIjBDQ0NDQ0NDJSEjOSQiMnMpSE8nZWBeKyIhIzw3JCQiME1NTU1NTU0lISM5JCIxKUdgOUh4UyR6ISM8NyQkIjBXV1dXV1dXJSEjOSQiMkM+dD4jcCJSZiYhIz03JCQiMGFhYWFhYWElISM5JCIxWEFoZEZmVEohIzw3JCQiMGxra2tra2slISM5JCIxTlguP1Y7RG8hIz03JCQiMHZ1dXV1dXUlISM5JCEwdlh2L3F0byIhIzs3JCQiMCZbW1tbW1tbISM5JCEwZnNuJUctJSlRISM7NyQkIjAmXFxcXFxcXCEjOSQhMiU9T28kb3lpJGUhIz03JCQiMDAwMDAwMDAmISM5JCEwPCcpUk0qKWVbKCEjOzckJCIwOjo6Ojo6OiYhIzkkITA+VUxUTnF5KSEjOzckJCIwRERERERERCYhIzkkITEkKXAqcElPY3EqISM8NyQkIjBOTk5OTk5OJiEjOSQhMiRlN2RHcic9LSIhIzw3JCQiMFhYWFhYWFgmISM5JCEyJilcelp4azguIiEjPDckJCIwY2JiYmJiYiYhIzkkITFLJT0/MzspKikqKiEjPDckJCIwbWxsbGxsbCYhIzkkITEvb11dPHNsIyohIzw3JCQiMHd2dnZ2dnYmISM5JCEwSGRLI0dqIj0pISM7NyQkIjAnZWVlZWVlZSEjOSQhMTdiVHM3dSZ6JyEjPDckJCIwJ2ZmZmZmZmYhIzkkITFOLmwhKmYjKnleISM8NyQkIjAxMTExMTExJyEjOSQhMFolPURPOTRNISM7NyQkIjA7Ozs7Ozs7JyEjOSQhMScqXGAiPUlpYyIhIzw3JCQiMEVFRUVFRUUnISM5JCIyJ2U6XGMuKlFzIyEjPjckJCIwT09PT09PTychIzkkIjElUWonUjcpWy4jISM8NyQkIjBZWVlZWVlZJyEjOSQiMTpjVkohPXJsJCEjPDckJCIwY2NjY2NjYychIzkkIjEtdEFIJykqUTMmISM8NyQkIjBubW1tbW1tJyEjOSQiMUpbQUQuX3BpISM8NyQkIjB4d3d3d3d3JyEjOSQiMSpcKT1tSyN6PCghIzw3JCQiMChvb29vb29vISM5JCIxKlxhVj8vRXkoISM8NyQkIjAocHBwcHBwcCEjOSQiMTtRJT1sVmsxKSEjPDckJCIwMjIyMjIyMighIzkkIjE4bylHTjA7LSkhIzw3JCQiMDw8PDw8PDwoISM5JCIxNj8lZnU1PGwoISM8NyQkIjBGRkZGRkZGKCEjOSQiMUB5I3kpcFJ3cCEjPDckJCIwUFBQUFBQUCghIzkkIjEnPm9SdXgyLichIzw3JCQiMFpaWlpaWlooISM5JCIyWUEzWiJRW2lbISM9NyQkIjBlZGRkZGRkKCEjOSQiMG5lWC1xel8kISM7NyQkIjBvbm5ubm5uKCEjOSQiMkF2JFEiZkYnKTMjISM9NyQkIjB5eHh4eHh4KCEjOSQiMWh4WTwnRzoyJyEjPTckJCIwKXl5eXl5eXkhIzkkITF3SD8ua0tiJikhIz03JCQiMCl6enp6enp6ISM5JCExMTlabmNuVUEhIzw3JCQiMDMzMzMzMzMpISM5JCEyWWZnIjM6Uy5OISM9NyQkIjA9PT09PT09KSEjOSQhMWFuXzVCJFFmJSEjPDckJCIwR0dHR0dHRykhIzkkITFcYS89UHV4YSEjPDckJCIwUVFRUVFRUSkhIzkkITJ2SGVbLCoqbzcnISM9NyQkIjBbW1tbW1tbKSEjOSQhMFcyayY+PUBsISM7NyQkIjBlZWVlZWVlKSEjOSQhMUxpTEpZaVttISM8NyQkIjBwb29vb29vKSEjOSQhMW0heTw1c2NdJyEjPDckJCIwenl5eXl5eSkhIzkkITFpWmY/WiYpKjQnISM8NyQkIjAqKSkpKSkpKSkpKSkpKSkhIzkkITFKKnlpYjc8WCYhIzw3JCQiMCoqKSopKikqKSopKikqKSEjOSQhMSdIJFt3JHpMZiUhIzw3JCQiMDQ0NDQ0NDQqISM5JCEyRFInZl8ncGhjJCEjPTckJCIwPj4+Pj4+PiohIzkkITFiTnZTdHM8QyEjPDckJCIwSEhISEhISCohIzkkITFvej43IikzKj4iISM8NyQkIjBSUlJSUlJSKiEjOSQiMWFLeWQjKSlmIVEhIz43JCQiMFxcXFxcXFwqISM5JCIyXyc9MHpUc1Y3ISM9NyQkIi8nZmZmZmZmKiEjOCQiMkRyWTNqeTlQIyEjPTckJCIvKHBwcHBwcCohIzgkIjJ2JSlSc0NhKXpMISM9NyQkIi8penp6enp6KiEjOCQiMkV6KT0yZElMVSEjPTckJCIvKiopKikqKSopKikqKSohIzgkIjFKYjhRIVFHIVwhIzw3JCQiJCsiISIiJCIxJ2U6MHUhUW1gISM8LSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJSdDVVJWRVNHNiU3VTckJCIiISEiIiQiMmxmLSEzUU8iUiQhIzw3JCQiIiMhIiIkIjFyeXcoPlZoJG8hIzs3JCQiIiUhIiIkIjFtIjQsWFFDIm8hIzs3JCQiIichIiIkIjFZeiplK1sjKkgpISM7NyQkIiIpISIiJCIxRVAjbylwOkJpISM7NyQkIiM1ISIiJCIxOUw4QXF3Tm0hIzs3JCQiIzchIiIkIjInKj4zJT14KUdGIiEjPDckJCIjOSEiIiQiMnNHIlFJYTpgNyEjPDckJCIjOyEiIiQhMTwoKkhvUSc9NCYhIzw3JCQiIz0hIiIkITFzNSx5KkhzIj4hIzs3JCQiIz8hIiIkITEvVWEpZjJXRSMhIzs3JCQiI0EhIiIkITImKT1GXCk9JVFKIiEjPDckJCIjQyEiIiQhMj1id01yJHooPiMhIzw3JCQiI0UhIiIkITJQbCM0ZzFuazwhIzw3JCQiI0chIiIkITJoZyQqR0AlZiMzIiEjPDckJCIjSSEiIiQhMS8kcC1ecSI9TSEjPDckJCIjSyEiIiQiMWRTaz5zJ1FeJCEjPDckJCIjTSEiIiQiMWQtOSZSOnpsKSEjPDckJCIjTyEiIiQiMl9yc1NZSmdCIiEjPDckJCIjUSEiIiQiMW00eT95IVtOIiEjOzckJCIjUyEiIiQiMiY9JylRWkJXRDkhIzw3JCQiI1UhIiIkIjFVWkIvJ2UjWycpISM8NyQkIiNXISIiJCIwKj0wdlRJbGQhIzs3JCQiI1khIiIkIjE8R2oqMyo+JUgjISM8NyQkIiNbISIiJCExJj5LUyJlMV5LISM8NyQkIiNdISIiJCExIylIJXkqZWVyIyohIzw3JCQiI18hIiIkITFfbFY2NEFJKiohIzw3JCQiI2EhIiIkITJhOGFnIzNDPzUhIzw3JCQiI2MhIiIkITFsYjhBKWYqNCUpISM8NyQkIiNlISIiJCExLWpRbTJbZlwhIzw3JCQiI2chIiIkITAnNDhjIik9WUchIzs3JCQiI2khIiIkITBbY3Z0YTBrJiEjPDckJCIjayEiIiQiMHMpR1xPIlxzIyEjOzckJCIjbSEiIiQiMT5tYWAvUyFcJSEjPDckJCIjbyEiIiQiMT0peiR5PDVLKSkhIzw3JCQiI3EhIiIkIjFzWm1GSEZ0dyEjPDckJCIjcyEiIiQiMS9nYSJRdSpbYyEjPDckJCIjdSEiIiQiMUJbK1tUQG4hKSEjPDckJCIjdyEiIiQiMiRwYixxWDtZQSEjPTckJCIjeSEiIiQiMk9vbi8oR3ErPyEjPjckJCIjISkhIiIkITIwXEVNPVAuWyIhIz03JCQiIyMpISIiJCExQ15HZUwjKkhyISM8NyQkIiMlKSEiIiQhMGY6JSpvPXpyKSEjOzckJCIjJykhIiIkITEtYmo6cSIzZyUhIzw3JCQiIykpISIiJCExaVlIJW9CLlgnISM8NyQkIiMhKiEiIiQhMSQqKm83Pms1XSchIzw3JCQiIyMqISIiJCEyX0p1YEQmRzs5ISM9NyQkIiMlKiEiIiQiMTFTdT9AJWUiKiohIz43JCQiIycqISIiJCIxOCYpZilSOCJmQyEjPDckJCIjKSohIiIkIjEnZS1lbnYxRiYhIzw3JCQiJCsiISIiJCIxN1glPlFKcjomISM8LSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJSZTVFlMRUc2IyUmUE9JTlRHLSUlVklFV0c2JDskIiIhISIiJCIkKyIhIiI7JCEyKHpIaSIzZClcQiEjPCQiMXgpNGluRzlPKSEjOy0mJSZfQVhJU0c2IyIiIjYmLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSYlJl9BWElTRzYjIiIjNiYtJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIS0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJStBWEVTTEFCRUxTRzYkLUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NVEidDYiLyUnZmFtaWx5R1EhNiIvJSVzaXplR1EjMTA2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EldHJ1ZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ2l0YWxpYzYiLUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NVEieTYiLyUnZmFtaWx5R1EhNiIvJSVzaXplR1EjMTA2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EldHJ1ZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ2l0YWxpYzYiLSUpX1ZJU0lCTEVHNiMiIiItJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiIhISIiLSUpQk9VTkRTX1lHNiMkIiIhISIiLSUtQk9VTkRTX1dJRFRIRzYjJCIlK10hIiItJS5CT1VORFNfSEVJR0hURzYjJCIlK10hIiItJSlDSElMRFJFTkc2Ii0lK0FOTk9UQVRJT05HNictJSlCT1VORFNfWEc2IyQiIiEhIiItJSlCT1VORFNfWUc2IyQiIiEhIiItJS1CT1VORFNfV0lEVEhHNiMkIiUrXSEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiUrXSEiIi0lKUNISUxEUkVORzYiRzYi

Figure 4 The noisy observations, best-fit solution (green) and unperturbed solution (red)

<Text-field style="Heading 1" layout="Heading 1">Example 2: Determining Chemical Rate Constants</Text-field>

In a recent on-line forum, estimating chemical rate-reaction constants LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JlEiYUYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnRjIvRjZRJ25vcm1hbEYn and LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JlEiYkYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnRjIvRjZRJ25vcm1hbEYn for the system governed by the initial value problem

q1 := diff(x(t),t) = -a*x(t)*y(t) + b*z(t); q2 := diff(y(t),t) = -a*x(t)*y(t) + b*z(t); q3 := diff(z(t),t) = a*x(t)*y(t) - b*z(t); ic := x(0)=10,y(0)=10,z(0)=0;

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSNxMUdGKC8tSSVkaWZmR0YmNiQtSSJ4R0YoNiNJInRHRihGNSwmKihJImFHRigiIiJGMkY5LUkieUdGKEY0RjkhIiIqJkkiYkdGKEY5LUkiekdGKEY0RjlGOTcjRi4=

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSNxMkdGKC8tSSVkaWZmR0YmNiQtSSJ5R0YoNiNJInRHRihGNSwmKihJImFHRigiIiItSSJ4R0YoRjRGOUYyRjkhIiIqJkkiYkdGKEY5LUkiekdGKEY0RjlGOTcjRi4=

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSNxM0dGKC8tSSVkaWZmR0YmNiQtSSJ6R0YoNiNJInRHRihGNSwmKihJImFHRigiIiItSSJ4R0YoRjRGOS1JInlHRihGNEY5RjkqJkkiYkdGKEY5RjJGOSEiIjcjRi4=

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSNpY0dGKDYlLy1JInhHRig2IyIiISIjNS8tSSJ5R0YoRjJGNC8tSSJ6R0YoRjJGMzclRi9GNUY4

was proposed. As observed data, the 26 equispaced measurements of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictRiM2Jy1GLDYmUSJ6RicvJSdpdGFsaWNHUSV0cnVlRicvJStleGVjdXRhYmxlR1EmZmFsc2VGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUTAmQXBwbHlGdW5jdGlvbjtGJy9GO1Enbm9ybWFsRicvJSZmZW5jZUdGOS8lKnNlcGFyYXRvckdGOS8lKXN0cmV0Y2h5R0Y5LyUqc3ltbWV0cmljR0Y5LyUobGFyZ2VvcEdGOS8lLm1vdmFibGVsaW1pdHNHRjkvJSdhY2NlbnRHRjkvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZTLUkobWZlbmNlZEdGJDYlLUYjNiUtRiw2JlEidEYnRjRGN0Y6RjdGQUY3RkFGN0ZBRitGN0ZB, captured in the vector

C := <0,4.981119825867072,6.619863704639184,7.424230767960948,7.895291258366989,8.200043893374819,8.410083696006469,8.561234239001111,8.673418582102649,8.75860469142499,8.824418203860477,8.875947188583103,8.916714710357871,8.949234971834077,8.975347266499856,8.996425002439704,9.01351139331151,9.027409882146703,9.038747104014094,9.048016122443334,9.055608432249834,9.061836841504377,9.066952814916256,9.071159284138764,9.074620916958418,9.07747159947345>:

were provided. Figure 5 displays these observations as black circles.

N := Dimension(C)-1: Pts := [seq([k-1,C[k]],k=1..N+1)]: p5 := plot(Pts, style=point, symbol=circle, color=black, labels=[t,"z"]): p5;

NiwtJSdDVVJWRVNHNiM3PDckJCIiISEiIiQiIiEhIiI3JCQiIzUhIiIkIjFzcSdlIyk+NilcISM6NyQkIiM/ISIiJCIxJT1SWXFqKT5tISM6NyQkIiNJISIiJCIxWzQnendJVVUoISM6NyQkIiNTISIiJCIxKilwT2U3SCYqeSEjOjckJCIjXSEiIiQiMT5bUCQqUS8rIykhIzo3JCQiI2chIiIkIjFwaysncCQzNSUpISM6NyQkIiNxISIiJCIxNjYrUlVCaCYpISM6NyQkIiMhKSEiIiQiMVxFNSNlPU1uKSEjOjckJCIjISohIiIkIjAqXFUicC8nZSgpISM5NyQkIiQrIiEiIiQiMXgvJ1E/PVcjKSkhIzo3JCQiJDUiISIiJCIxLkplKT1aZigpKSEjOjckJCIkPyIhIiIkIjFyeU41WnI7KikhIzo3JCQiJEkiISIiJCIxeFMkPShcQlwqKSEjOjckJCIkUyIhIiIkIjFjKSpcbXNNdiopISM6NyQkIiRdIiEiIiQiMS8oUkMrRGsqKikhIzo3JCQiJGciISIiJCIwXjZMUjZOLCohIzk3JCQiJHEiISIiJCIxLm45IykpNHUtKiEjOjckJCIkIT0hIiIkIjElNDlTNVooUSEqISM6NyQkIiQhPiEiIiQiMU1MV0FoLFshKiEjOjckJCIkKyMhIiIkIjFNKVxBVjNjMCohIzo3JCQiJDUjISIiJCIxeFZdVG8kPTEqISM6NyQkIiQ/IyEiIiQiMWNpIlwiRyZwMSohIzo3JCQiJEkjISIiJCIxayhRVEdmNjIqISM6NyQkIiRTIyEiIiQiMT0lZXAiNGl1ISohIzo3JCQiJF0jISIiJCIwWHQlKmZydTIqISM5LSYlJl9BWElTRzYjIiIiNiYtJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIS0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJiUmX0FYSVNHNiMiIiM2Ji0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqTElORVNUWUxFRzYjIiIhLSUqVEhJQ0tORVNTRzYjIiIhLSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0lK0FYRVNMQUJFTFNHNiQtSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1USJ0NiIvJSdmYW1pbHlHUSE2Ii8lJXNpemVHUSMxMDYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiJRIno2Ii0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUnU1lNQk9MRzYjJSdDSVJDTEVHLSUmU1RZTEVHNiMlJlBPSU5URy0lKV9WSVNJQkxFRzYjIiIiLSUlUk9PVEc2Jy0lKUJPVU5EU19YRzYjJCIiISEiIi0lKUJPVU5EU19ZRzYjJCIiISEiIi0lLUJPVU5EU19XSURUSEc2IyQiJStdISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJStdISIiLSUpQ0hJTERSRU5HNiItJStBTk5PVEFUSU9ORzYnLSUpQk9VTkRTX1hHNiMkIiIhISIiLSUpQk9VTkRTX1lHNiMkIiIhISIiLSUtQk9VTkRTX1dJRFRIRzYjJCIlK10hIiItJS5CT1VORFNfSEVJR0hURzYjJCIlK10hIiItJSlDSElMRFJFTkc2Ig==Ig==

Figure 5 Equispaced observations of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictRiM2Jy1GLDYmUSJ6RicvJSdpdGFsaWNHUSV0cnVlRicvJStleGVjdXRhYmxlR1EmZmFsc2VGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUTAmQXBwbHlGdW5jdGlvbjtGJy9GO1Enbm9ybWFsRicvJSZmZW5jZUdGOS8lKnNlcGFyYXRvckdGOS8lKXN0cmV0Y2h5R0Y5LyUqc3ltbWV0cmljR0Y5LyUobGFyZ2VvcEdGOS8lLm1vdmFibGVsaW1pdHNHRjkvJSdhY2NlbnRHRjkvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZTLUkobWZlbmNlZEdGJDYlLUYjNiUtRiw2JlEidEYnRjRGN0Y6RjdGQUY3RkFGN0ZBRitGN0ZB

Dr. Allan Wittkopf of Maplesoft showed that because of implied conservation laws, the given system can be reduced to the initial value problem

DE1 := diff(z(t),t)=a*(10-z(t))^2-b*z(t); IC := z(0) = 0;

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSRERTFHRigvLUklZGlmZkdGJjYkLUkiekdGKDYjSSJ0R0YoRjUsJiomSSJhR0YoIiIiLCYiIzVGOUYyISIiIiIjRjkqJkkiYkdGKEY5RjJGOUY8NyNGLg==

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSNJQ0dGKC8tSSJ6R0YoNiMiIiFGMjcjRi4=

In fact, he showed that, in general, the conservation laws

C1 := x(t)+z(t)=x[0]+z[0]; C2 := y(t)+z(t)=y[0]+z[0];

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSNDMUdGKC8sJi1JInhHRig2I0kidEdGKCIiIi1JInpHRihGMkY0LCYmRjE2IyIiIUY0JkY2RjlGNDcjRi4=

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSNDMkdGKC8sJi1JInlHRig2I0kidEdGKCIiIi1JInpHRihGMkY0LCYmRjE2IyIiIUY0JkY2RjlGNDcjRi4=

lead to the single equation

temp := rifsimp([q1,q2,q3,C1,C2],[x,y,z]): map(s -> collect(s,{a,b},factor), temp[Solved][1]);

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyMvLUklZGlmZkdGJjYkLUkiekdGKDYjSSJ0R0YoRjMsJiooLChGMCIiIiZGMTYjIiIhISIiJkkieUdGKEY5RjtGNywoRjBGN0Y4RjsmSSJ4R0YoRjlGO0Y3SSJhR0YoRjdGNyomSSJiR0YoRjdGMEY3RjtGKw==

Hence, the exact solution for LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictRiM2Jy1GLDYmUSJ6RicvJSdpdGFsaWNHUSV0cnVlRicvJStleGVjdXRhYmxlR1EmZmFsc2VGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUTAmQXBwbHlGdW5jdGlvbjtGJy9GO1Enbm9ybWFsRicvJSZmZW5jZUdGOS8lKnNlcGFyYXRvckdGOS8lKXN0cmV0Y2h5R0Y5LyUqc3ltbWV0cmljR0Y5LyUobGFyZ2VvcEdGOS8lLm1vdmFibGVsaW1pdHNHRjkvJSdhY2NlbnRHRjkvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZTLUkobWZlbmNlZEdGJDYlLUYjNiUtRiw2JlEidEYnRjRGN0Y6RjdGQUY3RkFGN0ZBRitGN0ZB is

Zt := rhs(dsolve({DE1,IC},z(t)));

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSNadEdGKCwkKiYsKComLUkldGFuaEdGJTYjLCYqJkkidEdGKCIiIiomSSJiR0YoRjgsJkkiYUdGKCIjU0Y6RjhGOCNGOCIiI0Y+LUkoYXJjdGFuaEdGJTYjKiYsJkY8IiM/RjpGOEY4RjkjISIiRj9GOEY4RjlGPkY4RjwhIz9GOkZHRjhGPEZHRkY3I0Yu

Unfortunately, for positive a and b, the argument of the arctanh function is greater than 1, and hence, the value of the function is complex. This leads to complications with the fitting command next employed, so we modify this solution as follows.

Ztr := simplify(evalc(Zt)) assuming a>0, b>0;

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSRadHJHRigsJCooLCYqJiwmSSJhR0YoIiM/SSJiR0YoIiIiRjYtSSVzaW5oR0YlNiMsKCooSSJ0R0YoRjZGNSNGNiIiIywmRjMiI1NGNUY2Rj1GPS1JI2xuR0YlNiMsKEYzRjRGNUY2KiZGNUY9Rj9GPUY2Rj0tRkI2IywoRjNGNEY1RjZGRSEiIiNGSUY+RjZGNiooLUklY29zaEdGJUY5RjZGNUY9Rj9GPUZJRjZGM0ZJRjdGSUY9NyNGLg==

This function can be fit to the given observations via the NonlinearFit command from the Statistics package. In particular, it provides

Params := NonlinearFit(Ztr,<seq(k,k=0..N)>, C, [t], parameternames=[a,b], output=parametervalues, parameterranges=[a=0..infinity, b=0..infinity]): A := eval(a, Params): B := eval(b, Params): a = A; b = B;

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyMvSSJhR0YoJCIwckxyMysrKyIhIzo3Iy9GLV0zRkI5OTk5OUJGMDYwRDYx

LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyMvSSJiR0YoJCIwMi0kKWY0NDQqISM8NyMvRi1dM0Y4MjlFNDEzQjUzMjYwNw==

as estimates of the rate constants LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JlEiYUYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnRjIvRjZRJ25vcm1hbEYn and LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JlEiYkYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnRjIvRjZRJ25vcm1hbEYn. Using these values, we obtain Figure 6, a graph of the numeric solution of the original IVP, along with the observed values of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictRiM2Jy1GLDYmUSJ6RicvJSdpdGFsaWNHUSV0cnVlRicvJStleGVjdXRhYmxlR1EmZmFsc2VGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUTAmQXBwbHlGdW5jdGlvbjtGJy9GO1Enbm9ybWFsRicvJSZmZW5jZUdGOS8lKnNlcGFyYXRvckdGOS8lKXN0cmV0Y2h5R0Y5LyUqc3ltbWV0cmljR0Y5LyUobGFyZ2VvcEdGOS8lLm1vdmFibGVsaW1pdHNHRjkvJSdhY2NlbnRHRjkvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZTLUkobWZlbmNlZEdGJDYlLUYjNiUtRiw2JlEidEYnRjRGN0Y6RjdGQUY3RkFGN0ZBRitGN0ZB.

M := dsolve(eval({q1,q2,q3, ic},[a=.1,b=.00909090090909]), {x(t),y(t),z(t)},numeric): p6 := odeplot(M,[t,z(t)],0..25): display([p5,p6]);

NiotJSdDVVJWRVNHNiY3PDckJCIiISEiIiQiIiEhIiI3JCQiIzUhIiIkIjFzcSdlIyk+NilcISM6NyQkIiM/ISIiJCIxJT1SWXFqKT5tISM6NyQkIiNJISIiJCIxWzQnendJVVUoISM6NyQkIiNTISIiJCIxKilwT2U3SCYqeSEjOjckJCIjXSEiIiQiMT5bUCQqUS8rIykhIzo3JCQiI2chIiIkIjFwaysncCQzNSUpISM6NyQkIiNxISIiJCIxNjYrUlVCaCYpISM6NyQkIiMhKSEiIiQiMVxFNSNlPU1uKSEjOjckJCIjISohIiIkIjAqXFUicC8nZSgpISM5NyQkIiQrIiEiIiQiMXgvJ1E/PVcjKSkhIzo3JCQiJDUiISIiJCIxLkplKT1aZigpKSEjOjckJCIkPyIhIiIkIjFyeU41WnI7KikhIzo3JCQiJEkiISIiJCIxeFMkPShcQlwqKSEjOjckJCIkUyIhIiIkIjFjKSpcbXNNdiopISM6NyQkIiRdIiEiIiQiMS8oUkMrRGsqKikhIzo3JCQiJGciISIiJCIwXjZMUjZOLCohIzk3JCQiJHEiISIiJCIxLm45IykpNHUtKiEjOjckJCIkIT0hIiIkIjElNDlTNVooUSEqISM6NyQkIiQhPiEiIiQiMU1MV0FoLFshKiEjOjckJCIkKyMhIiIkIjFNKVxBVjNjMCohIzo3JCQiJDUjISIiJCIxeFZdVG8kPTEqISM6NyQkIiQ/IyEiIiQiMWNpIlwiRyZwMSohIzo3JCQiJEkjISIiJCIxayhRVEdmNjIqISM6NyQkIiRTIyEiIiQiMT0lZXAiNGl1ISohIzo3JCQiJF0jISIiJCIwWHQlKmZydTIqISM5LSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSdTWU1CT0xHNiMlJ0NJUkNMRUctJSZTVFlMRUc2IyUmUE9JTlRHLSUnQ1VSVkVTRzYkN1Q3JCQiIiEhIiIkIiIhISIiNyQkIi9sSzszLy1eISM5JCIxO3I1OW9WckwhIzo3JCQiL2BFaiIzLy0iISM4JCIxV25QXWYzSl0hIzo3JCQiLikqW0M3MWAiISM3JCIyOnErWmppXCwnISM7NyQkIi8xYEVqIjMvIyEjOCQiMTYoKT5mO1ZqbSEjOjckJCIvSzszLy1eRCEjOCQiMVRQMVxPPEByISM6NyQkIi9meipbQzcxJCEjOCQiMVROPU02OWd1ISM6NyQkIi8nRzlkRzlkJCEjOCQiMTN2JlFvaSxzKCEjOjckJCIvNzFgRWoiMyUhIzgkIjJOKUhDRypwXSN6ISM7NyQkIi9RcE1uJD1mJSEjOCQiMTNLRWgoKSoqKjMpISM6NyQkIi9sSzszLy1eISM4JCIwUHFiSERdQSkhIzk3JCQiLyNmeipbQzdjISM4JCIwOnFzJHo2UCQpISM5NyQkIi89ZnoqW0M3JyEjOCQiMWBMbyZHXTdWKSEjOjckJCIvV0FoSWxLbSEjOCQiMV4qb11bbjVeKSEjOjckJCIvciZHOWRHOSghIzgkIjF4JkhMaC4keiYpISM6NyQkIi8pKltDNzFgdyEjOCQiMVduJyoqKilbIVEnKSEjOjckJCIvQzcxYEVqIikhIzgkIjElNEhZPEcqKW8pISM6NyQkIi5ieFFwTW4pISM3JCIxLiY+dnlES3QpISM6NyQkIi94UXBNbiQ9KiEjOCQiMVVmV2Ehbz54KSEjOjckJCIvLy1edihRcCohIzgkIjBySXNiJylmISkpISM5NyQkIi9gRWoiMy8tIiEjNyQiMV81SGI5J2YkKSkhIzo3JCQiLydHOWRHOTIiISM3JCIxJT5ob1xhQycpKSEjOjckJCIvPWZ6KltDNyIhIzckIjE9didHdUxmKSkpISM6NyQkIi9edihRcE08IiEjNyQiMURAKlI4I3oxKikhIzo3JCQiLyU9ZnoqW0M3ISM3JCIxRG5bUTtPRCopISM6NyQkIi87My8tXnY3ISM3JCIxO3JXcl4jPiUqKSEjOjckJCIvXEM3MWBFOCEjNyQiMEJOXjJDbiYqKSEjOTckJCIvIzMvLV52UCIhIzckIjFyWmgoeW0qcCopISM6NyQkIi85ZEc5ZEc5ISM3JCIxUXJTJ29LPSkqKSEjOjckJCIvWnRPPWZ6OSEjNyQiMCQqW2ImelojKiopISM5NyQkIi4pKltDNzFgIiEjNiQiMXQ2KyVvUT8rKiEjOjckJCIvNzFgRWoiZSIhIzckIjFGITQmeU1qNSEqISM6NyQkIi9YQWhJbEs7ISM3JCIxVSFlKXpwTz0hKiEjOjckJCIveFFwTW4kbyIhIzckIjFqOyQqRzRMRCEqISM6NyQkIi5edihRcE08ISM2JCIxdUY9S2lnSiEqISM6NyQkIi9WciZHOWR5IiEjNyQiMFpASCRbRVAhKiEjOTckJCIvdihRcE1uJD0hIzckIjF0IikzJT5xQi8qISM6NyQkIi8zLy1edigpPSEjNyQiMSxzQjEoeXAvKiEjOjckJCIvVD81YnhRPiEjNyQiMW52S2AzOV4hKiEjOjckJCIvdE89ZnoqKT4hIzckIi8oNHRWLFwwKiEjODckJCIvMWBFaiIzLyMhIzckIjE+UlBwLkllISohIzo3JCQiL1JwTW4kPTQjISM3JCIxJm9cJDNPUGghKiEjOjckJCIvciZHOWRHOSMhIzckIjExa3U4TDprISohIzo3JCQiLy8tXnYoUT4jISM3JCIxeCEpenUhb20xKiEjOjckJCIvUD1meipbQyMhIzckIjE3SW9gUCUqbyEqISM6NyQkIi9wTW4kPWZIIyEjNyQiMVgjZj1uLjUyKiEjOjckJCIvLV52KFFwTSMhIzckIjEpZShHUCdvRzIqISM6NyQkIi9ObiQ9ZnpSIyEjNyQiMCN6WUR0YnUhKiEjOTckJCIvbiQ9ZnoqW0MhIzckIjAoKSoqXHInM3chKiEjOTckJCIkXSMhIiIkIjFgeiQqUUBaeCEqISM6LSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJiUmX0FYSVNHNiMiIiI2Ji0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqTElORVNUWUxFRzYjIiIhLSUqVEhJQ0tORVNTRzYjIiIhLSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0mJSZfQVhJU0c2IyIiIzYmLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUrQVhFU0xBQkVMU0c2JS1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRInQ2Ii8lJ2ZhbWlseUdRKERFRkFVTFQ2Ii8lJXNpemVHUSMxMjYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiJRIno2Ii0lJUZPTlRHNiUlKERFRkFVTFRHJShERUZBVUxURyIjNy0lKV9WSVNJQkxFRzYjIiIiLSUlUk9PVEc2Jy0lKUJPVU5EU19YRzYjJCIiISEiIi0lKUJPVU5EU19ZRzYjJCIiISEiIi0lLUJPVU5EU19XSURUSEc2IyQiJStdISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJStdISIiLSUpQ0hJTERSRU5HNiItJStBTk5PVEFUSU9ORzYnLSUpQk9VTkRTX1hHNiMkIiIhISIiLSUpQk9VTkRTX1lHNiMkIiIhISIiLSUtQk9VTkRTX1dJRFRIRzYjJCIlK10hIiItJS5CT1VORFNfSEVJR0hURzYjJCIlK10hIiItJSlDSElMRFJFTkc2Ig==Ig==

Figure 6 Observations and numeric solution for LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictRiM2Jy1GLDYmUSJ6RicvJSdpdGFsaWNHUSV0cnVlRicvJStleGVjdXRhYmxlR1EmZmFsc2VGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUTAmQXBwbHlGdW5jdGlvbjtGJy9GO1Enbm9ybWFsRicvJSZmZW5jZUdGOS8lKnNlcGFyYXRvckdGOS8lKXN0cmV0Y2h5R0Y5LyUqc3ltbWV0cmljR0Y5LyUobGFyZ2VvcEdGOS8lLm1vdmFibGVsaW1pdHNHRjkvJSdhY2NlbnRHRjkvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZTLUkobWZlbmNlZEdGJDYlLUYjNiUtRiw2JlEidEYnRjRGN0Y6RjdGQUY3RkFGN0ZBRitGN0ZB

The parameter-estimation problem for differential equations is an important component of the modeling effort prevalent in the applied sciences. In particular, we point to the extensive work of the control-theory specialist H. T. Banks of North Carolina State University.

Parameters in a system of differential equations that can only be solved numerically can be estimated by the same technique used in Example 1. For each pair of parameteter valuesLUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictSShtZmVuY2VkR0YkNiUtRiM2Jy1GLDYmUSJhRicvJSdpdGFsaWNHUSV0cnVlRicvJStleGVjdXRhYmxlR1EmZmFsc2VGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYuUSIsRidGOi9GPlEnbm9ybWFsRicvJSZmZW5jZUdGPC8lKnNlcGFyYXRvckdGOS8lKXN0cmV0Y2h5R0Y8LyUqc3ltbWV0cmljR0Y8LyUobGFyZ2VvcEdGPC8lLm1vdmFibGVsaW1pdHNHRjwvJSdhY2NlbnRHRjwvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR1EsMC4zMzMzMzMzZW1GJy1GLDYmUSJiRidGN0Y6Rj1GOkZERjpGREYrRjpGRA==, the function

SS := proc(a,b) local F, V; if not type([a,b],['numeric','numeric']) then return 'SS'(a,b); elif a<0 or b<0 then return 1e100; end if; F := dsolve(eval({q1,q2,q3,ic},{:-a=a,:-b=b}),[x(t),y(t),z(t)],numeric, output=Array([seq(k,k=0..N)])); V := convert(Column(F[2,1],4),Vector); Norm(V-C,2); end proc:

integrates the given initial value problem numerically, then computes the LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictSSVtc3ViR0YkNiUtRiw2JlEiTEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW5HRiQ2JVEiMkYnRjgvRjxRJ25vcm1hbEYnLyUvc3Vic2NyaXB0c2hpZnRHUSIwRidGK0Y4RkI= norm of the difference between the vectors of computed and observed values of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictRiM2Jy1GLDYmUSJ6RicvJSdpdGFsaWNHUSV0cnVlRicvJStleGVjdXRhYmxlR1EmZmFsc2VGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUTAmQXBwbHlGdW5jdGlvbjtGJy9GO1Enbm9ybWFsRicvJSZmZW5jZUdGOS8lKnNlcGFyYXRvckdGOS8lKXN0cmV0Y2h5R0Y5LyUqc3ltbWV0cmljR0Y5LyUobGFyZ2VvcEdGOS8lLm1vdmFibGVsaW1pdHNHRjkvJSdhY2NlbnRHRjkvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZTLUkobWZlbmNlZEdGJDYlLUYjNiUtRiw2JlEidEYnRjRGN0Y6RjdGQUY3RkFGN0ZBRitGN0ZB. Hence, LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictRiM2Jy1GLDYmUSNTU0YnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2LVEwJkFwcGx5RnVuY3Rpb247RicvRjtRJ25vcm1hbEYnLyUmZmVuY2VHRjkvJSpzZXBhcmF0b3JHRjkvJSlzdHJldGNoeUdGOS8lKnN5bW1ldHJpY0dGOS8lKGxhcmdlb3BHRjkvJS5tb3ZhYmxlbGltaXRzR0Y5LyUnYWNjZW50R0Y5LyUnbHNwYWNlR1EmMC4wZW1GJy8lJ3JzcGFjZUdGUy1JKG1mZW5jZWRHRiQ2JS1GIzYnLUYsNiZRImFGJ0Y0RjdGOi1GPjYuUSIsRidGN0ZBRkMvRkZGNkZHRklGS0ZNRk9GUS9GVVEsMC4zMzMzMzMzZW1GJy1GLDYmUSJiRidGNEY3RjpGN0ZBRjdGQUY3RkFGK0Y3RkE=, from LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictSSVtc3VwR0YkNiUtRiw2JlEiUkYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW5HRiQ2JVEiMkYnRjgvRjxRJ25vcm1hbEYnLyUxc3VwZXJzY3JpcHRzaGlmdEdRIjBGJ0YrRjhGQg== to LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JlEiUkYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnRjIvRjZRJ25vcm1hbEYn, represents a sum of squares of deviations.

Seen in Figure 7, the surface determined by LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictRiM2Jy1GLDYmUSNTU0YnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2LVEwJkFwcGx5RnVuY3Rpb247RicvRjtRJ25vcm1hbEYnLyUmZmVuY2VHRjkvJSpzZXBhcmF0b3JHRjkvJSlzdHJldGNoeUdGOS8lKnN5bW1ldHJpY0dGOS8lKGxhcmdlb3BHRjkvJS5tb3ZhYmxlbGltaXRzR0Y5LyUnYWNjZW50R0Y5LyUnbHNwYWNlR1EmMC4wZW1GJy8lJ3JzcGFjZUdGUy1JKG1mZW5jZWRHRiQ2JS1GIzYnLUYsNiZRImFGJ0Y0RjdGOi1GPjYuUSIsRidGN0ZBRkMvRkZGNkZHRklGS0ZNRk9GUS9GVVEsMC4zMzMzMzMzZW1GJy1GLDYmUSJiRidGNEY3RjpGN0ZBRjdGQUY3RkFGK0Y3RkE= suggests the existence of at least one minimum.

plot3d(SS, .0995..0.1005, .009..0.0092, view=0..0.025, axes=box, labels=[a,b,"SS"], orientation=[-80,60]);

LSUnUExPVDNERzYuLSUlR1JJREc2JTskIiQmKiohIiUkIiUwNSEiJTskIiIqISIkJCIjIyohIiU3Ozc7JCIrLkpGMEEhIzYkIipWMmlAIyEjNSQiK0Epel5CIyEjNiQiK1w5KT5FIyEjNiQiKzw8TCdIIyEjNiQiK2s6KnlMIyEjNiQiKzsieWlRIyEjNiQiKyZRejVXIyEjNiQiK296Jz1dIyEjNiQiKyMqSEBvRCEjNiQiKiUqKm9SRSEjNSQiK0YjKSllciMhIzYkIittbVQneiMhIzYkIitlcSE0KUchIzYkIitPaSxwSCEjNiQiKihvVWdJISM1JCIrX3QlWzokISM2JCIrQjUsX0shIzYkIitmXG5eTCEjNiQiKyRbPU9YJCEjNiQiK086a2ROISM2JCIqLGpObSQhIzUkIisoMz43eCQhIzYkIitkO1khKVEhIzYkIis3bzoiKlIhIzY3OyQiK0x0Sj0/ISM2JCIrdUZfQD8hIzYkIitFMGlMPyEjNiQiK2ZpV2E/ISM2JCIrZjZ0JDMjISM2JCIrJSl5NkBAISM2JCIrSyw9bUAhIzYkIitZTlc9QSEjNiQiKzdmU3hBISM2JCIrT15iVUIhIzYkIis3VFE4QyEjNiQiK2hAUypbIyEjNiQiK2NKOXFEISM2JCIrTjI8YkUhIzYkIisuNjNXRiEjNiQiKyMqUl1PRyEjNiQiK29BNUtIISM2JCIrOzFkSUkhIzYkIitbT2pKSiEjNiQiKyVvVl1CJCEjNiQiKzAleTBNJCEjNiQiKywmUSFbTSEjNiQiKydSWHNiJCEjNiQiKjZSIW9PISM1JCIrJ0d3LXkkISM2NzskIislNF0uJT0hIzYkIit5IlxVJD0hIzYkIitAdih6JD0hIzYkIisteFlePSEjNiQiKi8sWCg9ISM1JCIrTFJzMT4hIzYkIis3Rm5aPiEjNiQiKyxLIW8qPiEjNiQiK081X2A/ISM2JCIrYCY0czYjISM2JCIrckpEKD0jISM2JCIramAwakEhIzYkIitcNzBXQiEjNiQiK0ZhckhDISM2JCIrX2ZjPkQhIzYkIis7YTs4RSEjNiQiKzMrNzVGISM2JCIrZXMyNUchIzYkIitNSXM3SCEjNiQiK24ieXgsJCEjNiQiKzNdKlw3JCEjNiQiK3BTOk1LISM2JCIrajMxWEwhIzYkIisxSWFkTSEjNiQiKyJmWjlkJCEjNjc7JCIrJnokPXU7ISM2JCIrblclb2wiISM2JCIrYSopSF07ISM2JCIrX2ptYTshIzYkIio4XylwOyEjNSQiK3lvYiZwIiEjNiQiK2pOSUo8ISM2JCIqdnhreCIhIzUkIituTVBJPSEjNiQiKzt5QiMqPSEjNiQiK28kNDgnPiEjNiQiKy95JW8uIyEjNiQiKzhXOj1AISM2JCIrN1dlL0EhIzYkIitwRmImSCMhIzYkIissYGAhUiMhIzYkIippbSEqWyMhIzUkIissaXQhZiMhIzYkIisjXCQ9JnAjISM2JCIrI28jNC1HISM2JCIrVnY9NkghIzYkIitzbEFBSSEjNiQiKmIpKlw4JCEjNSQiK3cmPSRcSyEjNiQiK2pYLWxMISM2NzskIitVQmlCOiEjNiQiK0Mhb0ZcIiEjNiQiK3hwY3Q5ISM2JCIrY29ZbTkhIzYkIiszRmpyOSEjNiQiKyIqbyMqKVsiISM2JCIqOUR6XiIhIzUkIitWWCd6YiIhIzYkIipXNiMzOyEjNSQiKzFZdG47ISM2JCIrejhkTjwhIzYkIis0JnoyIj0hIzYkIispZXNDKj0hIzYkIisnZlMpej4hIzYkIitgKGVAMiMhIzYkIiteIyl5b0AhIzYkIis+PjxwQSEjNiQiK29yI0dQIyEjNiQiKy15THpDISM2JCIrMGRNKWUjISM2JCIrbkhhKnAjISM2JCIrYlptN0chIzYkIisySVtGSCEjNiQiK3A0IVEvJCEjNiQiKz4mWzk7JCEjNjc7JCIqcGtPUiIhIzUkIit6J0hvTSIhIzYkIisjR1NASiIhIzYkIis4WmMhSCIhIzYkIitgS3YjRyIhIzYkIitUWSUqKUciISM2JCIrWilHKjM4ISM2JCIrRyl5P00iISM2JCIrTDBXKFEiISM2JCIrOksmUVciISM2JCIrIzRwK14iISM2JCIqZl1bZSIhIzUkIiphT3FtIiEjNSQiK3QqeWJ2IiEjNiQiK3RzYlw9ISM2JCIrLHU8Wz4hIzYkIipsaTIwIyEjNSQiKy0wdWNAISM2JCIqWkdjRSMhIzUkIis+MC14QiEjNiQiK0JbZCFcIyEjNiQiKjouZ2cjISM1JCIqJj0xQkYhIzUkIitAVmFURyEjNiQiK21aRmhIISM2NzskIitCdlohSCIhIzYkIit6KFxhQSIhIzYkIisyI3lBPCIhIzYkIipsP0U4IiEjNSQiKzBlIno1IiEjNiQiK1ApZSIqNCIhIzYkIitmXHIxNiEjNiQiK0dSQ0k2ISM2JCIrOClvKG82ISM2JCIrQyhlM0EiISM2JCIrJlJgW0ciISM2JCIqTGghZjghIzUkIispRycqPVciISM2JCIreilcPmAiISM2JCIrQmMsRzshIzYkIitUSzNIPCEjNiQiK1ciPVYkPSEjNiQiKzhlLlY+ISM2JCIraGpuYT8hIzYkIitAPHlvQCEjNiQiKyJHd1xHIyEjNiQiK0I0JkhTIyEjNiQiK0UsWEFEISM2JCIrRDtFVkUhIzYkIityIjNfdyMhIzY3OyQiK0M0ej83ISM2JCIrTDxXTzYhIzYkIitHRVdpNSEjNiQiK1tsMiw1ISM2JCIqJSlvd2EqISM2JCIxLCsrRiMpZmQjKiEjPSQiK2kkXG86KiEjNyQiKzAiXDlEKiEjNyQiKiMqSGFgKiEjNiQiK2wyXiMqKiohIzckIisqUT0rMSIhIzYkIitVQlRMNiEjNiQiKzopUXJAIiEjNiQiKz1kPzQ4ISM2JCIrSnUnelMiISM2JCIrQEk1NzohIzYkIis9eWM/OyEjNiQiK3ovYUs8ISM2JCIrVGBQWj0hIzYkIis9SmNrPiEjNiQiK2YpKnAkMyMhIzYkIisoKUhZL0EhIzYkIisuSWZFQiEjNiQiK2AmeilcQyEjNiQiK0gzOnVEISM2NzskIitKUVUhPiIhIzYkIipGJnAoMyIhIzUkIis3KnBiIyoqISM3JCIrLT0kUTIqISM3JCIrNC41XyQpISM3JCIrZElQJ3ooISM3JCIxKioqKioqM0pzVnUhIz0kIitPS01CdCEjNyQiK2F3SVl1ISM3JCIrOCJHNCF5ISM3JCIrQ0BmZCQpISM3JCIrVSlcIXohKiEjNyQiKy5tSEgqKiEjNyQiKzZTIXkzIiEjNiQiKzxEOyE+IiEjNiQiK0BbQSlIIiEjNiQiK14ubjU5ISM2JCIrKT5AbF8iISM2JCIrIj1cXWsiISM2JCIrd3BxbDwhIzYkIit6JnkhKSk9ISM2JCIrUFMlPSwjISM2JCIrWFZ2T0AhIzYkIit0TWhpQSEjNiQiKy1hRSpRIyEjNjc7JCIrdTlHLTchIzYkIitcaWMlMyIhIzYkIisjcHI9ciohIzckIitAK09RJykhIzckIitEKlwubCghIzckIitAYTImeSchIzckIisqKjR6JTQnISM3JCIyJyoqKioqXFozUWsmISM+JCIrR2g4IlwmISM3JCIqJCoqcGdjISM2JCIrbCFlYjcnISM3JCIrQGVVRG8hIzckIitUYkUncCghIzckIitaKlJtbykhIzckIitOI3krdyohIzckIitzOT4qMyIhIzYkIitwL2MxNyEjNiQiKyhwJCpwSyIhIzYkIisoZjwoXDkhIzYkIitNcz11OiEjNiQiKzYxLCs8ISM2JCIraSshcCM9ISM2JCIrRzBrYT4hIzYkIipab0kzIyEjNSQiKyRlZD9AIyEjNjc7JCIrdU00YjchIzYkIitzPU5GNiEjNiQiK1AyTiw1ISM2JCIraSdHI3koKSEjNyQiK3V5ZXp2ISM3JCIyJioqKioqPlFdeVYnISM+JCIqM0chKlEmISM2JCIrWmNEKVwlISM3JCIrVEdwdlEhIzckIiskZUcvbSQhIzckIitDc3c+UiEjNyQiK3NGTnRYISM3JCIqPnM7WyYhIzYkIitiNz9SbCEjNyQiMSoqKioqejsjUiVvKCEjPSQiKzs7QyQpKSkhIzckIitjWWw2NSEjNiQiK2F1SFA2ISM2JCIqJT5iazchIzUkIislKWYnSFIiISM2JCIrbmxCQTohIzYkIis9UTpfOyEjNiQiKj1vRHkiISM1JCIrLi9QOD4hIzYkIitQKnlXLyMhIzY3OyQiK146KFJNIiEjNiQiK1YjPjZAIiEjNiQiK0pya3k1ISM2JCIrPDMqcFkqISM3JCIrbEM2YiIpISM3JCIxKioqKio+OyYpXCZvISM9JCIyJioqKioqcC5xWGQmISM+JCIqYFQ1TCUhIzYkIihwdzskISIqJCIrLjMkUkAjISM3JCIqNEVBJD0hIzYkIioySUxLIyEjNiQiK2laeT5MISM3JCIreHAvKFwlISM3JCIrWUl3V2QhIzckIitDaEdEcSEjNyQiMSwrKzhASkIkKSEjPSQiKyozOjtqKiEjNyQiKzN4aiU0IiEjNiQiKS1hRTchIiokIisjNE0oZTghIzYkIiswQDgiXCIhIzYkIitLYG5COyEjNiQiKzdFS2M8ISM2JCIrTFYvKik9ISM2NzskIisyM0hpOSEjNiQiKzJsMUc4ISM2JCIreCJ6UT4iISM2JCIqekcoZjUhIzUkIisqPmBoRCohIzckIis5dVE6eiEjNyQiKydHISpcZCchIzckIitmOidcQiYhIzckIiskKTRJJipRISM3JCIrSCQzZ2IjISM3JCIpTjM8NyEjNSQiKyo0MltAIiEjOCQiK0dnbWY5ISM3JCIrTScqWyh6IyEjNyQiKzRvJVw4JSEjNyQiKyhwUD9aJiEjNyQiMiYqKioqKkhgaSgzbyEjPiQiK2w6N1giKSEjNyQiK1hdNiJbKiEjNyQiK0JWbiIzIiEjNiQiK1AxPzo3ISM2JCIrcy9wWzghIzYkIithUTkjWyIhIzYkIisyM2M6OyEjNiQiK2Y4JSpbPCEjNjc7JCIrQCoqWy47ISM2JCIrblpCcTkhIzYkIispKm8/UDghIzYkIitUI3BXPyIhIzYkIisuZDZzNSEjNiQiK2tuI0hTKiEjNyQiK1g0VCM0KSEjNyQiK2NuJlF6JyEjNyQiK2YtWDpiISM3JCIrMnQmXEYlISM3JCIrOCEqWzxKISM3JCIrbF4iKnlAISM3JCIrdHd5Sz0hIzckIis+cSdlTiMhIzckIitART1rTCEjNyQiK20iPl9hJSEjNyQiKUowJXomISM1JCIrJ3BsWTIoISM3JCIxLCsrSEVKcyQpISM9JCIxKioqKio0diUqKnonKiEjPSQiKyIqKSpSKjQiISM2JCIrayM9N0IiISM2JCIreUVLajghIzYkIisxd2kmXCIhIzYkIitQXzJHOyEjNjc7JCIrImZDP3ciISM2JCIrcmApNGoiISM2JCIrYCVSMF0iISM2JCIrT2AlM1AiISM2JCIrVmw3VTchIzYkIitZd3E5NiEjNiQiKzc3dyEqKSohIzckIis0MSQqZicpISM3JCIrNCkpKnBZKCEjNyQiK3dPPUxqISM3JCIrRnZJJ0gmISM3JCIrdU9yQ1chIzckIislKT5bS1EhIzckIislKVFeZE8hIzckIipDLF8mUiEjNiQiK2Q1L05ZISM3JCIrRmxGZWIhIzckIitwTG9CbSEjNyQiMi8rK2dLSkd4KCEjPiQiK01OV3QqKSEjNyQiKUxzPzUhIiokIiswbExZNiEjNiQiK3Q6XnQ3ISM2JCIpSiI9UyIhIiokIiszJVs0YCIhIzY3OyQiK2N5ZUw+ISM2JCIrTi1KMD0hIzYkIitQLyp6biIhIzYkIis+WSY9YiIhIzYkIisjSDFzVSIhIzYkIitIMFkvOCEjNiQiK1dsPiU9IiEjNiQiK1EoUnMxIiEjNiQiK2ZuIXphKiEjNyQiMSoqKioqXCRwRCdbKSEjPSQiK1lVJFJeKCEjNyQiK0tuKilwbSEjNyQiMiYqKioqKmYnbz8zZyEjPiQiK2B4byRmJiEjNyQiMi8rK11nQUNbJiEjPiQiK01TJz5wJiEjNyQiKzkoeiYqPSchIzckIituKkdKInAhIzckIittMSEpKnooISM3JCIrczBFKykpISM3JCIrJzRhKXopKiEjNyQiKy8xXyw2ISM2JCIrMEcxPjchIzYkIitebmFSOCEjNiQiKzRiQmk5ISM2NzskIis2IW9cNiMhIzYkIipGKWYqKT4hIzUkIitIZltsPSEjNiQiK3A1Kkd1IiEjNiQiKyRmW0BpIiEjNiQiKyIqW3AuOiEjNiQiKzdgNSlRIiEjNiQiK0RuOXc3ISM2JCIrJUddKW82ISM2JCIrKW8zdzEiISM2JCIrJ3pnSHUqISM3JCIrJypSKlEiKikhIzckIis7P0Q/IykhIzckIisqW04mKXAoISM3JCIrTjgyJlEoISM3JCIrJEg1bEkoISM3JCIrWTAycXUhIzckIitBaVlneSEjNyQiK0QnR2hXKSEjNyQiKkRUJyo9KiEjNiQiK28xZjA1ISM2JCIrbXJlLDYhIzYkIisicC5aPyIhIzYkIitaSUQ4OCEjNiQiKzBcKWZVIiEjNjc7JCIrJG8yUUkjISM2JCIqZSRHIj0jISM1JCIrPDNDZz8hIzYkIitvKVs0JT4hIzYkIitkSnVCPSEjNiQiK1xVLzQ8ISM2JCIrQkhRKGYiISM2JCIrUERWKlsiISM2JCIrPihcZ1EiISM2JCIrQElLKUciISM2JCIrLWJpKD4iISM2JCIrSCZmYzYiISM2JCIrSCdvVy8iISM2JCIrL2t0aikqISM3JCIrSFQheVYqISM3JCIrVkJKKT0qISM3JCIqancmSCIqISM2JCIrMXMxbCMqISM3JCIrWyYqUidlKiEjNyQiK3pqYzI1ISM2JCIrKkh0NDIiISM2JCIrIip5V1k2ISM2JCIrVD93SjchIzYkIismXDBdSyIhIzYkIitVPmhDOSEjNjc7JCIrZSN6JClcIyEjNiQiK1spUiZ5QiEjNiQiK14pUS5FIyEjNiQiKywuL1dAISM2JCIqJ1EnKkg/ISM1JCIrKD4pXD0+ISM2JCIrcnU2NT0hIzYkIitIN1MwPCEjNiQiKzx3MDA7ISM2JCIrPSlbKjQ6ISM2JCIrSmU2QDkhIzYkIis9bHpSOCEjNiQiKy9qVW43ISM2JCIrOGNnMDchIzYkIis7YSxjNiEjNiQiK1h6RT82ISM2JCIrJGYrKCo0IiEjNiQiK19iOiY0IiEjNiQiKyIpcCJvNSIhIzYkIiskPnNUOCIhIzYkIis7VTZ3NiEjNiQiKyhIVDZCIiEjNiQiK2smeXZIIiEjNiQiKik0d3Q4ISM1JCIrWj46ZTkhIzY3OyQiKy8oM3VwIyEjNiQiKjhhK2UjISM1JCIqc1RXWSMhIzUkIit4QCMzTiMhIzYkIismPSNcUkEhIzYkIiptLjM4IyEjNSQiKyhmdV4tIyEjNiQiK3I5NUI+ISM2JCIrXkc8RD0hIzYkIisqRyQzSzwhIzYkIismb1hZayIhIzYkIitMKil6ajohIzYkIitEZmchXCIhIzYkIitcXkJFOSEjNiQiKkQ9PlAiISM1JCIrZHYoKUc4ISM2JCIrb1BBKUgiISM2JCIrJTNOM0ciISM2JCIrRD5DeDchIzYkIitEJVt2RyIhIzYkIitsZFQ2OCEjNiQiK2ZJNls4ISM2JCIrKCkpPW1SIiEjNiQiK09MdWI5ISM2JCIrKnBUVV8iISM2NzskIiskelkqKipHISM2JCIrQCVwW3kjISM2JCIrJlsoZnJFISM2JCIrKnlrLmMjISM2JCIrX0dXXkMhIzYkIitCdDlYQiEjNiQiKyc9WD1DIyEjNiQiK1ZFJz45IyEjNiQiK0JMKmYvIyEjNiQiK2poXWE+ISM2JCIrZT06bz0hIzYkIitGaG0oeSIhIzYkIitVdCdRciIhIzYkIitIYGtaOyEjNiQiK015JCoqZSIhIzYkIitQPHBUOiEjNiQiK3ciM1FdIiEjNiQiK2BaMng5ISM2JCIrJz0kNGk5ISM2JCIqWDojZjkhIzUkIitcRl1vOSEjNiQiK3l4cipbIiEjNiQiK3lCTkE6ISM2JCIrTUxvbDohIzYkIit1SCUpPTshIzY3OyQiKyNveF81JCEjNiQiK2t1RiMqSCEjNiQiK0dnNiIpRyEjNiQiKkA0P3gjISM1JCIraFA/bEUhIzYkIis8PSk0YyMhIzYkIip3bCdmQyEjNSQiK0VMaWhCISM2JCIrKW9zc0UjISM2JCIrdnEzeEAhIzYkIitWIylmIjQjISM2JCIpeVI2PyEiKiQiK3ZdOFA+ISM2JCIrMSo0JnA9ISM2JCIrI2VlIzQ9ISM2JCIqVkpydiIhIzUkIitEOCdRciIhIzYkIislKVI3IW8iISM2JCIrNkpcYzshIzYkIitDZVJWOyEjNiQiK09hMlQ7ISM2JCIrcCZvJlw7ISM2JCIrI1EsKG87ISM2JCIrJ3AwIilwIiEjNiQiK2U7RFA8ISM2NzskIitKWSZHSiQhIzYkIisiKWZ1LEshIzYkIipwKFsjNCQhIzUkIitGekYmKUghIzYkIipGUy4pRyEjNSQiK2hvI3p4IyEjNiQiK0Q8S3lFISM2JCIrXlQlPWUjISM2JCIrNj4mKSlbIyEjNiQiK1VTdSpSIyEjNiQiK1lGJ1xKIyEjNiQiK0ZRKlxCIyEjNiQiK0ZcT2dAISM2JCIrZDFrIjQjISM2JCIrN1BUSD8hIzYkIiskNCJIdT4hIzYkIit0XChvIz4hIzYkIit5JFF4KT0hIzYkIithb1JkPSEjNiQiK1wleWkkPSEjNiQiKzZhcEM9ISM2JCIrNjojRyM9ISM2JCIrJmZ5MSQ9ISM2JCIrdlk4Wz0hIzYkIispKlEiXCg9ISM2NzskIitkWkRBTiEjNiQiK20qb0dUJCEjNiQiK3AlSGBJJCEjNiQiK2R4IikqPiQhIzYkIiswamAnNCQhIzYkIitDM3ImKkghIzYkIitvRGYoKkchIzYkIisqZmdDIUchIzYkIiopUWk1RiEjNSQiK3pIVUFFISM2JCIrODZCUUQhIzYkIio9YSVlQyEjNSQiK0ckSE5RIyEjNiQiK2k6I1JKIyEjNiQiK1J5Nl1BISM2JCIqPj1FPiMhIzUkIiswTiM+OSMhIzYkIis8K18pNCMhIzYkIisoKjMnRzEjISM2JCIrS2dNTj8hIzYkIis2O0k7PyEjNiQiKzs7J2YrIyEjNiQiK3lNWC8/ISM2JCIrOykqeTY/ISM2JCIra3UneS0jISM2NzskIit3cDlMUCEjNiQiK2RCTERPISM2JCIraylbJD5OISM2JCIrJ2ZpYFQkISM2JCIrKTNkTkokISM2JCIrTlw4OUshIzYkIis2JT50NiQhIzYkIitEZ05CSSEjNiQiKyIqUV5LSCEjNiQiKyV6JzNYRyEjNiQiK2NQUmhGISM2JCIrIykpeTxvIyEjNiQiKzstaDFFISM2JCIrJEh4aWAjISM2JCIreXE9ckMhIzYkIismKXp2NkMhIzYkIisieTYlZUIhIzYkIislSGo6SiMhIzYkIio3MztGIyEjNSQiKyh6MypRQSEjNiQiK2wteThBISM2JCIrRGRaJz4jISM2JCIrPlU8KD0jISM2JCIrVTcoZj0jISM2JCIrJlt0Rz4jISM2LSUlVklFV0c2JSUoREVGQVVMVEclKERFRkFVTFRHOyQiIiEhIiIkIiNEISIkLSYlJl9BWElTRzYjIiIiNiYtJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIS0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJiUmX0FYSVNHNiMiIiM2Ji0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqTElORVNUWUxFRzYjIiIhLSUqVEhJQ0tORVNTRzYjIiIhLSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0mJSZfQVhJU0c2IyIiJDYmLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUrQVhFU0xBQkVMU0c2JS1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRImE2Ii8lJ2ZhbWlseUdRITYiLyUlc2l6ZUdRIzEwNiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2Ii1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRImI2Ii8lJ2ZhbWlseUdRITYiLyUlc2l6ZUdRIzEwNiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2IlEjU1M2Ii0lK0xJR0hUTU9ERUxHNiNRKExJR0hUXzM2Ii0lK1BST0pFQ1RJT05HNiwkIig0Q28pISIpJCIoeCFbKSohIigkIil0JFFdIiEiKSQhKSdRUyNcISIpJCIpPVtPPCEiKSQhKCZvRyYpISIoJCEoYS1tKSEiKCQiIiEhIiIkIiImISIiJCIjNSEiIi0lKkFYRVNTVFlMRUc2IyUkQk9YRy0lKV9WSVNJQkxFRzYjIiIiLSUlUk9PVEc2Jy0lKUJPVU5EU19YRzYjJCIiISEiIi0lKUJPVU5EU19ZRzYjJCIiISEiIi0lLUJPVU5EU19XSURUSEc2IyQiJStxISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJStxISIiLSUpQ0hJTERSRU5HNiItJStBTk5PVEFUSU9ORzYnLSUpQk9VTkRTX1hHNiMkIiIhISIiLSUpQk9VTkRTX1lHNiMkIiIhISIiLSUtQk9VTkRTX1dJRFRIRzYjJCIlK3EhIiItJS5CT1VORFNfSEVJR0hURzYjJCIlK3EhIiItJSlDSElMRFJFTkc2Ig==Ig==

Figure 7 Surface determined by the function LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictRiM2Jy1GLDYmUSNTU0YnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2LVEwJkFwcGx5RnVuY3Rpb247RicvRjtRJ25vcm1hbEYnLyUmZmVuY2VHRjkvJSpzZXBhcmF0b3JHRjkvJSlzdHJldGNoeUdGOS8lKnN5bW1ldHJpY0dGOS8lKGxhcmdlb3BHRjkvJS5tb3ZhYmxlbGltaXRzR0Y5LyUnYWNjZW50R0Y5LyUnbHNwYWNlR1EmMC4wZW1GJy8lJ3JzcGFjZUdGUy1JKG1mZW5jZWRHRiQ2JS1GIzYnLUYsNiZRImFGJ0Y0RjdGOi1GPjYuUSIsRidGN0ZBRkMvRkZGNkZHRklGS0ZNRk9GUS9GVVEsMC4zMzMzMzMzZW1GJy1GLDYmUSJiRidGNEY3RjpGN0ZBRjdGQUY3RkFGK0Y3RkE=

To find the minimum shown in Figure 7, we again invoke the Nelder-Meade (simplex) algorithm via the NLPSolve command, and obtain

params := NLPSolve(SS(a,b), method=nonlinearsimplex, initialpoint=[a=.2, b=.1],evaluationlimit=200): An := eval(a, params[2]): Bn := eval(b, params[2]): a = An; b = Bn;

L0kiYUc2IiQiMyVIIio9eEIrKysiISM9

L0kiYkc2IiQiMyFIWiEzYVwiNDQqISM/

as estimates for LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JlEiYUYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnRjIvRjZRJ25vcm1hbEYn and LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JlEiYkYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnRjIvRjZRJ25vcm1hbEYn, with

params[1];

JCIzIyoqZnkleSVmPHkmISND

as the corresponding minimum value of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JlEjU1NGJy8lJ2l0YWxpY0dRJXRydWVGJy8lK2V4ZWN1dGFibGVHUSZmYWxzZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJ0YyL0Y2USdub3JtYWxGJw==. By way of comparison, the values of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JlEiYUYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnRjIvRjZRJ25vcm1hbEYn and LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JlEiYkYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnRjIvRjZRJ25vcm1hbEYn computed from the exact solution for LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYnLUkjbWlHRiQ2I1EhRictRiM2Jy1GLDYmUSJ6RicvJSdpdGFsaWNHUSV0cnVlRicvJStleGVjdXRhYmxlR1EmZmFsc2VGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUTAmQXBwbHlGdW5jdGlvbjtGJy9GO1Enbm9ybWFsRicvJSZmZW5jZUdGOS8lKnNlcGFyYXRvckdGOS8lKXN0cmV0Y2h5R0Y5LyUqc3ltbWV0cmljR0Y5LyUobGFyZ2VvcEdGOS8lLm1vdmFibGVsaW1pdHNHRjkvJSdhY2NlbnRHRjkvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZTLUkobWZlbmNlZEdGJDYlLUYjNiUtRiw2JlEidEYnRjRGN0Y6RjdGQUY3RkFGN0ZBRitGN0ZB differ from these new values by

An-A; Bn-B;

JCIkXiIhIzU=

JCIlZmAhIzc=

respectively, and the minimum of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JlEjU1NGJy8lJ2l0YWxpY0dRJXRydWVGJy8lK2V4ZWN1dGFibGVHUSZmYWxzZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJ0YyL0Y2USdub3JtYWxGJw== at the former point is

SS(A,B);

JCIzbl5pLCd5UUdnKiEjQw==

a number marginally higher than that determined by our numeric method.