runMc36sm.jl

include("mc36sm_ss.jl")
include("mc36sm.jl")
function runGj()
    println("starting..")
    vj = -[ zeros(1, 250) LinRange(0, 120, 750)' ];
    gj = zeros(length(vj))
    parAllFast = [ .15 28 24e-9 * 4 24e-9 * 4 / 10 800 50 .01 -1;       .1 36 140e-9 * 4 140e-9 * 4 / 10 1e10 1e10 .03 1 ];
    parAllSlow = [ .1 48 140e-9 * 4  1e-15 1e10 1e10 .03 -1;            .1 48 140e-9 * 4  1e-15 1e10 1e10 .03 1 ];
    par = [parAllFast; parAllSlow];
    pc1c2 = 0;
    pc2c1 = 0;
    pc1c2 = 0.01;
    pc2c1 = 0.001;
    ppp = zeros(1, 36);
    ppp[1] = 1;
    for i = 1:length(vj)
        gj[i] = MC36SM_Mindaugo_SS(vj[i], par, ppp, pc1c2, pc2c1);
    end
    println("done.")
    return vj, gj
end