Question

本人以及可以用共轭梯度法解不含参数的非线性方程，现在想在方程中加入迭代参数，这样可以写一个while语句

来求得满足一定精度的解。原题大体如下：
function[f,g]=y3(x)
%r=10;
f=exp((x(1)+x(2)+1)^2)+1000000*((exp(x(1)^2-x(2))+exp(x(2)-x(1)^2))/2-1);
g(1)=exp((x(1)+x(2)+1)^2)*2*(x(1)+x(2)+1)+1000000*x(1)*(exp(x(1)^2-x(2))-exp(x(2)-x(1)^2));
g(2)=exp((x(1)+x(2)+1)^2)*2*(x(1)+x(2)+1)+500000*(-exp(x(1)^2-x(2))+exp(x(2)-x(1)^2));
用共轭梯度即可求出数值解。改写为：
function[f,g]=y3(x)
%r=10;
f=exp((x(1)+x(2)+1)^2)+r*((exp(x(1)^2-x(2))+exp(x(2)-x(1)^2))/2-1);
g(1)=exp((x(1)+x(2)+1)^2)*2*(x(1)+x(2)+1)+r*x(1)*(exp(x(1)^2-x(2))-exp(x(2)-x(1)^2));
g(2)=exp((x(1)+x(2)+1)^2)*2*(x(1)+x(2)+1)+r*(-exp(x(1)^2-x(2))+exp(x(2)-x(1)^2));
现在想通过r的迭代来写一段代码，望高手不吝赐教。

Answer 1

加个 function【f,g】＝diedai（r） for a＝0:10:r； 再把第二个式子里的r换成a即可

追问

谢谢回答，还有疑惑：这个新的函数是写在新的m文件吗？怎么调用呢；如果不是，是放在y3.m文件中的哪里呢？

追答

你可以新建m脚本，之后的话你可以直接调用啊。。。。。调用整个函数时就就直接用【f，g】=diedai（r）。。。。。你个小气鬼，连分都不给

追问

不太会用，用其他方法的