企业的营销管理问题-毕业论文.doc

企业的营销管理问题
（物理系：段文瑞 数学系：张晔 韩莹）
摘要
本文属于企业的营销管理优化问题， 主要解决的是公司最大利润、营销部的最大收入等问题。通过对影响公司利润，营销部收入的因素如公司的生产成本、宣传费用、意向订单、客户资料等的分析以及对问题的整体把握，由于最大生产能力的制约，我们需要在满足已签约的销售合同的基础上，对意向签约量有选择的安排生产，使公司的利润达到最大，而营销部应尽力争取与意向客户的签订合同，争取多销售产品时，可使总收入极大化。我们先假设在理想情况下即客户意向为100%情况下求解出方案，然后一一考虑不同客户的意向需求量，根据不同商家的期望值，及产品的不同利润有选择的进行生产。本问题中我们使用非线性规划求最优解的方法并利用MATLAB和LinGo软件进行求解。得出问题（1）中公司的利润最大为 1209.685元 ；问题（2）营销部总收入极大为 405.9621元；问题（3）兼顾公司和营销部二者的利益为 元；问题（4）兼顾公司、营销部的利益以及客户的需求，尽量做到均衡销售时的利润最大为 元；问题（4）使营销部总收入最大的定价、生产及销售
关键字 利润 非线性规划 最优解 拟合

