摘要: clcclearx0=input('x0=');y0=input('y0=');x=input('x=');n=length(x0);for k=1:n F(k,1)=y0(k);endfor i=2:n for j=2:i F(i,j)=(F(i,j-1)-F(i-1,j... 阅读全文
posted @ 2016-01-11 19:34 薄樱 阅读(86) 评论(0) 推荐(0) 编辑
摘要: a=input('a=');b=input('b=');x0=input('x0');y=input('y');n=length(x0);h=(b-a)/n;x(1)=y(1)+y(n+1);x(2)=0;x(3)=0;for i=2:n xx=mod(i,2); if (xx==1) ... 阅读全文
posted @ 2016-01-11 19:33 薄樱 阅读(70) 评论(0) 推荐(0) 编辑
摘要: clcclearx0=input('x0=');y0=input('y0=');FPO=input('FPO=');FPN=input('FPN=');x=input('x=');n=length(x0);for i=1:n a(i)=y0(i);endfor i=1:n-1 h(i)=... 阅读全文
posted @ 2016-01-11 19:32 薄樱 阅读(337) 评论(0) 推荐(0) 编辑
摘要: function Romberg(a,b,n)h=b-a;R(1,1)=h*(f(a)+f(b))/2;for i=2:n m=0; for k=1:2^(i-2) m=m+f(a+(k-0.5)*h); end R(2,1)=1/2*(R(1,1)+h*m);... 阅读全文
posted @ 2016-01-11 19:31 薄樱 阅读(802) 评论(0) 推荐(0) 编辑