GoodLuckAC的博客

数值分析实验，站在巨人（百度）的肩膀上写了代码~~~~~。function [p,q]=chashang_newton(x,y)%%%%%%%%%%%%%%%%%%%%差商表的求法 n=length(x); p(:,1)=x; p(:,2)=y; for j=3:n+1 p(1:n+2-j,j)=diff(p(1:n+3-j,j-1)) ./(x(j-1:n)-x(1...

x=[0 10 20 30 40 80 90 95]';y=[68 67.1 66.4 65.6 64.6 61.8 61 60]';x1=ones(8,1);z=[x x1]\y;a=z(1);b=z(2);v=linspace(0,95,100);plot(x,y,'b-+',v,a*v+b,'k-');plot(x,y,'LineWidth',3,'Marker','dia...

实验内容计算圆周率π，要求误差小于10-8。 用复合梯形公式和复合辛普森求积公式计算。    ​复合梯度代码%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%function T_n=fht(a,b,n) h=(b-a)/n; for k=0:n x(k+1)=a+k*h; if...