数学实验程序
数学实验考试小程序 牛顿迭代法：
function x=polyzero(P,x0);dP=polyder(P);x=x0-polyval(P,x0)/polyval(dP,x0);while abs(polyval(P,x))>1e-6;x=x-polyval(P,x)/polyval(dP,x);end end 同一个图上画多条曲线：hold on x=-2:0.01:2 for a=1:3;y=x.^a;plot(x,y,"r")grid on end 选出一个数组内的最大数，并返回其位置：function [s,t]=mymax(A);[m,n]=size(A);s=A(1,1);t=[1,1];for i=1:m for j=1:n if A(i,j)>s s=A(i,j);t=[i,j];end end end end 画五角星：c=pi/2:4*pi/5:4.5*pi;x=cos(c);y=sin(c);plot(x,y,"m")axis equal grid on
判断一个数的尾数：function a=weishu(x,n)a=(fix(x/10.^(n-1))*10.^(n-1)-fix(x/10.^n)*10.^n)/10.^(n-1);end 判断一个数的位数：function y=weishu1(x);y=length(num2str(x));end 求log1+log2+……+logn 的算法：n=0;s=0;while s<=100 n=n+1;s=s+log(1+n);end m=n 画球：function drawball(x0,y0,z0,r,b)ca=linspace(0,2*pi,31)+b;fa=linspace(0,pi,16);[c,f]=meshgrid(ca,fa);x=r*sin(f).*cos(c)+x0;y=r*sin(f).*sin(c)+y0;z=r*cos(f)+z0;mesh(x,y,z);end 画 单位 圆：t=linspace(0,2*pi,100)x=cos(t);y=sin(t);plot(x,y)马鞍面（抛物双曲面）的画法：t=linspace(0,2*pi,100);x=sin(t);y=cos(t);z=cos(2*t);plot3(x,y,z)grid on axis equal
