kurs_ras4et2_1

clear all;
y=[15.82 16.83 15.89 16.18 17.17;
9.37 8.92 8.84 9.9 9.64;
10.98 10.27 10.96 10.62 11.24;
23.56 23.37 23.5 23 24.07];
yi=sum(y,2)/5;
S0i=var(y,0,2);
G0=max(S0i)/sum(S0i);
S0=sum(S0i)/4;
S0sr=S0/5;
b0=sum(yi)/4;
b1=(-yi(1)+yi(2)-yi(3)+yi(4))/4;
b2=(-yi(1)-yi(2)+yi(3)+yi(4))/4;
b12=(yi(1)-yi(2)-yi(3)+yi(4))/4;
Sb0=S0sr/4;
X1=-1:0.1:1;
X0=ones(length(X1));
X2=-1:0.1:1;
X12=[1;-1;-1;1];
[X1,X2]=meshgrid(X1,X2);
F=b0.*X0+b1.*X1+b2.*X2+b12.*X1.*X2;
S=Sb0*4;
Fn=F-2.*S.^0.5;
Fv=F+2.*S.^0.5;
meshc(X1,X2,F);
colormap(copper);
xlabel('x');
ylabel('y');
zlabel('z');
figure;
mesh(X1,X2,F);
colormap(winter);
xlabel('x');
ylabel('y');
zlabel('z');
hold on;
mesh(X1,X2,Fn);
mesh(X1,X2,Fv);