howardjp

Quantum Mechanics for Engineers

Problem 1

Find the fringe spacing in double-slit experiment where a particle is an electron at 1 MeV, and other values included below.

howardjp

import sepsis

sepsis.run_all(intensity=True, device='cuda')
sepsis.run_all(intensity=False, device='cuda')

howardjp

howardjp

howardjp

extern crate rand;

use rand::Rng;

fn montecarlo_pi(trials: u32) -> f64 {
    let mut count: u32 = 0;
    let mut rng = rand::thread_rng();

    for _i in 1..trials {
        let x: f64 = rng.gen::<f64>();

howardjp

### This is a mostly ggplot-ish implementation of the plot for dbchoice
### and sbchoice models in the DCchoice package.

ggplot.dbchoice <- function(x, ...) {
    minbid <- min(x$bid)
    maxbid <- max(x$bid)
    bidseq <- seq(minbid, maxbid, by = (maxbid - minbid) / 100)
    
    b <- x$coefficients
    npar <- length(b)

howardjp

-----BEGIN PGP PUBLIC KEY BLOCK-----
 
mDMEVoAa3xYJKwYBBAHaRw8BAQdAjmFdyF84L4mnsxlWQePC1nlWmCfvFa/AnA7L
cSFSrJG0GEphbWVzIFBhdHJpY2sgSG93YXJkLCBJSYh8BBMWCgAkAhsBBQsJCAcD
BRUKCQgLBRYCAwEAAh4BAheABQJYErV1AhkBAAoJEH4a3z4YNYeRVuwBAKFslfxr
/fWgttA4m6/AQ7FJ0APdOyMDRbm/b4ZJEJ8eAP478Bla7UMZgzYhQ5hQs9qFlNKd
hh5BfYMqVZnKKLnSA7QsSmFtZXMgUGF0cmljayBIb3dhcmQsIElJIDxqaEBqYW1l
c2hvd2FyZC51cz6IeQQTFgoAIQUCVoAcVwIbAQULCQgHAwUVCgkICwUWAgMBAAIe
AQIXgAAKCRB+Gt8+GDWHkUxPAQCPuEdmr9Dq/Y6z1iMWGk3v17PzcJQhzcCIw+bP
oC0SEwD7BRCwvIRETDDK226t1toDmUJ6EQ5yo5BmqD9BEQ1g+gO0L0phbWVzIFBh

howardjp

## Computing a bootsrap confidence interval
bootCI.parallel <- function (obj, nboot = 1000, CI = 0.95, individual = NULL, threads = 1) {

    ## Manage any potential incorrect inputs, error out as approrpiate
    if(!inherits(obj, c("sbchoice", "dbchoice", "oohbchoice")))
        stop("the object must be sbchoice, dbchoice, or oohbchoice class")
    if(CI > 1 | CI < 0)
        stop("CI must be between 0 and 1")
    if(!is.null(threads) && !is.numeric(threads))
        stop("threads must be an integer or NULL for autodetection")

howardjp

howardjp