Regression Analysis Daniel Skoog, Maplesoft Linear regression is one of the fundamental approaches for determining relationships between dependent and exploratory variables.This worksheet covers several introductory examples for regression analysis in Maple. The primary focus of the worksheet is the use of the Statistics:-Fit command to fit a model to data as well as generating reports on the fitted model. Additional examples include fitting a curve using the Optimization package.

The Statistics and CurveFitting packages contain a number of commands for fitting curves to data points including linear, exponential, polynomial, least squares, and spline. Statistics - Regression Curve Fitting ExponentialFitFitLinearFitLogarithmicFitLowessNonlinearFitOneWayANOVAPolynomialFitPowerFitPredictiveLeastSquaresRepeatedMedianEstimatorArrayInterpolationBSplineBSplineCurveInteractiveLeastSquaresLowessPolynomialInterpolationRationalInterpolationSplineThieleInterpolation JSFH

The Statistics:-Fit command is a general purpose command that allows for specification of linear or non-linear models.Syntax:Fit( Model, X::Independent Variable, Y::Dependent Variable, x::Name of independent variable in component function, options );The Fit command fits a model by minimizing the least-squares error.

with(Statistics): Start by specifying some data for an independent variable:X := Vector([1, 2, 3, 4, 5, 6], datatype = float);

LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEiWEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RKn4mIzg3ODg7fkYnL0YzUSdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGPS8lKXN0cmV0Y2h5R0Y9LyUqc3ltbWV0cmljR0Y9LyUobGFyZ2VvcEdGPS8lLm1vdmFibGVsaW1pdHNHRj0vJSdhY2NlbnRHRj0vJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZMLUkobWFjdGlvbkdGJDYlLUkobWZlbmNlZEdGJDYoLUYjNiYtSSdtdGFibGVHRiQ2Oi1JJG10ckdGJDYmLUkkbXRkR0YkNigtSSNtbkdGJDYkUSMxLkYnRjkvJSlyb3dhbGlnbkdRIUYnLyUsY29sdW1uYWxpZ25HRmBvLyUrZ3JvdXBhbGlnbkdGYG8vJShyb3dzcGFuR1EiMUYnLyUrY29sdW1uc3BhbkdGZ29GXm9GYW9GY28tRmVuNiYtRmhuNigtRltvNiRRIzIuRidGOUZeb0Zhb0Zjb0Zlb0Zob0Zeb0Zhb0Zjby1GZW42Ji1GaG42KC1GW282JFEjMy5GJ0Y5Rl5vRmFvRmNvRmVvRmhvRl5vRmFvRmNvLUZlbjYmLUZobjYoLUZbbzYkUSM0LkYnRjlGXm9GYW9GY29GZW9GaG9GXm9GYW9GY28tRmVuNiYtRmhuNigtRltvNiRRIzUuRidGOUZeb0Zhb0Zjb0Zlb0Zob0Zeb0Zhb0Zjby1GZW42Ji1GaG42KC1GW282JFEjNi5GJ0Y5Rl5vRmFvRmNvRmVvRmhvRl5vRmFvRmNvLyUmYWxpZ25HUSVheGlzRicvRl9vUSliYXNlbGluZUYnL0Zib1EnY2VudGVyRicvRmRvUSd8ZnJsZWZ0fGhyRicvJS9hbGlnbm1lbnRzY29wZUdGMS8lLGNvbHVtbndpZHRoR1ElYXV0b0YnLyUmd2lkdGhHRmpyLyUrcm93c3BhY2luZ0dRJjEuMGV4RicvJS5jb2x1bW5zcGFjaW5nR1EmMC44ZW1GJy8lKXJvd2xpbmVzR1Elbm9uZUYnLyUsY29sdW1ubGluZXNHRmVzLyUmZnJhbWVHRmVzLyUtZnJhbWVzcGFjaW5nR1EsMC40ZW1+MC41ZXhGJy8lKmVxdWFscm93c0dGPS8lLWVxdWFsY29sdW1uc0dGPS8lLWRpc3BsYXlzdHlsZUdGPS8lJXNpZGVHUSZyaWdodEYnLyUwbWlubGFiZWxzcGFjaW5nR0Zicy8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJSlyZWFkb25seUdGPUY5RjkvSSttc2VtYW50aWNzR0YkUSpDb2xWZWN0b3JGJy8lJW9wZW5HUSJbRicvJSZjbG9zZUdRIl1GJ0ZddS8lK2FjdGlvbnR5cGVHUS5ydGFibGVhZGRyZXNzRicvJSlydGFibGVpZEdRNTE4NDQ2NzQ0MDc0MzI5MDY0MDIyRidGaHRGW3VGOQ==

Next add data for a dependent variable:Y := Vector([2, 3, 4, 3.5, 5.8, 7], datatype = float);

LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEiWUYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RKn4mIzg3ODg7fkYnL0YzUSdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGPS8lKXN0cmV0Y2h5R0Y9LyUqc3ltbWV0cmljR0Y9LyUobGFyZ2VvcEdGPS8lLm1vdmFibGVsaW1pdHNHRj0vJSdhY2NlbnRHRj0vJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZMLUkobWFjdGlvbkdGJDYlLUkobWZlbmNlZEdGJDYoLUYjNiYtSSdtdGFibGVHRiQ2Oi1JJG10ckdGJDYmLUkkbXRkR0YkNigtSSNtbkdGJDYkUSMyLkYnRjkvJSlyb3dhbGlnbkdRIUYnLyUsY29sdW1uYWxpZ25HRmBvLyUrZ3JvdXBhbGlnbkdGYG8vJShyb3dzcGFuR1EiMUYnLyUrY29sdW1uc3BhbkdGZ29GXm9GYW9GY28tRmVuNiYtRmhuNigtRltvNiRRIzMuRidGOUZeb0Zhb0Zjb0Zlb0Zob0Zeb0Zhb0Zjby1GZW42Ji1GaG42KC1GW282JFEjNC5GJ0Y5Rl5vRmFvRmNvRmVvRmhvRl5vRmFvRmNvLUZlbjYmLUZobjYoLUZbbzYkUTEzLjUwMDAwMDAwMDAwMDAwRidGOUZeb0Zhb0Zjb0Zlb0Zob0Zeb0Zhb0Zjby1GZW42Ji1GaG42KC1GW282JFExNS44MDAwMDAwMDAwMDAwMEYnRjlGXm9GYW9GY29GZW9GaG9GXm9GYW9GY28tRmVuNiYtRmhuNigtRltvNiRRIzcuRidGOUZeb0Zhb0Zjb0Zlb0Zob0Zeb0Zhb0Zjby8lJmFsaWduR1ElYXhpc0YnL0Zfb1EpYmFzZWxpbmVGJy9GYm9RJ2NlbnRlckYnL0Zkb1EnfGZybGVmdHxockYnLyUvYWxpZ25tZW50c2NvcGVHRjEvJSxjb2x1bW53aWR0aEdRJWF1dG9GJy8lJndpZHRoR0Zqci8lK3Jvd3NwYWNpbmdHUSYxLjBleEYnLyUuY29sdW1uc3BhY2luZ0dRJjAuOGVtRicvJSlyb3dsaW5lc0dRJW5vbmVGJy8lLGNvbHVtbmxpbmVzR0Zlcy8lJmZyYW1lR0Zlcy8lLWZyYW1lc3BhY2luZ0dRLDAuNGVtfjAuNWV4RicvJSplcXVhbHJvd3NHRj0vJS1lcXVhbGNvbHVtbnNHRj0vJS1kaXNwbGF5c3R5bGVHRj0vJSVzaWRlR1EmcmlnaHRGJy8lMG1pbmxhYmVsc3BhY2luZ0dGYnMvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXUYnLyUpcmVhZG9ubHlHRj1GOUY5L0krbXNlbWFudGljc0dGJFEqQ29sVmVjdG9yRicvJSVvcGVuR1EiW0YnLyUmY2xvc2VHUSJdRidGXXUvJSthY3Rpb250eXBlR1EucnRhYmxlYWRkcmVzc0YnLyUpcnRhYmxlaWRHUTUxODQ0Njc0NDA3NDMyOTA2NTEwMkYnRmh0Rlt1Rjk=

Using the Fit command, you can fit a linear model to data and display summary of the regression model.Fit a least-squares model to the data:Fit(b*t+a, X, Y, t)

LCZdM0ZFREE3NDBEQTc0MERCOSIiIiomXTNGRUUxNDdBRTE0N0FFMTFGJEkidEc2IkYkRiQ=

Fit a polynomial of degree 2 to the data:Fit(c*t^2+b*t+a, X, Y, t)

LChdM0ZGRjVDMjhGNUMyOEY2QSIiIiomXTNGQzUxRUI4NTFFQjg0RjJGJEkidEc2IkYkRiQqJl0zRkJDNTdDNTdDNTdDNTg1RiQpRiciIiNGJEYk

The 'output' option is useful when you want to return specific values from the regression model. Available outputs are: AtkinsonTstatisticconfidenceintervalsCookDstatisticdegreesoffreedomexternallystandardizedresidualsinternallystandardizedresidualsleastsquaresfunctionleveragesparametervaluesparametervectorresidualsresidualmeansquareresidualstandarddeviationresidualsumofsquaresrsquaredrsquaredadjustedstandarderrorstprobabilitytvaluevariancecovariancematrixOutput the standard errors: Fit(c*t^2+b*t+a, X, Y, t, output=standarderrors)

LUkobWFjdGlvbkc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkobWZlbmNlZEdGJDYoLUknbXRhYmxlR0YkNjUtSSRtdHJHRiQ2KC1JJG10ZEdGJDYoLUkjbW5HRiQ2JFExMS4xNzE5OTA1NzM2NjU5OEYnLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGJy8lKXJvd2FsaWduR1EhRicvJSxjb2x1bW5hbGlnbkdGQC8lK2dyb3VwYWxpZ25HRkAvJShyb3dzcGFuR1EiMUYnLyUrY29sdW1uc3BhbkdGRy1GNTYoLUY4NiRRMjAuNzY2NzQ4MjU4MDA2ODAwRidGO0Y+RkFGQ0ZFRkgtRjU2KC1GODYkUTIwLjEwNzIyNjE1ODA5Mzk2NEYnRjtGPkZBRkNGRUZIRj5GQUZDLyUmYWxpZ25HUSVheGlzRicvRj9RKWJhc2VsaW5lRicvRkJRJ2NlbnRlckYnL0ZEUSd8ZnJsZWZ0fGhyRicvJS9hbGlnbm1lbnRzY29wZUdRJXRydWVGJy8lLGNvbHVtbndpZHRoR1ElYXV0b0YnLyUmd2lkdGhHRlxvLyUrcm93c3BhY2luZ0dRJjEuMGV4RicvJS5jb2x1bW5zcGFjaW5nR1EmMC44ZW1GJy8lKXJvd2xpbmVzR1Elbm9uZUYnLyUsY29sdW1ubGluZXNHRmdvLyUmZnJhbWVHRmdvLyUtZnJhbWVzcGFjaW5nR1EsMC40ZW1+MC41ZXhGJy8lKmVxdWFscm93c0dRJmZhbHNlRicvJS1lcXVhbGNvbHVtbnNHRmFwLyUtZGlzcGxheXN0eWxlR0ZhcC8lJXNpZGVHUSZyaWdodEYnLyUwbWlubGFiZWxzcGFjaW5nR0Zkb0Y7L0krbXNlbWFudGljc0dGJFEqUm93VmVjdG9yRicvJSVvcGVuR1EiW0YnLyUmY2xvc2VHUSJdRidGW3EvJSthY3Rpb250eXBlR1EucnRhYmxlYWRkcmVzc0YnLyUpcnRhYmxlaWRHUTUxODQ0Njc0NDA3NDM4MjM5MzQ2MkYn

Output the r-squared and adjusted r-squared value:Fit(c*t^2+b*t+a, X, Y, t, output=[rsquared,rsquaredadjusted])

NyRdM0ZFRDlBRkM1MkY3MzFGQl0zRkVDMDI0RjM0RjE1MzRE

It is often useful to display the results in summary form: Fit(c*t^2+b*t+a, X, Y, t, summarize=true):

Summary: ---------------- Model: 1.9600000+.16500000*t+.11071429*t^2 ---------------- Coefficients: Estimate Std. Error t-value P(>|t|) a 1.9600 1.1720 1.6724 0.1930 b 0.1650 0.7667 0.2152 0.8434 c 0.1107 0.1072 1.0325 0.3778 ---------------- R-squared: 0.9252, Adjusted R-squared: 0.8753

You can also display this summary as an embedded table, complete with more details on residuals.Fit(c*t^2+b*t+a, X, Y, t, summarize=embed): Model:LCgkIikrK2c+ISIoIiIiKiYkIikrK107ISIpRiYlInRHRiZGJiomJCIpSDkyNkYqRiYpRisiIiNGJkYmCoefficientsEstimateStandard Errort-valueP(>|t|)aJCInK2c+ISImJCInKj48IiEiJg==JCInUHM7ISImJCInWEk+ISInbJCInK107ISInJCInW253ISInJCInJT46IyEiJw==JCInOk0lKSEiJw==cJCInOTI2ISInJCInRXM1ISInJCInYEs1ISImJCIncHhQISInR-squared:JCIncF4jKiEiJw==Adjusted R-squared:JCInI0d2KSEiJw==

Residual Sum of SquaresResidual Mean SquareResidual Standard ErrorDegrees of FreedomJCIncihHIiEiJg==JCInUSNIJSEiJw==JCInal5sISInIiIk Five Point Summary MinimumFirst QuartileMedianThird QuartileMaximumJCEnSDkqKSEiJw==JCEnZC5IISInJCInOWQ6ISInJCInJmYhSCEiJw==JCInciZbJiEiJw==

In this example, you will use a random number generator to add noise to a signal, then fit a model to the resulting data. restart; with(Statistics):JSFH Define the noise as a random number generator with mean value 1 and standard deviation 0.5. Noise := Sample(RandomVariable(Normal(1, .5))): We use a sample size of 200 to create a vector of noise values. PureNoiseVector := Noise(200)

LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEwUHVyZU5vaXNlVmVjdG9yRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2LVEqfiYjODc4ODt+RicvRjNRJ25vcm1hbEYnLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0Y9LyUpc3RyZXRjaHlHRj0vJSpzeW1tZXRyaWNHRj0vJShsYXJnZW9wR0Y9LyUubW92YWJsZWxpbWl0c0dGPS8lJ2FjY2VudEdGPS8lJ2xzcGFjZUdRJjAuMGVtRicvJSdyc3BhY2VHRkwtSShtYWN0aW9uR0YkNiUtSShtZmVuY2VkR0YkNiYtSSdtdGFibGVHRiQ2OC1JJG10ckdGJDYmLUkkbXRkR0YkNigtRiM2Ji1GLDYlUSt+MX4uLn4yMDB+RidGL0YyLUklbXN1YkdGJDYlLUYsNiVRJ1ZlY3RvckYnRi9GMi1GIzYlLUYsNiVRJHJvd0YnRi9GMi8lK2ZvcmVncm91bmRHUShbMCwwLDBdRidGOS8lL3N1YnNjcmlwdHNoaWZ0R1EiMEYnRmhvRjkvJSlyb3dhbGlnbkdRIUYnLyUsY29sdW1uYWxpZ25HRmBwLyUrZ3JvdXBhbGlnbkdGYHAvJShyb3dzcGFuR1EiMUYnLyUrY29sdW1uc3BhbkdGZ3BGXnBGYXBGY3AtRlk2Ji1GZm42KC1GIzYmLUYsNiVRLERhdGF+VHlwZTp+RidGL0YyLUZebzYlLUYsNiVRJmZsb2F0RidGL0YyLUYjNiUtSSNtbkdGJDYkUSI4RidGOUZob0Y5RltwRmhvRjlGXnBGYXBGY3BGZXBGaHBGXnBGYXBGY3AtRlk2Ji1GZm42KC1GIzYmLUYsNiVRKlN0b3JhZ2U6fkYnRi9GMi1GLDYlUSxyZWN0YW5ndWxhckYnRi9GMkZob0Y5Rl5wRmFwRmNwRmVwRmhwRl5wRmFwRmNwLUZZNiYtRmZuNigtRiM2Ji1GLDYlUShPcmRlcjp+RidGL0YyLUYsNiVRLkZvcnRyYW5fb3JkZXJGJ0YvRjJGaG9GOUZecEZhcEZjcEZlcEZocEZecEZhcEZjcC8lJmFsaWduR1ElYXhpc0YnL0ZfcFEpYmFzZWxpbmVGJy9GYnBRJ2NlbnRlckYnL0ZkcFEnfGZybGVmdHxockYnLyUvYWxpZ25tZW50c2NvcGVHRjEvJSxjb2x1bW53aWR0aEdRJWF1dG9GJy8lJndpZHRoR0ZjdC8lK3Jvd3NwYWNpbmdHUSYxLjBleEYnLyUuY29sdW1uc3BhY2luZ0dRJjAuOGVtRicvJSlyb3dsaW5lc0dRJW5vbmVGJy8lLGNvbHVtbmxpbmVzR0ZedS8lJmZyYW1lR0ZedS8lLWZyYW1lc3BhY2luZ0dRLDAuNGVtfjAuNWV4RicvJSplcXVhbHJvd3NHRj0vJS1lcXVhbGNvbHVtbnNHRj0vJS1kaXNwbGF5c3R5bGVHRj0vJSVzaWRlR1EmcmlnaHRGJy8lMG1pbmxhYmVsc3BhY2luZ0dGW3VGOS8lJW9wZW5HUSJbRicvJSZjbG9zZUdRIl1GJy8lK2FjdGlvbnR5cGVHUS1icm93c2VydGFibGVGJy8lKXJ0YWJsZWlkR1E1MTg0NDY3NDQwNzQzMjkwNjQ2MjJGJ0Zoby8lKXJlYWRvbmx5R0Y9Rjk=

We can then extract entries from the vector. The vector indices go from 1 to 200. To select the LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYqLUkjbWlHRiQ2I1EhRictRiM2KC1GLDYnUSJpRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLyUqcmVmZXJlbmNlR1ElOTUwNEYnLyUyZGlzcGxheV9yZWZlcmVuY2VHUSU5NTAzRicvJStleGVjdXRhYmxlR1EmZmFsc2VGJy8lMGZvbnRfc3R5bGVfbmFtZUdRKTJEfklucHV0RicvRjhRJ25vcm1hbEYnL0Y7USYxODg4NEYnL0Y+USYxODg4M0YnRitGQEZDRkZGSEZKth entry, use the selection operator: PureNoiseVector[i].

SSJpRzYi

PureNoiseVector[1]

XTNGRERBRUIzOTg1OUNCQjY=

Next, we define a signal; the following creates a noisy signal by adding noise to the original signal: signal := x -> x*sin(x): Generate a data set using the seq command.noisyData := [seq([.1*(i-100), signal(.1*(i-100))+PureNoiseVector[i]], i = 1..200)]: Plot the signal and noisyData together.plot([noisyData, signal(x)], style=point, symbol=[cross, circle], symbolsize=15, color=["Black", "Blue"], legend=["Noisy Signal", "Signal"], view=[default,-5..5], axes=boxed)

