- 博客(3)
- 收藏
- 关注
原创 希望人没事3
jia ke bifunction 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') returnendx=zeros(n,1);k=0;while k<N for i = 1:n s = b(i); for j=1:
2021-06-07 01:03:54
94
原创 希望人没事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 tiaoX = [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)); endendfor i=1:n.
2021-06-07 00:34:51
54
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人