代码搜索结果

%Taylor法求解常微分方程 function y=Taylor(a,b,N,af); h=(b-a)/N; x(1)=a; y2(1)=af; y4(1)=af; jqj(1)=af; for i=2:N y2(i)=y2(i-1)+h*((1-h/2)*(x(i-1)-y2(i-1))+1);%二阶Taylor法 y4(i)=y4(i-1)+h*((1-h/

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

A=[5,-4,1,0;-4,6,-4,1;1,-4,6,-4;0,1,-4,5]; b=[2,-1,-1,2]'; n=length(b); w=10; D=diag(diag(A)); CL=-triu(A,1); CLZ=CL'; L=((D-w*CL)*D.^(1/2))/sqrt(w*(2-w)); M=L*L'; C=inv(D)*CL; u=[2,3,4,5]';

a=[4,-1,1;-1,4.25,2.75;1,2.75,3.5]; b=[8,4,10]; %a=[1,2,3;2,5,2;3,1,5]; %b=[14,18,20]; n=size(a,1); l(1,1)=sqrt(a(1,1)); for i=2:n l(i,1)=a(i,1)/l(1,1); end for j=2:n sum1=0; fo

a=0;b=1;N=10; h=(b-a)/N; T(1,1)=(b-a)*(ft(a)+ft(b))/2; %while (i>1) for i=2:10 sum=0; for k=1:2^(i-2) sum=sum+ft(a+(2*k-1)*(b-a)/2^(i-1)); end T(1,i

function y=fc(x); y=3*x^2-exp(x); %y=x*tan(x)-0.5;

function X=trisys(A,D,C,B); N=length(B); for k=2:N mult=A(k-1)/D(k-1); D(k)=D(k)-mult*C(k-1); B(k)=B(k)-mult*B(k-1); end X(N)=B(N)/D(N); for k=N-1:-1:1 X(k)=(B(k)-C(k)*X(k+1))

close all; clear all; a=0;b=2;c=0;d=1; n=6;m=5;TOL=1e-10; ITMAX=100; f=inline('x*exp(y)','x','y'); ga=inline('0','x','y');gb=inline('2*exp(y)','x','y'); gc=inline('x','x','y');gd=inline('exp(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

function y=fg(x) y=1./(1+25*x.^2);