Ni4tJSdDVVJWRVNHNiY3ZHc3JCQhIyoqISIiJCEyWDgqKkhUPGUxJSEjOzckJCEjKSohIiIkITJDYHRBU01ndiMhIzs3JCQhIygqISIiJCEyWypbWFdTaVc+ISM7NyQkISMnKiEiIiQhMVtOeHRzLmpjISM7NyQkISMmKiEiIiQiMmsyKXk6SF1MRiEjPDckJCEjJSohIiIkIjE7KilHVCcqSCg0IyEjOjckJCEjJCohIiIkIjJ3L0c9dygqKlI4ISM7NyQkISMjKiEiIiQiMi86eGVIdGgkUSEjOzckJCEjIiohIiIkIjJVPzB4SXUqKSlSISM7NyQkISMhKiEiIiQiMTBmUHJiSTdfISM6NyQkISMqKSEiIiQiMWIiekQnPlckZiYhIzo3JCQhIykpISIiJCIxNjFfdTVHOWEhIzo3JCQhIygpISIiJCIxXy0zdkwmKlxuISM6NyQkISMnKSEiIiQiMjpvby1PakRKKCEjOzckJCEjJikhIiIkIjEiUSQpZSYqWyZleCEjOjckJCEjJSkhIiIkIjFNdm8heUlqMykhIzo3JCQhIyQpISIiJCIxY2I9QVImKiplKSEjOjckJCEjIykhIiIkIjFQaylIbGFYbSkhIzo3JCQhIyIpISIiJCIxNCF6JT0lZVoiKSkhIzo3JCQhIyEpISIiJCIxJ0hSKjNrITMkKSkhIzo3JCQhI3ohIiIkIjFuUiQqXHNrcyopISM6NyQkISN5ISIiJCIxI3BqIlIoeTFYKSEjOjckJCEjeCEiIiQiMUxkKCp5bFwpUSkhIzo3JCQhI3chIiIkIjF2KEgnKlwqPi8pKSEjOjckJCEjdiEiIiQiMTcrcC04XmwjKSEjOjckJCEjdSEiIiQiMUJwOE4lRyUqRykhIzo3JCQhI3QhIiIkIjEiSCgqM0pxYUsoISM6NyQkISNzISIiJCIyJmVKTnQ4RidwJyEjOzckJCEjciEiIiQiMUZSQ0JoRF5mISM6NyQkISNxISIiJCIxV3dMWyc+XywmISM6NyQkISNwISIiJCIxeCRIPyRIRSNRJSEjOjckJCEjbyEiIiQiMmAob2AxW2kjeSQhIzs3JCQhI24hIiIkIjEqelQnSFRLblAhIzo3JCQhI20hIiIkIjFmJHonMyRvaXUjISM6NyQkISNsISIiJCIyOVlBYFViJD5GISM7NyQkISNrISIiJCIxeTNOXyVIMG0iISM6NyQkISNqISIiJCIwNGdbXzkhZnAhIzo3JCQhI2khIiIkIjIocGgiM0k4IWU2ISM7NyQkISNoISIiJCExRDF6Z2I1VWchIzs3JCQhI2chIiIkITJUIypcZFE3JWY2ISM7NyQkISNmISIiJCEyPihHIT49Jj0pKT0hIzs3JCQhI2UhIiIkITE5J3pAIkdgJ3ApISM7NyQkISNkISIiJCEySyIqUmkkbylbWSIhIzs3JCQhI2MhIiIkITEmSCQ0Tix5WkkhIzo3JCQhI2IhIiIkITJhL28wd3hidyMhIzs3JCQhI2EhIiIkITJ2dGZKJTRZSUMhIzs3JCQhI2AhIiIkITFOejByPjxoUyEjOjckJCEjXyEiIiQhMS8iPSJvKEg/eSQhIzo3JCQhI14hIiIkITEoUm07eWkycCQhIzo3JCQhI10hIiIkITIvW0glZnVtJFIkISM7NyQkISNcISIiJCEyKEhhRlBiK3pRISM7NyQkISNbISIiJCEySyR6cGN5TURQISM7NyQkISNaISIiJCExckU1OiFmV1EkISM6NyQkISNZISIiJCExJzNVb2VFOWUkISM6NyQkISNYISIiJCEyKCpwKVxUcDU5USEjOzckJCEjVyEiIiQhMSNmcStRdiQ9TCEjOjckJCEjViEiIiQhMS9TKSplT20qRyQhIzo3JCQhI1UhIiIkITJYaiRwOVAkR28iISM7NyQkISNUISIiJCExLzgiKVFLbiVbIyEjOjckJCEjUyEiIiQhMTo3RzdERG8/ISM6NyQkISNSISIiJCEyRV0obyVIKVx1NSEjOzckJCEjUSEiIiQhMilHKXo9KmVpOz0hIzs3JCQhI1AhIiIkIjFXcW5XTWdeXCEjPTckJCEjTyEiIiQhMU5pMTwvMTxgISM7NyQkISNOISIiJCEyKSopSFVVJ2VvNyIhIz03JCQhI00hIiIkITFEU0pkbGg1OCEjOzckJCEjTCEiIiQiMWQiUlkmKXBhImYhIzs3JCQhI0shIiIkIjA5MyZwJioqMzEnISM6NyQkISNKISIiJCIyJD5baEBKbVk2ISM7NyQkISNJISIiJCIxJioqUkoyO19UIyEjOjckJCEjSCEiIiQiMiZ5XlJcR14rRyEjOzckJCEjRyEiIiQiMVFNIlx3ImV5QiEjOjckJCEjRiEiIiQiMjFAIUctUk4oZiMhIzs3JCQhI0UhIiIkIjIqKWZYTGZxJTQ+ISM7NyQkISNEISIiJCIyTmJIJDNaKEcvIyEjOzckJCEjQyEiIiQiMjI5RDM5WUY9IyEjOzckJCEjQiEiIiQiMmIrIT5lblAxRSEjOzckJCEjQSEiIiQiMUthJjNdOTFoIyEjOjckJCEjQCEiIiQiMmFjNCFlelUoMyQhIzs3JCQhIz8hIiIkIjFNZShwaj8pPkEhIzo3JCQhIz4hIiIkIjFRLUIoKlxATUIhIzo3JCQhIz0hIiIkIjFpKSlHTC8mejMjISM6NyQkISM8ISIiJCIxNl5FaUBoP0AhIzo3JCQhIzshIiIkIjAsTi4qR0ROSyEjOTckJCEjOiEiIiQiMj1adiNbWil6diIhIzs3JCQhIzkhIiIkIjJsL291JzNLJz0kISM7NyQkISM4ISIiJCIyckQ7amdhKjQ+ISM7NyQkISM3ISIiJCIxQC5ZYExsZkUhIzo3JCQhIzYhIiIkIjJvMWdla2c+ZiIhIzs3JCQhIzUhIiIkIjF5RXZsLXRqOiEjOjckJCEiKiEiIiQiMUVUW1UsPlRBISM6NyQkISIpISIiJCIydng5Ok9oJjQ7ISM7NyQkISIoISIiJCIyQUFiKD00KilmOCEjOzckJCEiJyEiIiQiMiMpRzRlRGdJbCIhIzs3JCQhIiYhIiIkIjFoJlt6WXV5dykhIzs3JCQhIiUhIiIkIjElKSkqKmZucjJcKiEjOzckJCEiJCEiIiQiMkhJKWU3UCcpZj0hIzs3JCQhIiMhIiIkIjFgNFtDOnomSCohIzs3JCQhIiIhIiIkIjFDaylHSGBkKmUhIzs3JCQiIiEhIiIkIjFCTC4xSWAjPiIhIzo3JCQiIiIhIiIkIjJPRCZmJio+JikzNiEjOzckJCIiIyEiIiQiMTNzRnJzJ0hhJyEjOzckJCIiJCEiIiQiMihRTUc4LVFKNiEjOzckJCIiJSEiIiQiMiY0b19CJEdZYiIhIzs3JCQiIiYhIiIkIjFzSFthI2YrWighIzs3JCQiIichIiIkIjJWREF1clY5PCIhIzs3JCQiIighIiIkIjIkR1dlbClwW1ciISM7NyQkIiIpISIiJCIyQ0VYJCl6JClwOyMhIzs3JCQiIiohIiIkIjIpej1JRC0wbzohIzs3JCQiIzUhIiIkIjEuJ3oob3hNc0AhIzo3JCQiIzYhIiIkIjFzY1YnMygzVkQhIzo3JCQiIzchIiIkIjIpSHQ+Iykzd00/ISM7NyQkIiM4ISIiJCIyYjsxLXIjKUhkIyEjOzckJCIjOSEiIiQiMSZmVElJS3JcIyEjOjckJCIjOiEiIiQiMidcM2MzcipHRiMhIzs3JCQiIzshIiIkIjJCXmBMWVFOXSMhIzs3JCQiIzwhIiIkIjIyVTFVVigqbypIISM7NyQkIiM9ISIiJCIxMGo+KFxHJjNGISM6NyQkIiM+ISIiJCIyOWt1Pnk3YFsjISM7NyQkIiM/ISIiJCIxNE9iKm9lUipRISM6NyQkIiNAISIiJCIyJHA/REdwdiZ6IyEjOzckJCIjQSEiIiQiMXZ6TFoqW1teJCEjOjckJCIjQiEiIiQiMmtSWUpfTVIhRyEjOzckJCIjQyEiIiQiMk91QnBFUkRVIyEjOzckJCIjRCEiIiQiMkNJLyQ0KlEvciIhIzs3JCQiI0UhIiIkIjJsbVIlXD4qUVYiISM7NyQkIiNGISIiJCIyWGZBUHJaY1kiISM7NyQkIiNHISIiJCIyMHojR2UtUyF5IiEjOzckJCIjSCEiIiQiMil6JFxDRzdcPCIhIzs3JCQiI0khIiIkIjJmYzEkKXBtJCopPiEjOzckJCIjSiEiIiQiMTEyRWVxYyFHJSEjOzckJCIjSyEiIiQiMigpcGZtRFtzdSIhIzs3JCQiI0whIiIkITFGJHo2bzg6LCkhIzw3JCQiI00hIiIkITBLS0FPdDVoKCEjOjckJCIjTiEiIiQiMVJmKDRAZStIJiEjOzckJCIjTyEiIiQhMjwqZilbJVFIJ3kiISM7NyQkIiNQISIiJCEyI3pSUiQpRyQqKVsiISM7NyQkIiNRISIiJCEyaiQ9TjgiXCRbNyEjOzckJCIjUiEiIiQhMXotPWEoW14oeiEjOzckJCIjUyEiIiQhMU4yVHBRPXQ4ISM6NyQkIiNUISIiJCEwYCpcJ1JHblAkISM5NyQkIiNVISIiJCEyRXBaOUJnJykqRyEjOzckJCIjViEiIiQhMSdSQiRwJ1FbUCQhIzo3JCQiI1chIiIkITImKj43KHolSDNHJCEjOzckJCIjWCEiIiQhMSo0WClcLFVOTSEjOjckJCIjWSEiIiQhMnRoNid6YSJRUyQhIzs3JCQiI1ohIiIkITIwYmVPK0c6XCQhIzs3JCQiI1shIiIkITImb2giPiNwa0RLISM7NyQkIiNcISIiJCExVFpDIVIvVSpIISM6NyQkIiNdISIiJCExUnZaV3VhYUohIzo3JCQiI14hIiIkITJzLHVBTSdwMFAhIzs3JCQiI18hIiIkITJ2NiV5Z3okbyVIISM7NyQkIiNgISIiJCEyYU4yNmtjKj5KISM7NyQkIiNhISIiJCExWkxFIzQ/SyZHISM6NyQkIiNiISIiJCEyOip6QzVXd21IISM7NyQkIiNjISIiJCExOSwnb1teTjIjISM6NyQkIiNkISIiJCEyKSopXGokKjMxXDwhIzs3JCQiI2UhIiIkITJ2L1UzSC95KT4hIzs3JCQiI2YhIiIkITImenBseC5qIlsjISM7NyQkIiNnISIiJCExXjxZPW08VUQhIzs3JCQiI2ghIiIkIjJOIls2bT1cIlwiISM8NyQkIiNpISIiJCEyTSVHZWAwNy5AISM8NyQkIiNqISIiJCIxTlgzSjUwMXghIzs3JCQiI2shIiIkIjIoUSYqKmVjIz1sOCEjOzckJCIjbCEiIiQiMHI4KFtoeiEzIyEjOTckJCIjbSEiIiQiMjEnZXFtNSZcZyMhIzs3JCQiI24hIiIkIjEjSDkpPiJvKFxWISM6NyQkIiNvISIiJCIyREhhRmAib1BVISM7NyQkIiNwISIiJCIxJik0c2k9YXBfISM6NyQkIiNxISIiJCIxRWojKjNiKFxhJiEjOjckJCIjciEiIiQiMVx0JT4ocGNbYyEjOjckJCIjcyEiIiQiMFglKnA4W0ljJyEjOTckJCIjdCEiIiQiME4nRytrNVduISM5NyQkIiN1ISIiJCIxJCo9cFUyJkgkeiEjOjckJCIjdiEiIiQiMXJXJz0lb2FGeiEjOjckJCIjdyEiIiQiMFV6Jz5qdXQkKiEjOTckJCIjeCEiIiQiMUUhPlYpRzQreiEjOjckJCIjeSEiIiQiMT4iUmpqJCo9RSkhIzo3JCQiI3ohIiIkIjF5LFkqKTNHSygqISM6NyQkIiMhKSEiIiQiMWJpK0J5QlAqKSEjOjckJCIjIikhIiIkIjByaiw3XFVLKSEjOTckJCIjIykhIiIkIjFuWShRPyUpKnkhKiEjOjckJCIjJCkhIiIkIjFvV193NC1pKCkhIzo3JCQiIyUpISIiJCIxalFKUiYpZSc9KSEjOjckJCIjJikhIiIkIjBwQTlPKlxDeSEjOTckJCIjJykhIiIkIjFMOHBSKFsjenMhIzo3JCQiIygpISIiJCIxUHkjKnpzYkVsISM6NyQkIiMpKSEiIiQiMTEkZiFmUSRSIXAhIzo3JCQiIyopISIiJCIyJT0sJmUmR2Y5XyEjOzckJCIjISohIiIkIjEiKipHIlJYRlJdISM6NyQkIiMiKiEiIiQiMXAlXFUkPi4kKlshIzo3JCQiIyMqISIiJCIyanBdUWAqKVx6JCEjOzckJCIjJCohIiIkIjJadndUQnlOdyIhIzs3JCQiIyUqISIiJCIya0tcK2RuKVE4ISM7NyQkIiMmKiEiIiQhMkY/a14rWHYhSCEjPTckJCIjJyohIiIkITIjekArJkdCNjAiISM7NyQkIiMoKiEiIiQhMXl3OVlMUig9KSEjOzckJCIjKSohIiIkITIjPT11JSpcMW1FISM7NyQkIiMqKiEiIiQhMm0leTlsJ1wyNSQhIzs3JCQiJCsiISIiJCExOz44TydHTjIlISM6LSUnTEVHRU5ERzYjLSUpX1RZUEVTRVRHNiMtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjUtSSNtc0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjYjUS1Ob2lzeX5TaWduYWw2Ii8lJ2ZhbWlseUdRMFRpbWVzfk5ld35Sb21hbjYiLyUlc2l6ZUdRIzEyNiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJmZhbHNlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1Enbm9ybWFsNiItJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lJ1NZTUJPTEc2IyUmQ1JPU1NHLSUnQ1VSVkVTRzYmN2B4NyQkIS1oJUgmPSRHJyEjNiQhMVJfOzRUcSgqeSEjQjckJCExRURYIkc7ckAnISM6JCEyOXFDeFgrWTUlISM8NyQkITEtTFAnKSpIJ2ZoISM6JCExRDtMI3AsN2YoISM7NyQkITEpel9faiEpXDQnISM6JCExa1EpcClRTVM2ISM6NyQkITFVd0VZRiEqSGchIzokITI5TUc+JD4qNF4iISM7NyQkITFfKkc2O01eJ2YhIzokITEoZSc9X0pSbD0hIzo3JCQhMSZRXDxSJjMwZiEjOiQhMm49NjBMJCopekAhIzs3JCQhMEo8I3AkM0glZSEjOSQhMnVTZiRwIio9IVwjISM7NyQkITEpW3BOIVtneWQhIzokITF1QCIqUS1oJHojISM6NyQkITI6QXMtWTJYciYhIzskITJhRSEpNHRudDIkISM7NyQkITEyZzdrYWRbYyEjOiQhMSxBY2NNJClbTCEjOjckJCExKilHLkVEXSFmJiEjOiQhMm1yUmc3PCxkJCEjOzckJCExOiF5QTpFXl8mISM6JCEybGU1QywpZil6JCEjOzckJCExczcmb0wiW2ZhISM6JCExJCkpKjM1YFcwUyEjOjckJCExKjQqXCRSP2lSJiEjOiQhMWw4azFFI0g9JSEjOjckJCEyTmtHa1J0KFFgISM7JCExLz9WeCc+X0slISM6NyQkITFhJD4wT2ovRiYhIzokITF0ITN4cDQzWiUhIzo3JCQhMUU/JykqeiZmN18hIzokITF2X3hEaSdSZCUhIzo3JCQhMWtqLHE/R1heISM6JCEwKmVTcSxtcVkhIzk3JCQhMXQuQjghKXAmMyYhIzokITFpTFdbUVVOWiEjOjckJCEwO3IkZltLP10hIzkkITFGKyplKls8JXklISM6NyQkITFtP1YzJio+KilcISM6JCExQXlXWjBGKnolISM6NyQkITF0SFxkVDJlXCEjOiQhMXpxeG93PTRbISM6NyQkITFtQkZRaCQ9JVwhIzokITIwSnojKik9SzdbISM7NyQkITB3XiE+IilmRFwhIzkkITFpYyxCdDE5WyEjOjckJCExYTYkKSo0ZyQ0XCEjOiQhMWNgXGo9VjlbISM6NyQkITJ2YTUxM0FKKlshIzskITFYVnUoekJNIlshIzo3JCQhMTsmKlJSKilIalshIzokITE9OiYqPmctM1shIzo3JCQhMSVbKT0penZNJFshIzokITElMypmPVszKXolISM6NyQkITFhZlU2dkksWyEjOiQhMTp6eUsjUkJ5JSEjOjckJCExQU1tQyNSInBaISM6JCExRSNmUmNoOXclISM6NyQkITFsYUIxO0otWiEjOiQhMSlRRS0lRzItWiEjOjckJCExKjNVeid6OFdZISM6JCExalJ4LWJLTFkhIzo3JCQhMVdqJ1slKTM4ZSUhIzokITElZlUhUXIrVVghIzo3JCQhMSMpXCpHQCtrXiUhIzokITFQJClcU2skKkhXISM6NyQkITE6JXkoeigqKkdYJSEjOiQhMmwpRywqKj0jUUklISM7NyQkITEqcDFjJSlmOVIlISM6JCExVm1rNVJDblQhIzo3JCQhMUE1JFJGVEtLJSEjOiQhMXpoIylRVysrUyEjOjckJCExem84YFElPkUlISM6JCEyMGh1a1ZIbyRRISM7NyQkITFMbXdzZFwnPiUhIzokITEpUiFcJTNLLmwkISM6NyQkITEuUVglMyE+UFQhIzokITFzOVApXF05WiQhIzo3JCQhMVdHVy1MTnNTISM6JCExXiMqKXpgSWpFJCEjOjckJCExMyI+TSVwTTZTISM6JCEyOGN6SipHTmxJISM7NyQkITFvZjA/eGRaUiEjOiQhMkMlMyllJilHI1tHISM7NyQkITJDbmwjWztCJilRISM7JCEyeV8rSTMpPUlFISM7NyQkITFWVy4qM2gqPlEhIzokITJkIjR4ZGE2KFIjISM7NyQkITEiei9IeSg0ZFAhIzokITJMKj1ZKFJSI3BAISM7NyQkITFSKFFHKDQiR3AkISM6JCEyTyQpKjNHIz5TJD4hIzs3JCQhMVZ4LCFcYyFITyEjOiQhMmUkUkFqOiEpKnAiISM7NyQkITIxcCplVEZacU4hIzskITJ1LHkpKkgheiVbIiEjOzckJCExQyY9WjpITF0kISM6JCEyJ1trVCpHQylSNyEjOzckJCEyQ1cnPkleRlZNISM7JCExLl02Pi80QjUhIzo3JCQhMURmWDldQ3pMISM6JCExTXlKUk1aYnohIzs3JCQhMlhkTE8tZ3pKJCEjOyQhMVxaIT0lKTM6I2UhIzs3JCQhMVUoZXQiNCIqXEshIzokITFWZiNST25MXiQhIzs3JCQhMk1nTjlSJik0PiQhIzskITE8P3FuVVp2OiEjOzckJCEwamFDUi9UNyQhIzkkIjFvTT9hb0JqYSEjPDckJCEyY3YnW0EnXEoxJCEjOyQiMi0peiRlJVJPK0MhIzw3JCQhMihIcWMkb3VrKkghIzskIjIuVGs7Km88TFYhIzw3JCQhMicqXHhzW1gqUUghIzskIjFLSlRrMiw6ZiEjOzckJCEyJ2YkZnpXRU0oRyEjOyQiMUYiUWwxTk5oKCEjOzckJCExNzlfIUhRLCJHISM6JCIwZkUmZj5yVyIqISM6NyQkITFiY3F4OSpvdSMhIzokIjEqKSoqM2Qnb2kwIiEjOjckJCEyKEhMXVF0KFFvIyEjOyQiMk5DUDp0L2c9IiEjOzckJCEyUDNIIXkvTUJFISM7JCIyJVwqUSI9Wl0qSCIhIzs3JCQhMiZSJWUnRyYpKnliIyEjOyQiMnlvYSpRUm80OSEjOzckJCEyLkFsZGdGYlwjISM7JCIyRyQ+JTNNSUNdIiEjOzckJCEyTjR2ajkmKSlIQyEjOyQiMnUpNENtYC0oZSIhIzs3JCQhMkMoKipcOEhZcUIhIzskIjJiaGpBRyUzXzshIzs3JCQhMjpxK0FBIXovQiEjOyQiMkhXMyIqUj04ciIhIzs3JCQhMiNvLCJ5QSgpPUMjISM7JCIyUCtGXis/ZHYiISM7NyQkITJMYG81TVAiekAhIzskIjI7LCllXFQ+KXkiISM7NyQkITFBdTpUW2Y4QCEjOiQiMlpWJj50LyEpND0hIzs3JCQhMSpHRGcheUFgPyEjOiQiMmUhZlwhKT0qKT09ISM7NyQkITJobVZZcS45Kj4hIzskIjEpZi9GO055Ij0hIzo3JCQhMic9SWxdRDlCPiEjOyQiMkpBW3FZJipcIT0hIzs3JCQhMU54KHBtUDgnPSEjOiQiMiN5aC1HcUskeSIhIzs3JCQhMSd5SG5mSSIpeiIhIzokIjEyJ1FjKWUnPXYiISM6NyQkITI5XClRTTEmUXQiISM7JCIxai9IbEgmM3IiISM6NyQkITI8NzYwJVt3dTshIzskIjJ0d09pTUBkbSIhIzs3JCQhMi5henZIInk2OyEjOyQiMTpiKHl4Ri9oIiEjOjckJCEycygzWio9eCNcOiEjOyQiMk1AJikpKVs9Kls6ISM7NyQkITImSGxVQXNUI1siISM7JCIyN2JXJzM6anc5ISM7NyQkITEnPXZ5PjtNVSIhIzokIjJZYzt0QCYpelMiISM7NyQkIS9VZ2o3TWI4ISM4JCIyKEdaNilwL1NLIiEjOzckJCEyTnl5UCkpM1JIIiEjOyQiMioqUk8iSF5pVzchIzs3JCQhMmMqKT5ddkNKQiIhIzskIjFAPy4nWylbajYhIzo3JCQhMiM9Pl8qR2l4OyIhIzskIjI3Q2FUaSY+dTUhIzs3JCQhMmpieCQ+YTgtNiEjOyQiMHhkUDBsSCQpKiEjOjckJCExa2B0JW9DQy8iISM6JCIwJnlZJTMkby0hKiEjOjckJCExailvKW9gOSF6KiEjOyQiMTEuPyIpPkpEIikhIzs3JCQhMTNwQStWS3IiKiEjOyQiMVZYV3h4diFHKCEjOzckJCExKVtvTHQ6MF0pISM7JCIxcCllQEJnbFEnISM7NyQkITBcTyM9KXorI3ohIzokIjE0dypSeU9zaiYhIzs3JCQhMSJHV3kxR01EKCEjOyQiMWBIKD42bT0iWyEjOzckJCExbWpWdi51RW0hIzskIjElelh6Vl1wMiUhIzs3JCQhMUpecS1lNTFnISM7JCIwT1tCW0dWUiQhIzo3JCQhMVJXYDV2KFxRJiEjOyQiMWY5dUkicDt3IyEjOzckJCExdUxacGZZXlohIzskIjI7QWRAW1pPPCMhIzw3JCQhMWo7IipbViYqelMhIzskIjIkKjM7WypRISk9OyEjPDckJCExVDxhXFJDa00hIzskIjJjRiozJUhRaTwiISM8NyQkITJXMiVRQyhlayZHISM8JCIxRj1vVyhcKVshKSEjPDckJCEyTnJsOktZbj8jISM8JCIweCsleTtJSVshIzs3JCQhMnkiKWVlR1slZTohIzwkIjInUlx6QzolKj1DISM9NyQkITEmXHUnKipIaWwoKiEjPCQiMVI9JDREKWVAJiohIz03JCQhMi1dLTIpUXdzRyEjPSQiMTJPZkJzal4jKSEjPjckJCIyJG81KjMqNCdbKkchIz0kIjEpPmh1dFwheiQpISM+NyQkIjE0XHc3aHVTJyohIzwkIjEoZSsoNHcrIUcqISM9NyQkIjIjelAzJipSIXpnIiEjPCQiMT5hNFkmSFVkIyEjPDckJCIyLS13ZSVwdyM9IyEjPCQiMVIleiE9enNFWiEjPDckJCIxazUzZC9FSEchIzskIjEkejA9IipbJCkqeSEjPDckJCIyamlLcE1SK1skISM8JCIyO3loUENxbj0iISM8NyQkIjInSCY+JCk+RHg3JSEjPCQiMnY7b1U5UmVsIiEjPDckJCIyTT46QCpHQEdaISM8JCIxRi9eXHVBYEAhIzs3JCQiMV9mVjxKKSpcYCEjOyQiMihHaFckPUx3cyMhIzw3JCQiMDs5Ung9SSpmISM6JCIydCpcLkghZi9RJCEjPDckJCIxRG8pbz8jKlJqJyEjOyQiLyd5JEdOPyYzJSEjOTckJCIxcCpbJG9ASiRIKCEjOyQiMUZdXVNFMGdbISM7NyQkIjFcLEdcOi91eSEjOyQiMUc1KTMhZiUqeWIhIzs3JCQiMSgpKUdvRzB5XykhIzskIjEqZlJ0aChSQWshIzs3JCQiMUNqNFRNRCU9KiEjOyQiMWZBZG1rQilIKCEjOzckJCIxaSE9WShHJ28iKSohIzskIjFOKSpcKT0pM2kiKSEjOzckJCIyMEVLWEdMIlI1ISM7JCIxVkIrQG8ncCYqKSEjOzckJCIyKFw6Vz9MVzI2ISM7JCIxaWwyOEd0MSoqISM7NyQkIjIieSkpNCIpM0psNiEjOyQiMmJMb04rPDMyIiEjOzckJCIyMmFXNGhDRUIiISM7JCIxQiNbRj83RzsiISM6NyQkIjFKMHRuJzNBSCIhIzokIjJ0QyJvKUcoUVU3ISM7NyQkIjJXdSplQD1lZDghIzskIjIkRzE0Qi8lb0siISM7NyQkIjIzJHpZQkQkKT45ISM7JCIyI1ttbyY9JW8uOSEjOzckJCIyYk5dLmcleSVbIiEjOyQiMXlDY2FkSHo5ISM6NyQkIjIyVXNGKTNWVzohIzskIjF6NWk9VCpRYSIhIzo3JCQiMkZbKEhjdXczOyEjOyQiMldBMylIeWcyOyEjOzckJCIyJ1Jhc3VdZnY7ISM7JCIxdEcrXEBTbTshIzo3JCQiMTopPUlyb1B0IiEjOiQiMWMqUjg9Jno1PCEjOjckJCIxaFg0T3lmJ3oiISM6JCIyKDRMKipcNio0diIhIzs3JCQiMkIjZjFva11oPSEjOyQiMltrVkhyKVIkeSIhIzs3JCQiMi1cIz0scCtEPiEjOyQiMTYtKVxoLmIhPSEjOjckJCIyaT9hYCRvVycpPiEjOyQiMiQzKVs+YDN0Ij0hIzs3JCQiMVc/PkBoYj8/ISM6JCIydnBAPydcaD49ISM7NyQkIjI9KSlIcVNsWTAjISM7JCIyLSIqeUE7JHo9PSEjOzckJCIxYXBVPFRKJjMjISM6JCIybShSV2goR2AiPSEjOzckJCIxRVMjeSNHJ2Y2IyEjOiQiMkIoelIiKltCND0hIzs3JCQiMidwVT4zNFQiPSMhIzskIjEieTI/eD9zeSIhIzo3JCQiMjA1MmxmOzJDIyEjOyQiMTpeJjM6TGt2IiEjOjckJCIyJGYhPSZ5TGIwQiEjOyQiMTp5V2xmcTU8ISM6NyQkIjEoelR2dGZsTyMhIzokIjFKOWJVYihmbCIhIzo3JCQiMVBcITQnKkcuViMhIzokIjImKSo9XClIKVwnZSIhIzs3JCQiMjtCJnBLXW4jXCMhIzskIjJgPmQ0SiVRMTohIzs3JCQiMj1ZRT5mWHpiIyEjOyQiMkJKYDVVNCc0OSEjOzckJCIyTDZjISkqKTMzaSMhIzskIjJzaVYyMDVTSSIhIzs3JCQiMWtANzNkNCZvIyEjOiQiMlJYMHNiMk89IiEjOzckJCIxaEolND5dKVtGISM6JCIwIz5uJmZcPzAiISM5NyQkIjJMQHIkUlJWMkchIzskIjFyMVkpUm93PyohIzs3JCQiMFFVaV94WChHISM5JCIxYjddcFJuJWUoISM7NyQkIjFqV3ddOmpNSCEjOiQiMShlcmRsIUdJZyEjOzckJCIyeChcXW07bSkqSCEjOyQiMS9uQmIjUTlGJSEjOzckJCIyL0xGdGxZKmZJISM7JCIyIz0qbz4+XGJcIyEjPDckJCIyOTstT3cmKno3JCEjOyQiMF92LFFxSEQlISM7NyQkIjInKkhEJipHQHA9JCEjOyQhMiZ5ODYvTDRXOSEjPDckJCIyTUYxJilHLVFEJCEjOyQhMiVHKW8pKXBBTWskISM8NyQkIjIlW1RaZXF2OUwhIzskITFqMS8lKmVMNmQhIzs3JCQiMkwoUVIoKj5WIlEkISM7JCExZjIiM3p2Qy4pISM7NyQkIjJXUyRvJD5oKlFNISM7JCEyJDMhKlIyUGoyNSEjOzckJCIyVmEsSUIhWy9OISM7JCEyKCp5N0FxLlNDIiEjOzckJCIyOlxSL1JveGMkISM7JCEya0puK2kmKVtaIiEjOzckJCIyMERiSz86NWokISM7JCEyTU4kZWIlKipwcSIhIzs3JCQiMXV2WFUkSFNwJCEjOiQhMnAkRyFbJyoqW1E+ISM7NyQkIjFAPSRIP21YdiQhIzokITJrIilwLjo1KzsjISM7NyQkIjJhWS1COjMrI1EhIzskITJDSTVqJ1xHKFIjISM7NyQkIjJZbyY+diF6QilRISM7JCExJWZpOTomND9FISM6NyQkIjE2ZWVNOi1bUiEjOiQhMktrQDEsaShcRyEjOzckJCIyQyQ0WW5QVzJTISM7JCEyYEtfOnVkQTAkISM7NyQkIjEtLXdlazZ0UyEjOiQhMi5BdjA3KXpvSyEjOzckJCIxTzI6YCU+ZzglISM6JCEyTHNIMSU+JHlZJCEjOzckJCIxckIqKVIkcCgpPiUhIzokITFsPSplYjRxbCQhIzo3JCQiMSJbLilSPUprVSEjOiQhMWpKbzhuTVZRISM6NyQkIjE7YyRcKCl5WUslISM6JCExRko/NSMqby5TISM6NyQkIjFQc0p3SF0nUSUhIzokITFWTiYpPixqYlQhIzo3JCQiMSUpeUlJVHdhVyEjOiQhMWkoUWFPWnhJJSEjOjckJCIxckopUixwbF4lISM6JCExRT5JOS9ESVchIzo3JCQiMT42QiUzdyh6WCEjOiQhMjBvInAjb2smUlghIzs3JCQiMkVUc2wvY1NrJSEjOyQhMS0jPTMieUBMWSEjOjckJCIxJHlcLyU9OS5aISM6JCEuLCt3U0hxJSEjNzckJCIxbDhRJFFEaHclISM6JCExUkUoW3FYI2ZaISM6NyQkIjEocE51Vnh0eiUhIzokITElKm9KVz8xIXklISM6NyQkIjImRytcIlxIJ0dbISM7JCExOSIpZUtWLyd6JSEjOjckJCIyOj83XVpmPydbISM7JCExWVY7UUhxMlshIzo3JCQiMldQTSZlJSpbJipbISM7JCEybCI9TzFmbDhbISM7NyQkIjFaKyIzKCoqKVwjXCEjOiQhMiZcJCpmOmU1OVshIzs3JCQiMT5kMyRbIVxhXCEjOiQhMXpQSSg+KioqNFshIzo3JCQiMTc3QV16XykpXCEjOiQhMV5yJlEsUiYqeiUhIzo3JCQiMTFuTjxhY0FdISM6JCExeVZ5WXMpR3klISM6NyQkIjFAQD0oeigqUjMmISM6JCExdCpcZUknKXB0JSEjOjckJCIxNDUlZiM+eVdeISM6JCEyJil5OEwrJkdyWSEjOzckJCIydykqUTlSVyxAJiEjOyQhMS8mR2IvRnpkJSEjOjckJCIxW0xlaDd4dl8hIzokITFYdV5hK1VnVyEjOjckJCIxVGJBJyo+W05gISM6JCExd0YqKUhXI0dMJSEjOjckJCIyJj1TMldAKikpUiYhIzskITBAaDt5bGU8JSEjOTckJCIxOCNRNER1MlkmISM6JCExIjRpNCpIZixTISM6NyQkIjFjU2kyXiZ5XyYhIzokITElbz44ZDAmKnkkISM6NyQkIjFiczgqcCkqZWUmISM6JCEwdTwoKioqR3BlJCEjOTckJCIxd2s8dVFjX2MhIzokITE2MzkjMztJTCQhIzo3JCQiMW9zVFZFQjpkISM6JCExYjAieTpsVTIkISM6NyQkIjEkUiFwK2hIeGQhIzokITImM3A+KG8kZip6IyEjOzckJCIwWTIqSCozJVJlISM5JCExKm93Z01vcl0jISM6NyQkIjJtYzhTM2dGIWYhIzskITJjbVBDRyZ5Ij4jISM7NyQkIjInUShwZ0M2KnBmISM7JCExN0BJanR4Uj0hIzo3JCQiMHQxZ0cjW0pnISM5JCExVWBYWkc7LTohIzo3JCQiMSlcQSYzM0UjNCchIzokITI5JVxNIjMzaDoiISM7NyQkIjJOTC8pWz9CZGghIzskITFjMWdAYnVNeCEjOzckJCIxQV1QXz0xQWkhIzokITI4IWUnKnA5eCtRISM8NyQkIi1oJUgmPSRHJyEjNiQhMVJfOzRUcSgqeSEjQi0lJ0xFR0VOREc2Iy0lKV9UWVBFU0VURzYjLUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1LUkjbXNHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2I1EnU2lnbmFsNiIvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW42Ii8lJXNpemVHUSMxMjYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSZmYWxzZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ25vcm1hbDYiLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIzUhIiItJSdTWU1CT0xHNiMlJ0NJUkNMRUctJSVWSUVXRzYkJShERUZBVUxURzskISNdISIiJCIjXSEiIi0mJSZfQVhJU0c2IyIiIjYoLSUrX1RJQ0tNQVJLU0c2IyUpX1BJVElDS1NHLSUrX0dSSURMSU5FU0c2Iy0lKV9WSVNJQkxFRzYjIiIhLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiItJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSYlJl9BWElTRzYjIiIjNictJStfR1JJRExJTkVTRzYjLSUpX1ZJU0lCTEVHNiMiIiEtJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIi0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJStBWEVTTEFCRUxTRzYkLUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NVEieDYiLyUnZmFtaWx5R1EhNiIvJSVzaXplR1EjMTA2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EldHJ1ZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ2l0YWxpYzYiUSE2Ii0lJ1NZTUJPTEc2JCUoREVGQVVMVEciIzotJSZTVFlMRUc2IyUmUE9JTlRHLSUqQVhFU1NUWUxFRzYjJSRCT1hHLSUpX1ZJU0lCTEVHNiMiIiItJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiQ/JSEiIi0lKUJPVU5EU19ZRzYjJCIjISohIiItJS1CT1VORFNfV0lEVEhHNiMkIiVdTSEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiVnSiEiIi0lKUNISUxEUkVORzYiLSUrQU5OT1RBVElPTkc2Jy0lKUJPVU5EU19YRzYjJCIiISEiIi0lKUJPVU5EU19ZRzYjJCIiISEiIi0lLUJPVU5EU19XSURUSEc2IyQiJStTISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJStTISIiLSUpQ0hJTERSRU5HNiI=NiI=

