f.m

本文件提供了几个数值分析中的Matlab小程序。水平有限

function y=f(a,b,x0)
a=[1 0 0 1 0 0 1 0 1 0 0 0 0 1 1 0 0 0 0 1 1;
   0 1 0 1 0 0 1 0 0 0 1 1 0 0 0 0 1 1 0 0 1;
   1 0 1 0 0 0 1 0 0 1 0 0 0 1 0 1 0 1 0 0 1;
   0 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
   1 0 0 0 1 0 1 0 1 0 0 0 1 0 0 0 1 0 1 0 1;
   1 0 0 0 1 1 0 0 0 1 0 1 0 0 0 1 0 0 0 1 1;
   0 1 1 0 0 0 1 0 0 0 1 0 1 0 0 1 0 0 1 0 1;
   1 0 0 1 0 0 0 1 0 0 1 1 0 0 0 1 0 0 1 0 1;
   1 0 1 0 0 0 0 1 1 0 0 1 0 0 0 0 1 0 0 1 1;
   0 1 0 0 1 0 1 0 0 1 0 1 0 0 1 0 0 0 1 0 1;
   1 0 1 0 0 1 0 0 0 0 1 0 1 0 1 0 0 0 1 0 1;
   0 1 1 0 0 0 0 1 1 0 0 1 0 0 1 0 0 1 0 0 1;
   0 0 0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0;
   1 0 0 0 1 0 0 1 0 0 1 0 0 1 1 0 0 1 0 0 1;
   0 1 0 1 0 0 0 1 0 1 0 0 1 0 1 0 0 0 0 1 1;
   0 1 0 1 0 1 0 0 1 0 0 0 0 1 0 1 0 0 1 0 1;
   1 0 0 1 0 1 0 0 0 1 0 0 1 0 0 0 1 1 0 0 1;
   0 1 0 0 1 0 0 1 1 0 0 0 1 0 0 1 0 1 0 0 1;
   0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0;
   0 1 0 0 1 1 0 0 0 0 1 0 0 1 0 0 1 0 0 1 1;
   0 1 1 0 0 1 0 0 0 0 1 0 1 0 0 1 0 0 0 1 1]
b=[859.1 597.5 937.7 0 579.3 510.4 640.1 616.8 832.3 636.1 801.8 979.7 0 980.5 909.5 862.2 799.7 799.2 0 932.6 803.1]
x0=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]