代码搜索结果

%平方法计算相似矩阵的传递闭包 r1=[1,0.1,0.2;0.1,1,0.3;0.2,0.3,1]; n=size(r1,1); I=1; while (I

%Euler法求解常微分方程 function y=Euler1(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

%y''=p(x)y'+q(x)y+r(x),a

tic %其中a为对称正定矩阵 a=[5,-4,1,0;-4,6,-4,1;1,-4,6,-4;0,1,-4,5]; b=[2,-1,-1,2]; %a=[4,-1,1;-1,4.26,2.75;1,2.75,3.5]; %b=[8,4,10]; n=size(a,1); d(1)=a(1,1); for i=2:n for j=1:i-1 sum1=0

clear all %a1=[1,-1,2,-1;2,-2,3,-3;1,1,1,0;1,-1,4,3];%系数矩阵a1 %b=[-8,-20,-2,4]'; %a1=[0.012,0.01,0.167;1,0.8334,5.91;3200,1200,4.2]; %b=[0.6781,12.1,981]'; a1=[2,-2,5;2,3,1;-1,4,1]; b=[6,13,3]';

%f(x)=sin(x) %x=0.5:0.2:1.9; %y=[0.4794,0.6442,0.7833,0.8912,0.9636,0.9975,0.9917,0.9463]; %x=-1:0.01:1; %y=1./(1+25*x.^2); %n=length(x); %x1=-0.9:0.1:0.9; function s=fx(x,y,x1) n=length(x);

a(1)=1; t=2; h(1)=0.02; y(1)=f2(a); e=10^(-4); while(k>0) a(k+1)=a(k)+h(k); y(k+1)=f2(a(k+1)); if(y(k+1)

close all; clear all; 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+1); u(m)=0; for i=1:m-1 u(i)=f(i*h); end

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

a=[2,-1,0,0;-1,2,-1,0;0,-1,2,-1;0,0,-1,2]; b=[1,0,0,1]; n=size(a,1); l(1,1)=a(1,1); for i=1:n-1 l(i+1,i)=a(i+1,i); u(i,i+1)=a(i,i+1)/l(i,i); l(i+1,i+1)=a(i+1,i+1)-l(i+1,i)*u(i,i+1);