Fitting a Model to the data:To fit a functional model to the noisy data, use the Fit command from the Statistics package. The calling sequence is Fit(f, X, Y, v), where f is the function model, X and Y are the x and y data, and v is the name of the independent variable in the component function.Create two lists Xdata and Ydata from the data in the two columns of noisyData.Xdata := [seq(noisyData[i, 1], i = 1 .. numelems(noisyData))]: Ydata := [seq(noisyData[i, 2], i = 1 .. numelems(noisyData))]: Tip: It's often useful to do something for each element of a data structure. Here, we do this by i=1..numelems(noisyData). Recall, numelems is the number of elements.The function model is LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYuLUkjbWlHRiQ2J1EiZkYnLyUlc2l6ZUdRIzExRicvJSdpdGFsaWNHUSV0cnVlRicvJStleGVjdXRhYmxlR1EmZmFsc2VGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYvUSI9RidGL0Y1L0Y5USdub3JtYWxGJy8lJmZlbmNlR0Y3LyUqc2VwYXJhdG9yR0Y3LyUpc3RyZXRjaHlHRjcvJSpzeW1tZXRyaWNHRjcvJShsYXJnZW9wR0Y3LyUubW92YWJsZWxpbWl0c0dGNy8lJ2FjY2VudEdGNy8lJ2xzcGFjZUdRLDAuMjc3Nzc3OGVtRicvJSdyc3BhY2VHRlEtRiw2J1EiY0YnRi9GMkY1RjgtRjw2L1EnJnNkb3Q7RidGL0Y1Rj9GQUZDRkVGR0ZJRktGTS9GUFEmMC4wZW1GJy9GU0Zlbi1GLDYnUSJ4RidGL0YyRjVGOEZXLUYsNidRJHNpbkYnRi8vRjNGN0Y1Rj8tSShtZmVuY2VkR0YkNiYtRiM2Jy1GLDYnUSJhRidGL0YyRjVGOEZXRmduRjVGP0YvRjVGPy1GPDYvUSIrRidGL0Y1Rj9GQUZDRkVGR0ZJRktGTS9GUFEsMC4yMjIyMjIyZW1GJy9GU0Zqby1GLDYnUSJiRidGL0YyRjVGOEY1Rj8=.f := Fit(c*x*sin(a*x)+b, Xdata, Ydata, x)

LCYqKF0zRkVGQzU4RkRGRjNEQjYxIiIiSSJ4RzYiRiUtSSRzaW5HNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHRic2IywkKiZdM0ZFRkY0OUIwMUY5OEEwM0YlRiZGJUYlRiVGJV0zRkYwMjQwQTZDMUE2MEY0RiU=

Plot LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2J1EiZkYnLyUlc2l6ZUdRIzExRicvJSdpdGFsaWNHUSV0cnVlRicvJStleGVjdXRhYmxlR1EmZmFsc2VGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRidGNS9GOVEnbm9ybWFsRic= and the noisyData together to visualize the fit.plot([noisyData, f(x)], x = -10 .. 10, style=[point, line], symbol=cross, symbolsize=15, legend=["Noisy Signal", "Fitted Model"], axes=boxed, color=["Black","Red"], linestyle=dash, thickness=2)

