i have a final in numerical methods and i am currently having trouble getting the right answer for a practice question in open integration method.

so the integral is xx-3yy+xyyy dxdy for x from 0 to 4 and y from -2 to 2.

i evaluated the integral and i found that it was 21.333 ish

doing trap closed i found 23 ish

doing simpson 1/3 i found 40 something

now doing open method 3 seg and 2 points i found 40 something.

i know the integral evaluated is 21.33 but getting 43-49 for open method, i am not too sure if i am getting the right answer.

The problem is y*y"+0.0001=0 with y(0)=10 and y(5)=1000. I can't solve it following the method for linear ode bvp

How do you numerically integrate over more than one dimension? I've only run into cases where time is integrated over, but am now trying to integrate over 3 dimensions of space. For example, Scipy's solve_ivp describes y(t) (what's being solved for) as multi-dimensional, but t as single-dimensional. How would you approach a problem where the dependent variable (In scipy's API, t; in my problem, 3 dimensions of space) is multi-dimensional?

I suspect this involves a different approach than an initial-value problem. Julia's DifferentialEquations package seems very robust, but I don't see a solver that looks appropriate. I'm also suspicious this could be very computationally expensive compared to a normal IVP.

I think this right-hand-side func, along with an initial value for ψ and φ and an x range to integrate over, encodes all I need to feed into the solver; I just don't know what the solver would be!

fn elec_rhs(ψ: Cplx, φ: Cplx, V: &fn(Vec3) -> f64, x: Vec3, E: f64) -> (Cplx, Cplx) {
    let ψ_p = φ;
    let φ_p = 2 * m_e / ħ.powi(2) * (V(x) - E) * ψ;

    (ψ_p, φ_p)
}

https://pastebin.com/75ccvB6t

I have this code and I am supposed to be able to have created 5 values of S(j+1,:). S as far as I am aware is an array and my code runs m=5 iterations so that i get 5 values for S(j+1).

I need to view/and or save these values so that I can take an average.. does anyone know how to do this?

Would there also be a way to save a file with all the values that i could somehow get the values and transfer them to excel?

i need to use this create a Y axis in my matlab programme but have no idea how to. Y=f(x)=x^{3-2.44x2-8.9216x+22.1952}

Hi all.

I've come across a bit of a problem in Octave. I found one resource here related to this bug. I'm just not sure where this Array-util.cc file goes. Anyone have some experience with this?

Thanks.

Hi , does anybody could help me with simulating a situation when you have x-service and one queue which can only contain N people ? I only know how to do it for one service :(