blob: f36c1a590fa3e598206380dbea6e8fcb774fe026 (
plain)
1
2
3
4
5
6
7
8
9
10
11
|
clc
// compare the intensities of ordinary and extraordinary rays
//intensity of ordinary rays is given by Io=a^2 *(sin theta)^2
//where theta=30 degree
//we get Io=a^2/4
Io=1/4
//intensity of extraordinary ray is given by IE=(a*cos theta)^2
//we get IE=3*a^2/4
IE=3/4
I=IE/Io
disp("the intensities of ordinary and extraordinary rays is I="+string(I)+"unitless")
|