问题重述
任何一个企业都面临着将产品推向市场，销售出去。一般的，企业下设营销部，通过营销部对企业生产的产品进行销售。由于企业生产能力的制约，需要在满足已签约的销售合同量的基础上，对意向签约量有选择的安排生产。一方面，企业会组织安排生产，完成已签约的销售合同；另一方面，企业希望营销部门尽力争取与意向签约的客户签订正式销售合同；除此之外，企业还希望营销部门努力再多销售一些产品（计划外）。
对于所签约的销售合同和意向签约量（计划内），企业根据销售量向营销部发放经费（包括工资以及宣传费用等）；对于计划外销售的产品，实行承包制，营销部向企业缴纳利润。产品的生产费用由企业承担，与销售有关的费用（如产品的宣传费用等）由营销部承担。
宏宇电器公司2011年计划生产三类10种小家电，其中包括：热水壶（1.5升，1.8升，2升）、豆浆机（0.9升，1.1升，1.3升）、电饭煲（2升，2.5升，3升，3.5升）。
三类小家电的年最大生产能力分别为：热水壶：5万个；豆浆机：6.5万个；电饭煲6.2万个。
相关的各类数据见附录一（表1至表14）
分别按以下要求为宏宇电器公司制定相应的生产、销售方案：
（1）使公司的利润达到最大；
（2）使营销部的总收入极大化；
（3）兼顾公司和营销部二者的利益；
（4）兼顾公司、营销部的利益以及客户的需求，尽量做到均衡销售；
（5）公司为了调动营销部的积极性，为公司多创利润，规定对于计划外销售的产品，营销部可以自行定价。统计数据显示，自行定价后每单位产品的销售价格上涨或下降5％，则三类小家电的可能销售量相应的减少或增加10％。试确定使营销部总收入最大的定价、生产及销售方案。
问题分析
问题（1）（2）中要求根据给出的相关数据建立模型求出最佳的生产、销售方案，分别使公司的利润达到最大和营销部的总收入极大化。
通过对表分析，我们知道产品在计划内销售时带给公司的利润要大于将其在计划外销售所带来的利润益，那么必然要使计划内的订单订购量足够大（但不能超出最大生产能力），我们假设按订单的最大意向签约量签约有签约意向的客户，此时计划内的订够量最大。
对于问题（3）在兼顾公司和营销部二者的利益时，把两者看作一个整体建立模型求出最佳的生产、销售方案，那么此时就不用考虑营销部给公司上缴的那一部分利润，只考虑销售的总利润以及宣传费用
问题（4）当兼顾公司、营销部的利益以及客户的需求时，需要考虑顾客的要求，即要求客户的需求概率，这里用数学期望来解决
问题（5）产品的销量与产品单价有一定的关系，故提高或降低价格会减少或增加销量，利用MATLAB求出目标函数的最优解，即营销部的最大总收入，同时求出每种产品价格上涨率的具体值，从而确定营销部的定价方案。
模型假设
（1）假设客户的订购量就是生产量，生产量就是销售量；
（2）与销售有关的经费只考虑宣传费用，其他费用不考虑；
（3）营销部除了只给企业上缴利润和承担宣传费用外，无其他费用支出；
（4）假设问题（1）（2）（3）中，将意向客户的需求概率设为100%；
（5）该企业周围的各种因素都是稳定的
符号说明
xi 第i种产品的总销售量 i=1，2···10
yi 计划外第i种产品的销售量 i=1，2···10
zi 第i中产品的已签约量 i=1，2···10
pi 第i中产品的意向签约量 i=1，2···10
ki 第i种产品可能签订订购合同的最大合同量； i=1，2···10
Qi(xi) 第i种产品的宣传费函数 i=1，2···10
Si(xi) 第i种产品的销售额函数 i=1，2···10
Ui 第i种产品的的需求期望 i=1，2···10
Mi（yi） 计划外第i中产品的销售额 i=1，2···10
Ti(yi ) 计划外第i种产品的上缴利润函数 i=1，2···10
Fi(xi) 第i种产品的成本函数 i=1，2···10
Si(xi-yi) 计划内第i种产品的销售额函数 i=1，2···10
Fi(xi-yi) 计划内第i种产品的成本函数 i=1，2···10
Ni(xi-yi) 计划内第i种产品的销售部经费函数 i=1，2···10
Hi（yi(1-0.1S)） 计划外的第i种产品 价格改变后的销售额函数 i=1，2···10
Bi（yi(1-0.1S)） 价格改变后的第i种产品的上缴利润函数 i=1，2···10
模型的建立与求解
（一）问题（1）我们可以知道公司利润就是计划内的利润和计划外的利润之和，那公司的赢利可以表示为：
公司利润＝计划内销售额+计划外上缴利润－计划内成本－计划外营销部经费
我们分别分析每种产品产量与成本关系表（表2），每种小家电的销售额随订购量变化(表4），计划外销售部分营销部向企业缴纳利润表(表14 );
用Matlab软件（运用三次拟合）编写程序 ，可以得出 第i种产品的销售额函数 Si(xi-yi)；计划外第i种产品的上缴利润函数Ti(yi )；计划内第i种产品的成本函数Fi(xi-yi)；计划外销售部经费函数Ni(yi)
根据述各种准备可建立目标函数
M1=
约束条件
（1）因为各种产品的生产量不能超过其最大生产能力，所以可得:
热水壶：x1+x2+x3<=50
豆浆机：x4+x5+x6 <=65
电饭煲：x7+x8+x9+x10 <=62
（2）每种产品计划内的销售量应该介于已签约的量，已签约量和意向签约量之间:
xi-yi>=zi，xi-yi<=zi+pi
（3）计划外的每种产品的销售量应该小于可能签订订购合同的最大合同量:
yi<=ki
（4）每种产品总的销售额小于已签约量，意向签约量，计划外销售量三者之和：
xi<=zi+pi+ki
然后应用LinGo软件编程（程序见附录二）
可得结果公司的利润最大为1209.685元
（二）问题（2）与问题1相似，为了使销售部的利润达到极大化，同样需要建立一个目标最大化模型，：
营销部的收入＝经费+计划外销售额—计划外成本－上缴利润－宣传费
约束条件
（1）而各种产品的生产量不能超过其最大生产能力
热水壶：x1+x2+x3<=50
豆浆机：x4+x5+x6 <=65
电饭煲：x7+x8+x9+x10 <=62
（2）每种产品计划内的销售量应该介于已签约的量，已签约量和意向签约量之间:
xi-yi>=zi，xi-yi<=zi+pi
（3）计划外的每种产品的销售量应该小于可能签订订购合同的最大合同量:
yi<=ki
（3）每种产品总的销售额小于已签约量，意向签约量，计划外销售量三者之和：
xi<=zi+pi+ki
然后应用LinGo软件编程（程序见附录三）得出结果
营销部的总收入极大为 405.9621 元。
（三）问题（3）将公司和营销部看作一个整体来考虑,不用考虑营销部给公司上缴的那一部分利润，只考虑销售的总利润以及宣传费用，那么利润可以表示为：
公司和营销部兼顾时的利润为＝总销售额－总成本－计划外产品的宣传费
我们分别分析每种产品产量与成本关系表（表2），每种小家电的宣传费用随着销售量表(3)，每种小家电的销售额随订购量变化表（表4） 用Matlab软件（运用三次拟合）编写程序可得：第i种产品的宣传费函数 Qi(xi)； 第i种产品的销售额函数Si(xi) ；第i种产品的成本函数Fi(xi) ，从而可得目标函数：
M3=
约束条件
（1）因为各种产品的生产量不能超过其最大生产能力，所以有：
热水壶：x1+x2+x3<=50
豆浆机：x4+x5+x6 <=65
电饭煲：x7+x8+x9+x10 <=62
（2）每种产品计划内的销售量应该介于已签约的量，已签约量和意向签约量之间:
xi-yi>=zi，xi-yi<=zi+pi
（3）计划外的每种产品的销售量应该小于可能签订订购合同的最大合同量:
yi<=ki
（4）计划外的每种产品的总销售量总是不大于计划内每种产品的总销售额：
yi<=xi
然后应用LinGo软件编程（程序见附录四）得出结果：
兼顾公司和营销部的利益时，最大利润为 2155.776 元。
（四）问题（4）兼顾公司、营销部的利益以及客户的需求，尽量做到均衡销售时，我们要考虑到客户的需求，即要求客户的需求概率，这里用数学期望来解决，在问题（3）基础上，约束条件应该考虑客户期望，那么目标函数为：M4=
约束条件
（1）因为各种产品的生产量不能超过其最大生产能力，所以有：
热水壶：x1+x2+x3<=50
豆浆机：x4+x5+x6 <=65
电饭煲：x7+x8+x9+x10 <=62
（2）每种产品计划内的销售量应该介于已签约的量，已签约量和意向签约量之间:
xi-yi>=zi，xi-yi<=zi+pi
（3）通过对意向用户和计划外需求的概率表分析，根据期望计算公式
E(X) = X1*p(X1) + X2*p(X2) + …… + Xn*p(Xn)
X1,X2,X3,……,Xn为这几个数据，p(X1),p(X2),p(X3),……p(Xn)为这几个数据的概率函数
表 各产品的需求期望单位（千件）
产品
1
2
3
4
5
6
7
8
9
10
期望
3.650
3.250
3.050
3.950
4.750
1.800
4.300
8.000
3.700
3.900
所以 xi<=zi+ui
然后应用LinGo软件编程（程序见附录五）得出结果兼顾公司营销部的利益以及客户的需求，尽量做到均衡销售时，利润最大为 1025.885 元。
（五）问题（5）营销部在计划外自行定价，在问题（4）的基础上，计划内的订购量已经确定，只是计划外的销售量随定价的变化变化，用Matlab软件（运用三次拟合）编写程序可得
M5=
约束条件
（1）因为各种产品的生产量不能超过其最大生产能力，所以有：
热水壶：x1+x2+x3<=50
豆浆机：x4+x5+x6 <=65
电饭煲：x7+x8+x9+x10 <=62
（2）计划外的每种产品的总销售量总是不大于计划内每种产品的总销售额：
yi<=xi
（3）每种产品计划内的销售量应该等于已签约的量：xi-yi=zi
然后应用LinGo软件编程（程序见附录六）得出结果
模型解释
为了便于比较，我们分别对五个问题列表进行解释如下，（单位：千）

问题（1） 表格 1
第i种产品的总销售量xi
12
24
14
17.316
27.807
19.876
15.019
27.437
13.687
5.856
计划外第i种产品的销售量yi
4
4
4
2.316
2.807
2.876
0.019
0.437
0.687
0.856
问题（2） 表格2
第i种产品的总销售量xi
11
27
12
15
29
20
15
26
13
8
计划外第i种产品的销售量yi
0
4
0
0
4
3
4
5
3
3
问题（3） 表格 3
第i种产品的总销售量xi
8
27
15
16
29
20
11
32
14
5
计划外第i种产品的销售量yi
0
4
4
4
5
3
0
5
3
10
问题（4） 表格 4
第i种产品的总销售量xi
8
23.65
13.05
14.95
24.75
16.8
15.3
25
13.7
5.163
计划外第i种产品的销售量yi
0
0
0
0
0
0
0
0
0
0
问题（5） 表格 5
第i种产品的总销售量xi
12.687
27.312
10
10
22.11
32.8
11.0
17
29
5
计划外第i种产品的销售量yi
4.687
7.312
0
0
2.11
17.8
0
0
19
0
模型应用与推广
本问题是一个优化问题，通过对问题的分析建立了多个优化模型，可以推广到其他优化模型中，但是要根据实际情况进行假设，要注意条件的改变和现实情况的不同，建立真正适用的模型。
参考文献
（1）王连堂《数学建模》 陕西师范大学出版
（2）文波 《数学建模及其基础知识详解》
（3）陈怀申《MATLAB及其在理工课程中的应用指南》西安电子科技大学出版社

附录一
表1 10种小家电的“已签约合同量”以及“意向签约量”表
序 号
1
2
3
4
5
6
7
8
9
10
家电名称
1.5升热水壶
1.8升热水壶
2升热水壶
0.9升豆浆机
1.1升豆浆机
1.3升豆浆机
2升电饭煲
2．5升电饭煲
3升电饭煲
3.5升电饭煲
已签约量（万个）
0.8
2.0
1.0
1.0
2.0
1.5
1.1
1.7
1.0
0.5
意向签约量（万个）
0.3
0.3
0.2
0.5
0.5
0.2
0.4
1.0
0.3
0.4
表2 10种小家电产量与成本关系表
产量（千件）
家电
1
3
5
7
10
15
20
25
30
35
40
1
9
24.3
39.6
53.55
72.9
102.6
129.6
153
175.5
195.3
208.8
2
9.8
26.46
43.12
58.31
79.38
111.7
141.2
166.6
191.1
208.3
227.4
3
10.2
27.54
44.88
60.69
82.62
116.3
146.9
173.4
198.9
221.3
236.6
4
23
66.55
104.6
140.1
188.6
265.7
335.8
391
441.6
491.1
515.2
5
26.5
75.52
120.6
161.4
217.3
298.1
386.9
450.5
508.8
565.8
593.6
6
27.2
77.52
123.6
165.6
223.1
314.2
397.1
462.4
522.2
580.7
609.3
7
28
79.8
126
166.6
224
315
380.8
434
462
509.6
560
8
30
85.5
135
178.5
240
337.5
408
465
495
546
600
9
31.7
90.3
142.6
188.6
253.6
356.6
431.1
491.3
523.1
576.9
634
10
32.5
92.
146.2
193.4
260
365.6
442
503.7
536.2
591.5
650
说明：表2中的数据为每种家电在相应产量下的总成本（单位：万元）。

表3 每种小家电的宣传费用随着销售量变化表
产量（千件）
家电
1
3
5
7
10
15
20
25
30
35
40
1
0.88
2.33
3.63
4.76
6.28
8.69
10.98
12.96
14.74
16.41
17.54
2
0.96
2.54
3.96
5.19
6.83
9.46
11.96
14.11
16.05
17.49
19.1
3
1
2.64
4.12
5.39
7.11
9.85
12.44
14.69
16.71
18.59
19.87
4
2.25
6.38
9.59
12.45
16.23
22.51
28.45
33.12
37.09
41.25
43.27
5
2.6
7.24
11.06
14.34
18.71
25.25
32.77
38.16
42.74
47.52
49.86
6
2.67
7.43
11.33
14.72
19.21
26.61
33.63
39.16
43.86
48.78
51.18
7
2.74
7.65
11.56
14.81
19.28
26.67
32.26
36.76
38.81
42.81
47.04
8
2.94
8.19
12.38
15.87
20.66
28.59
34.56
39.38
41.58
45.86
50.4
9
3.11
8.66
13.08
16.76
21.83
30.21
36.51
41.62
43.94
48.46
53.26
10
3.18
8.82
13.41
17.19
22.38
30.97
37.45
42.66
45.04
49.69
54.6
说明：表3中的数据为每种小家电在相应销售量下的宣传费用（单位：万元）。

表4 每种小家电的销售额随订购量变化表
订购量（千件）
家电
1
2
3
4
5
6
7
8
9
10
11
1
12.6
22.9
33.3
43.1
51.9
60.1
68
75.7
82.1
89.7
94.3
2
13.72
25.3
36.3
46.6
56.5
65.1
74.1
81.1
88.7
97.6
102.6
3
14.28
26.5
37.7
48.4
58.8
68.2
77.1
83.8
92.5
101.6
108
4
32.2
61.8
91.2
114.7
137
157.8
177.9
194.8
210.9
231.9
246.2
5
37.1
71.9
103.5
131.6
158
181.9
204.9
224.7
245.1
267.3
283.3
6
38.08
74.4
106.2
134.9
161.9
186.9
210.3
230.4
249.9
274.4
292.7
7
39.2
75.3
109.3
137.9
165.1
188.7
211.6
230.8
251.8
275.5
294.6
8
42
80.2
117.1
147.8
176.8
202.3
226.7
248.7
270.3
295.2
314.9
9
44.38
84.8
123.7
156.1
186.8
213.6
239.5
252.9
275.2
311.9
331.8
10
45.5
86.9
126
159.6
191.5
219
245.6
269.4
292.1
319.8
341.3
说明：表4中的数据为每种小家电在相应订购量下的销售额（单位：万元）。
表5 八个签约订单的签约量表
家电
订单
1
2
3
4
5
6
7
8
9
10
1
0.45
0.8
0.1

0.35
0.5
0.25
0.08
2

0.5
0.35

3
0.35
0.45
0.2

0.45
0.45
0.35
0.07
4

0.25
0.35
0.2

5

0.35

0.4
0.3
0.35
0.1
0.15
6

0.45
0.4
0.25

0.4
0.3
0.2
7

0.65
0.35

8

0.95
0.5

说明：表5中的数据为已签约的每个订单中对应小家电的数量（单位：万个）。
表6 五家意向签约客户的最大意向签约量表（万个）
家电
意向客户
1
2
3
4
5
6
7
8
9
10
1
0.15
0.1
0.1

0.1
0.25
0.2

2
0.15
0.2
0.1

0.1
0.4
0.1
0.2
3

0.3
0.15
0.1

4

0.2
0.35
0.1

5

0.2
0.35

0.2
表7 意向客户1可能签订订购合同量的概率表
家电
签约量（万个）
1
2
3
4
5
6
7
8
9
10
0
0.3
0.6
0.55

0.3
0.1

0.1
0.3
0.4
0.45

0.7
0.1
0.3

0.15
0.4

0.2

0.2

0.4
0.7

0.25

0.2

表8 意向客户2可能签订订购合同量的概率表
家电
签约量（万个）
1
2
3
4
5
6
7
8
9
10
0
0.2
0.3
0.5

0.4
0.1
0.2
0.1
0.1
0.3
0.2
0.5

0.6
0.1
0.8
0.3
0.15
0.5
0.3

0.3

0.4
0.2

0.2

0.2

0.2
0.3

0.2

0.4

0.1

表9 意向客户3可能签订订购合同量的概率表
家电
签约量（万个）
1
2
3
4
5
6
7
8
9
10
0

0.3
0.4
0.5

0.1

0.1
0.2
0.5

0.15

0.3
0.4

0.2

0.2

0.3

0.1

表10 意向客户4可能签订订购合同量的概率表
家电
签约量（万个）
1
2
3
4
5
6
7
8
9
10
0

0.4
0
0.7

0.1

0.3
0.1
0.3

0.15

0.2
0.3

0.2

0.1
0.3

0.3

0.1

0.35

0.2

表11 意向客户5可能签订订购合同量的概率表
家电
签约量（万个）
1
2
3
4
5
6
7
8
9
10
0

0.5
0.2

0.2
0.1

0.3
0.2

0.5
0.2

0.2
0.3

0.3
0.3

0.2

0.35

0.1

表12 计划外可能签订订购合同量的概率表
产品
合同量
1
2
3
4
5
6
7
8
9
10
0
0.1
0.1
0.2
0.1
0.1
0.4
0.1
0.05
0.2
0.2
0.1
0.4
0.3
0.1
0.3
0.3
0.3
0.1
0.1
0.5
0.2
0.2
0.2
0.4
0.2
0.3
0.4
0.2
0.3
0.25
0.2
0.5
0.3
0.2
0.1
0.4
0.2
0.1
0.1
0.4
0.3
0.1
0.1
0.4
0.1
0.1
0.1
0.1
0.1

0.1
0.2

0.5

0.1

表13 企业向营销部发放计划内销售产品的经费（万元）表
销售量（千件）
家电
1
3
5
7
10
15
20
25
30
35
40
1
1.872
4.68
6.396
7.514
8.736
11.18
14.144
16.692
18.252
20.332
21.736
2
2.0384
5.1168
6.9576
8.2108
9.4744
12.22
15.392
18.2
19.864
21.632
23.608
3
2.1216
5.2832
7.2384
8.5332
9.8696
12.688
16.068
18.928
20.696
23.036
24.596
4
4.784
12.818
16.848
19.656
22.516
29.016
36.712
42.692
45.916
51.064
53.56
5
5.512
14.549
19.448
22.62
26
32.552
42.224
49.192
52.936
58.812
61.724
6
5.6576
14.913
19.916
23.244
26.676
34.32
43.368
50.492
54.288
60.372
63.336
7
5.824
15.34
20.332
23.4
26.78
34.372
41.6
47.372
48.048
52.988
58.24
8
6.24
16.432
21.736
25.064
28.704
36.868
44.564
50.752
51.48
56.784
62.4
9
6.5936
17.368
22.984
26.468
30.316
38.948
47.06
53.664
54.392
60.008
65.936
10
6.76
17.68
23.556
27.144
31.096
39.936
48.256
55.016
55.744
61.516
67.6
表14 计划外销售部分营销部向企业缴纳利润表
销售量（千件）
家电
1
2
3
4
5
6
7
8
9
10
11
1
2.16
3.75
5.4
6.69
7.38
8.115
8.67
9.42
9.69
10.08
10.476
2
2.352
4.302
5.904
7.086
8.028
8.631
9.474
9.762
10.164
10.932
11.253
3
2.448
4.578
6.096
7.314
8.352
9.09
9.246
9.48
10.314
11.388
12.218
4
5.52
10.215
14.79
17.475
19.44
21.27
22.68
23.124
24.282
25.980
26.508
5
6.36
12.534
16.788
20.124
22.44
24.54
26.1
26.802
27.858
30
31.104
6
6.528
13.224
17.208
20.604
22.98
25.38
26.82
27.384
29.382
30.780
32.028
7
6.72
12.84
17.7
21
23.46
25.44
27
28.122
29.364
30.9
32.640
8
7.2
13.47
18.96
22.53
25.08
27.33
28.92
29.82
30.51
33.12
34.440
9
7.608
14.28
20.04
23.79
26.52
28.8
30.54
31.584
31.962
33.78
35.760
10
7.8
14.79
20.4
24.3
27.18
29.52
31.32
32.28
32.58
34.68
37.308

附录二
问题一的程序
Max=
(0.0003*(x1-y1)*(x1-y1)*(x1-y1)-0.2720*(x1-y1)*(x1-y1)+11.4429*(x1-y1)+1.3576
+0.0078*(x2-y2)*(x2-y2)*(x2-y2)-0.4139*(x2-y2)*(x2-y2)+12.8769*(x2-y2)+1.2185
+0.0165*(x3-y3)*(x3-y3)*(x3-y3)-0.5603*(x3-y3)*(x3-y3)+13.9414*(x3-y3)+0.7558
+0.0682*(x4-y4)*(x4-y4)*(x4-y4)-1.9654*(x4-y4)*(x4-y4)+36.0104*(x4-y4)-2.0803
+0.0584*(x5-y5)*(x5-y5)*(x5-y5)-1.9019*(x5-y5)*(x5-y5)+39.7564*(x5-y5)-0.5955
+0.0783*(x6-y6)*(x6-y6)*(x6-y6)-2.2395*(x6-y6)*(x6-y6)+41.9655*(x6-y6)-1.5942
+0.1045*(x7-y7)*(x7-y7)*(x7-y7)-2.7425*(x7-y7)*(x7-y7)+44.6484*(x7-y7)-3.1439
+0.0985*(x8-y8)*(x8-y8)*(x8-y8)-2.7237*(x8-y8)*(x8-y8)+46.9763*(x8-y8)-2.6833
+0.1870*(x9-y9)*(x9-y9)*(x9-y9)-4.3020*(x9-y9)*(x9-y9)+55.7898*(x9-y9)-9.0291
+0.1048*(x10-y10)*(x10-y10)*(x10-y10)-2.9051*(x10-y10)*(x10-y10)+50.6082*(x10-y10)-2.7303)-
(0.0002*(x1-y1)*(x1-y1)*(x1-y1)-0.0735*(x1-y1)*(x1-y1)+7.7904*(x1-y1)+1.8542+0.0005*(x2-y2)*(x2-y2)*(x2-y2)-0.0949*(x2-y2)*(x2-y2)+8.7093*(x2-y2)+1.4070+0.0002*(x3-y3)*(x3-y3)*(x3-y3)-0.0834*(x3-y3)*(x3-y3)+8.8314*(x3-y3)+2.0939+0.0004*(x4-y4)*(x4-y4)*(x4-y4)-0.2036*(x4-y4)*(x4-y4)+20.3812*(x4-y4)+5.7324-0.2010*(x5-y5)*(x5-y5)+22.8334*(x5-y5)+8.0202+0.0005*(x6-y6)*(x6-y6)*(x6-y6)-0.2436*(x6-y6)*(x6-y6)+24.1736*(x6-y6)+6.2363+0.0065*(x7-y7)*(x7-y7)*(x7-y7)-0.6369*(x7-y7)*(x7-y7)+29.0853*(x7-y7)-3.4501+0.007*(x8-y8)*(x8-y8)*(x8-y8)-0.6824*(x8-y8)*(x8-y8)+31.1628*(x8-y8)-3.6965+0.0073*(x9-y9)*(x9-y9)*(x9-y9)-0.7210*(x9-y9)*(x9-y9)+32.9269*(x9-y9)-3.9249+0.0076*(x10-y10)*(x10-y10)*(x10-y10)-0.7409*(x10-y10)*(x10-y10)+33.7974*(x10-y10)-4.2823)+
(0.007*y1*y1*y1-0.2024*y1*y1+2.3288*y1-0.0140+0.0101*y2*y2*y2-0.2572*y2*y2+2.6463*y2-0.0346+0.0189*y3*y3*y3-0.3986*y3*y3+3.2644*y3-0.4915+0.0293*y4*y4*y4-0.7339*y4*y4+7.0365*y4-0.8632+0.0367*y5*y5*y5-0.8843*y5*y5+8.2061*y5-0.8670+0.0354*y6*y6*y6-0.8631*y6*y6+8.1944*y6-0.5351+0.0387*y7*y7*y7-0.9269*y7*y7+8.5683*y7-0.9010+0.0435*y8*y8*y8-1.03