内容简介作者简介 - 科学与工程计算系

2 2 6 2 4 4 2 6
x~ + y - 1/6 x~ - 1/2 y x~ - 1/2 y x~ - 1/6 y
5.7.3.2 运行 MAPLE 程序
【例 5.7.3.2-1】目标:设计求取一般隐函数 f ( x,
y)
= 0 的导数 y′ (x)
解析解的程序,并要求:
该程序能象 MAPLE 原有函数一样可被永久调用。
(1)
[DYDZZY.src]
DYDZZY:=proc(f)
# DYDZZY(f) is used to get the derivate of
# an implicit function
local Eq,deq,imderiv;
Eq:='Eq';
Eq:=f;
deq:=diff(Eq,x);
readlib(isolate);
imderiv:=isolate(deq,diff(y(x),x));
end;
(2)
procread('DYDZZY.src')
ans =
DYDZZY := proc (f) local Eq, deq, imderiv; Eq := 'Eq'; Eq := f; deq :=
diff(Eq,x); readlib(isolate); imderiv := isolate(deq,diff(y(x),x)) end
(3)
s1=maple('DYDZZY(x=log(x+y(x)));')
s2=maple('DYDZZY(x^2*y(x)-exp(2*x)=sin(y(x)))')
s3=maple('DYDZZY','cos(x+sin(y(x)))=sin(y(x))')
s1 =
diff(y(x),x) = x+y(x)-1
s2 =
diff(y(x),x) = (-2*x*y(x)+2*exp(2*x))/(x^2-cos(y(x)))
s3 =
diff(y(x),x) =
sin(x+sin(y(x)))/(-sin(x+sin(y(x)))*cos(y(x))-cos(y(x)))
(4)
clear maplemex
procread('DYDZZY.src');
maple('save(`DYDZZY.m`)');
(5)
maple('read','`DYDZZY.m`');
ss2=maple('DYDZZY(x^2*y(x)-exp(2*x)=sin(y(x)))')
ss2 =
diff(y(x),x) = (-2*x*y(x)+2*exp(2*x))/(x^2-cos(y(x)))
5.7.3.3 数值、符号计算集成 M 文件的编写
5.8 可视化数学分析界面
17