Ni8tJSdDVVJWRVNHNiY3ZHc3JCQhIyoqISIiJCEyWDgqKkhUPGUxJSEjOzckJCEjKSohIiIkITJDYHRBU01ndiMhIzs3JCQhIygqISIiJCEyWypbWFdTaVc+ISM7NyQkISMnKiEiIiQhMVtOeHRzLmpjISM7NyQkISMmKiEiIiQiMmsyKXk6SF1MRiEjPDckJCEjJSohIiIkIjE7KilHVCcqSCg0IyEjOjckJCEjJCohIiIkIjJ3L0c9dygqKlI4ISM7NyQkISMjKiEiIiQiMi86eGVIdGgkUSEjOzckJCEjIiohIiIkIjJVPzB4SXUqKSlSISM7NyQkISMhKiEiIiQiMTBmUHJiSTdfISM6NyQkISMqKSEiIiQiMWIiekQnPlckZiYhIzo3JCQhIykpISIiJCIxNjFfdTVHOWEhIzo3JCQhIygpISIiJCIxXy0zdkwmKlxuISM6NyQkISMnKSEiIiQiMjpvby1PakRKKCEjOzckJCEjJikhIiIkIjEiUSQpZSYqWyZleCEjOjckJCEjJSkhIiIkIjFNdm8heUlqMykhIzo3JCQhIyQpISIiJCIxY2I9QVImKiplKSEjOjckJCEjIykhIiIkIjFQaylIbGFYbSkhIzo3JCQhIyIpISIiJCIxNCF6JT0lZVoiKSkhIzo3JCQhIyEpISIiJCIxJ0hSKjNrITMkKSkhIzo3JCQhI3ohIiIkIjFuUiQqXHNrcyopISM6NyQkISN5ISIiJCIxI3BqIlIoeTFYKSEjOjckJCEjeCEiIiQiMUxkKCp5bFwpUSkhIzo3JCQhI3chIiIkIjF2KEgnKlwqPi8pKSEjOjckJCEjdiEiIiQiMTcrcC04XmwjKSEjOjckJCEjdSEiIiQiMUJwOE4lRyUqRykhIzo3JCQhI3QhIiIkIjEiSCgqM0pxYUsoISM6NyQkISNzISIiJCIyJmVKTnQ4RidwJyEjOzckJCEjciEiIiQiMUZSQ0JoRF5mISM6NyQkISNxISIiJCIxV3dMWyc+XywmISM6NyQkISNwISIiJCIxeCRIPyRIRSNRJSEjOjckJCEjbyEiIiQiMmAob2AxW2kjeSQhIzs3JCQhI24hIiIkIjEqelQnSFRLblAhIzo3JCQhI20hIiIkIjFmJHonMyRvaXUjISM6NyQkISNsISIiJCIyOVlBYFViJD5GISM7NyQkISNrISIiJCIxeTNOXyVIMG0iISM6NyQkISNqISIiJCIwNGdbXzkhZnAhIzo3JCQhI2khIiIkIjIocGgiM0k4IWU2ISM7NyQkISNoISIiJCExRDF6Z2I1VWchIzs3JCQhI2chIiIkITJUIypcZFE3JWY2ISM7NyQkISNmISIiJCEyPihHIT49Jj0pKT0hIzs3JCQhI2UhIiIkITE5J3pAIkdgJ3ApISM7NyQkISNkISIiJCEySyIqUmkkbylbWSIhIzs3JCQhI2MhIiIkITEmSCQ0Tix5WkkhIzo3JCQhI2IhIiIkITJhL28wd3hidyMhIzs3JCQhI2EhIiIkITJ2dGZKJTRZSUMhIzs3JCQhI2AhIiIkITFOejByPjxoUyEjOjckJCEjXyEiIiQhMS8iPSJvKEg/eSQhIzo3JCQhI14hIiIkITEoUm07eWkycCQhIzo3JCQhI10hIiIkITIvW0glZnVtJFIkISM7NyQkISNcISIiJCEyKEhhRlBiK3pRISM7NyQkISNbISIiJCEySyR6cGN5TURQISM7NyQkISNaISIiJCExckU1OiFmV1EkISM6NyQkISNZISIiJCExJzNVb2VFOWUkISM6NyQkISNYISIiJCEyKCpwKVxUcDU5USEjOzckJCEjVyEiIiQhMSNmcStRdiQ9TCEjOjckJCEjViEiIiQhMS9TKSplT20qRyQhIzo3JCQhI1UhIiIkITJYaiRwOVAkR28iISM7NyQkISNUISIiJCExLzgiKVFLbiVbIyEjOjckJCEjUyEiIiQhMTo3RzdERG8/ISM6NyQkISNSISIiJCEyRV0obyVIKVx1NSEjOzckJCEjUSEiIiQhMilHKXo9KmVpOz0hIzs3JCQhI1AhIiIkIjFXcW5XTWdeXCEjPTckJCEjTyEiIiQhMU5pMTwvMTxgISM7NyQkISNOISIiJCEyKSopSFVVJ2VvNyIhIz03JCQhI00hIiIkITFEU0pkbGg1OCEjOzckJCEjTCEiIiQiMWQiUlkmKXBhImYhIzs3JCQhI0shIiIkIjA5MyZwJioqMzEnISM6NyQkISNKISIiJCIyJD5baEBKbVk2ISM7NyQkISNJISIiJCIxJioqUkoyO19UIyEjOjckJCEjSCEiIiQiMiZ5XlJcR14rRyEjOzckJCEjRyEiIiQiMVFNIlx3ImV5QiEjOjckJCEjRiEiIiQiMjFAIUctUk4oZiMhIzs3JCQhI0UhIiIkIjIqKWZYTGZxJTQ+ISM7NyQkISNEISIiJCIyTmJIJDNaKEcvIyEjOzckJCEjQyEiIiQiMjI5RDM5WUY9IyEjOzckJCEjQiEiIiQiMmIrIT5lblAxRSEjOzckJCEjQSEiIiQiMUthJjNdOTFoIyEjOjckJCEjQCEiIiQiMmFjNCFlelUoMyQhIzs3JCQhIz8hIiIkIjFNZShwaj8pPkEhIzo3JCQhIz4hIiIkIjFRLUIoKlxATUIhIzo3JCQhIz0hIiIkIjFpKSlHTC8mejMjISM6NyQkISM8ISIiJCIxNl5FaUBoP0AhIzo3JCQhIzshIiIkIjAsTi4qR0ROSyEjOTckJCEjOiEiIiQiMj1adiNbWil6diIhIzs3JCQhIzkhIiIkIjJsL291JzNLJz0kISM7NyQkISM4ISIiJCIyckQ7amdhKjQ+ISM7NyQkISM3ISIiJCIxQC5ZYExsZkUhIzo3JCQhIzYhIiIkIjJvMWdla2c+ZiIhIzs3JCQhIzUhIiIkIjF5RXZsLXRqOiEjOjckJCEiKiEiIiQiMUVUW1UsPlRBISM6NyQkISIpISIiJCIydng5Ok9oJjQ7ISM7NyQkISIoISIiJCIyQUFiKD00KilmOCEjOzckJCEiJyEiIiQiMiMpRzRlRGdJbCIhIzs3JCQhIiYhIiIkIjFoJlt6WXV5dykhIzs3JCQhIiUhIiIkIjElKSkqKmZucjJcKiEjOzckJCEiJCEiIiQiMkhJKWU3UCcpZj0hIzs3JCQhIiMhIiIkIjFgNFtDOnomSCohIzs3JCQhIiIhIiIkIjFDaylHSGBkKmUhIzs3JCQiIiEhIiIkIjFCTC4xSWAjPiIhIzo3JCQiIiIhIiIkIjJPRCZmJio+JikzNiEjOzckJCIiIyEiIiQiMTNzRnJzJ0hhJyEjOzckJCIiJCEiIiQiMihRTUc4LVFKNiEjOzckJCIiJSEiIiQiMiY0b19CJEdZYiIhIzs3JCQiIiYhIiIkIjFzSFthI2YrWighIzs3JCQiIichIiIkIjJWREF1clY5PCIhIzs3JCQiIighIiIkIjIkR1dlbClwW1ciISM7NyQkIiIpISIiJCIyQ0VYJCl6JClwOyMhIzs3JCQiIiohIiIkIjIpej1JRC0wbzohIzs3JCQiIzUhIiIkIjEuJ3oob3hNc0AhIzo3JCQiIzYhIiIkIjFzY1YnMygzVkQhIzo3JCQiIzchIiIkIjIpSHQ+Iykzd00/ISM7NyQkIiM4ISIiJCIyYjsxLXIjKUhkIyEjOzckJCIjOSEiIiQiMSZmVElJS3JcIyEjOjckJCIjOiEiIiQiMidcM2MzcipHRiMhIzs3JCQiIzshIiIkIjJCXmBMWVFOXSMhIzs3JCQiIzwhIiIkIjIyVTFVVigqbypIISM7NyQkIiM9ISIiJCIxMGo+KFxHJjNGISM6NyQkIiM+ISIiJCIyOWt1Pnk3YFsjISM7NyQkIiM/ISIiJCIxNE9iKm9lUipRISM6NyQkIiNAISIiJCIyJHA/REdwdiZ6IyEjOzckJCIjQSEiIiQiMXZ6TFoqW1teJCEjOjckJCIjQiEiIiQiMmtSWUpfTVIhRyEjOzckJCIjQyEiIiQiMk91QnBFUkRVIyEjOzckJCIjRCEiIiQiMkNJLyQ0KlEvciIhIzs3JCQiI0UhIiIkIjJsbVIlXD4qUVYiISM7NyQkIiNGISIiJCIyWGZBUHJaY1kiISM7NyQkIiNHISIiJCIyMHojR2UtUyF5IiEjOzckJCIjSCEiIiQiMil6JFxDRzdcPCIhIzs3JCQiI0khIiIkIjJmYzEkKXBtJCopPiEjOzckJCIjSiEiIiQiMTEyRWVxYyFHJSEjOzckJCIjSyEiIiQiMigpcGZtRFtzdSIhIzs3JCQiI0whIiIkITFGJHo2bzg6LCkhIzw3JCQiI00hIiIkITBLS0FPdDVoKCEjOjckJCIjTiEiIiQiMVJmKDRAZStIJiEjOzckJCIjTyEiIiQhMjwqZilbJVFIJ3kiISM7NyQkIiNQISIiJCEyI3pSUiQpRyQqKVsiISM7NyQkIiNRISIiJCEyaiQ9TjgiXCRbNyEjOzckJCIjUiEiIiQhMXotPWEoW14oeiEjOzckJCIjUyEiIiQhMU4yVHBRPXQ4ISM6NyQkIiNUISIiJCEwYCpcJ1JHblAkISM5NyQkIiNVISIiJCEyRXBaOUJnJykqRyEjOzckJCIjViEiIiQhMSdSQiRwJ1FbUCQhIzo3JCQiI1chIiIkITImKj43KHolSDNHJCEjOzckJCIjWCEiIiQhMSo0WClcLFVOTSEjOjckJCIjWSEiIiQhMnRoNid6YSJRUyQhIzs3JCQiI1ohIiIkITIwYmVPK0c6XCQhIzs3JCQiI1shIiIkITImb2giPiNwa0RLISM7NyQkIiNcISIiJCExVFpDIVIvVSpIISM6NyQkIiNdISIiJCExUnZaV3VhYUohIzo3JCQiI14hIiIkITJzLHVBTSdwMFAhIzs3JCQiI18hIiIkITJ2NiV5Z3okbyVIISM7NyQkIiNgISIiJCEyYU4yNmtjKj5KISM7NyQkIiNhISIiJCExWkxFIzQ/SyZHISM6NyQkIiNiISIiJCEyOip6QzVXd21IISM7NyQkIiNjISIiJCExOSwnb1teTjIjISM6NyQkIiNkISIiJCEyKSopXGokKjMxXDwhIzs3JCQiI2UhIiIkITJ2L1UzSC95KT4hIzs3JCQiI2YhIiIkITImenBseC5qIlsjISM7NyQkIiNnISIiJCExXjxZPW08VUQhIzs3JCQiI2ghIiIkIjJOIls2bT1cIlwiISM8NyQkIiNpISIiJCEyTSVHZWAwNy5AISM8NyQkIiNqISIiJCIxTlgzSjUwMXghIzs3JCQiI2shIiIkIjIoUSYqKmVjIz1sOCEjOzckJCIjbCEiIiQiMHI4KFtoeiEzIyEjOTckJCIjbSEiIiQiMjEnZXFtNSZcZyMhIzs3JCQiI24hIiIkIjEjSDkpPiJvKFxWISM6NyQkIiNvISIiJCIyREhhRmAib1BVISM7NyQkIiNwISIiJCIxJik0c2k9YXBfISM6NyQkIiNxISIiJCIxRWojKjNiKFxhJiEjOjckJCIjciEiIiQiMVx0JT4ocGNbYyEjOjckJCIjcyEiIiQiMFglKnA4W0ljJyEjOTckJCIjdCEiIiQiME4nRytrNVduISM5NyQkIiN1ISIiJCIxJCo9cFUyJkgkeiEjOjckJCIjdiEiIiQiMXJXJz0lb2FGeiEjOjckJCIjdyEiIiQiMFV6Jz5qdXQkKiEjOTckJCIjeCEiIiQiMUUhPlYpRzQreiEjOjckJCIjeSEiIiQiMT4iUmpqJCo9RSkhIzo3JCQiI3ohIiIkIjF5LFkqKTNHSygqISM6NyQkIiMhKSEiIiQiMWJpK0J5QlAqKSEjOjckJCIjIikhIiIkIjByaiw3XFVLKSEjOTckJCIjIykhIiIkIjFuWShRPyUpKnkhKiEjOjckJCIjJCkhIiIkIjFvV193NC1pKCkhIzo3JCQiIyUpISIiJCIxalFKUiYpZSc9KSEjOjckJCIjJikhIiIkIjBwQTlPKlxDeSEjOTckJCIjJykhIiIkIjFMOHBSKFsjenMhIzo3JCQiIygpISIiJCIxUHkjKnpzYkVsISM6NyQkIiMpKSEiIiQiMTEkZiFmUSRSIXAhIzo3JCQiIyopISIiJCIyJT0sJmUmR2Y5XyEjOzckJCIjISohIiIkIjEiKipHIlJYRlJdISM6NyQkIiMiKiEiIiQiMXAlXFUkPi4kKlshIzo3JCQiIyMqISIiJCIyanBdUWAqKVx6JCEjOzckJCIjJCohIiIkIjJadndUQnlOdyIhIzs3JCQiIyUqISIiJCIya0tcK2RuKVE4ISM7NyQkIiMmKiEiIiQhMkY/a14rWHYhSCEjPTckJCIjJyohIiIkITIjekArJkdCNjAiISM7NyQkIiMoKiEiIiQhMXl3OVlMUig9KSEjOzckJCIjKSohIiIkITIjPT11JSpcMW1FISM7NyQkIiMqKiEiIiQhMm0leTlsJ1wyNSQhIzs3JCQiJCsiISIiJCExOz44TydHTjIlISM6LSUnTEVHRU5ERzYjLSUpX1RZUEVTRVRHNiMtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjUtSSNtc0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjYjUS1Ob2lzeX5TaWduYWw2Ii8lJ2ZhbWlseUdRMFRpbWVzfk5ld35Sb21hbjYiLyUlc2l6ZUdRIzEyNiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJmZhbHNlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1Enbm9ybWFsNiItJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lJlNUWUxFRzYjJSZQT0lOVEctJSdDVVJWRVNHNiY3ank3JCQhJCsiISIiJCExPV4iUSIqM2lGJSEjOjckJCEpeiVbKikqISIoJCExaVJTJ1xCLEskISM6NyQkISs+Y04uKSohIiokITJqI1J3XjQsbEMhIzs3JCQhKyIpSFkrKCohIiokITIncFVVPTdoKFsiISM7NyQkISsuJCkpb2YqISIqJCExXFAsJCkqUmcoXCEjOzckJCErKSplIVFcKiEiKiQiMXRPUXV5IylHWyEjOzckJCErOGBCKVIqISIqJCIycyhvMlRwTnk4ISM7NyQkISszdkYqSCohIiokIjBCb3lZJClHRyMhIzk3JCQhMSoqKioqXD9OcD4qISM6JCIxdUEiPik+KFw9JCEjOjckJCErOTYjXDQqISIqJCIwJSopRydbTzovJSEjOTckJCErJDMoKSoqKSopISIqJCIxalgnW2knW3BbISM6NyQkIStKNmMoKikpISIqJCIxJyk9JzRCT3dhJiEjOjckJCErLDleJHopISIqJCIxWiRwVSVbQ1ppISM6NyQkIStLVy4qbykhIiokIjEpR05LbzpjKG8hIzo3JCQhKzg6TillKSEiKiQiMT0xLiFvaWVTKCEjOjckJCEreDwjcFwpISIqJCIyJm85SHUqKikqPnkhIzs3JCQhK3lIPylRKSEiKiQiMHlrcyQ+c0QjKSEjOTckJCErJipRNSdIKSEiKiQiMSVvN2w5eFFcKSEjOjckJCErbTIoKik9KSEiKiQiMV5yS0dbTjwoKSEjOjckJCEqR2I6OSkhIikkIjFgISlHRUZ6JnkpISM6NyQkISslelJUNCkhIiokIjFjQHNARGVOKSkhIzo3JCQhLERrRyJvISkhIzUkIjFha2NwbytiKSkhIzo3JCQhMSoqKioqNFw8QC8pISM6JCIxd3FeR3QnKW8pKSEjOjckJCExKioqKlxSajU7ISkhIzokIjF0JFtcdyk9eCkpISM6NyQkISspPSY0ISp6ISIqJCIxc0JeVCoqKip6KSkhIzo3JCQhMSoqKlwyVUVgJ3ohIzokIjA3dHYoW2V4KSkhIzk3JCQhLE5sZDAleiEjNSQiMWpCMyRHTS0oKSkhIzo3JCQhMjArXWkpKSl5OnohIzskIjE+XSxERSl6JikpISM6NyQkISs+LC0iKnkhIiokIjFRKlxaVW4zJSkpISM6NyQkISwmKSpHTFJ5ISM1JCIxZWJtaihvJyp5KSEjOjckJCEreWNrKHkoISIqJCIxTUxuYzEpeXIpISM6NyQkIStyY3IjcCghIiokIjFwJHBoOCMpUmApISM6NyQkIStMM0shZighIiokIjA3NCc+NFNqIykhIzk3JCQhK3Y4J1JbKCEiKiQiMSxvSzRcNzN6ISM6NyQkISthXlAiUighIiokIjBpSWhnIjNVdiEjOTckJCEyJioqKioqelx6OEgoISM7JCIxUDArQ2cpRzQoISM6NyQkITEqKioqKj4/dSEpPSghIzokIjE2USRHJGVdd2whIzo3JCQhK0MsLCgzKCEiKiQiMSc0LSpRa3hFZyEjOjckJCEyJioqKioqNFZEIyopcCEjOyQiMSkzJTQ0SHBmYSEjOjckJCErLUFsISlvISIqJCIxQzo/aHhtKHolISM6NyQkIStFVjQkeSchIiokIjEyTFMoKVE5Ij0lISM6NyQkITEqKioqKjRaSSp5bSEjOiQiMlpUTSM0LD0zTiEjOzckJCErT0RhJWUnISIqJCIyOlVgJD4zKyMqRyEjOzckJCExKioqKipIZl44WychIzokIjJGTmAkKWZqKD1BISM7NyQkIStjcEQlUSchIiokIjImUXIicHVaQmYiISM7NyQkISsmM2xGRychIiokIjB2ajRkZG1eKiEjOjckJCEqQlFOPSchIikkIjFcX2xnbjpiTSEjOzckJCEyMCsrITNwbHpnISM7JCExLXUhM28pUjdFISM7NyQkIStcb2d6ZiEiKiQhMTwxOXYqUi44KSEjOzckJCErODdIeGUhIiokITI3OHVdIiopKSpROCEjOzckJCErLEcjZXgmISIqJCEyJ29SQlghenciPSEjOzckJCEqKFFlI28mISIpJCEyYi4vN0NKa0AjISM7NyQkIStLOnN2YiEiKiQhMT4lZXBYIWVARSEjOjckJCErJGVVLFsmISIqJCEyYCNwSjYkeVgkSCEjOzckJCErKltOI3lgISIqJCExS1IqSEssXUAkISM6NyQkISsnUihwIUcmISIqJCEyOysycHIocElNISM7NyQkITInKioqKioqelFSc14hIzskITE9NzV3dzk0TyEjOjckJCErT0NdRF4hIiokITFjJXB3aiJbbU8hIzo3JCQhKyMqNGh5XSEiKiQhMk1xI0cvWiU9ciQhIzs3JCQhKyV5KVFEXSEiKiQhMTpJPycpSCUqW1AhIzo3JCQhK3dsO3NcISIqJCExSj53W1AiNHgkISM6NyQkITEsK11WTSJ6JVwhIzokITFKcWh2didmeCQhIzo3JCQhMSwrKzYubUJcISM6JCEyeFFIc2FdengkISM7NyQkITEsK115clMqKlshIzokITFIQEIoMyopb3gkISM6NyQkITIwKytnL2FeKFshIzskITFZUCVRYTZHeCQhIzo3JCQhLGIiZTRBWyEjNSQhMjFoa3BZQ052JCEjOzckJCEyJioqKioqXGVQIXBaISM7JCEwWWN4SShIP1AhIzk3JCQhKzB0REJaISIqJCEyJ3lxVjk5cyFvJCEjOzckJCErRHFaeFkhIiokITImcDFITiNHOGokISM7NyQkIStOLT90WCEiKiQhMm5NdyE0dnckWyQhIzs3JCQhKy9TWnNXISIqJCExYmwhXCJwJHpIJCEjOjckJCErdk4iPVAlISIqJCExYWQyXDppdEkhIzo3JCQhK2RNX3JVISIqJCEybmhjIlteImYiRyEjOzckJCErXmc8dlQhIiokITJudmFEYSdHU0QhIzs3JCQhK2E+LXJTISIqJCExJEhhR3EzaUAjISM6NyQkIStaYnZyUiEiKiQhMCd5KyYqZSdvKT0hIzk3JCQhK0ZCR25RISIqJCEyR1YiUThMJFJfIiEjOzckJCErMSozRngkISIqJCEyeXZQLDAlXCY9IiEjOzckJCErc3Y9b08hIiokITFZP0dpXWtiISkhIzs3JCQhK1lWMm9OISIqJCEybE1BOF9IbFMlISM8NyQkISlbP29NISIoJCExJG9NS1FwVS0pISM8NyQkIStyMSpRTyQhIiokIjJZYicpNGp2cChHISM8NyQkISpiOHlFJCEiKSQiMTYxODJpJipcaCEjOzckJCErUHVUcEohIiokIjFcNVsiKikqPV4kKiEjOzckJCErT2t4Z0khIiokIjI3ZDdvJTMtbjchIzs3JCQhKzM0VGlIISIqJCIyRnY4YSZ5I1xhIiEjOzckJCErSlAiPSdHISIqJCIyb09HS0gnSC89ISM7NyQkISskKiozJmZGISIqJCIyJjQjeTxgMCwvIyEjOzckJCErSjVabEUhIiokIjFgMXdma1pJQSEjOjckJCErNCdIX2MjISIqJCIxJjNoYzBKV1MjISM6NyQkIStKN3ZsQyEiKiQiMmAraUYwKXpZRCEjOzckJCErKkdTJGZCISIqJCIydG9DW2tHYW0jISM7NyQkIStUc1ZsQSEiKiQiMihITWtRJzQ5dSMhIzs3JCQhLCZIXEU2QSEjNSQiMVhWUUteP3RGISM6NyQkISs9RTRkQCEiKiQiMiRvV1ciKnBQJ3ojISM7NyQkISxORzEjM0AhIzUkIjAoW0lJQy81RyEjOTckJCErXCo+JGY/ISIqJCIyMitGNChHKzxHISM7NyQkISt3JVw0LCMhIiokIjJWKVtvQUFaPEchIzs3JCQhKy4hekQnPiEiKiQiMHdSejRlPCJHISM5NyQkITItKytiP2wwIj4hIzskIjF2aD9KVScqKXojISM6NyQkISszOWJlPSEiKiQiMiY+Xl9oIW8mekYhIzs3JCQhKC5UdiIhIickIjJXejR4VyEpPnMjISM7NyQkISs5LzJmOyEiKiQiMjxaRSoqUV4qXEUhIzs3JCQhK0grOmU6ISIqJCIxL29vZ1RtYkQhIzo3JCQhKyhmaCdmOSEiKiQiMiVmLjhrInkoW0MhIzs3JCQhKyg0KipHTiIhIiokIjIoM04nUjxPKD5CISM7NyQkIStfKT4wRSIhIiokIjEwKD5eaCQqKSo+IyEjOjckJCEqJCo9VzoiISIpJCIyeDt3RkcyajAjISM7NyQkISsnW3lZMCIhIiokIjJqPjUoKjMlej0+ISM7NyQkISpFOSFmJiohIiokIjJHY2xdRThMeSIhIzs3JCQhKil5WHEmKSEiKiQiMkJfVCk+U09eOyEjOzckJCEqJUg+aXYhIiokIjIpUipwQUVBTV8iISM7NyQkIShcTVwnISIoJCIyLUVCYjlxIilSIiEjOzckJCEqWTpOXiYhIiokIjJpQEcoUi8/JkgiISM7NyQkISpCJj5ZWCEiKiQiMUUmZSEpPl5uPyIhIzo3JCQhKiNlOTdOISIqJCIxelsvPiUzJ0c2ISM6NyQkISpYWi5bIyEiKiQiMHJGMml1InA1ISM5NyQkISo5WlViIiEiKiQiMiJcaTY3VGxLNSEjOzckJCEpYTlzWCEiKiQiMVY1LE00KDMsIiEjOjckJCIoOnRnJSEiKSQiMid6TV4/SCE0LCIhIzs3JCQiKUNQTTohIikkIjImWzFDKy4wSzUhIzs3JCQiKShlIWZEISIpJCIxKG8tOSFbLXQ1ISM6NyQkIik6KVJaJCEiKSQiMl8iSCV6OnJnNyIhIzs3JCQiKXohSF0lISIpJCIyJVtyVDF4Ny43ISM7NyQkIilabFFiISIpJCIySzRkTiE+byhIIiEjOzckJCIpKHklcGwhIikkIjIsX2N5I1ttMTkhIzs3JCQiKVk9RHYhIikkIjJYayYpR20qKSk9OiEjOzckJCIpRnc5JikhIikkIjFaRSNHKWY0VzshIzo3JCQiKWQ9USYqISIpJCIyMioqR29gJlsheSIhIzs3JCQiKm1LZTAiISIpJCIyODgrPDkoUT8+ISM7NyQkIipwbTI7IiEiKSQiMm0lUidHVzZdMSMhIzs3JCQiKmsjPmA3ISIpJCIyZCxUP1BWLD4jISM7NyQkIipRVXNOIiEiKSQiMlZrJGUmZjFfSyMhIzs3JCQiKk0+PFkiISIpJCIyYiZlIyl6PTleQyEjOzckJCIqRS1DYyIhIikkIjI4bnlnPXAqZkQhIzs3JCQiKio+JFFsIiEiKSQiMUBgOnJZW1hFISM6NyQkIikzYmk8ISIoJCIxWy5MaFZeRkYhIzo3JCQiKigpXFkmPSEiKSQiMk9iVD9vYnl4IyEjOzckJCIqVzsjMz4hIikkIjEiMyIzKyVHIyl6IyEjOjckJCIqLCR5aD4hIikkIjElNF5EYjg7IkchIzo3JCQiKlwpPjQ/ISIpJCIybjc8NCk+UDxHISM7NyQkIiooUmhjPyEiKSQiMVhYLDlmPjxHISM6NyQkIipHTyczQCEiKSQiMXgvcjo8Jio0RyEjOjckJCIqZmUxOyMhIikkIjEjeXdSQjheeiMhIzo3JCQiKjcnPjVBISIpJCIyRlZMUThYUHgjISM7NyQkIipsTChmQSEiKSQiMV5xRU8kb151IyEjOjckJCIpIjNKTyMhIigkIjJEOjApcDUjPW0jISM7NyQkIio2USFlQyEiKSQiMkU1I0g5bWNjRCEjOzckJCIqJUhWZ0QhIikkIjJkeUAlPjkqPlQjISM7NyQkIilDem1FISIoJCIyMFcpW2I8KXpBIyEjOzckJCIqaXkkZkYhIikkIjJOKFJ3XXNRUz8hIzs3JCQiKkd1JGZHISIpJCIwPShSeF5ENT0hIzk3JCQiKmR6RSdIISIpJCIyRlQ2blwrVWEiISM7NyQkIipsVlAxJCEiKSQiMiY9I3BARipHZTchIzs3JCQiKk5HOjskISIpJCIxQjZ6JlJvLWcqISM7NyQkIipkLCxGJCEiKSQiMSZ6NjJBdU4yJyEjOzckJCIqV2Z3TyQhIikkIjFuIipIbXUyWUYhIzs3JCQiKUwjPVokISIoJCExbjpYXXMjKT0kKiEjPDckJCIqQzZpYyQhIikkITEiUSIzQWQjKVFWISM7NyQkIio9LSVwTyEiKSQhMXV6I1woXCYqKjQpISM7NyQkIiojb1xtUCEiKSQhMnR0OFNgUEk7IiEjOzckJCIqcCkpeidRISIpJCEyT1xbOExKa18iISM7NyQkIiphOnMnUiEiKSQhMl8/ZSZlUlNyPSEjOzckJCIqKG80clMhIikkITJDSmwoSC5YO0EhIzs3JCQiKiNwOXJUISIpJCEvLF1cRkFHRCEjODckJCIqY2lNRiUhIikkITFAQiE+TCU+QEchIzo3JCQiKig0JFxQJSEiKSQhMjoxdCdScTUiMyQhIzs3JCQiKSpwIm9XISIoJCEwWVxqZDEiKkckISM5NyQkIipDS11kJSEiKSQhMmtsWSwnXHcnWyQhIzs3JCQiKj42MW4lISIpJCExVCc9JlI9NEJPISM6NyQkIip1azpzJSEiKSQhMTBmci0kbyF6TyEjOjckJCIqSD1EeCUhIikkITEnXChHYG8qR3MkISM6NyQkIitOdEdAWyEiKiQhMTstPz09N2BQISM6NyQkIipRYysoWyEiKSQhMSUqUTk0UGRyUCEjOjckJCIxLCtdZEE4KCpbISM6JCEyUHhBdk1NbXgkISM7NyQkIjEsKytOIjNVI1whIzokITI8KSpwaFlSengkISM7NyQkIixELCVHXlwhIzUkITIxcyEqKkcmXGF4JCEjOzckJCIqKilmJHlcISIpJCExaXdVISp6N3BQISM6NyQkIitOOEREXSEiKiQhMmtFKVtCIj4hXFAhIzs3JCQiKnlVQDImISIpJCEybEBPKFJzOzxQISM7NyQkIioqXE9EXiEiKSQhMkxjZSpmP2ptTyEjOzckJCIpc2V5XiEiKCQhMmJVakFrdzFnJCEjOzckJCIqcypmdl8hIikkITJ5Sk4vLlMwVyQhIzs3JCQiKj47PFEmISIpJCExRSlHUTxZaj8kISM6NyQkIip2d0taJiEiKSQhMW55Kilwb0BiSCEjOjckJCIqYGB2ZCYhIikkITFAKTNZdz1eaCMhIzo3JCQiKnd6I3ljISIpJCEyRF8uNyczJVFCIyEjOzckJCIpLSUqeWQhIigkITJqYjEwTVxPIT0hIzs3JCQiKktJI3plISIpJCEyJ1E+a3c0VEg4ISM7NyQkIipzeGIoZiEiKSQhMXorOURINlgkKSEjOzckJCIqIj10emchIikkITFfPEQ8NjkzRSEjOzckJCIqQikqKnloISIpJCIyYSc9bS9kZSQ9JCEjPDckJCIqV3JNRychIikkIjFrcU5RdGZnJiohIzs3JCQiKihbL3lqISIpJCIyVydlZSQ9X0ViIiEjOzckJCIpaWMjWychIigkIjFUV0dWcGxFQSEjOjckJCIqVXpFZSchIikkIjI6InlsaC0jKXpHISM7NyQkIiooKltEbychIikkIjJGV0pkaUM8YCQhIzs3JCQiKUonb3knISIoJCIxMVtFdFhEMFUhIzo3JCQiKkJTSClvISIpJCIxVi9XXDAiPiJbISM6NyQkIipNTzgpcCEiKSQiMVkiKipbN09FVCYhIzo3JCQiKkx4KiozKCEiKSQiMWIkKnpPdVpWZyEjOjckJCIqJ0dNKT0oISIpJCIxWzhaNjsiemQnISM6NyQkIiowUyopRyghIikkIjJOLSgpUUZ5NzMoISM7NyQkIip5VzdSKCEiKSQiMEdZdTpKOmEoISM5NyQkIip1I0cmWyghIikkIjImcGdEPCV5SCJ6ISM7NyQkIio8Q2JlKCEiKSQiMTptPmBVJypbIykhIzo3JCQiKmEtXW8oISIpJCIxS1dkYnM6OyYpISM6NyQkIipbODl6KCEiKSQiMShvW1ZyKnpCKCkhIzo3JCQiKzBdT1F5ISIqJCIwYyY9bEFeKXkpISM5NyQkIipgO2ApeSEiKSQiMTcqenlgVGkkKSkhIzo3JCQiKyZvLUMieiEiKiQiMVkrSSEqKUdmJikpISM6NyQkIiolKSlbUnohIikkIjBeXWtfMSlwKSkhIzk3JCQiKyYqXGRteiEiKiQiMT4jKik+LUR5KCkpISM6NyQkIio6aE8qeiEiKSQiMTohXCI9JFIqeikpISM6NyQkIipLLyI9ISkhIikkIjAwNnRMWG4oKSkhIzk3JCQiKlxaRC8pISIpJCIxIypHTShIJG9vKSkhIzo3JCQiKm0hKnAxKSEiKSQiMSc0JFx6JUhkJikpISM6NyQkIiokUVYiNCkhIikkIjBRZ1I/ankkKSkhIzk3JCQiKVYhKVIiKSEiKCQiMWwhXHFiaXp5KSEjOjckJCIqeHUiKT0pISIpJCIxbm1Jb3psPSgpISM6NyQkIipQLUFIKSEiKSQiMTZKVyEzKm8uJikhIzo3JCQiKnhdbVIpISIpJCIxNEdeV0RqKD4pISM6NyQkIipPJG8iXCkhIikkIjElel1TSEA8JXkhIzo3JCQiKnYuRWYpISIpJCIxbEBQcyE+XVEoISM6NyQkIio9IzQicCkhIikkIjB6IzNBQipSJ28hIzk3JCQiKm9heXopISIpJCIxdltUJilSWj5pISM6NyQkIiojUkIhKikpISIpJCIxI0dPKVJnPSpmJiEjOjckJCIqJVtMJyoqKSEiKSQiMTxIb0dgM0BbISM6NyQkIipHdmc0KiEiKSQiMiYpZl1PJVw2S1MhIzs3JCQiKk5fWz4qISIpJCIyVkA5d3g6SD8kISM7NyQkIioqenEkSCohIikkIjIwZVMlMy0kSEwjISM7NyQkIipbTVhSKiEiKSQiMm9bOkhNYEVUIiEjOzckJCIqKCkzOV0qISIpJCIxa3daazxsNFQhIzs3JCQiKkItJSpmKiEiKSQhMW53JGVYJz47XyEjOzckJCIqRE1ocCohIikkITJtQCR5bTtHWTkhIzs3JCQiKj5SJip6KiEiKSQhMiVRRmVPWSoqR0MhIzs3JCQiKi0+RiEqKiEiKSQhMWFnJSozdCFHUiQhIzo3JCQiJCsiISIiJCExPV4iUSIqM2lGJSEjOi0lJ0xFR0VOREc2Iy0lKV9UWVBFU0VURzYjLUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1LUkjbXNHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2I1EtRml0dGVkfk1vZGVsNiIvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW42Ii8lJXNpemVHUSMxMjYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSZmYWxzZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ25vcm1hbDYiLSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJSZTVFlMRUc2IyUlTElORUctJSVWSUVXRzYkOyQhJCsiISIiJCIkKyIhIiIlKERFRkFVTFRHLSYlJl9BWElTRzYjIiIiNictJStfR1JJRExJTkVTRzYjLSUpX1ZJU0lCTEVHNiMiIiEtJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIS0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJiUmX0FYSVNHNiMiIiM2Jy0lK19HUklETElORVNHNiMtJSlfVklTSUJMRUc2IyIiIS0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqTElORVNUWUxFRzYjIiIhLSUqVEhJQ0tORVNTRzYjIiIhLSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0lK0FYRVNMQUJFTFNHNiQtSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1USJ4NiIvJSdmYW1pbHlHUSE2Ii8lJXNpemVHUSMxMDYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiJRITYiLSUqTElORVNUWUxFRzYjIiIkLSUnU1lNQk9MRzYkJSZDUk9TU0ciIzotJSpUSElDS05FU1NHNiMiIiMtJSpBWEVTU1RZTEVHNiMlJEJPWEctJSlfVklTSUJMRUc2IyIiIi0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJD8lISIiLSUpQk9VTkRTX1lHNiMkIiRJIiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJWdNISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJWdKISIiLSUpQ0hJTERSRU5HNiItJStBTk5PVEFUSU9ORzYnLSUpQk9VTkRTX1hHNiMkIiIhISIiLSUpQk9VTkRTX1lHNiMkIiIhISIiLSUtQk9VTkRTX1dJRFRIRzYjJCIlK1MhIiItJS5CT1VORFNfSEVJR0hURzYjJCIlK1MhIiItJSlDSElMRFJFTkc2Ig==Ig==

