|
4 ~) o. E: x+ q5 N4 o1 x: ` K = DSolve[{-V/(x*f''[x])*Sqrt[1 + f'[x]^2] == 2*V, f[-L] == 0,2 p( W. j& T( m+ J( T) z f'[-L] == 0}, {f[x]}, x]; 6 f4 \0 W; o* YM = K /. {V -> 1, L -> 100};7 t5 g: c! w! ]0 c Manipulate[( D0 [; ]# U4 G! }" l) b& ` Plot[M[[2, 1, 2]], {x, -100, 0},% ?. v' C6 S$ E: w L9 k9 p Epilog -> {PointSize[0.01], Black,/ l8 z1 u7 X: f: U9 K4 A Point[{{a, M[[2, 1, 2]] /. x -> a}, {0, 6 `& V5 A# }& b, E7 n& V6 d' |M[[2, 1, 2]] - x*D[M[[2, 1, 2]], x] /. x -> a}}], Red, Dashed, 0 X+ _# D- |$ T- \Line[{{a, M[[2, 1, 2]] /. x -> a}, {0,; z! Y- p5 F2 B' f/ ]* _ W M[[2, 1, 2]] - x*D[M[[2, 1, 2]], x] /. x -> a}}]}], {a, -100,2 K* t: a# @: u0 d. X* {+ P 0}] 9 [. u( t9 `. }# U% W8 s. ^* QExport["1.gif",Z3 x3 ^7 v! W* @6 d Table[Plot[M[[2, 1, 2]], {x, -100, 0},. J- `5 B$ z8 {. j o( q7 B Epilog -> {PointSize[0.01], Black,{4 M. J* Y0 g1 ?" `9 b# q1 f5 q Point[{{a, M[[2, 1, 2]] /. x -> a}, {0, ; n7 Y) `- L: V8 _. I& i/ YM[[2, 1, 2]] - x*D[M[[2, 1, 2]], x] /. x -> a}}], Red, Dashed, $ Y& S/ t. Q5 \0 g. Y9 _8 @Line[{{a, M[[2, 1, 2]] /. x -> a}, {0,0 z7 G, A% u$ M: k7 t M[[2, 1, 2]] - x*D[M[[2, 1, 2]], x] /. x -> a}}]}], {a,& r. \. n& a2 z$ w, P: d Join[Range[-100, -1], {-0.0001}]}]] 3 k6 G9 z% X7 ?! j& J% S+ E# F! L2 p6 s) [7 f, K
|
|