// Example 7.10.1  page 7.53

clc;
clear;

lamda=0.85d-6;
h=6.626d-34;    //plank's constant
c=3d8;      //speed of light
q=1.6d-19;  //charge of electron
eta=75/100; //quantum efficiency
P0=0.6d-6;  //incident optical power
Im=15d2;   //avalanche gain

R= eta*q*lamda/(h*c);   //computing responsivity
Ip=10^8*P0*R;       //computing photocurrent
Ip=floor(Ip);
M=Im/Ip;        //computing multiplication factor
printf("\nMultiplication factor is %d.",M);