Exercise: Repeat the noisy signal example using a random number generator based on a Weibull distribution.

The following example adjusts a fitted model by adding or removing variables in order to find better adjusted r-squared values.restart; with(Statistics): The Import command reads a datafile into a new DataFrame.ExperimentalData := Import(FileTools:-JoinPath(["Excel", "ExperimentalData.xls"], base = datadir));

LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVExRXhwZXJpbWVudGFsRGF0YUYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RKn4mIzg3ODg7fkYnL0YzUSdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGPS8lKXN0cmV0Y2h5R0Y9LyUqc3ltbWV0cmljR0Y9LyUobGFyZ2VvcEdGPS8lLm1vdmFibGVsaW1pdHNHRj0vJSdhY2NlbnRHRj0vJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZMLUkobWFjdGlvbkdGJDYlLUkobWZlbmNlZEdGJDYoLUYjNiYtSSdtdGFibGVHRiQ2Pi1JJG10ckdGJDYtLUkkbXRkR0YkNigtRiw2I1EhRicvJSlyb3dhbGlnbkdGXG8vJSxjb2x1bW5hbGlnbkdGXG8vJStncm91cGFsaWduR0Zcby8lKHJvd3NwYW5HUSIxRicvJStjb2x1bW5zcGFuR0Zlby1GaG42KC1GLDYlUSVUaW1lRidGL0YyRl1vRl9vRmFvRmNvRmZvLUZobjYoLUYsNiVRMENvbmNlbnRyYXRpb25+QUYnRi9GMkZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GLDYlUTJDb25jZW50cmF0aW9ufkNPMkYnRi9GMkZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GLDYlUTFDb25jZW50cmF0aW9ufkFURidGL0YyRl1vRl9vRmFvRmNvRmZvLUZobjYoLUYsNiVRMENvbmNlbnRyYXRpb25+QkYnRi9GMkZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GLDYlUTBDb25jZW50cmF0aW9ufkNGJ0YvRjJGXW9GX29GYW9GY29GZm8tRmhuNigtRiw2JVEyQ29uY2VudHJhdGlvbn5BLkJGJ0YvRjJGXW9GX29GYW9GY29GZm9GXW9GX29GYW8tRmVuNi0tRmhuNigtSSNtbkdGJDYkRmVvRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRIzAuRidGOUZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GYHI2JFEoMC41NDM5OUYnRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRKDAuMDAxMjFGJ0Y5Rl1vRl9vRmFvRmNvRmZvLUZobjYoLUZgcjYkUScwLjAwNDhGJ0Y5Rl1vRl9vRmFvRmNvRmZvLUZobjYoLUZgcjYkUSgwLjAwNzIyRidGOUZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GYHI2JFEoMC4wMDA3OUYnRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRKDAuMzQ0MjNGJ0Y5Rl1vRl9vRmFvRmNvRmZvRl1vRl9vRmFvLUZlbjYtLUZobjYoLUZgcjYkUSIyRidGOUZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GYHI2JFElMjAuMEYnRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRKDAuMjc4MzdGJ0Y5Rl1vRl9vRmFvRmNvRmZvLUZobjYoLUZgcjYkUSgwLjAwMDgyRidGOUZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GYHI2JFEoMC4xMDMzOUYnRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRJzAuMDAzN0YnRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRKDAuMDE0MzRGJ0Y5Rl1vRl9vRmFvRmNvRmZvLUZobjYoLUZgcjYkUSgwLjA5MTYyRidGOUZdb0Zfb0Zhb0Zjb0Zmb0Zdb0Zfb0Zhby1GZW42LS1GaG42KC1GYHI2JFEiM0YnRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRJTQwLjBGJ0Y5Rl1vRl9vRmFvRmNvRmZvLUZobjYoLUZgcjYkUSgwLjIzODEyRidGOUZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GYHI2JFEoMC4wMDEwOEYnRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRKDAuMTI1MDNGJ0Y5Rl1vRl9vRmFvRmNvRmZvLUZobjYoLUZgcjYkUSgwLjAwMTI0RidGOUZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GYHI2JFEoMC4wMTk2M0YnRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRKDAuMDM4MDRGJ0Y5Rl1vRl9vRmFvRmNvRmZvRl1vRl9vRmFvLUZlbjYtLUZobjYoLUZgcjYkUSI0RidGOUZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GYHI2JFElNjAuMEYnRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRKDAuMjUxMjNGJ0Y5Rl1vRl9vRmFvRmNvRmZvLUZobjYoLUZgcjYkUSgwLjAwMTE2RidGOUZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GYHI2JFEoMC4xMzk5N0YnRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRKDAuMDAwODdGJ0Y5Rl1vRl9vRmFvRmNvRmZvLUZobjYoLUZgcjYkUTUwLjAyMTMwNjk3MTE5ODAwMzk1OEYnRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRKDAuMDE4ODZGJ0Y5Rl1vRl9vRmFvRmNvRmZvRl1vRl9vRmFvLUZlbjYtLUZobjYoLUZgcjYkUSI1RidGOUZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GYHI2JFElODAuMEYnRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRKDAuMTc1MTZGJ0Y5Rl1vRl9vRmFvRmNvRmZvRmp6LUZobjYoLUZgcjYkUSgwLjE1MDE5RidGOUZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GYHI2JFEoMC4wMDA4MUYnRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRNTAuMDIxOTM2NDA3NzU1MDA1ODM0RidGOUZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GYHI2JFEoMC4wMDc0OEYnRjlGXW9GX29GYW9GY29GZm9GXW9GX29GYW8tRmVuNi0tRmhuNigtRmByNiRRIjZGJ0Y5Rl1vRl9vRmFvRmNvRmZvLUZobjYoLUZgcjYkUSYxMDAuMEYnRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRKDAuMTY0MTdGJ0Y5Rl1vRl9vRmFvRmNvRmZvLUZobjYoLUZgcjYkUSgwLjAwMTE3RidGOUZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GYHI2JFEnMC4xNTE5RidGOUZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GYHI2JFEoMC4wMDAxNEYnRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRMzAuMDIyMTg3NDA5ODM2MzUxNEYnRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRKDAuMDAxMTVGJ0Y5Rl1vRl9vRmFvRmNvRmZvRl1vRl9vRmFvLUZlbjYtLUZobjYoLUZgcjYkUSI3RidGOUZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GYHI2JFEmMTIwLjBGJ0Y5Rl1vRl9vRmFvRmNvRmZvLUZobjYoLUZgcjYkUSgwLjE1NjAyRidGOUZdb0Zfb0Zhb0Zjb0Zmb0ZpX2wtRmhuNigtRmByNiRRKDAuMTYzMDVGJ0Y5Rl1vRl9vRmFvRmNvRmZvLUZobjYoLUZgcjYkUSgwLjAwMDIwRidGOUZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GYHI2JFE0MC4wMjIzNzg4NDEyNDExNDY1MUYnRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRKDAuMDA5NTRGJ0Y5Rl1vRl9vRmFvRmNvRmZvRl1vRl9vRmFvLUZlbjYtLUZobjYoLUZgcjYkUSI4RidGOUZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GYHI2JFEmMTQwLjBGJ0Y5Rl1vRl9vRmFvRmNvRmZvLUZobjYoLUZgcjYkUSgwLjEyOTU5RidGOUZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GYHI2JFEoMC4wMDEyMkYnRjlGXW9GX29GYW9GY29GZm8tRmhuNigtRmByNiRRKDAuMTYwMTVGJ0Y5Rl1vRl9vRmFvRmNvRmZvLUZobjYoLUZgcjYkUSgwLjAwMDU3RidGOUZdb0Zfb0Zhb0Zjb0Zmby1GaG42KC1GYHI2JFE1MC4wMjI1MDY4NTIxNzY1NzQxMjhGJ0Y5Rl1vRl9vRmFvRmNvRmZvLUZobjYoLUZgcjYkUSgwLjAwMTAxRidGOUZdb0Zfb0Zhb0Zjb0Zmb0Zdb0Zfb0Zhby1GZW42LS1GaG42KC1GLDYlUSQuLi5GJ0YvRjJGXW9GX29GYW9GY29GZm9GY2ZsRmNmbEZjZmxGY2ZsRmNmbEZjZmxGY2ZsRl1vRl9vRmFvLyUmYWxpZ25HUSVheGlzRicvRl5vUSliYXNlbGluZUYnL0Zgb1EnY2VudGVyRicvRmJvUSd8ZnJsZWZ0fGhyRicvJS9hbGlnbm1lbnRzY29wZUdGMS8lLGNvbHVtbndpZHRoR1ElYXV0b0YnLyUmd2lkdGhHRmVnbC8lK3Jvd3NwYWNpbmdHUSYxLjBleEYnLyUuY29sdW1uc3BhY2luZ0dRJjAuOGVtRicvJSlyb3dsaW5lc0dRJW5vbmVGJy8lLGNvbHVtbmxpbmVzR0ZgaGwvJSZmcmFtZUdGYGhsLyUtZnJhbWVzcGFjaW5nR1EsMC40ZW1+MC41ZXhGJy8lKmVxdWFscm93c0dGPS8lLWVxdWFsY29sdW1uc0dGPS8lLWRpc3BsYXlzdHlsZUdGPS8lJXNpZGVHUSZyaWdodEYnLyUwbWlubGFiZWxzcGFjaW5nR0ZdaGwvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXUYnLyUpcmVhZG9ubHlHRj1GOUY5L0krbXNlbWFudGljc0dGJFEnTWF0cml4RicvJSVvcGVuR1EiW0YnLyUmY2xvc2VHUSJdRidGaGlsLyUrYWN0aW9udHlwZUdRLnJ0YWJsZWFkZHJlc3NGJy8lKXJ0YWJsZWlkR1E1MTg0NDY3NDQwNzQzMzM0Mjc3MDJGJ0ZjaWxGZmlsRjk=

