dltn = 1e-3；

V = 2.4

lam MDA = 1550 e-9；

C=pi。 /lammda。 *( 1- 1/(V.^2)).* dltn

R=tanh(C.*L)。 ^2;

Figure (1);

plot(L，R)；

Hold on;

L=linspace(0，4.5e-3， 1000)；