restart; clamped_spline := proc(xx::list,y::list,fp0::numeric,fpn::numeric,S::name) local N, I, X, A, AA, FP0, FPN, M, H, XA, XL, XU, XZ, C, J, B, D, SS, L, P; if nops(xx)<>nops(y) then ERROR("Both lists must be the same size"); fi; N:=nops(xx)-1; for I from 0 to N do X[I]:= xx[I+1]; A[I]:=y[I+1]; od; FP0 := fp0; FPN := fpn; STEP 1 M := N - 1; for I from 0 to M do H[I] := X[I+1] - X[I]; od; STEP 2 XA[0] := 3.0 * (A[1] - A[0]) / H[0] - 3.0 * FP0; XA[N] := 3.0 * FPN - 3.0 * (A[N] - A[N-1]) / H[N-1]; STEP 3 for I from 1 to M do XA[I] := 3.0*(A[I+1]*H[I-1]-A[I]*(X[I+1]-X[I-1])+A[I-1]*H[I])/(H[I]*H[I-1]); od; STEP 4 XL[0] := 2.0 * H[0]; XU[0] := 0.5; XZ[0] := XA[0] / XL[0]; STEP 5 for I from 1 to M do XL[I] := 2.0 * (X[I+1] - X[I-1]) - H[I-1] * XU[I-1]; XU[I] := H[I] / XL[I]; XZ[I] := (XA[I] - H[I-1] * XZ[I-1]) / XL[I]; od; STEP 6 XL[N] := H[N-1] * (2.0 - XU[N-1]); XZ[N] := (XA[N] - H[N-1] * XZ[N-1]) / XL[N]; C[N] := XZ[N]; STEP 7 for I from 1 to N do J := N - I; C[J] := XZ[J] - XU[J] * C[J+1]; B[J] := (A[J+1] - A[J]) / H[J] - H[J] * (C[J+1] + 2.0 * C[J]) / 3.0; D[J] := (C[J+1] - C[J]) / (3.0 * H[J]); od; STEP 8 printf(`\n\nThe coefficients of the clamped cubic spline S[i]\non the subintervals are:\n`); printf(` a(i) b(i) c(i) d(i)\n`); for I from 0 to M do printf(` %13.8f %13.8f %13.8f %13.8f \n`, A[I], B[I], C[I], D[I]); od; printf(`\n\n`); for I from 0 to M do SS[I]:=A[I]+B[I]*(x-X[I])+C[I]*(x-X[I])^2+D[I]*(x-X[I])^3; od; L:=[seq([x<X[I+1],SS[I]],I=0..M)]; P:=[seq(op(L[I]),I=1..M+1)]; P := subsop(2*M+1=NULL,P); S:=simplify(piecewise(op(P))); end;

Warning, imaginary unit `I` used as a local variable in procedure clamped_spline

NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5US9jbGFtcGVkX3NwbGluZUYoLyUnZmFtaWx5R1EmVGltZXNGKC8lJXNpemVHUSMxMkYoLyUlYm9sZEdRJmZhbHNlRigvJSdpdGFsaWNHUSV0cnVlRigvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUShbMCwwLDBdRigvJSdvcGFxdWVHRjgvJStleGVjdXRhYmxlR0Y4LyUpcmVhZG9ubHlHRjgvJSljb21wb3NlZEdGOC8lKmNvbnZlcnRlZEdGOC8lK2ltc2VsZWN0ZWRHRjgvJSxwbGFjZWhvbGRlckdGOC8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUqbWF0aGNvbG9yR0ZELyUvbWF0aGJhY2tncm91bmRHRkcvJStmb250ZmFtaWx5R0YyLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGKC8lKW1hdGhzaXplR0Y1LUkjbW9HRiU2M1EjOj1GKC8lJWZvcm1HUSZpbmZpeEYoLyUmZmVuY2VHRjgvJSpzZXBhcmF0b3JHRjgvJSdsc3BhY2VHUS90aGlja21hdGhzcGFjZUYoLyUncnNwYWNlR0ZbcC8lKXN0cmV0Y2h5R0Y4LyUqc3ltbWV0cmljR0Y4LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOC8lLm1vdmFibGVsaW1pdHNHRjgvJSdhY2NlbnRHRjgvJTBmb250X3N0eWxlX25hbWVHRlgvJSVzaXplR0Y1LyUrZm9yZWdyb3VuZEdGRC8lK2JhY2tncm91bmRHRkctRiQ2KS1GX282M1ElcHJvY0YoL0Zjb1EhRihGZW9GZ28vRmpvUSQwZW1GKC9GXXBGXnJGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2JS1GX282M1EiKEYoL0Zjb1EncHJlZml4RigvRmZvRjtGZ28vRmpvUS50aGlubWF0aHNwYWNlRigvRl1wRmlyL0ZfcEY7RmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNistRiQ2JS1GLTY5USN4eEYoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvLUZfbzYzUSM6OkYoRltyRmVvRmdvRl1yRl9yRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYtNjlRJWxpc3RGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcby1GX282M1EiLEYoRmJvRmVvL0Zob0Y7Rl1yL0ZdcFEzdmVyeXRoaWNrbWF0aHNwYWNlRihGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2JS1GLTY5USJ5RihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GY3NGZnNGaXMtRiQ2JS1GLTY5USRmcDBGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Zjcy1GLTY5UShudW1lcmljRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GaXMtRiQ2JS1GLTY5USRmcG5GKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Zjc0ZpdEZpcy1GJDYlLUYtNjlRIlNGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Zjcy1GLTY5USVuYW1lRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG8tRl9vNjNRIilGKC9GY29RKHBvc3RmaXhGKEZnckZnb0Zoci9GXXBRMnZlcnl0aGlubWF0aHNwYWNlRihGW3NGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtSSdtc3BhY2VHRiU2Ji8lJ2hlaWdodEdRJzAuMH5leEYoLyUmd2lkdGhHUScwLjN+ZW1GKC8lJmRlcHRoR0Zldi8lKmxpbmVicmVha0dRMWZpcnN0cHJvY25ld2xpbmVGKC1GJDYmLUZfbzYzUSZsb2NhbEYoRltyRmVvRmdvRl1yL0ZdcFEwbWVkaXVtbWF0aHNwYWNlRihGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2SS1GLTY5USJORihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GaXMtRi02OVEiSUYoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvRmlzLUYtNjlRIlhGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Zpcy1GLTY5USJBRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GaXMtRi02OVEjQUFGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Zpcy1GLTY5USRGUDBGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Zpcy1GLTY5USRGUE5GKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Zpcy1GLTY5USJNRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GaXMtRi02OVEiSEYoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvRmlzLUYtNjlRI1hBRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GaXMtRi02OVEjWExGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Zpcy1GLTY5USNYVUYoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvRmlzLUYtNjlRI1haRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GaXMtRi02OVEiQ0YoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvRmlzLUYtNjlRIkpGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Zpcy1GLTY5USJCRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GaXMtRi02OVEiREYoRjBGM0Y2L0Y6RjhGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduL0ZqblEnbm9ybWFsRihGXG9GaXMtRi02OVEjU1NGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Zpcy1GLTY5USJMRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GaXMtRi02OVEiUEYoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvLUZfbzYzUSI7RihGYm9GZW9GXHRGXXJGXHBGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHFGYHYtRiQ2PS1GJDYlLUYkNiQtRiQ2KC1GX282M1EjaWZGKEZbckZlb0Znb0ZdckZjd0ZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYlLUYkNiUtRi02OVElbm9wc0YoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvLUZfbzYzUTAmQXBwbHlGdW5jdGlvbjtGKEZib0Zlb0Znb0ZdckZfckZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYlRmJyLUYkNiNGYHNGaXUtRl9vNjNRIzw+RihGYm9GZW9GZ28vRmpvRl92Rl52Rl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiVGaFxsRltdbC1GJDYlRmJyLUYkNiNGYXRGaXUtRl9vNjNRJXRoZW5GKEZbckZlb0Znby9Gam9GZHdGY3dGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRmF2NiZGY3ZGZnZGaXYvRlx3UTZpbmNyZWFzZWluZGVudG5ld2xpbmVGKC1GJDYlLUYtNjlRJkVSUk9SRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GW11sLUYkNiVGYnItRiQ2Iy1JI21zR0YlNiNRQUJvdGh+bGlzdHN+bXVzdH5iZX50aGV+c2FtZX5zaXplRihGaXUtRmF2NiZGY3YvRmd2UScwLjB+ZW1GKEZpdi9GXHdRNmRlY3JlYXNlaW5kZW50bmV3bGluZUYoLUZfbzYzUSdlbmR+aWZGKEZbckZlb0Znb0ZfXmxGX3JGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHFGZltsRmB2LUYkNiUtRiQ2JUZnd0Zeby1GJDYlRmZcbC1GX282M1EqfiZtaW51czt+RihGW3JGZW9GZ29GXXJGX3JGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtSSNtbkdGJTY5RmdwRjBGM0Y2Rmp6RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZbW2xGXG9GZltsRmB2LUYkNiUtRiQ2Jy1GJDYpLUZfbzYzUSRmb3JGKEZbckZlb0Znb0ZdckZjd0ZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcUZqdy1GX282M1ElZnJvbUYoRltyRmVvRmdvRl9ebEZjd0ZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GZGBsNjlRIjBGKEYwRjNGNkZqekY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GW1tsRlxvLUZfbzYzUSN0b0YoRltyRmVvRmdvRl9ebEZjd0ZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcUZndy1GX282M1EjZG9GKEZbckZlb0Znb0ZfXmxGY3dGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHFGYF5sLUYkNiQtRiQ2JS1GJDYlLUYkNiZGXXgtRl9vNjNRIltGKEZlckZnckZnb0ZockZqckZbc0ZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYjRmp3LUZfbzYzUSJdRihGXHZGZ3JGZ29GaHJGXnZGW3NGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHFGXm8tRiQ2JkZgc0ZjYmwtRiQ2Iy1GJDYlRmp3LUZfbzYzUSIrRihGYm9GZW9GZ29GX15sRmN3Rl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxRmNgbEZoYmxGZltsRmB2LUYkNiUtRiQ2JkZgeEZjYmxGZmJsRmhibEZeby1GJDYmRmF0RmNibEZdY2xGaGJsRmFfbC1GX282M1EnZW5kfmRvRihGW3JGZW9GZ29GX15sRl9yRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxRmZbbEZgdi1GJDYlLUYkNiVGZnhGXm9GZnRGZltsRmB2LUYkNiUtRiQ2JUZpeEZeb0ZedUZmW2xGYHYtRiQ2JS1GJDYlRlx5Rl5vLUYkNiVGZ3dGYGBsRmNgbEZmW2xGYHYtRiQ2JS1GJDYnLUYkNilGXGFsRmp3Rl9hbEZiYWxGZWFsRlx5RmhhbEZgXmwtRiQ2JS1GJDYmRl95RmNibEZmYmxGaGJsRl5vLUYkNiUtRiQ2JkZdeEZjYmxGXWNsRmhibEZgYGxGYWJsRmFfbEZqY2xGZltsRmB2LUYkNiUtRiQ2JS1GJDYmRmJ5RmNibC1GJDYjRmJhbEZoYmxGXm8tRiQ2JS1GJDYlLUYkNiVGYnItRiQ2JS1GZGBsNjlRJDMuMEYoRjBGM0Y2Rmp6RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZbW2xGXG8tRl9vNjNRIipGKEZib0Zlb0Znb0ZockZqckZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYlRmJyLUYkNiUtRiQ2JkZgeEZjYmwtRiQ2I0ZjYGxGaGJsRmBgbC1GJDYmRmB4RmNibEZfZmxGaGJsRml1Rml1LUZfbzYzUSIvRihGYm9GZW9GZ29GaHJGanJGW3NGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2Iy1GJDYmRl95RmNibEZfZmxGaGJsRmBgbC1GJDYlRmlmbEZcZ2xGZnhGZltsRmB2LUYkNiUtRiQ2JS1GJDYmRmJ5RmNibC1GJDYjRmd3RmhibEZeby1GJDYlLUYkNiVGaWZsRlxnbEZpeEZgYGwtRiQ2JS1GJDYlRmJyLUYkNiVGaWZsRlxnbC1GJDYlRmJyLUYkNiUtRiQ2JkZgeEZjYmxGaGhsRmhibEZgYGwtRiQ2JkZgeEZjYmwtRiQ2I0ZpZGxGaGJsRml1Rml1RmlnbC1GJDYjLUYkNiZGX3lGY2JsRlxqbEZoYmxGZltsRmB2LUYkNiUtRiQ2Jy1GJDYnRlxhbEZqd0ZlYWxGXHlGaGFsRmBebC1GJDYlLUYkNiZGYnlGY2JsRmZibEZoYmxGXm8tRiQ2JS1GJDYlRmJyLUYkNiVGaWZsRlxnbC1GJDYlRmJyLUYkNictRiQ2JS1GJDYmRmB4RmNibEZdY2xGaGJsRlxnbC1GJDYmRl95RmNibC1GJDYjLUYkNiVGandGYGBsRmNgbEZoYmxGYGBsLUYkNiVGZmNsRlxnbC1GJDYlRmJyLUYkNiVGZ2VsRmBgbC1GJDYmRl14RmNibEZcXG1GaGJsRml1RmFjbC1GJDYlLUYkNiZGYHhGY2JsRlxcbUZoYmxGXGdsRmNlbEZpdUZpdUZpZ2wtRiQ2JUZici1GJDYjLUYkNiVGY2VsRlxnbEZqW21GaXVGYV9sRmpjbEZmW2xGYHYtRiQ2JS1GJDYlLUYkNiZGZXlGY2JsRl9mbEZoYmxGXm8tRiQ2JS1GZGBsNjlRJDIuMEYoRjBGM0Y2Rmp6RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZbW2xGXG9GXGdsRl5obEZmW2xGYHYtRiQ2JS1GJDYlLUYkNiZGaHlGY2JsRl9mbEZoYmxGXm8tRmRgbDY5USQwLjVGKEYwRjNGNkZqekY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GW1tsRlxvRmZbbEZgdi1GJDYlLUYkNiUtRiQ2JkZbekZjYmxGX2ZsRmhibEZeby1GJDYlLUYkNiNGXWZsRmlnbC1GJDYjRmZdbUZmW2xGYHYtRiQ2JS1GJDYnRmZqbEZgXmwtRiQ2JS1GJDYlLUYkNiUtRiQ2JkZleUZjYmxGZmJsRmhibEZeby1GJDYlLUYkNiVGal1tRlxnbEZiXG1GYGBsLUYkNiVGalttRlxnbC1GJDYmRmh5RmNibEZcXG1GaGJsRmZbbEZgdi1GJDYlLUYkNiUtRiQ2JkZoeUZjYmxGZmJsRmhibEZeby1GJDYlLUYkNiNGY2VsRmlnbC1GJDYjRlxgbUZmW2xGYHYtRiQ2JS1GJDYmRlt6RmNibEZmYmxGaGJsRl5vLUYkNiUtRiQ2JUZici1GJDYjLUYkNiVGampsRmBgbC1GJDYlRmpbbUZcZ2wtRiQ2JkZbekZjYmxGXFxtRmhibEZpdUZpZ2xGYGFtRmFfbEZqY2xGZltsRmB2LUYkNiUtRiQ2JS1GJDYmRmV5RmNibEZoaGxGaGJsRl5vLUYkNiVGYGpsRlxnbC1GJDYlRmJyLUYkNiVGal1tRmBgbC1GJDYmRmh5RmNibEZcamxGaGJsRml1RmZbbEZgdi1GJDYlLUYkNiUtRiQ2JkZbekZjYmxGaGhsRmhibEZeby1GJDYlLUYkNiVGYnItRiQ2Iy1GJDYlRmZobEZgYGwtRiQ2JUZgamxGXGdsLUYkNiZGW3pGY2JsRlxqbEZoYmxGaXVGaWdsLUYkNiNGZmJtRmZbbEZgdi1GJDYlLUYkNiUtRiQ2JkZeekZjYmxGaGhsRmhibEZeb0ZkY21GZltsRmB2LUYkNiUtRiQ2Jy1GJDYnRlxhbEZqd0ZlYWxGZ3dGaGFsRmBebC1GJDYmLUYkNiUtRiQ2JUZhekZeby1GJDYlRmd3RmBgbEZqd0ZmW2xGYHYtRiQ2JS1GJDYlLUYkNiZGXnpGY2JsLUYkNiNGYXpGaGJsRl5vLUYkNiUtRiQ2JkZbekZjYmxGXmZtRmhibEZgYGwtRiQ2JS1GJDYmRmh5RmNibEZeZm1GaGJsRlxnbC1GJDYmRl56RmNibC1GJDYjLUYkNiVGYXpGYWNsRmNgbEZoYmxGZltsRmB2LUYkNiUtRiQ2JS1GJDYmRmR6RmNibEZeZm1GaGJsRl5vLUYkNiUtRiQ2JS1GJDYlRmJyLUYkNiMtRiQ2JS1GJDYmRmB4RmNibEZqZm1GaGJsRmBgbC1GJDYmRmB4RmNibEZeZm1GaGJsRml1RmlnbC1GJDYjLUYkNiZGX3lGY2JsRl5mbUZoYmxGYGBsLUYkNiUtRiQ2JUZici1GJDYlRmRobUZcZ2wtRiQ2JUZici1GJDYlRmhmbUZhY2wtRiQ2JUZqXW1GXGdsRlxmbUZpdUZpdUZpZ2wtRiQ2I0ZpZmxGZltsRmB2LUYkNiUtRiQ2JkZnekZjYmxGXmZtRmhibEZeby1GJDYlLUYkNiVGYnItRiQ2Iy1GJDYlRmhmbUZgYGxGXGZtRml1RmlnbC1GJDYlRmJyLUYkNiVGaWZsRlxnbEZkaG1GaXVGYV9sRmpjbEZmW2xGYHYtRiQ2JS1GJDYlLUYtNjlRJ3ByaW50ZkYoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvRltdbC1GJDYlRmJyLUYkNiMtRi02OVFqb2B8K3wrVGhlfmNvZWZmaWNpZW50c35vZn50aGV+Y2xhbXBlZH5jdWJpY35zcGxpbmV+U1tpXXwrb25+dGhlfnN1YmludGVydmFsc35hcmU6fCtgRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GaXVGZltsRmB2LUYkNiUtRiQ2JUZoam1GW11sLUYkNiVGYnItRiQ2Iy1GLTY5UWZuYH5+fn5+fn5+fn5hKGkpfn5+fn5+fn5+fmIoaSl+fn5+fn5+fn5+YyhpKX5+fn5+fn5+fn5kKGkpfCtgRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GaXVGZltsRmB2LUYkNiUtRiQ2J0ZfZWxGYF5sLUYkNiVGaGptRltdbC1GJDYlRmJyLUYkNistRi02OVFBYH4lMTMuOGZ+JTEzLjhmfiUxMy44Zn4lMTMuOGZ+fCtgRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GaXNGZmNsRmlzLUYkNiZGZHpGY2JsRmZibEZoYmxGaXMtRiQ2JkZeekZjYmxGZmJsRmhibEZpcy1GJDYmRmd6RmNibEZmYmxGaGJsRml1RmFfbEZqY2xGZltsRmB2LUYkNiUtRiQ2JUZoam1GW11sLUYkNiVGYnItRiQ2Iy1GLTY5USN8K3wrRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GaXVGZltsRmB2LUYkNiUtRiQ2J0ZfZWxGYF5sLUYkNiUtRiQ2JkZdW2xGY2JsRmZibEZoYmxGXm8tRiQ2KUZmY2xGYWNsLUYkNiVGalxuRlxnbC1GJDYlRmJyLUYkNiUtRi02OVEieEYoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvRmBgbEZhYmxGaXVGYWNsLUYkNiVGXF1uRlxnbC1GJDYlRmdebi1GX282M1EiXkYoRmJvRmVvRmdvRmVdbEZedkZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GZGBsNjlRIjJGKEYwRjNGNkZqekY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GW1tsRlxvRmFjbC1GJDYlRl5dbkZcZ2wtRiQ2JUZnXm5GYl9uLUZkYGw2OVEiM0YoRjBGM0Y2Rmp6RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZbW2xGXG9GYV9sRmpjbEZmW2xGYHYtRiQ2JS1GJDYlRmBbbEZeby1GJDYlRmNibC1GJDYjLUYkNiUtRi02OVEkc2VxRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GW11sLUYkNiVGYnItRiQ2JS1GJDYlRmNibC1GJDYlLUYkNiVGW19uLUZfbzYzUSI8RihGYm9GZW9GZ29GaW9GXHBGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHFGZ2VsRmlzRmFebkZoYmxGaXMtRiQ2JUZqdy1GX282M1EiPUYoRmJvRmVvRmdvRmlvRlxwRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiVGYmFsLUZfbzYzUSMuLkYoRlx2RmVvRmdvRl9ebEZfckZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcUZceUZpdUZoYmxGZltsRmB2LUYkNiUtRiQ2JUZjW2xGXm8tRiQ2JUZjYmwtRiQ2Iy1GJDYlRmlgbkZbXWwtRiQ2JUZici1GJDYlLUYkNiUtRi02OVEjb3BGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0ZbXWwtRiQ2JUZici1GJDYjLUYkNiZGYFtsRmNibEZmYmxGaGJsRml1RmlzLUYkNiVGandGW2JuLUYkNiVGY2BsRmBibi1GJDYlRlx5RmFjbEZjYGxGaXVGaGJsRmZbbEZgdi1GJDYlLUYkNiVGY1tsRl5vLUYkNiUtRi02OVEnc3Vic29wRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GW11sLUYkNiVGYnItRiQ2JS1GJDYlLUYkNiUtRiQ2JUZlX25GXGdsRlx5RmFjbEZjYGxGW2JuLUYtNjlRJU5VTExGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Zpc0ZjW2xGaXVGZltsRmB2LUYkNiVGY3VGXm8tRiQ2JS1GLTY5USlzaW1wbGlmeUYoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvRltdbC1GJDYlRmJyLUYkNiMtRiQ2JS1GLTY5USpwaWVjZXdpc2VGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0ZbXWwtRiQ2JUZici1GJDYjLUYkNiVGY2NuRltdbC1GJDYlRmJyLUYkNiNGY1tsRml1Rml1Rml1RmFfbC1GX282M1EpZW5kfnByb2NGKEZbckZlb0Znb0ZfXmxGX3JGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHE3Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3Iz5JL2NsYW1wZWRfc3BsaW5lR0YoZio2JydJI3h4R0YoSSVsaXN0RyUqcHJvdGVjdGVkRydJInlHRihGYWhuJ0kkZnAwR0YoSShudW1lcmljR0ZiaG4nSSRmcG5HRihGZ2huJ0kiU0dGKEklbmFtZUdGYmhuNjZJIk5HRihJIklHRihJIlhHRihJIkFHRihJI0FBR0YoSSRGUDBHRihJJEZQTkdGKEkiTUdGKEkiSEdGKEkjWEFHRihJI1hMR0YoSSNYVUdGKEkjWFpHRihJIkNHRihJIkpHRihJIkJHRihJIkRHRihJI1NTR0YoSSJMR0YoSSJQR0YoRihGKEM9QCQwLUklbm9wc0dGYmhuNiNGYGhuLUZmam42I0ZkaG4tSSZFUlJPUkdGYmhuRl9fbD5GXmluLCZGZWpuIiIiRl5bbyEiIj8oRl9pbiIiIUZeW29GXmluSSV0cnVlR0ZiaG5DJD4mRmBpbjYjRl9pbiZGYGhuNiMsJkZfaW5GXltvRl5bb0ZeW28+JkZhaW5GZltvJkZkaG5GaFtvPkZjaW5GZmhuPkZkaW5GaWhuPkZlaW4sJkZeaW5GXltvRl5bb0ZfW28/KEZfaW5GYVtvRl5bb0ZlaW5GYltvPiZGZmluRmZbbywmJkZgaW5GaFtvRl5bb0ZlW29GX1tvPiZGZ2luNiNGYVtvLCYqKCQiI0lGX1tvRl5bbywmJkZhaW42I0ZeW29GXltvJkZhaW5GaFxvRl9bb0ZeW28mRmZpbkZoXG9GX1tvRl5bbyomRltdb0ZeW29GY2luRl5bb0ZfW28+JkZnaW42I0ZeaW4sJiomRltdb0ZeW29GZGluRl5bb0ZeW28qKEZbXW9GXltvLCYmRmFpbkZlXW9GXltvJkZhaW42I0ZgXG9GX1tvRl5bbyZGZmluRlxeb0ZfW29GX1tvPyhGX2luRl5bb0ZeW29GZWluRmJbbz4mRmdpbkZmW28qKEZbXW9GXltvLCgqJiZGYWluRmhbb0ZeW28mRmZpbjYjLCZGX2luRl5bb0ZeW29GX1tvRl5bb0ZeW28qJkZbXG9GXltvLCZGZVxvRl5bbyZGYGluRmZeb0ZfW29GXltvRl9bbyomJkZhaW5GZl5vRl5bb0ZjXG9GXltvRl5bb0ZeW28qJkZjXG9GXltvRmVeb0ZeW29GX1tvPiZGaGluRmhcbyomJCIjP0ZfW29GXltvRmFdb0ZeW28+JkZpaW5GaFxvJCIiJkZfW28+JkZqaW5GaFxvKiZGZ1xvRl5bb0ZfX29GX1tvPyhGX2luRl5bb0ZeW29GZWluRmJbb0MlPiZGaGluRmZbbywmKiZGYV9vRl5bb0ZpXm9GXltvRl5bbyomRmVeb0ZeW28mRmlpbkZmXm9GXltvRl9bbz4mRmlpbkZmW28qJkZjXG9GXltvRl1gb0ZfW28+JkZqaW5GZltvKiYsJkZgXm9GXltvKiZGZV5vRl5bbyZGamluRmZeb0ZeW29GX1tvRl5bb0ZdYG9GX1tvPiZGaGluRmVdbyomRl1eb0ZeW28sJkZhX29GXltvJkZpaW5GXF5vRl9bb0ZeW28+JkZqaW5GZV1vKiYsJkZkXW9GXltvKiZGXV5vRl5bbyZGamluRlxeb0ZeW29GX1tvRl5bb0ZcYW9GX1tvPiZGW2puRmVdb0ZhYW8/KEZfaW5GXltvRl5bb0ZeaW5GYltvQyY+RlxqbiwmRl5pbkZeW29GX2luRl9bbz4mRltqbjYjRlxqbiwmJkZqaW5GXmJvRl5bbyomJkZpaW5GXmJvRl5bbyZGW2puNiMsJkZcam5GXltvRl5bb0ZeW29GXltvRl9bbz4mRl1qbkZeYm8sJiomLCYmRmFpbkZkYm9GXltvJkZhaW5GXmJvRl9bb0ZeW28mRmZpbkZeYm9GX1tvRl5bbyooRl1jb0ZeW28sJkZjYm9GXltvKiZGYV9vRl5bb0ZdYm9GXltvRl5bb0ZeW29GW11vRl9bb0ZfW28+JkZeam5GXmJvKigsJkZjYm9GXltvRl1ib0ZfW29GXltvRltdb0ZfW29GXWNvRl9bby1JJ3ByaW50ZkdGKDYjSWhvfCt8K1RoZX5jb2VmZmljaWVudHN+b2Z+dGhlfmNsYW1wZWR+Y3ViaWN+c3BsaW5lflNbaV18K29ufnRoZX5zdWJpbnRlcnZhbHN+YXJlOnwrR0YoLUZmY282I0lafn5+fn5+fn5+fmEoaSl+fn5+fn5+fn5+YihpKX5+fn5+fn5+fn5jKGkpfn5+fn5+fn5+fmQoaSl8K0dGKD8oRl9pbkZhW29GXltvRmVpbkZiW28tRmZjbzYnST9+JTEzLjhmfiUxMy44Zn4lMTMuOGZ+JTEzLjhmfnwrR0YoRltcbyZGXWpuRmZbbyZGW2puRmZbbyZGXmpuRmZbby1GZmNvNiNJI3wrfCtHRig/KEZfaW5GYVtvRl5bb0ZlaW5GYltvPiZGX2puRmZbbywqRltcb0ZeW28qJkZgZG9GXltvLCZJInhHRihGXltvRmVbb0ZfW29GXltvRl5bbyomRmFkb0ZeW28pRltlbyIiI0ZeW29GXltvKiZGYmRvRl5bbylGW2VvIiIkRl5bb0ZeW28+RmBqbjcjLUkkc2VxR0ZiaG42JDckMkZcZW9GZVxvRmhkby9GX2luO0ZhW29GZWluPkZham43Iy1GZmVvNiQtSSNvcEdGYmhuNiMmRmBqbkZmW28vRl9pbjtGXltvLCZGZWluRl5bb0ZeW29GXltvPkZham4tSSdzdWJzb3BHRmJobjYkLywmKiZGX2VvRl5bb0ZlaW5GXltvRl5bb0ZeW29GXltvSSVOVUxMR0ZiaG5GYWpuPkZbaW4tSSlzaW1wbGlmeUdGKDYjLUkqcGllY2V3aXNlR0ZiaG42Iy1GYWZvNiNGYWpuRihGKEYoNyNGXWhu

