\( (*METODA\ DYFUZYJNA\ zastosowana\ \((bez\ sukcesu)\)\ \ do\ pewnej\ sumy\ funkcji\ Gaussa . \ Fukcja\ V \((x, y)\)\ ma
\ 
        dwa\ \(minima : \ jedno\ plytsze\ \((po\ lewej)\)\), \ 
      drugie\ glebsze\ \((globalne, \ 
              po\ prawej)\) . \ 
Metoda\ dyfuzyjna\ polega\ na\ odksztalceniu\ funkcji\ V -> Vt\ takim, \ 
      ze\ Vt\ 
ma\ juz\ tylko\ jedno\ minimum, \ 
      ktore\ z\ latwoscia\ znajdujemy . \ 
            Potem\ nastepuje\ sekwencja\ par\ \(operacji : \ 
zmniejszenie\ odksztalcenia\ i\ minimalizacja\ az\ do\ uzyskania\ ksztaltu\ Vt\) = 
        V . \ Wtedy\ 
"powinnismy byc"\ w\ minimum\ globalnym . \ Tymczasem\ tu\ tak\ nie\ jest . \ 
            Minimum\ globalne\ ma\ studnie\ tak
waska, \ 
      ze\ odksztalcenie\ funkcji\ powoduje\ jego\ znikniecie\ i\ metoda\ wskazuje\ na\ minimum\ lokalne\ \((lewe)\)*) 
V = \(-0.8\)\ E\^\(\(-0.1\)\ \((x + 5)\)\^2 - 0.2\ \((y - 4)\)\^2\) - 
        3\ E\^\(\(-\((x - 5)\)\^2\) - \((y + 4)\)\^2\)\)
Plot3D[V, {x, -8, 8}, {y, -8, 8}, 
      PlotPoints -> 30]
ContourPlot[V, {x, -8, 8}, {y, -8, 8}]
t = 20
Vt = -((0.8 (E^(-((0.1 (((x + 5))^2))/(1 + 4 0.1 t)) - ((0.2 (((y - 4))^2))/(1 + 4 0.2 t)))))/(Sqrt[(1 + 4 0.1 t) (1 + 4 0.2 t)])) - ((3 (E^(-((1.  (((x - 5))^2))/(1 + 4 t)) - ((1.  (((y + 4))^2))/(1 + 4 t)))))/(Sqrt[(1 + 4 t) (1 + 4 t)]))
Plot3D[Vt, {x, -8, 8}, {y, -8, 8}, 
      PlotPoints -> 30]
ContourPlot[Vt, {x, -8, 8}, {y, -8, 8}]
t = 10
Vt = -((0.8 (E^(-((0.1 (((x + 5))^2))/(1 + 4 0.1 t)) - ((0.2 (((y - 4))^2))/(1 + 4 0.2 t)))))/(Sqrt[(1 + 4 0.1 t) (1 + 4 0.2 t)])) - ((3 (E^(-((1.  (((x - 5))^2))/(1 + 4 t)) - ((1.  (((y + 4))^2))/(1 + 4 t)))))/(Sqrt[(1 + 4 t) (1 + 4 t)]))
Plot3D[Vt, {x, -8, 8}, {y, -8, 8}, 
      PlotPoints -> 30]
ContourPlot[Vt, {x, -8, 8}, {y, -8, 8}]
t = 2.5
Vt = -((0.8 (E^(-((0.1 (((x + 5))^2))/(1 + 4 0.1 t)) - ((0.2 (((y - 4))^2))/(1 + 4 0.2 t)))))/(Sqrt[(1 + 4 0.1 t) (1 + 4 0.2 t)])) - ((3 (E^(-((1.  (((x - 5))^2))/(1 + 4 t)) - ((1.  (((y + 4))^2))/(1 + 4 t)))))/(Sqrt[(1 + 4 t) (1 + 4 t)]))
Plot3D[Vt, {x, -8, 8}, {y, -8, 8}, 
      PlotPoints -> 30]
ContourPlot[Vt, {x, -8, 8}, {y, -8, 8}]
t = 0.5
Vt = -((0.8 (E^(-((0.1 (((x + 5))^2))/(1 + 4 0.1 t)) - ((0.2 (((y - 4))^2))/(1 + 4 0.2 t)))))/(Sqrt[(1 + 4 0.1 t) (1 + 4 0.2 t)])) - ((3 (E^(-((1.  (((x - 5))^2))/(1 + 4 t)) - ((1.  (((y + 4))^2))/(1 + 4 t)))))/(Sqrt[(1 + 4 t) (1 + 4 t)]))
Plot3D[Vt, {x, -8, 8}, {y, -8, 8}, 
      PlotPoints -> 30]
ContourPlot[Vt, {x, -8, 8}, {y, -8, 8}]
t = 0. 
Vt = -((0.8 (E^(-((0.1 (((x + 5))^2))/(1 + 4 0.1 t)) - ((0.2 (((y - 4))^2))/(1 + 4 0.2 t)))))/(Sqrt[(1 + 4 0.1 t) (1 + 4 0.2 t)])) - ((3 (E^(-((1.  (((x - 5))^2))/(1 + 4 t)) - ((1.  (((y + 4))^2))/(1 + 4 t)))))/(Sqrt[(1 + 4 t) (1 + 4 t)]))
Plot3D[Vt, {x, -8, 8}, {y, -8, 8}, 
      PlotPoints -> 30]
ContourPlot[Vt, {x, -8, 8}, {y, -8, 8}]