CODE:

function delar

data=xlsread('delar.xlsx');

xdata=data(:,1);  %Im

ydata=data(:,2);  %

x0=[100,250,6,0.5];    %wp2,gama,T,S

options= optimset('MaxFunEvals',5000,'Maxiter',1000);

lb=[10,20,0,0];

%ub=[10000,10000,100,1000];

ub=[];

[x,resnorm]=lsqcurvefit(@myfun,x0,xdata,ydata,lb,ub,options)

y=myfun(x,xdata);

plot(xdata,ydata,'g',xdata,y,'r')

end

function F=myfun(x,xdata)  %x=[]

wl=734;

wt=533;

e8=5.5;

C=0.48;

eq1=(wl^2-wt^2)./(wt^2-xdata.^2-1i*x(3)*xdata);

eq2=x(1)./(xdata.^2+1i*xdata*x(2));

epslong=e8*(1+eq1-eq2);

delta=x(1)*x(2)*((wt^2-xdata.^2).^2+(xdata*x(3)).^2)  +  xdata.^2*x(3)*(wl^2-wt^2).*(xdata.^2+x(2)^2);

eq3=2*C*wt^2* (x(1)*x(2)*(wt^2-xdata.^2)   -  xdata.^2*x(3).*(xdata.^2+x(2)^2-x(1)))   ./delta;

eq4=C^2*(wt^4./delta)  .*(x(1)*(x(2)*(wl^2-wt^2)+x(3)*(x(1)-2*xdata.^2))  +xdata.^2*x(3).*(xdata.^2+x(2)^2))  ./(wl^2-wt^2);

A=1+eq3+eq4;

F=x(4)*A.*imag(-1./epslong);

end