-0.0000 -2.2361 请输入正数a=8 请输入正数b=3 x = 1.7016 r = 1.7016 -1.6180 r = 1.7016 -4.7016 s = 0.0 -6.4031 请输入正数a=10 请输入正数b=0.1 x = 3.1127 r = 3.1127 -4.7016 r = 3.1127 -3.2127 s = -0.0000 -6.3254
4. 已知
?f1?1?f?0?2??f3?1??fn?fn?1?2fn?2?fn?3求f1~f100中：
(1) 最大值、最小值、各数之和。 (2) 正数、零、负数的个数。 解：M—文件
n?1n?2n?3n?3
以下是运算结果： max(f)=437763282635 min(f)=-899412113528 sum(f)=-742745601951 c1=49 c2=2 c3=49
5. 若两个连续自然数的乘积减1是素数，则称这两个边疆自然数是亲密数对，该素数是亲密素数。例如，2×3-1=5，由于5是素数，所以2和3是亲密数，5是亲密素数。求[2,50]区间内：
(1) 亲密数对的对数。
(2) 与上述亲密数对对应的所有亲密素数之和。 解：
M文件：
运算结果为： j =
29 s =
23615
实验五 函数文件
1. 定义一个函数文件，求给定复数的指数、对数、正弦和余弦，并在命令文件中调用该函数文件。
解：M文件如下： 函数fushu.M文件： function [e,l,s,c] = fushu(z) %fushu 复数的指数，对数，正弦，余弦的计算 %e 复数的指数函数值 %l 复数的对数函数值 %s 复数的正弦函数值 %c 复数的余弦函数值 e=exp(z); l=log(z); s=sin(z); c=cos(z); 命令文件M： z=input('请输入一个复数z='); [a,b,c,d]=fushu(z) 运算结果如下： z=input('请输入一个复数z='); [a,b,c,d]=fushu(z) 请输入一个复数z=1+i a = 1.4687 + 2.2874i b = 0.3466 + 0.7854i c =
1.2985 + 0.6350i d = 0.8337 - 0.9889i
2. 一物理系统可用下列方程组来表示：
?m1cos??msin??1?0??0?m10m20?sin?cos??sin??cos?0??a1??0??a??mg?0???2???1? 0??N1??0??????1??N2??m2g?从键盘输入m1、m2和θ的值，求a1、a2、N1和N2的值。其中g取9.8，输入θ时以角度为单位。
要求：定义一个求解线性方程组AX=B的函数文件，然后在命令文件中调用该函数文件。
解： M文件
函数fc.M文件: function X= fc(A,B) ü fc是求解线性方程的函数 %A A是未知矩阵的系数矩阵 X=A\\B； 命令M文件： clc; m1=input('输入m1='); m2=input('输入m2='); theta=input('输入theta='); x=theta*pi/180; g=9.8; A=[m1*cos(x) -m1 -sin(x) 0 m1*sin(x) 0 cos(x) 0 0 m2 -sin(x) 0 0 0 -cos(x) 1]; B=[0;m1*g;0;m2*g]; X=fc(A,B) 运算结果： 输入m1=1 输入m2=1 输入theta=30 X = 7.8400