1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
|
/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 17.07.0 ] */
/* [wxMaxima: title start ]
A manual testbench for plotdf
[wxMaxima: title end ] */
/* [wxMaxima: comment start ]
Maxima cannot check if the output of this testbench matches the expected results. But it at least can issue commands that test most of plotdf's features.
[wxMaxima: comment end ] */
/* [wxMaxima: input start ] */
plotdf(exp(-x)+y,[trajectory_at,2,-0.1])$
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
plotdf(x-y^2,[xfun,"sqrt(x);-sqrt(x)"],
[trajectory_at,-1,3], [direction,forward],
[y,-5,5], [x,-4,16])$
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
plotdf([v,-k*z/m], [z,v], [parameters,"m=2,k=2"],
[sliders,"m=1:5"], [trajectory_at,6,0])$
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
plotdf([y,-(k*x + c*y + b*x^3)/m],
[parameters,"k=-1,m=1.0,c=0,b=1"],
[sliders,"k=-2:2,m=-1:1"],[tstep,0.1])$
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
plotdf([w,-g*sin(a)/l - b*w/m/l], [a,w],
[parameters,"g=9.8,l=0.5,m=0.3,b=0.05"],
[trajectory_at,1.05,-9],[tstep,0.01],
[a,-10,2], [w,-14,14], [direction,forward],
[nsteps,300], [sliders,"m=0.1:1"], [versus_t,1])$
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
V: 900/((x+1)^2+y^2)^(1/2)-900/((x-1)^2+y^2)^(1/2)$
ploteq(V,[x,-2,2],[y,-2,2],[fieldlines,"blue"])$
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
plotdf([-y,x],[trajectory_at,5,0])$
/* [wxMaxima: input end ] */
/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$
|
