Numerical calculation of the supercritical BNLS asymptotic (S) profile.

These matlab scripts and functions calculate the S profile, as given in:
G. Baruch, G. Fibich, "Singular Solutions of the L^2 supercritical biharmonic nonlinear Schroedinger Equation" mansucript.

The numerical method is a variant of the Petviashvili (SLSR) method, and is
given in the manuscript's appendix.


In order to calculate the d=1, sigma=6 profile:
under matlab, type:
    W = calc_S_main(1, 6, 1.037, 0.361870, 4, 160, 32000)

In order to calculate the d=2, digma=3 profile: 
under matlab, type:
    W = calc_S_main(2, 3, 0.909, 0.228251, 4, 160, 32000)

The results are given in the W struct.

You may also want to plot or verify the WKB approximations. 
In order to do that, under matlab, type:
    plot_solns

For questions, please contact Guy Baruch:
guy.baruch@math.tau.ac.il
    
