I am not sure which plots are the ones you are getting and which ones are the ones you want, but I have followed the Intusoft method here:
http://www.intusoft.com/nlhtm/nl78.htm#ref2
and get some sensible results. It didn't seem as straightforward as it should be though!
With my simulator I set up the following in addition to what you can see in the schematic:
.param imax=0.33
.param vmax=21
.param is=1e-11
.param n={vmax*38.6/log(imax/is)}
.model DTEST d(is={is} n={n} xti={3*n} eg={1.11*n})
This is equivalent of what they are doing with the text at the side but in a "standard" syntax (I think).
The simulation I ran was with a fixed I1 current using:
.DC Rload 1 100 198.396793587174m
This sweeps the load from 1 ohm to 100 ohms. The strange step size is simply the simulator calculating the step size to make it 100 steps in the sweep.
Then, the optimum power curve comes from plotting power in the load on the Y axis against output voltage on the X axis. The current curve is from plotting output current against output voltage. Both of these plots require you to plot graphs against something other than what you swept in the DC analysis. You may not be used to doing this as it is not the normal way of doing things.
You need to set up a diode model as described to get the correct results, although it should really be based on your own solar panel values.
So, I think your problems are probably one or more of the following:
1. not creating the diode model,
2. not sweeping the right parameter (Rload)
3. not plotting the right parameters (in particular, the X axis is Vout).
So, while the model is very simple, the implementation of it is a bit tricky!
I hope that helps.
Keith.