summaryrefslogtreecommitdiff
path: root/2939/CH1/EX1.8/Ex1_8.sce
blob: f21b9ee28c62dc81225d1b04c082caeeeb9380b4 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
//Ex1_8

clc;

//Given:
e1=1; // in MeV
e2=2; // in MeV
ma=4; // in u(amu)
md=2; // in u(amu)
mp=1; // in u(amu)

// 1u = 1.6*10^-27 Kg

//solution: part a)For alpha particles

v1a=((2*e1*10^6*1.6*10^-19)/(ma*1.6605*10^-27))^.5;
printf("\n The velocity of alpha particles for 1 MeV in m/s is = %f ",v1a)// For 1 MeV

v2a=((2*e2*10^6*1.6*10^-19)/(ma*1.6605*10^-27))^.5;
printf("\n The velocity of alpha particles for 2 MeV in m/s is = %f ",v2a)// For 2 MeV

//solution: part b)For deuteron particles

v1b=((2*e1*10^6*1.6*10^-19)/(md*1.6605*10^-27))^.5;
printf("\n The velocity of deuteron particles for 1 MeV in m/s is = %f ",v1b) // For 1 MeV


v2b=((2*e2*10^6*1.6*10^-19)/(md*1.6605*10^-27))^.5;
printf("\n The velocity of deuteron particles for 2 MeV in m/s is = %f ",v2b) // For 2 MeV

//solution: part c)For proton particles

v1p=((2*e1*10^6*1.6*10^-19)/(mp*1.6605*10^-27))^.5;
printf("\n The velocity of proton particles for 1 MeV in m/s is = %f ",v1p) // For 1 MeV


v2p=((2*e2*10^6*1.6*10^-19)/(mp*1.6605*10^-27))^.5;
printf("\n The velocity of proton particles for 2 MeV in m/s is = %f ",v2p) // For 2 MeV