function?[?I?]?=?F_GS（?I1I2k?）
%UNtitleD4?Summary?of?this?function?goes?here
%???Detailed?explanation?goes?here
%???I1為輸入光場的振幅
%???I2為輸出光場的振幅
%???k為算法程序的迭代次數
%???I為輸出的光場，為復數，包括相位和振幅

[mn]=size（I1）;
I1=double（I1）;
I2=double（I2）;

I1re=rand（mn）;
I1lm=rand（mn）;
I1x=I1re+i*I1lm;
I1x=I1x./abs（I1x）;%隨機生成輸入面的初始相位

I2q=I2.^2;
sumG=sum（sum（I2q））;
sse=zeros（1k）;
l=0;

for?l=1:k
????????