线性代数/投入产出分析 M 文件

# Octave 命令 用于 _线性代数_ 由 Jim Hefferon, # 主题: leontif.tex a=[(25448-5395)/25448 -2664/30346; -48/25448 (30346-9030)/30346]; b=[17589; 21243]; ans=a \ b; printf("第一个系统的答案是 s=%0.0f 和 a=%0.0f\n",ans(1),ans(2)); b=[17489; 21243]; ans=a \ b; printf("第二个系统的答案是 s=%0.0f 和 a=%0.0f\n",ans(1),ans(2)); # 问题 1 b=[17789; 21243]; ans=a \ b; printf("问题 (1a) 的答案是 s=%0.0f 和 a=%0.0f\n",ans(1),ans(2)); b=[17689; 21443]; ans=a \ b; printf("问题 (1b) 的答案是 s=%0.0f 和 a=%0.0f\n",ans(1),ans(2)); b=[17789; 21443]; ans=a \ b; printf("问题 (1c) 的答案是 s=%0.0f 和 a=%0.0f\n",ans(1),ans(2)); # 问题 2 printf("汽车对钢铁使用的比率是 %0.4f\n",2664/30346); a=[(25448-5395)/25448 -0.0500; -48/25448 (30346-9030)/30346]; b=[17589; 21243]; ans=a \ b; printf("2(a) 的答案是 s=%0.0f 和 a=%0.0f\n",ans(1),ans(2)); b=[17589; 21500]; ans=a \ b; printf("2(b) 的答案是 s=%0.0f 和 a=%0.0f\n",ans(1),ans(2)); # 问题 3 printf("其他人使用的钢铁价值是 %0.2f\n",18.69-(6.90+1.28)); printf("其他人使用的汽车价值是 %0.2f\n",14.27-(0+4.40)); a=[(18.69-6.90)/18.69 -1.28/14.27; -0/18.69 (14.27-4.40)/14.27]; b=[1.10*(18.69-(6.90+1.28)); 1.15*(14.27-(0+4.40))]; ans=a \ b; printf("3(a) 的答案是 s=%0.2f 和 a=%0.2f\n",ans(1),ans(2)); printf("1947 年钢铁对钢铁使用的比率是 %0.2f\n",(18.69-6.90)/18.69); printf("1947 年钢铁对汽车使用的比率是 %0.2f\n",1.28/14.27); printf("1947 年汽车对钢铁使用的比率是 %0.2f\n",0/18.69); printf("1947 年汽车对汽车使用的比率是 %0.2f\n",(14.27-4.40)/14.27); printf("1958 年钢铁对钢铁使用的比率是 %0.2f\n",(25448-5395)/25448); printf("1958 年钢铁对汽车使用的比率是 %0.2f\n",2664/30346); printf("1958 年汽车对钢铁使用的比率是 %0.2f\n",48/25448); printf("1958 年汽车对汽车使用的比率是 %0.2f\n",(30346-9030)/30346); b=[17.598/1.30; 21.243/1.30]; ans=a \ b; newans=1.30 * ans; printf("3(c) 的答案是 (以 1947 年十亿美元为单位) s=%0.2f 和 a=%0.2f\n 以及以 1958 年十亿美元为单位它分别是 s=%0.2f 和 a=%0.2f\n",ans(1),ans(2),newans(1),newans(2));