1.改动rho_design，变成a1*int1(x,y)+a2*int2(x,y)，没有问题，a1=0.25,a2=0.25
2.改动几何，变成100*100的，rho_design变成int0(x,y)+a1*int1(x,y)，a1没有变化，警告：无法最小化非线性
2.1.F_load变成100，没用；
2.2.a1初值变成0.1，a1最后结果是-8.544；初值变成-0.1，a1最后结果是1.792——SNOPT算法
2.3.LM算法不支持，错误：不支持的Levenberg-Marquardt求解器对象混合贡献，请禁用非最小二乘对象或使用SNOPT，并重试
2.4.Monte Carlo计算也会出错，弹窗：优化计算器错误
2.5.NM算法和坐标查找算法出错，弹窗：For input string：“reg param”
