29 实 验 Mathcad (1)编 程 举 例
1 1 100 3 7,编 写 程 序 计 算 － 之 中 可 被 和 整 除 的 整 数 之 和
解 求 解 程 序 为,
div s 0←
s s k+← mod k 3,( ) 0= mod k 7,( ) 0=∨if
k 1 100..∈for
s
:=
div 2208=
d s 0← s1 0← s2 0←( )
s s k+← mod k 3,( ) 0=if
s1 s1 k+← mod k 7,( ) 0=if
s2 s2 k+← mod k 21,( ) 0=if
k 1 100..∈for
s s1+ s2?
:=
d 2208=
2 1990 12 3,12,我 国 年 人 口 为 亿 千 万 假 如 以 后 按 千 分 之 的 增 长 率 增 长 问 需 要 经
15,过 多 少 年 将 达 到 亿 届 时 全 国 人 口 数 是 多 少？ 编 写 程 序 求 解 。
, r 解 求 解 程 序 为 以 增 长 率 为 自 变 量
speed r( ) k 0←
p 1.23 109?←
p p 1 r+( )?←
k k 1+←
p 1.5 109?<while
"Years"
k
"Populatin"
p

:=
令 r 0.012= 可 以 得 到,
speed 0.012( ) "Years"17 "Populatin"1506516262.011= 17 15.07,经 过 年 人 口 将 达 到 亿
3 (Fibonacci) 编 写 程 序 输 出 菲 波 那 契 数 列 的 n第 项
,解 求 解 程 序 为
F n( ) g 1← n 1≤if
h 1←
g 1←
tmp h g+←
h g←
g tmp←
k 2 n..∈for
otherwise
g
:=
F 14( ) 610= 14第 项
k 1 15..:= fibk F k( ):= 15前 项
fibT
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
0 0 1 2 3 5 8 13 21 34 55 89 144 233 377 610 987
=
4 X,,,设 为 如 下 形 式 的 一 个 分 组 的 次 数 分 布 表 求 加 权 平 均 标 准 差 以 及 变 异 系 数 比 如
X,2,的 第 一 列 为 某 工 厂 各 等 级 工 资 标 准 值 第 列 为 对 应 各 级 别 工 资 的 人 数 求 该 厂 职
,,,.工 的 工 资 总 额 平 均 工 资 标 准 差 变 异 系 数 等
W X( ) sf 0← s 0← s2 0←( )
n rows X( ) 1?←
s s X 0( )k X 1( )k?+←
s2 s2 X 0( )k
2
X 1( )k?+←
sf sf X 1( )k+←
k 0 n..∈for
μ ssf←
σ s2sf μ 2?←
" "工 资 总 额
s
" "平 均 工 资
μ
" "工 资 标 准 差
σ
" "变 异 系 数
σ
μ
100?

:=
X
435.5
485.5
535.5
585.5
635.5
685.5
725.5
12
25
58
125
86
48
10

:=
W X( ) " "工 资 总 额216422 " "平 均 工 资594.57 " "工 资 标 准 差64.75 " "变 异 系 数10.89=
,:这 个 程 序 如 利 用 向 量 运 算 将 可 以 使 得 程 序 结 构 更 加 简 洁
W1 X( ) s X 0 X 1←
s2 X 0( )
2
X 1( )?←
sf X 1∑←
" "工 资 总 额
s
" "平 均 工 资
s
sf
" "工 资 标 准 差
s2
sf
s
sf

2
" "变 异 系 数
s2 sf? s2?
s 100?

:=
W1 X( ) " "工 资 总 额216422 " "平 均 工 资594.57 " "工 资 标 准 差4192.53 " "变 异 系 数10.89=