代码搜索结果

a=[1,2,3;2,5,2;3,1,5]; b=[14,18,20]; n=size(a,1); u=zeros(size(a)); l=zeros(size(a)); u(1,:)=a(1,:); for i=2:n l(i,1)=a(i,1)/u(1,1); end for r=2:n for i=r:1:n sum1=0; f

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))

%a=[0.012,0.01,0.167;1,0.8334,5.91;3200,1200,4.2]; %b=[0.6781,12.1,981]'; 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.003000,59.14;5.291,-6.130];%系数矩阵a1 %b=[59.17,46.7

x=1:8; y=[15.3,20.5,27.4,36.6,49.1,65.6,87.8,117.6]; y1=log(y); n=1;%表示用一次函数拟和 m=length(x); A=zeros(m,n+1); for i=1:2 A(:,i)=x.^(i-1)'; end B=A'*A; C=A'*y1'; re=B\C; a=exp(re(1)); b=r

function [x,i]=erfen(a,b,e) for i=1:1000 x=a+(b-a)/2; %a与b中间的值 if (fc(x)==0|(b-a)/2

a=0;b=1;m=10;T=0.5;N=50;af=1; f=inline('sin(pi*x)','x'); h=(b-a)/m; k=T/N; lmd=af^2*k/h^2; x=linspace(a,b,m+1); x=x(2:m); i=1:m-1; u=f(i.*h); for j=1:N t=j*k; u=trisys(-lmd*ones(m-2

function fe=fitfun() %连续函数的最佳逼近 %取基{1, x, ...} %fitfun.m %默认算例为课本：P60，例3.1 %原函数f(x)=x^(1/2)， 定义域 [1/4, 1] %结果：P(x) = 10/27 + 88/135x 平方误差=0.00010803 fs=input('\n输入原函数f(x)：[直接回车表

X=0:pi/20:pi; Y=sin(X); x1=0.5:0.2:1.9; y1=[0.4794,0.6442,0.7833,0.8912,0.9636,0.9975,0.9917,0.9463]; x=0.6:0.2:1.8; L=Langrange(x1,y1,x); s=fx(x1,y1,x); sm=sanci2(x1,y1,x) plot(X,Y,'r'); hol

close all; clear all; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%具体实例%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% a=0;b=2;c=0;d=1; n=6;m=5;TOL=1e-10; ITMAX=200; f=inline('x*exp(y)','x','y'); ga=inline('0','x','y

function [x1,n]=newton(x0) x1=x0-fc(x0)/df(x0); n=1; while(abs(x1-x0)>=1.0e-6)&(n