blob: 40bfbcdb8ba22a14f54f83d79d886f4cc0b694b5 (
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
|
// Given:-
p1 = 3*(10**5) // initial pressure of gas in pascal
v1 = 0.1 // initial volumme of gas in meter^3
v2 = 0.2 // final volume of gas in meter^3
// calculations
// Part (a) i.e. n=1.5
//constant = p1*(v1**n) // p*(v^n) = constant
constant1 = p1*(v1**1.5)
constant2 = p1*(v1**1)
constant3 = p1*(v1**0)
// function p
function v = p1(v)
v = constant1/(v^1.5)
endfunction
function v = p2(v)
v = constant2/(v^1)
endfunction
function v = p3(v)
v = constant3/(v^0)
endfunction
work1 = intg(v1,v2,p1) // integrating pdv from initial to final volume
w1 = work1(1)/1000 // divided by 1000 to convert to KJ
printf( 'The work done for n=1.5 in KJ is %.2f',w1)
//part(b) i.e. n = 1
work2 = intg(v1,v2,p2)
w2 = work2(1)/1000
printf( 'The work done for n=1 in KJ is %.2f',w2)
//part(c) i.e. n=0
work3 = intg(v1,v2,p3)
w3 = work3(1)/1000
printf( 'The work done for n=0 in KJ is %.2f',w3)
|
