Can't solve Coulomb Force differential equation in Mathematica -



Can't solve Coulomb Force differential equation in Mathematica -

i'm trying solve problem: have find trajectory of electron in graphene lattice using mathematica. i've tried solve coulomb forcefulness equation ndsolve , plot result each direction, obtain white plot. help me please? give thanks in advance. here's code x direction:

coordx = {0.6327, 1.88058, 3.03927, 4.28716, 5.44584, 6.69373, 7.85241, 9.10029, 1.9728, 3.22069, 4.37937, 5.62726, 6.78594, 8.03382, 9.19251, 10.4404, 3.3129, 4.56079, 5.71947, 6.96736, 8.12604, 9.37393, 10.53261, 11.7805, 4.653, 5.90089, 7.05956, 8.30746, 9.46614, 10.71403, 11.87271, 13.1206}; me = 9.01*10^-31; pi = 3.14159; epsilon0 = 8.854*10^-12; q = -1.6*10^-19; q = 1.6*10^-19; step = 0.01; forzax[p_, r_] := sum[(q*q)/(4 pi*epsilon0*norm[r - p[[i]]]^2), {i, length[p]}] forzax[coordx, {x[t]}]; ndsolve[{x''[t] == forzax[coordx, {x[t]}]/me , x[0] == 0, x'[0] == 0}, {x[t]}, {t, 0, 1500}] show[parametricplot[evaluate[{x[t]} /. %], {t, 0, 1500}, plotrange -> all]]

i don't know you're trying plot, these few modifications seem plot function.

sol = ndsolve[{x''[t] == forzax[coordx, {x[t]}]/me, x[0] == 0, x'[0] == 0}, {x}, {t, 0, 1500}]; f = sol[[1, 1, 2]]; plot[f[t], {t, 0, 1500}, plotrange -> all]

wolfram-mathematica

Comments

Post a Comment

Popular posts from this blog

php - Edges appear in image after resizing -

php - Edges appear in image after resizing - i have implemented next function: private static function generatepicture($sizekey, $src, $initialwidth, $initialheight, $imagetype) { $destination = adimage::addsizekey($sizekey, $src); if (file_exists($destination)) { return; } $finalsize = adimage::$sizes[$sizekey]; $initialratio = $initialwidth / $initialheight; $finalratio = $finalsize["w"] / $finalsize["h"]; $newimage = imagecreatetruecolor($finalsize["w"], $finalsize["h"]); imagealphablending( $newimage, false ); imagesavealpha( $newimage, true ); $oldimage = imagecreatefromstring(file_get_contents($src)); if ($initialratio >= $finalratio) { $scale = $finalsize["h"] / $initialheight; $surpluss = abs($initialwidth - ($initialheight * $finalratio)); imagecopyresized($newimage, $oldimage, 0, 0, $surpluss / 2, 0, $finalsize["w"], $finalsize[&q
Read more

ios8 - iOS custom keyboard - preserve state between appearances -

ios8 - iOS custom keyboard - preserve state between appearances - i'm having performance problems custom keyboard i'm working on. loading words in spell correction tree takes quite bit of time. seems done each time keyboard appears - there way preserve state of keyboard? apps can suspend / resume etc - can't find documentation on how extensions, or if there's mechanism doing att all. thanks! you loading dictionary net or disk? first case nsuserdefaults, encoding (nscopying) options store info. can load in background. according docs, os typically kills extension process, so, think there no way prevent keyboard objects deallocation (you may utilize nskeydarchiver). ios8 ios-app-extension custom-keyboard
Read more

Delphi change the assembly code of a running process -

Delphi change the assembly code of a running process - i'm trying alter address 00741fa5 has push test.009e721c . alter push test.009e71c8 . procedure callback; asm force $9e71c8 end; procedure tform2.btn1click(sender: tobject); var ppid :dword; pprocess : cardinal; begin getwindowthreadprocessid(findwindow(nil,pchar('test')), @ppid); pprocess := openprocess(process_all_access, false, ppid); injectasm(pprocess, $741fa5, 10, integer(@callback), 0); end; here injectasm function (with error checking) function injectasm(process: longword; injectaddress, injectsize, codeaddress, codesize : integer): pointer; var csize, rsize : integer; replaced : array of byte; jmp : array [0..4] of byte; jmpaddress : integer; nopv : byte; : integer; nbr: ulong_ptr; begin //injectsize must equal or greater 5, because need space our //(far) jmp wwxxyyzz instruction //if there's no space, inject in place of few instructions if injectsize &
Read more