The convert command converts the values in the DataFrame to a Matrix of values.ExMat := convert( ExperimentalData, Matrix )[..,2..7];

LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEmRXhNYXRGJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRictSSNtb0dGJDYtUSp+JiM4Nzg4O35GJy9GM1Enbm9ybWFsRicvJSZmZW5jZUdRJmZhbHNlRicvJSpzZXBhcmF0b3JHRj0vJSlzdHJldGNoeUdGPS8lKnN5bW1ldHJpY0dGPS8lKGxhcmdlb3BHRj0vJS5tb3ZhYmxlbGltaXRzR0Y9LyUnYWNjZW50R0Y9LyUnbHNwYWNlR1EmMC4wZW1GJy8lJ3JzcGFjZUdGTC1JKG1hY3Rpb25HRiQ2JS1JKG1mZW5jZWRHRiQ2Ji1JJ210YWJsZUdGJDY4LUkkbXRyR0YkNiYtSSRtdGRHRiQ2KC1GIzYmLUYsNiVRKX4xMX54fjZ+RidGL0YyLUYsNiVRJ01hdHJpeEYnRi9GMi8lK2ZvcmVncm91bmRHUShbMCwwLDBdRidGOS8lKXJvd2FsaWduR1EhRicvJSxjb2x1bW5hbGlnbkdGZW8vJStncm91cGFsaWduR0Zlby8lKHJvd3NwYW5HUSIxRicvJStjb2x1bW5zcGFuR0ZccEZjb0Zmb0Zoby1GWTYmLUZmbjYoLUYjNiYtRiw2JVEsRGF0YX5UeXBlOn5GJ0YvRjItRiw2JVEpYW55dGhpbmdGJ0YvRjJGYG9GOUZjb0Zmb0Zob0Zqb0ZdcEZjb0Zmb0Zoby1GWTYmLUZmbjYoLUYjNiYtRiw2JVEqU3RvcmFnZTp+RidGL0YyLUYsNiVRLHJlY3Rhbmd1bGFyRidGL0YyRmBvRjlGY29GZm9GaG9Gam9GXXBGY29GZm9GaG8tRlk2Ji1GZm42KC1GIzYmLUYsNiVRKE9yZGVyOn5GJ0YvRjItRiw2JVEuRm9ydHJhbl9vcmRlckYnRi9GMkZgb0Y5RmNvRmZvRmhvRmpvRl1wRmNvRmZvRmhvLyUmYWxpZ25HUSVheGlzRicvRmRvUSliYXNlbGluZUYnL0Znb1EnY2VudGVyRicvRmlvUSd8ZnJsZWZ0fGhyRicvJS9hbGlnbm1lbnRzY29wZUdGMS8lLGNvbHVtbndpZHRoR1ElYXV0b0YnLyUmd2lkdGhHRmBzLyUrcm93c3BhY2luZ0dRJjEuMGV4RicvJS5jb2x1bW5zcGFjaW5nR1EmMC44ZW1GJy8lKXJvd2xpbmVzR1Elbm9uZUYnLyUsY29sdW1ubGluZXNHRlt0LyUmZnJhbWVHRlt0LyUtZnJhbWVzcGFjaW5nR1EsMC40ZW1+MC41ZXhGJy8lKmVxdWFscm93c0dGPS8lLWVxdWFsY29sdW1uc0dGPS8lLWRpc3BsYXlzdHlsZUdGPS8lJXNpZGVHUSZyaWdodEYnLyUwbWlubGFiZWxzcGFjaW5nR0Zoc0Y5LyUlb3BlbkdRIltGJy8lJmNsb3NlR1EiXUYnLyUrYWN0aW9udHlwZUdRLWJyb3dzZXJ0YWJsZUYnLyUpcnRhYmxlaWRHUTUxODQ0Njc0NDA3NDMzMzQyODMwMkYnRmBvLyUpcmVhZG9ubHlHRj1GOQ==

We start by fitting a model that includes predicting variables for each of the columns of data:Fit( C + C2*v + C3*w + C4*x + C5*y + C6*z, ExMat[..,2..6], ExMat[..,1], [v,w,x,y,z], summarize=embed ): Model:LC4kIikpKSpvbiMhIikhIiIqJiQiKV07eT0hIiYiIiIlInZHRitGJiomJCIpI1siZVMhIihGKyUid0dGK0YmKiYkIilNOEc+ISInRislInhHRitGKyomJCIpcjV6Y0Y1RislInlHRitGKyomJCIpInokUURGMEYrJSJ6R0YrRis=CoefficientsEstimateStandard Errort-valueP(>|t|)CJCEnIXBuIyEiJw==JCInbUxoISInJCEnR2tWISInJCInRzJvISInC2JCEnPHk9ISIkJCInPFs7ISIkJCEnYVI2ISImJCInM2hJISInC3JCEnOmVTISImJCIneSUqKikhIic=JCEnbjZYISImJCInO0pqISIpC4JCInOEc+ISIlJCInUl5NISIlJCInYSdlJiEiJw==JCInIlwrJyEiJw==C5JCInNnpjISIlJCInOy1OISIlJCInZ0A7ISImJCInOmU7ISInC6JCInUVFEISImJCInKT5sIiEiJg==JCInZE86ISImJCInLV09ISInR-squared:JCInSncpKiEiJw==Adjusted R-squared:JCInaV8oKiEiJw==

Residual Sum of SquaresResidual Mean SquareResidual Standard ErrorDegrees of FreedomJCInRms9ISIpJCInYkdQISIqJCInJTQkPiEiKA==IiIm Five Point Summary MinimumFirst QuartileMedianThird QuartileMaximumJCEnYkI+ISIoJCEneTk1ISIoJCInPF04ISIpJCInVlFwISIpJCInPDZJISIo

From the above, we can observe that both the r-squared and adjusted r-squared are reasonably high, however only one of the coefficient values has a significant p-value (C3).Let's try to fit the data again, this time keeping the two coefficients with the lowest p-values and the intercept.Fit( C + C3*v + C5*w, ExMat[..,[3,5]], ExMat[..,1], [v,w], summarize=embed ): Model:LCgkIilIbSpcJiEiKSIiIiomJCIpanZMTSEiKEYmJSJ2R0YmISIiKiYkIilGdk1uRipGJiUid0dGJkYmCoefficientsEstimateStandard Errort-valueP(>|t|)CJCInbSpcJiEiJw==JCInVHY9ISIoJCInXktIISIlJCIndCEpPiEjOQ==C3JCEnd0xNISImJCInV2UiKSEiJw==JCEnJSkzVSEiJg==JCInYmdIISIpC5JCIndk1uISImJCInJ2UjZiEiJg==JCInXU82ISImJCInVScpRyEiJw==R-squared:JCInRy8pKiEiJw==Adjusted R-squared:JCInTmIoKiEiJw==

Residual Sum of SquaresResidual Mean SquareResidual Standard ErrorDegrees of FreedomJCInIipcSCEiKQ==JCInUShvJCEiKg==JCInRT8+ISIoIiIp Five Point Summary MinimumFirst QuartileMedianThird QuartileMaximumJCEnJFE/IyEiKA==JCEnY2I4ISIoJCEnRiFvIyEiKQ==JCInSCZcKCEiKQ==JCInJSpRUSEiKA==

From the above, we can see that the r-squared value does go down, however the adjusted r-squared goes up! Let's fit the model one last time to see if removing C5 increases or decreases the adjusted r-squared.Fit( C + C3*v, ExMat[..,3], ExMat[..,1], [v], summarize=embed ): Model:LCYkIikjM0xiJiEiKSIiIiomJCIpdW08RCEiKEYmJSJ2R0YmISIiCoefficientsEstimateStandard Errort-valueP(>|t|)CJCInSmBiISInJCInQFc9ISIoJCInQDZJISIlJCInWC1DISM6C3JCEnbjxEISImJCInJSp6NyEiJw==JCEnLW4+ISIlJCInTV41ISM4R-squared:JCInb3MoKiEiJw==Adjusted R-squared:JCInVVooKiEiJw==

Residual Sum of SquaresResidual Mean SquareResidual Standard ErrorDegrees of FreedomJCInPUVNISIpJCInKG8hUSEiKg==JCInN14+ISIoIiIq Five Point Summary MinimumFirst QuartileMedianThird QuartileMaximumJCEna2BBISIoJCEnU2ciKiEiKQ==JCEnTVU/ISIpJCInJjNgIyEiKQ==JCInckhbISIo

We can see that the final adjusted r-squared value is lower than the previous two, so we are probably better to keep the additional C5 coefficient value.

There are more ways to fit curves or surfaces to data. The Optimization package is a collection of commands for numerically solving optimization problems, which involve finding the minimum or maximum of an objective function possibly subject to constraints. The package takes advantage of built-in library routines provided by the Numerical Algorithms Group (NAG).The package solves linear programs (LPs), quadratic programs (QPs), nonlinear programs (NLPs), and both linear and nonlinear least-squares problems. Both constrained and unconstrained problems are accepted.

restart; with(Optimization);

NypJKkltcG9ydE1QU0c2IkksSW50ZXJhY3RpdmVHRiRJKExQU29sdmVHRiRJKExTU29sdmVHRiRJKU1heGltaXplR0YkSSlNaW5pbWl6ZUdGJEkpTkxQU29sdmVHRiRJKFFQU29sdmVHRiQ=

We can use the least squares command to do curve fitting. For example, if we have the following data points, (LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYqLUklbXN1YkdGJDYlLUkjbWlHRiQ2KFEieEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUrZXhlY3V0YWJsZUdRJmZhbHNlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLyUqcmVmZXJlbmNlR1ElOTUwMkYnLyUyZGlzcGxheV9yZWZlcmVuY2VHUSU5NTAxRictRiM2KC1GLzYoUSJpRidGMkY1RjhGO0Y+RjUvJTBmb250X3N0eWxlX25hbWVHUSpfcHN0eWxlNDFGJy9GOVEnbm9ybWFsRicvRjxRJjE4ODk4RicvRj9RJjE4ODk3RicvJS9zdWJzY3JpcHRzaGlmdEdRIjBGJy1JI21vR0YkNjBRIixGJ0Y1RklGO0Y+LyUmZmVuY2VHRjcvJSpzZXBhcmF0b3JHRjQvJSlzdHJldGNoeUdGNy8lKnN5bW1ldHJpY0dGNy8lKGxhcmdlb3BHRjcvJS5tb3ZhYmxlbGltaXRzR0Y3LyUnYWNjZW50R0Y3LyUnbHNwYWNlR1EmMC4wZW1GJy8lJ3JzcGFjZUdRLDAuMzMzMzMzM2VtRictRiw2JS1GLzYoUSJ5RidGMkY1RjhGO0Y+RkFGT0Y1RkZGSUZLRk0=):data := [[1, 1.5], [2, 3.5], [2.5, 1.9], [3.1, 4.5], [4.3, 1.9], [4.7, 2.4], [5.8, 3], [6, 5.7], [6.6, 3.4], [6.7, .5], [7.1, 4], [7.4, 1.7]];

Ny43JCIiIiQiIzohIiI3JCIiIyQiI05GJzckJCIjREYnJCIjPkYnNyQkIiNKRickIiNYRic3JCQiI1ZGJ0YvNyQkIiNaRickIiNDRic3JCQiI2VGJyIiJDckIiInJCIjZEYnNyQkIiNtRickIiNNRic3JCQiI25GJyQiIiZGJzckJCIjckYnIiIlNyQkIiN1RickIiM8Ric=

