· function varargout=Yxian(varargin) fx=inline('exp(x)*(sin(x)-2*cos(x))'); px=inline('0'); qx=inline('1'); a=0;b=pi;a0=0;b0=0;h=pi/10; [u,x]=YY(fx,px,qx,a,b,a0,b0,h); uu=u;n=length(u); for i=1:n...

· function varargout=liu(varargin) a=0;b=2;c=0;d=1;h1=1/16;h2=1/16; f=inline('(pi^2-1)*exp(x)*sin(pi*y)','x','y'); g1x=inline('0'); g2x=inline('0'); g1y=inline('sin(pi*y)'); g2y=inline('exp(2)*sin...

· function varargout=liu(varargin) T=1;a=1;h=1/10;k=1/400; f=inline('0','x','t'); fx=inline('exp(x)'); ft1=inline('exp(t)'); ft2=inline('exp(1+t)'); [X,Y,Z]=chfenmethed(f,fx,ft1,ft2,a,T,h,k); mes...

· function varargout=liu(varargin) a=1;T=1;a=1;b=0.5;h=1/20;k=1/40; f=inline('0','x','t'); fx1=inline('exp(x)'); fx2=inline('exp(x)'); ft1=inline('exp(t)'); ft2=inline('exp(1+t)'); [X,Y,Z]=chfenm...

· function varargout=saxplaxliu(varargin) A=[3 3 5;3 5 9;5 9 17]; b=[10 16 30]'; x=f(A,b) function x=f(A,b) n=length(b); L=zeros(n,n); for i=1:n s=0; if i>1 for k=1:i-1 s=s+L(i,k)^2;...

0 Bytes
2014-01-12
0
· function varargout=Yxian(varargin) xj=[0 1 2]; yj=[5 4 7]; mj=[1 2 4]; x=1.5; y=f(xj,yj,mj,x) function y=f(xj,yj,mj,x) n=length(xj); y=0; for i=1:n s=1; for j=1:n if i~=j s=s*(x-xj(j...

· function varargout=Yxian(varargin) qx=inline('1'); fx=inline('exp(x)*(sin(x)-2*cos(x))'); a=0;b=pi;a0=0;b0=0;h=pi/10; [u,x]=YY(fx,qx,a,b,a0,b0,h); uu=u;n=length(u); for i=1:n uu(i)=exp(x(i))*si...

· function varargout=liu(varargin) a=1;T=1;b=0.5;h=1/10;k=1/20; f=inline('0','x','t'); fx1=inline('exp(x)'); fx2=inline('exp(x)'); ft1=inline('exp(t)'); ft2=inline('exp(1+t)'); [X,Y,Z]=chfenmethe...

· function varargout=saxplax(varargin) clear,clc x0=0;xn=1;y0=1;z0=1;h=0.01; [y,x]=lgkt4j(x0,xn,y0,z0,h); plot(x,y) n1=length(x); fprintf(' x y\n'); for i=1:n1 fprintf('%4.2f %4.2f\n',x(i),y(i)...

0 Bytes
2014-01-12
0
· % MATLAB M-file fractal.m iterations=25000; Mat1=[0 0;0 0.16]; Mat2=[0.85 0.04;-0.04 0.85]; Mat3=[0.2 -0.26;0.23 0.22]; Mat4=[-0.15 0.28;0.26 0.24]; Vector1=[0;0]; Vector2=[0;1.6]; Vector3=[0;...

· function varargout=liu(varargin) %测试问题 来自文献 孙志忠 偏微分数值解法 C=1;a1=0;a2=1;b1=0;b2=1;h1=1/20;h2=1/800; fx=inline('exp(x)'); gy1=inline('exp(y)'); gy2=inline('exp(1+y)'); f0=inline('...

· function varargout=Yxian(varargin) %算法及测试问题 来自文献 孙志忠 偏微分数值解法 qx=inline('1'); fx=inline('exp(x)*(sin(x)-2*cos(x))'); a=0;b=pi;a0=-1;b0=-exp(pi);r1=0;r2=0;h=pi/10; [u,x]=YY(fx...

· function varargout=Yxian(varargin) %算法及测试问题 来自文献 孙志忠 偏微分数值解法 qx=inline('1'); fx=inline('exp(x)*(sin(x)-2*cos(x))'); a=0;b=pi;a0=-1;b0=-exp(pi);r1=0;r2=0;h=pi/160; [u,x]=YY(f...

· function varargout=liu(varargin) %测试问题 来自文献 孙志忠 偏微分数值解法 C=1;a1=0;a2=1;b1=0;b2=1;h1=1/20;h2=1/800; fx=inline('exp(x)'); gy1=inline('exp(y)'); gy2=inline('exp(1+y)'); [X,Y,U]=Hea...

· function varargout=SORliu(varargin) clc A=ones(4,4); for i=1:4 A(i,i)=-4; end b=[1;1;1;1]; x0=[0;0;0;0];w=1.3; x=SOR(A,b,x0,w); for i=1:4 fprintf('%4.2f ',x(i)); end fprintf('\n'); fu...

· function varargout=Yxian(varargin) %算法及测试问题 来自文献 孙志忠 偏微分数值解法 1.1节 qx=inline('1'); fx=inline('exp(x)*(sin(x)-2*cos(x))'); a=0;b=pi;a0=0;b0=0;h=pi/10; [u,x]=YY(fx,qx,a,b,a0,b0...

0 Bytes
2014-01-12
0
· function varargout=clliu(varargin) A=[1 2;2 5]; B=[5 4;4 8]; y=Clo1(A,B) function y=Clo1(A,B) [a1,a2]=size(A); [b1,b2]=size(B); y=zeros(a1*b1,a2*b2); for i=1:a1 for j=1:a2 y((i-1)*b1+1:i*b...

· function varargout=newtonliu(varargin) clear,clc %测试 取自 李庆扬等 数值分析 x=[0.40 0.55 0.65 0.80 0.90 1.05]; fx=[0.41075 0.57815 0.69675 0.88811 1.02652 1.25382]; newtonchzh(x,fx); function...

0 Bytes
2014-01-12
0
· function varargout=saxplax1(varargin) clear,clc xx=[1 2 3]; fchaf(xx); function fchaf(xx) %线性有理插值法 epsliu=0.00001; fprintf(' k x(k) f(x(k))\n '); for i=1:3 if f(xx(i))<0 fprintf('%2...