clamped_spline_dir:=proc() printf(`clamped_spline returns the piecewise cubic spline.\n\n`); printf(`The output also includes a table of values\n`); printf(`for the cubic polynomials\n\n`); printf(`S[i](x)=a[i]+b[i](x-x[i])+c[i](x-x[i])^2+d[i](x-x[i])^3\n\n`); printf(`for i=0..n-1 making up the spline.\n\n`); printf(`The arguments for clamped_spline are:\n`); printf(`(1)the list of x-values\n`); printf(`(2)the list of f(x) or y-values\n`); printf(`(3)f'(x[0])\n`); printf(`(4)f'(x[n])\n`); printf(`(5)the variable for returning the piecewise cubic spline\n\n`); printf(`If assigning the result to a variable, have the\n`); printf(`variable and the 5th argument the same.\n\n`); printf(`If S is the variable for returning the spline\n`); printf(`and has already been given a value,\n`); printf(`the procedure should be preceded by the statement:\n`); printf(`S:='S'`); end;

NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5UTNjbGFtcGVkX3NwbGluZV9kaXJGKC8lJ2ZhbWlseUdRJlRpbWVzRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR1EldHJ1ZUYoLyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EoWzAsMCwwXUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y4LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGRC8lL21hdGhiYWNrZ3JvdW5kR0ZHLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1EnaXRhbGljRigvJSltYXRoc2l6ZUdGNS1JI21vR0YlNjNRIzo9RigvJSVmb3JtR1EmaW5maXhGKC8lJmZlbmNlR0Y4LyUqc2VwYXJhdG9yR0Y4LyUnbHNwYWNlR1EvdGhpY2ttYXRoc3BhY2VGKC8lJ3JzcGFjZUdGW3AvJSlzdHJldGNoeUdGOC8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR0ZYLyUlc2l6ZUdGNS8lK2ZvcmVncm91bmRHRkQvJStiYWNrZ3JvdW5kR0ZHLUYkNigtRl9vNjNRJXByb2NGKC9GY29RIUYoRmVvRmdvL0Zqb1EkMGVtRigvRl1wRl5yRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiUtRl9vNjNRIihGKC9GY29RJ3ByZWZpeEYoL0Zmb0Y7RmdvL0Zqb1EudGhpbm1hdGhzcGFjZUYoL0ZdcEZpci9GX3BGO0ZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYjLUYtNjlGXHJGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG8tRl9vNjNRIilGKC9GY29RKHBvc3RmaXhGKEZnckZnb0Zoci9GXXBRMnZlcnl0aGlubWF0aHNwYWNlRihGW3NGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtSSdtc3BhY2VHRiU2Ji8lJ2hlaWdodEdRJzAuMH5leEYoLyUmd2lkdGhHUScwLjN+ZW1GKC8lJmRlcHRoR0ZcdC8lKmxpbmVicmVha0dRMWZpcnN0cHJvY25ld2xpbmVGKC1GJDYzLUYkNiUtRiQ2JS1GLTY5USdwcmludGZGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcby1GX282M1EwJkFwcGx5RnVuY3Rpb247RihGYm9GZW9GZ29GXXJGX3JGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2JUZici1GJDYjLUYtNjlRV2BjbGFtcGVkX3NwbGluZX5yZXR1cm5zfnRoZX5waWVjZXdpc2V+Y3ViaWN+c3BsaW5lLnwrfCtgRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GYHMtRl9vNjNRIjtGKEZib0Zlby9GaG9GO0ZdckZccEZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcUZncy1GJDYlLUYkNiVGW3VGXnUtRiQ2JUZici1GJDYjLUYtNjlRTmBUaGV+b3V0cHV0fmFsc29+aW5jbHVkZXN+YX50YWJsZX5vZn52YWx1ZXN8K2BGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Zgc0ZodUZncy1GJDYlLUYkNiVGW3VGXnUtRiQ2JUZici1GJDYjLUYtNjlRPmBmb3J+dGhlfmN1YmljfnBvbHlub21pYWxzfCt8K2BGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Zgc0ZodUZncy1GJDYlLUYkNiVGW3VGXnUtRiQ2JUZici1GJDYjLUYtNjlRZm5gU1tpXSh4KT1hW2ldK2JbaV0oeC14W2ldKStjW2ldKHgteFtpXSleMitkW2ldKHgteFtpXSleM3wrfCtgRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GYHNGaHVGZ3MtRiQ2JS1GJDYlRlt1Rl51LUYkNiVGYnItRiQ2Iy1GLTY5UUdgZm9yfmk9MC4ubi0xfm1ha2luZ351cH50aGV+c3BsaW5lLnwrfCtgRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GYHNGaHVGZ3MtRiQ2JS1GJDYlRlt1Rl51LUYkNiVGYnItRiQ2Iy1GLTY5UUlgVGhlfmFyZ3VtZW50c35mb3J+Y2xhbXBlZF9zcGxpbmV+YXJlOnwrYEYoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvRmBzRmh1RmdzLUYkNiUtRiQ2JUZbdUZedS1GJDYlRmJyLUYkNiMtRi02OVE7YCgxKXRoZX5saXN0fm9mfngtdmFsdWVzfCtgRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GYHNGaHVGZ3MtRiQ2JS1GJDYlRlt1Rl51LUYkNiVGYnItRiQ2Iy1GLTY5UUNgKDIpdGhlfmxpc3R+b2Z+Zih4KX5vcn55LXZhbHVlc3wrYEYoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvRmBzRmh1RmdzLUYkNiUtRiQ2JUZbdUZedS1GJDYlRmJyLUYkNiMtRi02OVEvYCgzKWYnKHhbMF0pfCtgRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GYHNGaHVGZ3MtRiQ2JS1GJDYlRlt1Rl51LUYkNiVGYnItRiQ2Iy1GLTY5US9gKDQpZicoeFtuXSl8K2BGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Zgc0ZodUZncy1GJDYlLUYkNiVGW3VGXnUtRiQ2JUZici1GJDYjLUYtNjlRZ25gKDUpdGhlfnZhcmlhYmxlfmZvcn5yZXR1cm5pbmd+dGhlfnBpZWNld2lzZX5jdWJpY35zcGxpbmV8K3wrYEYoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvRmBzRmh1RmdzLUYkNiUtRiQ2JUZbdUZedS1GJDYlRmJyLUYkNiMtRi02OVFTYElmfmFzc2lnbmluZ350aGV+cmVzdWx0fnRvfmF+dmFyaWFibGUsfmhhdmV+dGhlfCtgRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GYHNGaHVGZ3MtRiQ2JS1GJDYlRlt1Rl51LUYkNiVGYnItRiQ2Iy1GLTY5UUxgdmFyaWFibGV+YW5kfnRoZX41dGh+YXJndW1lbnR+dGhlfnNhbWUufCt8K2BGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Zgc0ZodUZncy1GJDYlLUYkNiVGW3VGXnUtRiQ2JUZici1GJDYjLUYtNjlRUWBJZn5TfmlzfnRoZX52YXJpYWJsZX5mb3J+cmV0dXJuaW5nfnRoZX5zcGxpbmV8K2BGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Zgc0ZodUZncy1GJDYlLUYkNiVGW3VGXnUtRiQ2JUZici1GJDYjLUYtNjlRR2BhbmR+aGFzfmFscmVhZHl+YmVlbn5naXZlbn5hfnZhbHVlLHwrYEYoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvRmBzRmh1RmdzLUYkNiUtRiQ2JUZbdUZedS1GJDYlRmJyLUYkNiMtRi02OVFWYHRoZX5wcm9jZWR1cmV+c2hvdWxkfmJlfnByZWNlZGVkfmJ5fnRoZX5zdGF0ZW1lbnQ6fCtgRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GYHNGaHVGZ3MtRiQ2JUZbdUZedS1GJDYlRmJyLUYkNiMtRi02OVEpYFM6PSdTJ2BGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Zgcy1GaHM2JkZqcy9GXnRRJzAuMH5lbUYoRmB0L0ZjdFE2ZGVjcmVhc2VpbmRlbnRuZXdsaW5lRigtRl9vNjNRKWVuZH5wcm9jRihGW3JGZW9GZ28vRmpvUTBtZWRpdW1tYXRoc3BhY2VGKEZfckZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcTcjLV9GKUksbXByaW50c2xhc2hHRig2JDcjPkkzY2xhbXBlZF9zcGxpbmVfZGlyR0YoZipGKEYoRihGKEMzLUkncHJpbnRmR0YoNiNJVWNsYW1wZWRfc3BsaW5lfnJldHVybnN+dGhlfnBpZWNld2lzZX5jdWJpY35zcGxpbmUufCt8K0dGKC1GYGJsNiNJTFRoZX5vdXRwdXR+YWxzb35pbmNsdWRlc35hfnRhYmxlfm9mfnZhbHVlc3wrR0YoLUZgYmw2I0k8Zm9yfnRoZX5jdWJpY35wb2x5bm9taWFsc3wrfCtHRigtRmBibDYjSVpTW2ldKHgpPWFbaV0rYltpXSh4LXhbaV0pK2NbaV0oeC14W2ldKV4yK2RbaV0oeC14W2ldKV4zfCt8K0dGKC1GYGJsNiNJRWZvcn5pPTAuLm4tMX5tYWtpbmd+dXB+dGhlfnNwbGluZS58K3wrR0YoLUZgYmw2I0lHVGhlfmFyZ3VtZW50c35mb3J+Y2xhbXBlZF9zcGxpbmV+YXJlOnwrR0YoLUZgYmw2I0k5KDEpdGhlfmxpc3R+b2Z+eC12YWx1ZXN8K0dGKC1GYGJsNiNJQSgyKXRoZX5saXN0fm9mfmYoeCl+b3J+eS12YWx1ZXN8K0dGKC1GYGJsNiNJLSgzKWYnKHhbMF0pfCtHRigtRmBibDYjSS0oNClmJyh4W25dKXwrR0YoLUZgYmw2I0llbig1KXRoZX52YXJpYWJsZX5mb3J+cmV0dXJuaW5nfnRoZX5waWVjZXdpc2V+Y3ViaWN+c3BsaW5lfCt8K0dGKC1GYGJsNiNJUUlmfmFzc2lnbmluZ350aGV+cmVzdWx0fnRvfmF+dmFyaWFibGUsfmhhdmV+dGhlfCtHRigtRmBibDYjSUp2YXJpYWJsZX5hbmR+dGhlfjV0aH5hcmd1bWVudH50aGV+c2FtZS58K3wrR0YoLUZgYmw2I0lPSWZ+U35pc350aGV+dmFyaWFibGV+Zm9yfnJldHVybmluZ350aGV+c3BsaW5lfCtHRigtRmBibDYjSUVhbmR+aGFzfmFscmVhZHl+YmVlbn5naXZlbn5hfnZhbHVlLHwrR0YoLUZgYmw2I0lUdGhlfnByb2NlZHVyZX5zaG91bGR+YmV+cHJlY2VkZWR+Ynl+dGhlfnN0YXRlbWVudDp8K0dGKC1GYGJsNiNJJ1M6PSdTJ0dGKEYoRihGKDcjRl1ibA==