- 博客(3)
- 收藏
- 关注
原创 希望人没事3
jia ke bi function x=jacobi(A,b,n,x0,eps,N) D = diag(diag(A)); L = -tril(A,-1); U = -triu(A,1); B = D\(L+U); if(max(abs(eig(B))>1)) disp('i dont know') return end x=zeros(n,1); k=0; while k<N for i = 1:n s = b(i); for j=1:
2021-06-07 01:03:54 88
原创 希望人没事2
two fen 法 function [outputArg1] = division(fun,left,right) mid=(left+right)/2; while (abs(left-right)>10^(-3)) if fun(mid)*fun(right)<=0 left=mid; else right=mid; end mid=(left+right)/2; end outputArg1 = mid;
2021-06-07 00:51:17 111
原创 希望人没事
%san ci yang tiao X = [0 1 3 4]; Y = [0 1 81 196]; x = 0:0.01:4; n=length(X); A=zeros(n,n); A(:,1)=Y'; D=zeros(n-2,n-2); g=zeros(n-2,1); s=zeros(size(x)); for j=2:n for i=j:n A(i,j)=(A(i,j-1)- A(i-1,j-1))/(X(i)-X(i-j+1)); end end for i=1:n.
2021-06-07 00:34:51 53
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人