How to do PAT image reconstruction

Hello I have to reconstruct image by model based inversion for a circular medium in 2D system. I have witten the following code for solving

Figure a is true image and figure b is reconstructed image. I want output like figure b.

for iter=1:200

Pc=W*p02;%%p02 is the randomly initialize value of initial pressure

delta_P=abs(Pc)-abs(P);

delta_p02=(J1'*J1+lambda*I)\(J1'*delta_P);

p02=p02-0.8*delta_p02;

delta_P1=angle(Pc)-angle(P);

delta_p002=(J2'*J2+lambda*I)\(J2'*delta_P1);

p02=p02-0.8*delta_p002;

[t1,t]=find(p02<0);p02(t1,t)=0;

max_err_re(iter)=norm(delta_P);

max_err_im(iter)=norm(delta_P1);

if (abs(delta_P)<0.6000)

sol=p0_2;

break

end

sol=p0_2;

end

end

but the result is not coming.