program5_3.m

此为matlab的代码。包括数值o45算法及o23算法

%MATLAB Program5_3
A=[2 1 0 0;0 2 0 0;0 0 -1 0;0 0 0 -1];
B=[0;1;1;1];
syms f1 f2 f3 f4;
F=[f1 f2 f3 f4];
P_A=poly(A-B*F)
P_1x=poly([-2;-2;-2;-1])
P_2x=poly([-2;-2;-2;-2])
a_x14=P_1x(5);
a_x13=P_1x(4);
a_x12=P_1x(3);
a_x11=P_1x(1);

a_x24=P_2x(5);
a_x23=P_2x(4);
a_x22=P_2x(3);
a_x21=P_2x(1);
4+f1-2*f2+4*f3+4*f4=8;
4-3*f2+2*f1=20;
-3-3*f4+3*f3+f1=18;
-2+f4+f3+f2=7;

4+f1-2*f2+4*f3+4*f4=16;
4-3*f2+2*f1=32;
-3-3*f4+3*f3+f1=24;
-2+f4+f3+f2=8;