We want to fit the data with a quintic polynomial as follows:LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzY6LUkjbWlHRiQ2I1EhRictRiM2PEYrLUYjNjotRiw2OFEicEYnLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRicvJSVzaXplR1EjMTJGJy8lJWJvbGRHUSZmYWxzZUYnLyUnaXRhbGljR1EldHJ1ZUYnLyUqdW5kZXJsaW5lR0Y+LyUqc3Vic2NyaXB0R0Y+LyUsc3VwZXJzY3JpcHRHRj4vJStmb3JlZ3JvdW5kR1EoWzAsMCwwXUYnLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV1GJy8lJ29wYXF1ZUdGPi8lK2V4ZWN1dGFibGVHRj4vJSlyZWFkb25seUdGPi8lKWNvbXBvc2VkR0Y+LyUqY29udmVydGVkR0Y+LyUraW1zZWxlY3RlZEdGPi8lLHBsYWNlaG9sZGVyR0Y+LyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR0Y+LyUwZm9udF9zdHlsZV9uYW1lR1EqX3BzdHlsZTU0RicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLyUqcmVmZXJlbmNlR1EmMTg5MDBGJy8lMmRpc3BsYXlfcmVmZXJlbmNlR1EmMTg4OTlGJy1JI21vR0YkNi1RMCZBcHBseUZ1bmN0aW9uO0YnL0Zcb1Enbm9ybWFsRicvJSZmZW5jZUdGPi8lKnNlcGFyYXRvckdGPi8lKXN0cmV0Y2h5R0Y+LyUqc3ltbWV0cmljR0Y+LyUobGFyZ2VvcEdGPi8lLm1vdmFibGVsaW1pdHNHRj4vJSdhY2NlbnRHRj4vJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZqcC1JKG1mZW5jZWRHRiQ2OC1GIzY4LUYsNjhRInhGJ0Y2RjlGPEY/RkJGREZGRkhGS0ZORlBGUkZURlZGWEZaRmZuRmhuRltvRl5vRmFvRjZGOUY8L0ZARj5GQkZERkZGSEZLRk5GUEZSRlRGVkZYRlpGZm5GaG5GaG9GXm9GYW9GNkY5RjxGZXFGQkZERkZGSEZLRk5GUEZSRlRGVkZYRlpGZm5GaG5GaG9GXm9GYW9GNkY5RjxGZXFGQkZERkZGSEZLRk5GUEZSRlRGVkZYRlpGZm5GaG5GaG9GXm9GYW8tRmVvNkFRIj1GJ0Y2RjlGPEZlcUZCRkRGRkZIRktGTkZQRlJGVEZWRlhGWkZmbkZobkZob0Zeb0Zhb0Zqb0ZccEZecEZgcEZicEZkcEZmcC9GaXBRLDAuMjc3Nzc3OGVtRicvRlxxRmpxLUYjNkNGKy1GIzY7LUYsNjhRImFGJ0Y2RjlGPEY/RkJGREZGRkhGS0ZORlBGUkZURlZGWEZaRmZuRmhuRltvRl5vRmFvLUZlbzYtUTEmSW52aXNpYmxlVGltZXM7RidGaG9Gam9GXHBGXnBGYHBGYnBGZHBGZnBGaHBGW3EtRiM2OC1JJW1zdXBHRiQ2JUZicS1JI21uR0YkNjhRIjVGJ0Y2RjlGPEZlcUZCRkRGRkZIRktGTkZQRlJGVEZWRlhGWkZmbkZobkZob0Zeb0Zhby8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRidGNkY5RjxGZXFGQkZERkZGSEZLRk5GUEZSRlRGVkZYRlpGZm5GaG5GaG9GXm9GYW9GK0Y2RjlGPEZlcUZCRkRGRkZIRktGTkZQRlJGVEZWRlhGWkZmbkZobkZob0Zeb0Zhby1GZW82QVEiK0YnRjZGOUY8RmVxRkJGREZGRkhGS0ZORlBGUkZURlZGWEZaRmZuRmhuRmhvRl5vRmFvRmpvRlxwRl5wRmBwRmJwRmRwRmZwL0ZpcFEsMC4yMjIyMjIyZW1GJy9GXHFGZnMtRiM2Oy1GLDY4USJiRidGNkY5RjxGP0ZCRkRGRkZIRktGTkZQRlJGVEZWRlhGWkZmbkZobkZbb0Zeb0Zhb0Zjci1GIzY4LUZpcjYlRmJxLUZcczY4USI0RidGNkY5RjxGZXFGQkZERkZGSEZLRk5GUEZSRlRGVkZYRlpGZm5GaG5GaG9GXm9GYW9GX3NGNkY5RjxGZXFGQkZERkZGSEZLRk5GUEZSRlRGVkZYRlpGZm5GaG5GaG9GXm9GYW9GK0Y2RjlGPEZlcUZCRkRGRkZIRktGTkZQRlJGVEZWRlhGWkZmbkZobkZob0Zeb0Zhb0Zicy1GIzY7LUYsNjhRImNGJ0Y2RjlGPEY/RkJGREZGRkhGS0ZORlBGUkZURlZGWEZaRmZuRmhuRltvRl5vRmFvRmNyLUYjNjgtRmlyNiVGYnEtRlxzNjhRIjNGJ0Y2RjlGPEZlcUZCRkRGRkZIRktGTkZQRlJGVEZWRlhGWkZmbkZobkZob0Zeb0Zhb0Zfc0Y2RjlGPEZlcUZCRkRGRkZIRktGTkZQRlJGVEZWRlhGWkZmbkZobkZob0Zeb0Zhb0YrRjZGOUY8RmVxRkJGREZGRkhGS0ZORlBGUkZURlZGWEZaRmZuRmhuRmhvRl5vRmFvRmJzLUYjNjstRiw2OFEiZEYnRjZGOUY8Rj9GQkZERkZGSEZLRk5GUEZSRlRGVkZYRlpGZm5GaG5GW29GXm9GYW9GY3ItRiM2OC1GaXI2JUZicS1GXHM2OFEiMkYnRjZGOUY8RmVxRkJGREZGRkhGS0ZORlBGUkZURlZGWEZaRmZuRmhuRmhvRl5vRmFvRl9zRjZGOUY8RmVxRkJGREZGRkhGS0ZORlBGUkZURlZGWEZaRmZuRmhuRmhvRl5vRmFvRitGNkY5RjxGZXFGQkZERkZGSEZLRk5GUEZSRlRGVkZYRlpGZm5GaG5GaG9GXm9GYW9GYnMtRiM2Oi1GLDY4USJlRidGNkY5RjxGP0ZCRkRGRkZIRktGTkZQRlJGVEZWRlhGWkZmbkZobkZbb0Zeb0Zhb0ZjckZicUY2RjlGPEZlcUZCRkRGRkZIRktGTkZQRlJGVEZWRlhGWkZmbkZobkZob0Zeb0Zhb0Zicy1GLDY4USJmRidGNkY5RjxGP0ZCRkRGRkZIRktGTkZQRlJGVEZWRlhGWkZmbkZobkZbb0Zeb0Zhb0Y2RjlGPEZlcUZCRkRGRkZIRktGTkZQRlJGVEZWRlhGWkZmbkZobkZob0Zeb0Zhb0YrRjZGOUY8RmVxRkJGREZGRkhGS0ZORlBGUkZURlZGWEZaRmZuRmhuRmhvRl5vRmFvRitGNkY5RjxGZXFGQkZERkZGSEZLRk5GUEZSRlRGVkZYRlpGZm5GaG5GaG9GXm9GYW8=.Now we need to find a, b, c, d, e, and, f such thatLUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzY6LUkmbWZyYWNHRiQ2Ki1JI21uR0YkNjhRIjFGJy8lJ2ZhbWlseUdRMFRpbWVzfk5ld35Sb21hbkYnLyUlc2l6ZUdRIzEyRicvJSVib2xkR1EmZmFsc2VGJy8lJ2l0YWxpY0dGOi8lKnVuZGVybGluZUdGOi8lKnN1YnNjcmlwdEdGOi8lLHN1cGVyc2NyaXB0R0Y6LyUrZm9yZWdyb3VuZEdRKFswLDAsMF1GJy8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRicvJSdvcGFxdWVHRjovJStleGVjdXRhYmxlR0Y6LyUpcmVhZG9ubHlHRjovJSljb21wb3NlZEdGOi8lKmNvbnZlcnRlZEdGOi8lK2ltc2VsZWN0ZWRHRjovJSxwbGFjZWhvbGRlckdGOi8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdGOi8lMGZvbnRfc3R5bGVfbmFtZUdRKl9jc3R5bGUyNkYnLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGJy8lKnJlZmVyZW5jZUdRJjE4OTAyRicvJTJkaXNwbGF5X3JlZmVyZW5jZUdRJjE4OTAxRictRiM2OC1GLzY4USIyRidGMkY1RjhGO0Y9Rj9GQUZDRkZGSUZLRk1GT0ZRRlNGVUZXRllGZm5GaW5GXG9GMkY1RjhGO0Y9Rj9GQUZDRkZGSUZLRk1GT0ZRRlNGVUZXRllGZm5GaW5GXG8vJS5saW5ldGhpY2tuZXNzR0YxLyUrZGVub21hbGlnbkdRJ2NlbnRlckYnLyUpbnVtYWxpZ25HRmhvLyUpYmV2ZWxsZWRHRjpGQ0ZGLUkrbXVuZGVyb3ZlckdGJDYnLUkjbW9HRiQ2QVEmJlN1bTtGJ0YyRjVGOEY7Rj1GP0ZBRkNGRkZJRktGTUZPRlFGU0ZVRldGWUZmbkZpbkZcby8lJmZlbmNlR1EmdW5zZXRGJy8lKnNlcGFyYXRvckdGZnAvJSlzdHJldGNoeUdRJXRydWVGJy8lKnN5bW1ldHJpY0dGZnAvJShsYXJnZW9wR0ZbcS8lLm1vdmFibGVsaW1pdHNHRltxLyUnYWNjZW50R0ZmcC8lJ2xzcGFjZUdRJjAuMGVtRicvJSdyc3BhY2VHUSwwLjE2NjY2NjdlbUYnLUYjNjotSSNtaUdGJDY4USJ4RidGMkY1RjgvRjxGW3FGPUY/RkFGQ0ZGRklGS0ZNRk9GUUZTRlVGV0ZZL0ZnblEnaXRhbGljRidGaW5GXG8tRmFwNkFRIj1GJ0YyRjVGOEY7Rj1GP0ZBRkNGRkZJRktGTUZPRlFGU0ZVRldGWUZmbkZpbkZcby9GZXBGOi9GaHBGOi9GanBGOi9GXXFGOi9GX3FGOi9GYXFGOi9GY3FGOi9GZXFRLDAuMjc3Nzc3OGVtRicvRmhxRl5zRi5GMkY1RjhGO0Y9Rj9GQUZDRkZGSUZLRk1GT0ZRRlNGVUZXRllGZm5GaW5GXG8tRiM2OC1GLzY4RjdGMkY1RjhGO0Y9Rj9GQUZDRkZGSUZLRk1GT0ZRRlNGVUZXRllGZm5GaW5GXG9GMkY1RjhGO0Y9Rj9GQUZDRkZGSUZLRk1GT0ZRRlNGVUZXRllGZm5GaW5GXG9GXHMvJSxhY2NlbnR1bmRlckdGOi1JJW1zdXBHRiQ2JS1JKG1mZW5jZWRHRiQ2OC1GIzY7LUZdcjY4USJwRidGMkY1RjhGYHJGPUY/RkFGQ0ZGRklGS0ZNRk9GUUZTRlVGV0ZZRmFyRmluRlxvLUZqczY4LUYjNjgtSSVtc3ViR0YkNiZGXHItRiM2OC1GXXI2OFEiaUYnRjJGNUY4RmByRj1GP0ZBRkNGRkZJRktGTUZPRlFGU0ZVRldGWUZhckZpbkZcb0YyRjVGOEY7Rj1GP0ZBRkNGRkZJRktGTUZPRlFGU0ZVRldGWUZmbkZpbkZcby8lL3N1YnNjcmlwdHNoaWZ0R1EiMEYnRlVGMkY1RjhGO0Y9Rj9GQUZDRkZGSUZLRk1GT0ZRRlNGVUZXRllGZm5GaW5GXG9GMkY1RjhGO0Y9Rj9GQUZDRkZGSUZLRk1GT0ZRRlNGVUZXRllGZm5GaW5GXG8tRmFwNkFRKCZtaW51cztGJ0YyRjVGOEY7Rj1GP0ZBRkNGRkZJRktGTUZPRlFGU0ZVRldGWUZmbkZpbkZcb0ZmckZnckZockZpckZqckZbc0Zccy9GZXFRLDAuMjIyMjIyMmVtRicvRmhxRmR1LUZmdDYmLUZdcjY4USJ5RidGMkY1RjhGYHJGPUY/RkFGQ0ZGRklGS0ZNRk9GUUZTRlVGV0ZZRmFyRmluRlxvRmh0Rl11RlVGMkY1RjhGO0Y9Rj9GQUZDRkZGSUZLRk1GT0ZRRlNGVUZXRllGZm5GaW5GXG9GMkY1RjhGO0Y9Rj9GQUZDRkZGSUZLRk1GT0ZRRlNGVUZXRllGZm5GaW5GXG9GX28vJTFzdXBlcnNjcmlwdHNoaWZ0R0ZfdUYyRjVGOEY7Rj1GP0ZBRkNGRkZJRktGTUZPRlFGU0ZVRldGWUZmbkZpbkZcbw==is a minimum.LSSolve takes as input a list of the residues LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzY6LUkjbWlHRiQ2I1EhRictRiM2O0YrLUYjNjotRiw2OFEicEYnLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRicvJSVzaXplR1EjMTJGJy8lJWJvbGRHUSZmYWxzZUYnLyUnaXRhbGljR1EldHJ1ZUYnLyUqdW5kZXJsaW5lR0Y+LyUqc3Vic2NyaXB0R0Y+LyUsc3VwZXJzY3JpcHRHRj4vJStmb3JlZ3JvdW5kR1EoWzAsMCwwXUYnLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV1GJy8lJ29wYXF1ZUdGPi8lK2V4ZWN1dGFibGVHRj4vJSlyZWFkb25seUdGPi8lKWNvbXBvc2VkR0Y+LyUqY29udmVydGVkR0Y+LyUraW1zZWxlY3RlZEdGPi8lLHBsYWNlaG9sZGVyR0Y+LyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR0Y+LyUwZm9udF9zdHlsZV9uYW1lR1EqX3BzdHlsZTQxRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLyUqcmVmZXJlbmNlR1EmMTg5MDRGJy8lMmRpc3BsYXlfcmVmZXJlbmNlR1EmMTg5MDNGJy1JI21vR0YkNi1RMCZBcHBseUZ1bmN0aW9uO0YnL0Zcb1Enbm9ybWFsRicvJSZmZW5jZUdGPi8lKnNlcGFyYXRvckdGPi8lKXN0cmV0Y2h5R0Y+LyUqc3ltbWV0cmljR0Y+LyUobGFyZ2VvcEdGPi8lLm1vdmFibGVsaW1pdHNHRj4vJSdhY2NlbnRHRj4vJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZqcC1JKG1mZW5jZWRHRiQ2OC1GIzY4LUklbXN1YkdGJDYmLUYsNjhRInhGJ0Y2RjlGPEY/RkJGREZGRkhGS0ZORlBGUkZURlZGWEZaRmZuRmhuRltvRl5vRmFvLUYjNjgtRiw2OFEiaUYnRjZGOUY8Rj9GQkZERkZGSEZLRk5GUEZSRlRGVkZYRlpGZm5GaG5GW29GXm9GYW9GNkY5RjwvRkBGPkZCRkRGRkZIRktGTkZQRlJGVEZWRlhGWkZmbkZobkZob0Zeb0Zhby8lL3N1YnNjcmlwdHNoaWZ0R1EiMEYnRlpGNkY5RjxGXXJGQkZERkZGSEZLRk5GUEZSRlRGVkZYRlpGZm5GaG5GaG9GXm9GYW9GNkY5RjxGXXJGQkZERkZGSEZLRk5GUEZSRlRGVkZYRlpGZm5GaG5GaG9GXm9GYW9GNkY5RjxGXXJGQkZERkZGSEZLRk5GUEZSRlRGVkZYRlpGZm5GaG5GaG9GXm9GYW8tRmVvNkFRKCZtaW51cztGJ0Y2RjlGPEZdckZCRkRGRkZIRktGTkZQRlJGVEZWRlhGWkZmbkZobkZob0Zeb0Zhb0Zqb0ZccEZecEZgcEZicEZkcEZmcC9GaXBRLDAuMjIyMjIyMmVtRicvRlxxRmVyLUZjcTYmLUYsNjhRInlGJ0Y2RjlGPEY/RkJGREZGRkhGS0ZORlBGUkZURlZGWEZaRmZuRmhuRltvRl5vRmFvRmhxRl5yRlpGNkY5RjxGXXJGQkZERkZGSEZLRk5GUEZSRlRGVkZYRlpGZm5GaG5GaG9GXm9GYW9GK0Y2RjlGPEZdckZCRkRGRkZIRktGTkZQRlJGVEZWRlhGWkZmbkZobkZob0Zeb0Zhbw==. p := a*x^5+b*x^4+c*x^3+d*x^2+e*x+f;

LC4qJkkiYUc2IiIiIilJInhHRiUiIiZGJkYmKiZJImJHRiVGJilGKCIiJUYmRiYqJkkiY0dGJUYmKUYoIiIkRiZGJiomSSJkR0YlRiYpRigiIiNGJkYmKiZJImVHRiVGJkYoRiZGJkkiZkdGJUYm

residues := map( d-> eval(p, x = d[1])-d[2], data);

Ny4sMEkiYUc2IiIiIkkiYkdGJUYmSSJjR0YlRiZJImRHRiVGJkkiZUdGJUYmSSJmR0YlRiYkIiM6ISIiRi4sMComIiNLRiZGJEYmRiYqJiIjO0YmRidGJkYmKiYiIilGJkYoRiZGJiomIiIlRiZGKUYmRiYqJiIiI0YmRipGJkYmRitGJiQiI05GLkYuLDAqJiQiKERjdyohIiZGJkYkRiZGJiomJCInRDFSISIlRiZGJ0YmRiYqJiQiJkRjIiEiJEYmRihGJkYmKiYkIiREJyEiI0YmRilGJkYmKiYkIiNERi5GJkYqRiZGJkYrRiYkIiM+Ri5GLiwwKiYkIileIkgnR0ZARiZGJEYmRiYqJiQiJ0BOIypGREYmRidGJkYmKiYkIiYiekhGSEYmRihGJkYmKiYkIiRoKkZMRiZGKUYmRiYqJiQiI0pGLkYmRipGJkYmRitGJiQiI1hGLkYuLDAqJiQiKlYlM3E5RkBGJkYkRiZGJiomJCIoLCk9TUZERiZGJ0YmRiYqJiQiJjImekZIRiZGKEYmRiYqJiQiJVw9RkxGJkYpRiZGJiomJCIjVkYuRiZGKkYmRiZGK0YmRlBGLiwwKiYkIioyXU1II0ZARiZGJEYmRiYqJiQiKCJveltGREYmRidGJkYmKiYkIidCUTVGSEYmRihGJkYmKiYkIiU0QUZMRiZGKUYmRiYqJiQiI1pGLkYmRipGJkYmRitGJiQiI0NGLkYuLDAqJiQiKm9uTmMnRkBGJkYkRiZGJiomJCIpJ1w7OCJGREYmRidGJkYmKiYkIic3Xj5GSEYmRihGJkYmKiYkIiVrTEZMRiZGKUYmRiYqJiQiI2VGLkYmRipGJkYmRitGJiIiJEYuLDAqJiIld3hGJkYkRiZGJiomIiUnSCJGJkYnRiZGJiomIiQ7I0YmRihGJkYmKiYiI09GJkYpRiZGJiomIiInRiZGKkYmRiZGK0YmJCIjZEYuRi4sMComJCIrd0RMXzdGQEYmRiRGJkYmKiYkIilPWigqPUZERiZGJ0YmRiYqJiQiJydcKEdGSEYmRihGJkYmKiYkIiVjVkZMRiZGKUYmRiYqJiQiI21GLkYmRipGJkYmRitGJiQiI01GLkYuLDAqJiQiKzJeN104RkBGJkYkRiZGJiomJCIpQDY6P0ZERiZGJ0YmRiYqJiQiJ2oySUZIRiZGKEYmRiYqJiQiJSpbJUZMRiZGKUYmRiYqJiQiI25GLkYmRipGJkYmRitGJiQiIiZGLkYuLDAqJiQiK14kSFUhPUZARiZGJEYmRiYqJiQiKSJvNmEjRkRGJkYnRiZGJiomJCInNnpORkhGJkYoRiZGJiomJCIlVF1GTEYmRilGJkYmKiYkIiNyRi5GJkYqRiZGJkYrRiZGN0YuLDAqJiQiK0NtKz5BRkBGJkYkRiZGJiomJCIpd2wpKkhGREYmRidGJkYmKiYkIidDX1NGSEYmRihGJkYmKiYkIiV3YUZMRiZGKUYmRiYqJiQiI3VGLkYmRipGJkYmRitGJiQiIzxGLkYu

Now we can call LSSolve with the residues.sol := LSSolve(residues);

NyQkIjMnNCd6SDNsVCp6KiEjPDcoL0kiYUc2Il0zRjc0QzI1RUIyNTE5RjhBL0kiYkdGKV1CRkMzMDkyREVGMEU4QTcxL0kiY0dGKV0zRkY4OTg3RDE5MUUyMEI4L0kiZEdGKV1DMDFDMkRGMjlDQ0MzNjVCL0kiZUdGKV00MDJDN0MzNEM5NDNEMzIyL0kiZkdGKV1DMDFDN0M0RThGRjg3NDM0

Our polynomial becomes:poly := eval(p, sol[2]);

LC4qJl0zRjc0QzI1RUIyNTE5RjhBIiIiKUkieEc2IiIiJkYlRiUqJl0zRkMzMDkyREVGMEU4QTcxRiUpRiciIiVGJSEiIiomXTNGRjg5ODdEMTkxRTIwQjhGJSlGJyIiJEYlRiUqJl00MDFDMkRGMjlDQ0MzNjVCRiUpRiciIiNGJUYuKiZdNDAyQzdDMzRDOTQzRDMyMkYlRidGJUYlXTQwMUM3QzRFOEZGODc0MzRGLg==

We can now plot the data and the curvep1 := Statistics:-ScatterPlot(data, symbolsize=15, color="Black", legend="Data"): p2 := plot(poly, x = .9 .. 8, thickness=4, color="Red", linestyle=dash, legend="Fitted Model"): plots:-display(p1, p2, axes=boxed);

