jmarangola / dynamical-coulomb-blockade-model Goto Github PK

Write P(E) to file or .csv so it can be used by the second program to perform I(v) plots separately.

must keep vars symbolic

get fullSumsEwithRealValues working with physical constants.

Summing two integrals over [a, b-delta] and [b+delta, c] (where a<<0, b=0, c>>0), the function P(E) appears to oscillate around an exact number as a and c get proportionally smaller and larger respectively using vpaintegral. There appears to be no error bound (which could either mean it is not a displayed output (or only outputted in an overloaded function call unlike integral and integral2 or it could be a really good sign meaning there is actually no error using vpaintegral with arbitrarily large a and c). If vpaintegral is to approximate P(E) or J(t) or any combination of the two, a and c must be small enough so that it is able to sequentially compute integrals and large enough so that P(E) and J(t) are accurate.

determine if the dimensions of P(E) plot are properly designated as that determines the output of the approximation

Currently spits back int functions due to the combination of nested infinite (or semi-infinite bounds) and complex integrands. Potential solution is breaking up each integral into a series of smaller integrals with nested loops --> eats up computation time and resources and may not be accurate however