From: parul agrawal on
syms z
Ps = max_signal_power*(10^(-total_loss_s/(10*z)));
Pp = max_pumping_power*(10^(-total_loss_p/(10*z)));
Ps_desh = Ps/Ps_saturation;
Pp_desh = Pp/Pp_saturation;
A = Nt*(Pp_desh + Ps_desh);
B = (1+(2*Ps_desh)+Pp_desh);
N2 = A/B;
N1 = Nt-N2;

L = input('enter the length in meter');
L_desh = [0:1:L];
l = length(L_desh);

F1 = int(N2,z,0,0);
F2 = int(N1,z,0,0);
F = ((sigma_signal_emission*F1)-(sigma_signal_absorption*F2));
G = gama_s*exp(F);

for i =L_desh(1,2):L_desh(1,l);
F1 = int(N2,z,0,i);
F1 = double(F1);
F2 = int(N1,z,0,i);
F2 = double(F2);
F = ((sigma_signal_emission*F1)-(sigma_signal_absorption*F2));
F = double(F);
G = [G,gama_s*exp(F)];
G = double(G);
end
plot(L_desh,G);
grid;




its giving first taht
explicit integral is not found ...when i tried to debug it....i found that problem is in reshape function(inside coding of int function)...when i checked the value of' ps' and tried to change it 'double' it is giving :

>>double(Ps)
??? Error using ==> reshape
To RESHAPE the number of elements must not change.


help plz