|
|
1.1 root 1: .G1
2: frame ht 2.5 wid 2.5
3: coord x 0,1 y 0,1
4: label bot "Direction field is $y sup prime = x sup 2 / y$"
5: label left "$y = sqrt {(2x sup 3 +1)/3}$" right .3
6: ticks left in 0 at 0,1
7: ticks bot in 0 at 0,1
8: len = .04
9: for tx from .01 to .91 by .1 do {
10: for ty from .01 to .91 by .1 do {
11: deriv = tx*tx/ty
12: scale = len/sqrt(1+deriv*deriv)
13: line from tx,ty to tx+scale,ty+scale*deriv
14: }
15: }
16: draw solid
17: for tx = 0 to 1 by .05 do {
18: next at tx, sqrt((2*tx*tx*tx+1)/3)
19: }
20: .G2
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.