代码搜索结果

%x=-1:0.01:1; %y1=5000./(1+25.*x.^2).^3.*x.^2-50./(1+25.*x.^2).^2; function sm=sanci2(x,y,x1) y1=0; y11=0; 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

%波动方程有限差分法 f=inline('sin(pi*x)','x'); g=inline('0'); l=1;T=0.5;m=10;N=10;af=1; xx=[];tt=[];uu=[]; h=l/m; k=T/N; lmd=k*af/h; u(1,2:N+1)=0; u(m+1,2:N+1)=0; u(1,1)=f(0); u(m+1,1)=f(l); for i=

function y=fh(x) y=sin(x);

function y=jacobi(a,b,x0) D=diag(diag(a)); U=-triu(a,1); L=-tril(a,-1); B=D\(L+U); f=D\b; y=B*x0+f;n=1; while norm(y-x0)>=1.0e-6 x0=y; y=B*x0+f;n=n+1; end n

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

function s=trapr1(f,a,b,n) h=(b-a)/n; s=0; for k=1:n-1 x=a+h*k; s=s+feval(f,x); end s=h*(feval(f,a)+feval(f,b))/2+h*s;

%四阶RK法求解常微分方程 function y=RK(a,b,N,af); h=(b-a)/N; x(1)=a; y(1)=af; jqj(1)=af; for i=2:N+1 K1=f(x(i-1),y(i-1)); K2=f(x(i-1)+h/2,y(i-1)+h*K1/2); K3=f(x(i-1)+h/2,y(i-1)+h*K2/2);

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

function simp=sp(a,b) h=(b-a)/500; f0=f3(a)+f3(b); f1=0; f2=0; for i=1:499 x=a+i*h; if(rem(i,2)==0) f2=f2+f3(x); else f1=f1+f3(x); end end simp=h*(f0+2*f

a=0; b=1; TOL=10^(-4); N=200; %自适应求积法 % f=inline('(100*sin(10/x))/x^2','x'); f=inline('4./(1+x^2)','x'); APP=0; i=1; TOL(i)=10*TOL; a(i)=a; h(i)=(b-a)/2; FA(i)=feval(f,a); FC(i)=feval(f,a