代码搜索结果

%Newton插值法 function y=Newton(x1,y1,x); m=length(x); n=length(x1); for i=2:n for j=n:-1:i y1(j)=(y1(j)-y1(j-1))/(x1(j)-x1(j-i+1)); end end for j=1:m p(j)=y1(1); for k

function y=df(x); y=diff('fc',1);

x=-1:0.01:1; y1=-50./(1+25.*x.^2).^2.*x; n=length(x); h(1)=x(2)-x(1); for i=2:n-1 h(i)=x(i+1)-x(i); lm(i)=h(i)/(h(i-1)+h(i)); mu(i)=1-lm(i); c1(i)=3*(lm(i)*(fg(x(i))-fg(x(i-1))

x1=[0,1,3]; y1=[1,3,2]; x=0:0.5:3; Newton(x1,y1,x);

%x=[-1.0,-0.5,0.0,0.5,1.0,1.5,2.0]; %y=[-4.447,-0.425,0.551,0.048,-0.447,0.549,4.552]; %x=-3:3; %y=[4,2,3,0,-1,-2,-5]; x=[0.1,0.2,0.3,0.4,0.5]; y=[5.1234,5.3057,5.5687,5.9378,6.4370]; n=2; m=le

%Langrange插值法 function y=Langrange(x1,y1,x); m=length(x); n=length(x1); for i=1:m y(i)=0; for j=1:n l=1; for k=1:n if (k~=j) l=(x(i)-x1(k)

function f=fun1(x); f=x.^2./(1+sin(x)+x.^2); end

x0=1.5; TOL=10^-2; N=10; i=1; while(i

%Euler法求解常微分方程 function y=Euler(a,b,N,af); h=(b-a)/N; x(1)=a; y(1)=af; yg(1)=af; yh(1)=af; jqj(1)=af; for i=2:N+1 y(i)=y(i-1)+h*f(x(i-1),y(i-1));%Euler法 yh(i)=yh(i-1)+(h/4)*(f(x(i-1)

%a=[0.003000,59.14;5.291,-6.130];%系数矩阵a1 %b=[59.17,46.78]'; %a=[1,-1,2,-1;2,-2,3,-3;1,1,1,0;1,-1,4,3];%系数矩阵a %b=[-8,-20,-2,4]'; a=[0.012,0.01,0.167;1,0.8334,5.91;3200,1200,4.2]; b=[0.6781,12.1,98