NiwtJSdDVVJWRVNHNic3LjckJCIjNSEiIiQiIzohIiI3JCQiIz8hIiIkIiNOISIiNyQkIiNEISIiJCIjPiEiIjckJCIjSiEiIiQiI1ghIiI3JCQiI1YhIiIkIiM+ISIiNyQkIiNaISIiJCIjQyEiIjckJCIjZSEiIiQiI0khIiI3JCQiI2chIiIkIiNkISIiNyQkIiNtISIiJCIjTSEiIjckJCIjbiEiIiQiIiYhIiI3JCQiI3IhIiIkIiNTISIiNyQkIiN1ISIiJCIjPCEiIi0lJ0xFR0VOREc2Iy0lKV9UWVBFU0VURzYjLUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1LUkjbXNHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2I1ElRGF0YTYiLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuNiIvJSVzaXplR1EjMTI2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EmZmFsc2U2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdub3JtYWw2Ii0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUnU1lNQk9MRzYkJShERUZBVUxURyIjOi0lJlNUWUxFRzYjJSZQT0lOVEctJSdDVVJWRVNHNic3ZHc3JCQiIiohIiIkIjJ0XSZvJT1ebSwiISM7NyQkIjAjUnljKipHdCQqISM6JCIyJzQ4N2g4KlE+IiEjOzckJCIxYjVAX3YzKXAqISM7JCIyTVMwREtTLU0iISM7NyQkIjJYKil5PGtOaisiISM7JCIxNkNlXT9OJ1wiISM6NyQkIjJrRF5BbC9KLyIhIzskIjI+JSpcJjNSc1c7ISM7NyQkIjJlPE5nISopcHo1ISM7JCIydlEsdSYqUVJ5IiEjOzckJCIxPE1vamtpODYhIzokIjIiUmd0alRzMD4hIzs3JCQiMm1MbldbYyhbNiEjOyQiMWF2eT47b0M/ISM6NyQkIjIuMztBKykzJj0iISM7JCIxKTNGOjo7LjkjISM6NyQkIjFMbUsvSUlANyEjOiQiMlFSeldCZSRbQSEjOzckJCIyWiRwUSZlYSZlNyEjOyQiMl9BIypSOVJBTiMhIzs3JCQiMmxJaCMpemw4SCIhIzskIjIwcm9oKipleVYjISM7NyQkIjJqQFZFWC4kRzghIzskIjFvJj0qSGYqeV8jISM6NyQkIjJZIkhld0ZSbDghIzskIjIxWUNyTSN5NkUhIzs3JCQiMnJUJG83XzgsOSEjOyQiMTFUVS1ybidvIyEjOjckJCIyaUNcKW9GZkw5ISM7JCIwM0MqUkIjKVxGISM5NyQkIjJBWSNcVXo9czkhIzskIjIwIWUnZlw3Ij5HISM7NyQkIjJ4XzByNiQpW10iISM7JCIybSdcQmlHMnRHISM7NyQkIjIkcFF4eGAiSGEiISM7JCIyTUQhPiQqKilmSUghIzs3JCQiMic+UnlzLmV3OiEjOyQiMTdjKVxwLnEoSCEjOjckJCIxQlkjej87TmgiISM6JCIyVjI2I2ZmQ0JJISM7NyQkIjJOclV2JnlvWzshIzskIjFRJ1s9MEhIMSQhIzo3JCQiMlUmMzwleSZRJm8iISM7JCIyJEdtMzA4LCtKISM7NyQkIjI9TW9tJGYzPjwhIzskIjJgYmtCJ3lNSUohIzs3JCQiMm5NcFE1T2F2IiEjOyQiMVIjeXBcaiNmSiEjOjckJCIxOUdjNlA+JHoiISM6JCIyQylIR3FRSiY9JCEjOzckJCIyakZiKyM9MUU9ISM7JCIxaCRSLkpXWz8kISM6NyQkIjEjUnl3Rmc6Jz0hIzokIjFOY0kxJzRHQSQhIzo3JCQiMiNbJ0h6bEwjKSo9ISM7JCIxKT4oekojRyJRSyEjOjckJCIyUXplZFM2VCQ+ISM7JCIxS2kqUjMmNF1LISM6NyQkIjJaIkhlclwjKW8+ISM7JCIwKlFiLFgpKmVLISM5NyQkIjIlWydIellvdCsjISM7JCIyd2okRyhbSmZFJCEjOzckJCIxVSUpbzg6K1U/ISM6JCIxOyIqNGQ7cXBLISM6NyQkIjFlOkojb3oqeT8hIzokIjFpOG44PklySyEjOjckJCIxLjE3KlIoWzdAISM6JCIxMysoR3I+MkYkISM6NyQkIjJOcVNUPT8iXEAhIzskIjE9InAmUU4sb0shIzo3JCQiMlgjXCkqeigpZSQ9IyEjOyQiMndCUDNGQE9FJCEjOzckJCIxJGU7VlI9Jz5BISM6JCIyJyo+XHNEXHNEJCEjOzckJCIyJylwUnojUiVbRCMhIzskIjFEIls8QiZRXEshIzo3JCQiMkViNXJ6QDxIIyEjOyQiMiUpKmZoU0xiUkshIzs3JCQiMlUhMztvJlJzSyMhIzskIjE6QHNiOW1HSyEjOjckJCIyYC0wIT47Y2pCISM7JCIxJzNaY0wtaUAkISM6NyQkIjFbJ0hcIUhlKlIjISM6JCIwIip5cyJlay1LISM5NyQkIjI8TW9jcyNvS0MhIzskIjFyJlxyQEwjKj0kISM6NyQkIjFSeGFiKT0xWiMhIzokIjEkNFpGW1xHPCQhIzo3JCQiMidlPE4jUVxYXSMhIzskIjI5Jlw9Tyd5dDokISM7NyQkIjJzVyopM1NFMmEjISM7JCIyeCUzOy5hOFNKISM7NyQkIjFvTnJCQ052RCEjOiQiMmxqZTJeK0k3JCEjOzckJCIyRV0rQDwrUWgjISM7JCIyWS5pLiF6Ti5KISM7NyQkIjJNcVNyNCQ0WkUhIzskIjJtKj4kKWUpNGYzJCEjOzckJCIxWSNcWyczKVtvIyEjOiQiMSY0QUBmMmQxJCEjOjckJCIxc1YoR0o/JD5GISM6JCIydnA4Uy8zcS8kISM7NyQkIjEiPlFtbCIqcHYjISM6JCIyeHE6SFBIai0kISM7NyQkIjFNb09iY1wqeSMhIzokIjJuejhsb24kM0khIzs3JCQiMSYqKil6OlJeRUchIzokIjFLIioqKT5KJnkpSCEjOjckJCIybUhmeXNyQSdHISM7JCIyJG9DN1Y1MG9IISM7NyQkIjAuMS09MSEpKkchIzkkIjJsQCR5OGtNW0ghIzs3JCQiMVwoXD48NE8kSCEjOiQiMShvZiZmWicpR0ghIzo3JCQiMV8vNCxEInknSCEjOiQiMkZpVjhOXC4iSCEjOzckJCIyQlwpcDBzeS9JISM7JCIyTjlYN0U8MSpHISM7NyQkIjApZj4heUUrLyQhIzkkIjEvdSxKSDlzRyEjOjckJCIyVShbKEh4OXIyJCEjOyQiMmx0enlnM0omRyEjOzckJCIya0xudUwpbzVKISM7JCIxWT0obyZ6R09HISM6NyQkIjEjUXc3VCR6WkohIzokIjJ2dFU3ai0jPUchIzs3JCQiMClmPjFPTCQ9JCEjOSQiMUQpRzYqUlUsRyEjOjckJCIxS2tHJkgoeT1LISM6JCIyUHF5W2ZnX3kjISM7NyQkIjF1WiY0Iik9ZUQkISM6JCIxKTRsTmVNIXBGISM6NyQkIjFXKW8oKT1FKipHJCEjOiQiMlAiUmxyLnNhRiEjOzckJCIxRmEzNG8mW0skISM6JCIxZjYnPVhAMnUjISM6NyQkIjE2QVdrVlVqTCEjOiQiMm1mIjQnXHpncyMhIzs3JCQiMidvUHZFVE0pUiQhIzskIjBgPltvImY4RiEjOTckJCIxUnhhQ2gwTU0hIzokIjJjb3JTejA7cSMhIzs3JCQiMmQ4RmtJdS5aJCEjOyQiMlpKaCUzeUQhcCMhIzs3JCQiMVokcEd5ZFBdJCEjOiQiMlYpeVEtbmYhbyMhIzs3JCQiMkdmPXhbViRSTiEjOyQiMVc/LzQzOHJFISM6NyQkIjJuTXA9SmVZZCQhIzskIjJQWmpYeS9FbSMhIzs3JCQiMk1tS2w+TUNoJCEjOyQiMXVcYXotWWFFISM6NyQkIjJ0XSwkeShwZGskISM7JCIydlYpKXBYQyJbRSEjOzckJCIydWAyJj5AQiVvJCEjOyQiMlE+cEooKT49ayMhIzs3JCQiMUdjNzw5JSo9UCEjOiQiMlVeOWpibHFqIyEjOzckJCIybUhmR1glR2BQISM7JCIxMHYuckBDTEUhIzo3JCQiMmxQdiFcVUAheiQhIzskIjJiVmBzXjksaiMhIzs3JCQiMiRbJ0hSViRIRlEhIzskIjJrOSV5LlQrR0UhIzs3JCQiMXRZJD4rSTUnUSEjOiQiMm5JZCdSQilwaSMhIzs3JCQiMlgiSGVRbiZvKlEhIzskIjIocCZwOUtNb2kjISM7NyQkIitLLCM9JFIhIiokIjFlIlIsWj13aSMhIzo3JCQiMndnQCQ+M3NwUiEjOyQiMXMsO1ksXUhFISM6NyQkIjElKW9QXmFeLVMhIzokIjIobzxyeiEqKT5qIyEjOzckJCIxUnhheTc9U1MhIzokIjF3JHpNekNlaiMhIzo3JCQiMUppQ08iKmV2UyEjOiQiMTBUYUNNUFNFISM6NyQkIjJsRWBFKlxsNVQhIzskIjIneWVLQ3V3WEUhIzs3JCQiMWM1QFcoW2Q5JSEjOiQiMSlIOSJcIlI/bCMhIzo3JCQiMCc+UlxAYSI9JSEjOSQiMiU9bWRoN0tmRSEjOzckJCIxYzdEZkRbPlUhIzokIjI5WF9Dcygqem0jISM7NyQkIjEnPlIpKj5xVUQlISM6JCIxTHh6RTIhb24jISM6NyQkIjFFXy9vK2gpRyUhIzokIjEwIT56LW5pbyMhIzo3JCQiMlgoXCo+Iyk9YEslISM7JCIyOGROcXc8c3AjISM7NyQkIjApZj5rdyU+TyUhIzkkIjEoXHEhemsoKjNGISM6NyQkIjEjR2NfQUNbUiUhIzokIjJZWCdmJkg7LXMjISM7NyQkIjE8TG0jKSlvUFYlISM6JCIyZHlUPSozTU1GISM7NyQkIi8sLW1mTm1XISM4JCIxKlwkZT4uIm91IyEjOjckJCIxYzhGPy1aL1ghIzokIjInW1kxYGE2aUYhIzs3JCQiMSI0PU9lWTNhJSEjOiQiMG81KT5eVHhGISM5NyQkIjJPITM7VmpLdFghIzskIjJFMzhqYT47eiMhIzs3JCQiMSRmPShISyYpNFkhIzokIjJYalgnZVc8M0chIzs3JCQiMWI0PlNBaVlZISM6JCIyY25NPChSVURHISM7NyQkIjF1WyhSXDtLbyUhIzokIjI6NEE8bk5KJUchIzs3JCQiMTtKaV5TOTxaISM6JCIyLyM+KDQxMisnRyEjOzckJCIxTnFTc1NGX1ohIzokIjJ2XC0lPnUqeShHISM7NyQkIjF6ZDohZjAnKXklISM6JCIxcjJoUDUibypHISM6NyQkIjFLakUjZj9bI1shIzokIjJ1Pjxga1FnIkghIzs3JCQiMUxtS3RAMmlbISM6JCIwKmVYTDk7T0ghIzk3JCQiMTA1PydRJClbKlshIzokIjIwZjNZUlBUJkghIzs3JCQiMTpIZVM1Iz0kXCEjOiQiMWZYNVREaHVIISM6NyQkIjE4RV9rLiIqb1whIzokIjFMLDNLJHlgKkghIzo3JCQiMTtKaStHbC9dISM6JCIyaitrQ2RRYiwkISM7NyQkIjFXKil5Yy42UF0hIzokIjEjPidHJSlIJFIuJCEjOjckJCIxaEBWSWJxdl0hIzokIjFkLEY/SidlMCQhIzo3JCQiMUVfLzAyUzNeISM6JCIyJGUscEF0V3VJISM7NyQkIjFwTnJsSFZZXiEjOiQiMkJGTmVTP2c0JCEjOzckJCIxPU9zZ3o0IT0mISM6JCIyMUAvISpmSV02JCEjOzckJCIxQVYnZXpMcUAmISM6JCIxKD4iMzRNdU5KISM6NyQkIjE3Q1tYYT9fXyEjOiQiMTZdbylwd186JCEjOjckJCIyRGI1QFAuKilHJiEjOyQiMidlIkdYIWZTdkohIzs3JCQiMjAvM1lfLkVLJiEjOyQiMiZlVmtQbmgkPiQhIzs3JCQiMVkhND1wYCplYCEjOiQiMlssIUd1SiJIQCQhIzs3JCQiMTdEXSpINm5SJiEjOiQiMiUpPnlPXDxEQiQhIzs3JCQiMXZcKnpTeiZIYSEjOiQiMjhkPyU+KzxcSyEjOzckJCIwND1jJ3kybGEhIzkkIjEoUU4qKVFzbUUkISM6NyQkIjFaJHBlQ148XSYhIzokIjFAJypbbz88JUckISM6NyQkIjEkXClwJCoqR3dgJiEjOiQiMmRIJT5zS20rTCEjOzckJCIxOEVfZkRNc2IhIzokIjI4eUxQRHVmSiQhIzs3JCQiMVokcGUwJykzaCYhIzokIjEkejciKlxvQEwkISM6NyQkIjFUIkc7NT5iayYhIzokIjFqbVJNLyVmTSQhIzo3JCQiMWQ3RHFzXCNvJiEjOiQiMmJraFxzcShmTCEjOzckJCIxLS4xKClcKztkISM6JCIxXi9FUEBackwhIzo3JCQiMSwvM3N4al9kISM6JCIxKkhFOCwvTFEkISM6NyQkIjFDWSN6TzFyeSYhIzokIjJsdnA7Qm5NUiQhIzs3JCQiMSNHY0EpZjhCZSEjOiQiMV9NMjUjSElTJCEjOjckJCIybFJ6ZV5oJGVlISM7JCIxNl5IKmZ0N1QkISM6NyQkIjFfLTAmUVJfKmUhIzokIjFKKVEhel1uPU0hIzo3JCQiMS4wNWNydklmISM6JCIyLlY7I2VvYkNNISM7NyQkIjFDWiVwP3pxJ2YhIzokIjJfISp5R2ZaI0hNISM7NyQkIjFaJHBHXCtKKychIzokIjImRzZUYSY9RFYkISM7NyQkIjAvM09KK2kuJyEjOSQiMSZlK2FzaFVWJCEjOjckJCIxUHVbVms4dWchIzokIjFtQEx5RnJNTSEjOjckJCIxZTlIcXAxM2ghIzokIjJFcU8kb09tTE0hIzs3JCQiMVgiSCkpKVJDV2ghIzokIjI7YWE0IXooNFYkISM7NyQkIjFuS2w2KygpeWghIzokIjEkUVpzd1JvVSQhIzo3JCQiMSwtL2d4SjxpISM6JCIya2FUPydmUT9NISM7NyQkIjEtLzMmbzUxRCchIzokIjEwQUF2Yjs4TSEjOjckJCIxWCopeV8lKVIpRychIzokIjFoeiZSJGYyLk0hIzo3JCQiMDI5MyF6JEdLJyEjOSQiMWZNYGhMMiNSJCEjOjckJCIxKil5ZFcjNDBPJyEjOiQiMigqPitDajAheUwhIzs3JCQiMUxsSVZLLCRSJyEjOiQiMiR5KTMmKXo4VE8kISM7NyQkIjElcFFQXUorVichIzokIjFLKSpmPStFWUwhIzo3JCQiMSYqKil6OiQqeWxrISM6JCIyYVkmUScqZSFwSyQhIzs3JCQiMUhkOW9QXyxsISM6JCIxQjFccjtXMEwhIzo3JCQiMVolKikpZm43UGwhIzokIjJMJVIpXGclKj1HJCEjOzckJCIxXywuKjNJOGQnISM6JCIyVkRDPCJSPmRLISM7NyQkIjEjPlFPei8kM20hIzokIjFecVYkKXo5R0shIzo3JCQiMXpjOG9WYVZtISM6JCIxOmgiPltcIik+JCEjOjckJCIxdFgiNE9LMW8nISM6JCIxa2x2MVMza0ohIzo3JCQiMU5xU0RmPzluISM6JCIxTEhSRnEqNDgkISM6NyQkIjEiMzsjKio0Sl5uISM6JCIyQl5VTiVSKj00JCEjOzckJCIxemM4JT5eb3knISM6JCIyTCpRdF1vKj0wJCEjOzckJCIxSmhBJClbSUFvISM6JCIxcCNcQUpzJTRJISM6NyQkIjFzVyopKVJPJGZvISM6JCIyJz1NWFAnPkMnSCEjOzckJCIxViYzbnhWTSpvISM6JCIyPFx5NDJrbCJIISM7NyQkIjJiN0RxUnUkR3AhIzskIjJOLTdNaGVxJ0chIzs3JCQiMCI+UV8+JXAncCEjOSQiMUwvd0RGTzRHISM6NyQkIjFvTXA5PCc9KyghIzokIjJrUWYzJ1FMYUYhIzs3JCQiMVB1WzdQZFBxISM6JCIybyFwbkZURiZwIyEjOzckJCIyWCRvTyUqPSpRMighIzskIjIkRzBfYG5HS0UhIzs3JCQiMVghNDNQdnM1KCEjOiQiMTEqKm9YT3ZyRCEjOjckJCIxI0hlYzJoRzkoISM6JCIxUHJ0WEtUL0QhIzo3JCQiMmIvNCkqKmU8eXIhIzskIjJjQEoqNDZvTUMhIzs3JCQiMWhCWiV5XmZAKCEjOiQiMmFuciUqZVxvTiMhIzs3JCQiMTE3Q210R1xzISM6JCIyPUlvTzFiYEcjISM7NyQkIjFQc1cyKFx4RyghIzokIjFVLzBLYWIqPiMhIzo3JCQiMUZgMTAhZkNLKCEjOiQiMkR5J1FKVzA+QCEjOzckJCIxJioqKXpTPyFvTighIzokIjJsKW8kUWYxbC4jISM7NyQkIjF3XSxQPXQkUighIzokIjJtPEFbLypcVz4hIzs3JCQiMVokcD0tNjNWKCEjOiQiMnIiW048OXJbPSEjOzckJCIxc1YoKSplWlhZKCEjOiQiMndfI2YqUnAmZTwhIzs3JCQiMThFX0VWUCt2ISM6JCIxOidbO0UtKGY7ISM6NyQkIjEpcFIqPnhMTnYhIzokIjIoZWNIdCJ6K2MiISM7NyQkIjExOEUyJVFLZCghIzokIjE4NyVIIlFjWzkhIzo3JCQiMSRlOyRSSS4xdyEjOiQiMjEtUlRxKTRcOCEjOzckJCIxUHVbbSkpcFZ3ISM6JCIyXyVmYl0hUjlCIiEjOzckJCIwJGY9Q241enchIzkkIjI8akg/LCVbPDYhIzs3JCQiMWxIZiFlc1RyKCEjOiQiMWQuSl0hPjkrIiEjOjckJCIyT3ZdQExtI1x4ISM7JCIwbnRWIXpfPykpISM6NyQkIjFlO0xQKGZdeSghIzokIjEtI1IxeDYqcHYhIzs3JCQiMWI0PlosK0J5ISM6JCIxb2c9P2drMmkhIzs3JCQiMSYqKXl4eSh5ZHkhIzokIjJCa3RYaCpSRFwhIzw3JCQiMURcKWZsRkAqeSEjOiQiMSIqbyRbQyZcR08hIzs3JCQiMXVZJCo0ayQpR3ohIzokIjFXN0YtaSh5PyMhIzs3JCQiMXljOF9fWWx6ISM6JCIxYT8jKVtoIzNiKCEjPDckJCIjISkhIiIkITFpKSlIYDYiKm9rISM8LSUnTEVHRU5ERzYjLSUpX1RZUEVTRVRHNiMtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjUtSSNtc0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjYjUS1GaXR0ZWR+TW9kZWw2Ii8lJ2ZhbWlseUdRMFRpbWVzfk5ld35Sb21hbjYiLyUlc2l6ZUdRIzEyNiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJmZhbHNlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1Enbm9ybWFsNiItJSZDT0xPUkc2JiUkUkdCRyQiIzUhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiJC0lKlRISUNLTkVTU0c2IyIiJS0lJVZJRVdHNiQ7JCIiKiEiIiQiIyEpISIiJShERUZBVUxURy0mJSZfQVhJU0c2IyIiIjYnLSUrX0dSSURMSU5FU0c2Iy0lKV9WSVNJQkxFRzYjIiIhLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSYlJl9BWElTRzYjIiIjNictJStfR1JJRExJTkVTRzYjLSUpX1ZJU0lCTEVHNiMiIiEtJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIS0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJStBWEVTTEFCRUxTRzYlUSE2IlEhNiItJSVGT05URzYlJShERUZBVUxURyUoREVGQVVMVEciIzctJSpBWEVTU1RZTEVHNiMlJEJPWEctJSlfVklTSUJMRUc2IyIiIi0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJD8kISIiLSUpQk9VTkRTX1lHNiMkIiRJIiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJStPISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJVNKISIiLSUpQ0hJTERSRU5HNiItJStBTk5PVEFUSU9ORzYnLSUpQk9VTkRTX1hHNiMkIiIhISIiLSUpQk9VTkRTX1lHNiMkIiIhISIiLSUtQk9VTkRTX1dJRFRIRzYjJCIlK1MhIiItJS5CT1VORFNfSEVJR0hURzYjJCIlK1MhIiItJSlDSElMRFJFTkc2Ig==Ig==

Additional Example: Bivariate Polynomial Regression

Legal Notice: \302\251 Maplesoft, a division of Waterloo Maple Inc. 2017. Maplesoft and Maple are trademarks of Waterloo Maple Inc. This application may contain errors and Maplesoft is not liable for any damages resulting from the use of this material. This application is intended for non-commercial, non-profit use only. Contact Maplesoft for permission if you wish to use this application in for-profit activities.