生成一些数据测试你写的函数。 function [v]=myvander(x) N=length(x); for k=1:N
v(k,:)=x.^(k-1); end
三、思考题
编写程序,用如下迭代公式求
a
,a 的值分别为:3,17,113。迭代的终止条件为
5n 1n 10x x -+≤-,迭代初值0.1x 0=,迭代次数不超过100次。分别对迭代结果和
准确值进行比较,并统计迭代次数。
n
n 1n 222x 2a 2x x a x x 2a x +=
?+=?=+
clc clear
a=input('请输入a='); x0=1; for n=1:100
x1=x0/2+a/(2*x0); if (abs(x1-x0)<=1e-5) break end x0=x1; end
disp('sqrt(a)数值解为:') x1
disp('迭代次数:') n
disp('sqrt(a)的准确解为:') sqrt(a)
请输入a=3
sqrt(a)数值解为: x1 =
1.7321 迭代次数:5
sqrt(a)的准确解为: ans =
1.7321
实验六 MATLAB 的绘图