detamyfun_h.m

来自「同伦方程组求非线性方程组所有解」· M 代码 · 共 27 行

function y=detamyfun_H(x)
% find the deta of function myfun! 
%clear; 
%x=[-0.15535199070086;0.34341369822925;-0.75632689194207;  -1.18890170828280 ; -0.31367935040340 ; -10.93106207068521];
%x=[0;0.2];
[n0,n1]=size(x);h0=0.01;h1=norm(myfun_H(x),inf);
eps1=min(h0,h1);detaf=[];
format long
 for j=1:n0
     x(j)=x(j)-eps1;A11=myfun_H(x);
     x(j)=x(j)+2*eps1;A12=myfun_H(x);
     x(j)=x(j)-eps1;  
     for i=1:n0
       w1=A11(i,1);
       w2=A12(i,1);
       deta(i,j)=(w2-w1)/(2*eps1);
     end
end
y=deta;
%__________ 2004.3.6 __________compiled______ok!_____