文档介绍：附录一 图
附录二 程序代码
代码（1）
clc,clear
u=1;
v=;
w=;
a=10000;
xx=u*a/((-log())^)
yy=v*a/((-log())^)
zz=w*a/((-log())^)
x=linspace(0,60000);
y1=1-exp(-(1/xx*x).^2);
y2=1-exp(-(1/yy*x).^2);
y3=1-exp(-(1/zz*x).^2);
plot(x,y1,'b-',x,y2,'b--',x,y3,'r')
legend('y1企业单位','y2机关单位','y3事业单位 ')
title('企业、机关和事业单位满意度函数')
代码（2）
clc,clear
a=[0. 0. 0. 0. 0. 0. 0. 1. 1. 1. 0. 0. 0. 1. 1. 1. 1. 1. 2. 2. 0. 0. 1. 1. 1. 1. 1. 2. 3. 3.];
s=10000*a(1,:);
d=10000*a(2,:);
f=10000*a(3,:);
a=10000*a;
g=a';
bar(g,1)
legend('企业单位平均养老金','事业单位平均养老金','机关单位平均养老金')
title('企业、机关与事业单位平均养老金比较图')
colormap([1 0 0;0 1 0;0 0 1])
代码（3）
clc,clear
x=[1 2 3 4 5 6 7 8 9 10];
y=[0. 0. 0. 0. 0. 0. 0. 1. 1. 1.];
y=10000*y;
z=[0. 0. 0. 1. 1. 1. 1. 1. 2. 2.];
z=10000*z;
w=[0. 0. 1. 1. 1. 1. 1. 2. 3. 3.];
w=10000*w;
plot(x,y,'bo');
hold on
plot(x,z,'b*')
hold on
plot(x,w,'ro')
hold on
代码（4）
clc,clear
x=[1 2 3 4 5 6 7 8 9 10];
y=[ 0. 0. 0. 0. 0. 0. 0. 1. 1. 1.];
y=10000*y;
z=[ 9550 10227 11547 13188 15169 17730 20973 24186 26909 30587];
plot(x,y,'ro');
hold on
plot(x,z,'b*')
hold on
代码（5）
x=[1 2 3 4 5 6 7 8 9 10];
y=[12028 16356 18772 20315 21603 22742 24489 25963 27881 29208];
z=[22439 23919 25522 27301 28445 29239 30605 31588 33834 34671];
w=[20836 22203 23684 25329 26384 27114 28375 29281