《坐标轮换法c程序汇总.docx》由会员分享，可在线阅读，更多相关《坐标轮换法c程序汇总.docx（3页珍藏版）》请在三一文库上搜索。
1、#includeusing namespace std;#includedouble function(double x2)double z;/ z=10*(x0+x1-5)*(x0+x1-5)+(x0-x1)*(x0-x1);z=x0*x0*x0*x0-2*x0*x0*x1+x0*x0+2*x1*x1-2*x0*x1+4.5*x0-4*x1+4;return z;double pointju(double x1口,double x2,int n)double z;double sum=0;for(int i=0;in;i+)sum=(x1i-x2i)*(x1i-x2i)+sum;z=sqrt
2、(sum);return z;double * zuobiaolunhuan(int n,double x02,double h,double s口,double e)double a2;double b2;double *p;double result2;int i, k=1;double x12,x22,x32,f1,f2,f3;for(i=0;in;i+)x1i=x0i;x2i=x0i+k*h*si;f1=function(x1);f2=function(x2);if(f1f2) do k=k-1; for(i=0;in;i+) x3i=x0i+(k-1)*h*si; f3=functi
3、on(x3);x1i=x0i+h*(k+1)*si; f1=function(x1); while(f3f1); for(i=0;in;i+) ai=x3i; bi=x3i+2*h*si; else /仙+pY? do pk=k+1; for(i=0;in;i+) x3i=x0i+k*h*si; f3=function(x3); x2i=x0i+h*(k-1)*si;f2=function(x2); while(f3f2); for(i=0;in;i+) ai=x2i-h*si; bi=x3i; /coutx0 ?c? 依?Ta0,b0Tendl;coutx1 ? ?c? 依?Ta1,b1Te
4、ndl; double result2; do for(i=0;in;i+) x1i=ai+0.382*(bi-ai); x2i=ai+0.618*(bi-ai); f1=function(x1); f2=function(x2); if(f1f2) for(i=0;in;i+) bi=x2i; else for(i=0;in;i+) ai=x1i; cout？D?a0Tb0Te); for(i=0;in;i+) resulti=0.5*(ai+bi);x1i=resulti;coutx10=x10endlx11=x11endl;p=x1;return p;void main()double
5、x12,x22;int n=2;double x0=0,0;float h=1;double s12=1,0;double s22=0,1;double e=0.0001;double *r;/ double *r=shengli(n,x0,h,s1,e);/ coutr0endlr1endl;do(r=zuobiaolunhuan(n,x0,h,s1,e);for(int i=0;in;i+)(x1i=ri;x0i=x1i;zuobiaolunhuan(n,x0,h,s2,e);for(i=0;in;i+)(x2i=ri;for(i=0;ie);cout?o ey p?u?x?o ?a?a命endlx0=x20endlx1=x2 1endl; cout? bey p?x?D?(i?a塔endlf(x)= function(x2)endl;