2 views (last 30 days)
Mustafa Batuhan Turaç on 5 Apr 2023
Edited: VBBV on 5 Apr 2023
I could not integrate the "K" value in the code I wrote into my function. In other words, my function has 6 different T values, 6 different K values. The V value changes at every T value, I integrated it with numel codes, but I cannot change the K value at every T value. If I write a single K value in the function, the code works, but I have to write 6 different codes. I would be glad if you help. Thanks a lot
Must-have K format and values
K = [1.250161 1.118431 1.06921 1.031865 1.004945 0.999667]
%Plot P-V diagram of a given substance using peng robinson equation of state
%Substance properties are defined by the peng robinson constants and
%the critical properties.
clc; clear; close all;
Tc = 765.62; % R
Pc = 550.60; % psi
Vc = 4.0860; % ft3/lbmole
R = 10.732 ; % psi ft3/(lbmole-R)
% peng robinson equation of state constants
% for n-butane
a = 56065.52286;
b = 1.161014;
% alpha
K = 1.250161;
V = linspace(b*1.2,40*Vc,1000); % vector of volume
% temperature in F
T = [60 180 230 270 300 306];
T = T + 460; % temperature in R
%peng robinson equation
fprEOSp = @(Tx,Vx)(R*Tx./(Vx-b)-a*K./((Vx.^2+Vx*b)+(Vx*b-b.^2)));
P = zeros(numel(T),numel(V));
for i= 1:numel(T)
Tx = T(i);
P(i,:) = fprEOSp(Tx,V);
end
% plot(V,P); xlim([0 100]);
semilogx(V,P); xlim([1 100]);
ylim([0 800]);
xlabel('Volume, ft^3');
ylabel('Pressure, psi');

MarKf on 5 Apr 2023
By "I cannot change the K value at every T value" I imagine you need K to change for every iteration and therefore have sitll 6 results. You can also have 36 results instead (6*6 T*K parameters). Otherwise it's not very clear what you need by your description. I'll assume the former with still 6 results, but the other case is also trivial to do.
%Plot P-V diagram of a given substance using peng robinson equation of state
Tc = 765.62; % R
Pc = 550.60; % psi
Vc = 4.0860; % ft3/lbmole
R = 10.732 ; % psi ft3/(lbmole-R)
a = 56065.52286; b = 1.161014;
K = [1.250161 1.118431 1.06921 1.031865 1.004945 0.999667];
V = linspace(b*1.2,40*Vc,1000); % vector of volume
T = [60 180 230 270 300 306]; % temperature in F
T = T + 460; % temperature in R
fprEOSp = @(Tx,Vx,Kx)(R*Tx./(Vx-b)-a*Kx./((Vx.^2+Vx*b)+(Vx*b-b.^2))); %just add Kx if you need to change K as well
P = zeros(numel(T),numel(V));
for i= 1:numel(T)
Kx = K(i); %maybe I did not get what you mean since you already got this far with Tx below
Tx = T(i);
P(i,:) = fprEOSp(Tx,V,Kx);
end
semilogx(V,P); xlim([1 100]);
ylim([-3100 800]); % I wanna see the whole thing
xlabel('Volume, ft^3'); ylabel('Pressure, psi');

VBBV on 5 Apr 2023
%Plot P-V diagram of a given substance using peng robinson equation of state
%Substance properties are defined by the peng robinson constants and
%the critical properties.
clc; clear; close all;
Tc = 765.62; % R
Pc = 550.60; % psi
Vc = 4.0860; % ft3/lbmole
R = 10.732 ; % psi ft3/(lbmole-R)
% peng robinson equation of state constants
% for n-butane
a = 56065.52286;
b = 1.161014;
% alpha
K = [1.250161 1.118431 1.06921 1.031865 1.004945 0.999667];
V = linspace(b*1.2,40*Vc,1000); % vector of volume
% temperature in F
T = [60 180 230 270 300 306];
T = T + 460; % temperature in R
%peng robinson equation
P = zeros(numel(K),numel(T),numel(V));
for k = 1:length(K)
fprEOSp = @(Tx,Vx)(R*Tx./(Vx-b)-a*K(k)./((Vx.^2+Vx*b)+(Vx*b-b.^2)));
for i= 1:numel(T)
Tx = T(i);
P(k,i,:) = fprEOSp(Tx,V);
end
end
P
P =
P(:,:,1) = 1.0e+04 * 0.5724 1.1270 1.3581 1.5430 1.6817 1.7094 0.7654 1.3200 1.5511 1.7359 1.8746 1.9023 0.8374 1.3921 1.6232 1.8080 1.9467 1.9744 0.8921 1.4468 1.6778 1.8627 2.0014 2.0291 0.9316 1.4862 1.7173 1.9021 2.0408 2.0685 0.9393 1.4939 1.7250 1.9099 2.0485 2.0763 P(:,:,2) = 1.0e+03 * -0.8174 2.4478 3.8083 4.8967 5.7130 5.8763 0.7597 4.0249 5.3854 6.4738 7.2901 7.4534 1.3489 4.6141 5.9746 7.0630 7.8794 8.0426 1.7960 5.0612 6.4217 7.5101 8.3264 8.4897 2.1183 5.3835 6.7440 7.8324 8.6487 8.8120 2.1815 5.4467 6.8072 7.8956 8.7119 8.8752 P(:,:,3) = 1.0e+03 * -2.5111 -0.1974 0.7666 1.5378 2.1162 2.2319 -1.1901 1.1236 2.0876 2.8589 3.4373 3.5530 -0.6965 1.6172 2.5812 3.3525 3.9309 4.0466 -0.3220 1.9917 2.9557 3.7270 4.3054 4.4211 -0.0520 2.2617 3.2257 3.9969 4.5754 4.6910 0.0009 2.3146 3.2786 4.0499 4.6283 4.7440 P(:,:,4) = 1.0e+03 * -2.9359 -1.1443 -0.3978 0.1993 0.6472 0.7368 -1.8085 -0.0169 0.7296 1.3268 1.7746 1.8642 -1.3873 0.4043 1.1508 1.7480 2.1959 2.2855 -1.0676 0.7239 1.4704 2.0676 2.5155 2.6051 -0.8372 0.9543 1.7008 2.2980 2.7459 2.8355 -0.7921 0.9995 1.7460 2.3432 2.7911 2.8807 P(:,:,5) = 1.0e+03 * -2.9326 -1.4708 -0.8618 -0.3745 -0.0091 0.0640 -1.9561 -0.4944 0.1147 0.6019 0.9673 1.0404 -1.5913 -0.1295 0.4795 0.9668 1.3322 1.4053 -1.3145 0.1473 0.7563 1.2436 1.6090 1.6821 -1.1149 0.3468 0.9559 1.4431 1.8085 1.8816 -1.0758 0.3859 0.9950 1.4822 1.8477 1.9208 P(:,:,6) = 1.0e+03 * -2.7733 -1.5388 -1.0245 -0.6130 -0.3044 -0.2426 -1.9174 -0.6829 -0.1686 0.2429 0.5515 0.6132 -1.5976 -0.3631 0.1512 0.5627 0.8713 0.9330 -1.3549 -0.1205 0.3939 0.8053 1.1140 1.1757 -1.1800 0.0544 0.5688 0.9803 1.2889 1.3506 -1.1457 0.0887 0.6031 1.0145 1.3232 1.3849 P(:,:,7) = 1.0e+03 * -2.5614 -1.4930 -1.0479 -0.6918 -0.4247 -0.3713 -1.8037 -0.7353 -0.2902 0.0659 0.3330 0.3864 -1.5205 -0.4522 -0.0071 0.3491 0.6161 0.6696 -1.3057 -0.2374 0.2077 0.5639 0.8309 0.8844 -1.1509 -0.0826 0.3626 0.7187 0.9858 1.0392 -1.1205 -0.0522 0.3929 0.7491 1.0161 1.0696 P(:,:,8) = 1.0e+03 * -2.3394 -1.3978 -1.0054 -0.6915 -0.4561 -0.4090 -1.6629 -0.7213 -0.3289 -0.0151 0.2203 0.2674 -1.4102 -0.4685 -0.0762 0.2377 0.4731 0.5202 -1.2184 -0.2768 0.1156 0.4295 0.6649 0.7120 -1.0802 -0.1385 0.2538 0.5677 0.8031 0.8502 -1.0531 -0.1114 0.2809 0.5948 0.8302 0.8773 P(:,:,9) = 1.0e+03 * -2.1252 -1.2834 -0.9326 -0.6520 -0.4416 -0.3995 -1.5169 -0.6751 -0.3243 -0.0437 0.1667 0.2088 -1.2896 -0.4478 -0.0971 0.1835 0.3940 0.4361 -1.1171 -0.2753 0.0754 0.3560 0.5664 0.6085 -0.9928 -0.1510 0.1997 0.4803 0.6908 0.7328 -0.9685 -0.1267 0.2241 0.5047 0.7151 0.7572 P(:,:,10) = 1.0e+03 * -1.9259 -1.1648 -0.8477 -0.5940 -0.4037 -0.3656 -1.3754 -0.6143 -0.2972 -0.0435 0.1468 0.1848 -1.1698 -0.4087 -0.0915 0.1622 0.3524 0.3905 -1.0137 -0.2526 0.0645 0.3182 0.5085 0.5465 -0.9012 -0.1401 0.1770 0.4307 0.6210 0.6590 -0.8792 -0.1181 0.1991 0.4528 0.6430 0.6811 P(:,:,11) = 1.0e+03 * -1.7437 -1.0492 -0.7598 -0.5283 -0.3547 -0.3200 -1.2429 -0.5483 -0.2590 -0.0275 0.1462 0.1809 -1.0557 -0.3612 -0.0718 0.1597 0.3333 0.3680 -0.9137 -0.2192 0.0702 0.3017 0.4753 0.5100 -0.8114 -0.1169 0.1725 0.4040 0.5777 0.6124 -0.7913 -0.0968 0.1926 0.4241 0.5977 0.6325 P(:,:,12) = 1.0e+03 * -1.5787 -0.9401 -0.6740 -0.4611 -0.3014 -0.2695 -1.1208 -0.4821 -0.2160 -0.0031 0.1565 0.1885 -0.9496 -0.3110 -0.0449 0.1680 0.3277 0.3596 -0.8198 -0.1812 0.0849 0.2978 0.4575 0.4894 -0.7262 -0.0876 0.1785 0.3914 0.5511 0.5830 -0.7079 -0.0692 0.1969 0.4098 0.5694 0.6014 P(:,:,13) = 1.0e+03 * -1.4299 -0.8388 -0.5925 -0.3955 -0.2477 -0.2181 -1.0093 -0.4182 -0.1719 0.0251 0.1729 0.2024 -0.8522 -0.2611 -0.0148 0.1822 0.3300 0.3596 -0.7330 -0.1419 0.1044 0.3015 0.4492 0.4788 -0.6470 -0.0559 0.1904 0.3874 0.5352 0.5648 -0.6302 -0.0391 0.2072 0.4043 0.5520 0.5816 P(:,:,14) = 1.0e+03 * -1.2959 -0.7458 -0.5165 -0.3332 -0.1956 -0.1681 -0.9082 -0.3580 -0.1288 0.0546 0.1921 0.2196 -0.7633 -0.2131 0.0161 0.1995 0.3370 0.3645 -0.6534 -0.1032 0.1260 0.3094 0.4469 0.4744 -0.5741 -0.0240 0.2053 0.3886 0.5262 0.5537 -0.5586 -0.0084 0.2208 0.4042 0.5417 0.5692 P(:,:,15) = 1.0e+03 * -1.1753 -0.6608 -0.4464 -0.2749 -0.1463 -0.1206 -0.8165 -0.3020 -0.0877 0.0838 0.2125 0.2382 -0.6825 -0.1680 0.0464 0.2179 0.3465 0.3722 -0.5808 -0.0663 0.1481 0.3196 0.4482 0.4739 -0.5074 0.0070 0.2214 0.3929 0.5215 0.5473 -0.4931 0.0214 0.2358 0.4073 0.5359 0.5616 P(:,:,16) = 1.0e+03 * -1.0666 -0.5834 -0.3821 -0.2210 -0.1002 -0.0761 -0.7336 -0.2504 -0.0491 0.1120 0.2328 0.2569 -0.6092 -0.1260 0.0753 0.2364 0.3572 0.3814 -0.5148 -0.0316 0.1697 0.3308 0.4516 0.4758 -0.4467 0.0365 0.2378 0.3989 0.5197 0.5438 -0.4334 0.0498 0.2511 0.4122 0.5330 0.5572 P(:,:,17) = -968.5825 -513.1205 -323.3447 -171.5241 -57.6586 -34.8855 -658.5563 -203.0943 -13.3185 138.5022 252.3677 275.1408 -542.7147 -87.2528 102.5231 254.3437 368.2092 390.9823 -454.8234 0.6386 190.4144 342.2351 456.1006 478.8737 -391.4672 63.9948 253.7706 405.5912 519.4567 542.2298 -379.0455 76.4165 266.1923 418.0130 531.8785 554.6516 P(:,:,18) = -880.0152 -449.2643 -269.7847 -126.2011 -18.5133 3.0242 -590.6045 -159.8536 19.6260 163.2096 270.8973 292.4349 -482.4660 -51.7151 127.7645 271.3481 379.0358 400.5734 -400.4191 30.3318 209.8114 353.3950 461.0828 482.6203 -341.2759 89.4750 268.9546 412.5382 520.2260 541.7635 -329.6801 101.0708 280.5504 424.1340 531.8217 553.3593 P(:,:,19) = -799.8539 -391.2707 -221.0277 -84.8332 17.3126 37.7417 -529.0117 -120.4285 49.8145 186.0090 288.1548 308.5839 -427.8114 -19.2281 151.0149 287.2093 389.3551 409.7843 -351.0286 57.5547 227.7977 363.9921 466.1380 486.5671 -295.6800 112.9033 283.1463 419.3407 521.4866 541.9157 -284.8282 123.7551 293.9981 430.1925 532.3383 552.7675 P(:,:,20) = -727.1665 -338.5809 -176.6703 -47.1418 50.0046 69.4339 -473.1145 -84.5289 77.3817 206.9103 304.0566 323.4859 -378.1878 10.3978 172.3084 301.8370 398.9833 418.4126 -306.1649 82.4207 244.3313 373.8598 471.0062 490.4355 -254.2475 134.3381 296.2487 425.7772 522.9236 542.3529 -244.0684 144.5171 306.4278 435.9563 533.1027 552.5320 P(:,:,21) = -661.1320 -290.6780 -136.3221 -12.8375 79.7761 98.2988 -422.3170 -51.8630 102.4929 225.9775 318.5911 337.1138 -333.0836 37.3704 191.7262 315.2109 407.8244 426.3471 -265.3804 105.0736 259.4295 382.9142 475.5277 494.0504 -216.5768 153.8772 308.2331 431.7178 524.3313 542.8540 -207.0082 163.4458 317.8016 441.2863 533.8998 552.4225 P(:,:,22) = -601.0294 -247.0902 -99.6156 18.3641 106.8489 124.5459 -376.0882 -22.1491 125.3256 243.3053 331.7901 349.4870 -292.0388 61.9004 209.3750 327.3547 415.8395 433.5365 -228.2687 125.6704 273.1450 391.1248 479.6095 497.3065 -182.3003 171.6388 319.1134 437.0932 525.5779 543.2749 -173.2877 180.6515 328.1261 446.1058 534.5906 552.2876 P(:,:,23) = -546.2251 -207.3913 -66.2105 46.7341 131.4426 148.3843 -333.9559 4.8780 146.0588 259.0034 343.7119 360.6536 -254.6413 84.1926 225.3733 338.3180 423.0264 439.9681 -194.4637 144.3702 285.5510 398.4956 483.2041 500.1457 -151.0849 187.7490 328.9298 441.8744 526.5828 543.5245 -142.5799 196.2539 337.4347 450.3793 535.0878 552.0295 P(:,:,24) = -496.1629 -171.1977 -35.7956 72.5261 153.7674 170.0157 -295.5012 29.4639 164.8661 273.1878 354.4291 370.6773 -220.5238 104.4413 239.8434 348.1652 429.4064 445.6547 -163.6370 161.3282 296.7303 405.0521 486.2933 502.5416 -122.6303 202.3349 337.7370 446.0587 527.3000 543.5483 -114.5904 210.3748 345.7769 454.0986 535.3399 551.5882 P(:,:,25) = -450.3531 -138.1660 -8.0881 95.9743 174.0210 189.6304 -260.3530 51.8341 181.9121 285.9744 364.0212 379.6306 -189.3593 122.8278 252.9057 356.9681 435.0149 450.6242 -135.4949 176.6922 306.7701 410.8325 488.8793 504.4886 -96.6670 215.5201 345.5981 449.6604 527.7072 543.3166 -89.0543 223.1328 353.2108 457.2731 535.3199 550.9293 P(:,:,26) = -408.3648 -107.9889 17.1677 117.2930 192.3870 207.4058 -228.1821 72.1939 197.3505 297.4758 372.5698 387.5886 -160.8566 139.5193 264.6759 364.8012 439.8952 454.9140 -109.7754 190.6005 315.7571 415.8824 490.9764 505.9952 -72.9538 227.4222 352.5788 452.7041 527.7981 542.8169 -65.7344 234.6415 359.7981 459.9234 535.0174 550.0362 P(:,:,27) = -369.8176 -80.3917 40.2024 136.6777 209.0342 223.5054 -198.6964 90.7295 211.3236 307.7989 380.1554 394.6266 -134.7569 154.6690 275.2631 371.7384 444.0949 458.5662 -86.2446 203.1813 323.7754 420.2507 492.6072 507.0785 -51.2747 238.1512 358.7453 455.2206 527.5771 542.0483 -44.4184 245.0074 365.6016 462.0769 534.4333 548.9046 P(:,:,28) = -334.3752 -55.1291 61.2235 154.3055 224.1170 238.0793 -171.6367 107.6094 223.9619 317.0440 386.8555 400.8178 -110.8294 168.4167 284.7693 377.8513 447.6629 461.6252 -64.6936 214.5525 330.9051 423.9871 493.7987 507.7610 -31.4368 247.8094 364.1619 457.2439 527.0555 541.0178 -24.9164 254.3298 370.6823 463.7644 533.5759 547.5382 P(:,:,29) = -301.7398 -31.9817 80.4176 170.3369 237.7765 251.2644 -146.7724 122.9857 235.3850 325.3043 392.7439 406.2318 -88.8687 180.8894 293.2886 383.2080 450.6475 464.1354 -44.9360 224.8221 337.2214 427.1407 494.5803 508.0682 -13.2672 256.4909 368.8901 458.8095 526.2490 539.7369 -7.0582 262.6999 375.0991 465.0185 532.4580 545.9460 P(:,:,30) = -271.6470 -10.7533 97.9524 184.9170 250.1404 263.1851 -123.8980 136.9957 245.7014 332.6660 397.8894 410.9341 -68.6915 192.2022 300.9079 387.8725 453.0959 466.1406 -26.8052 234.0885 342.7942 429.7588 494.9822 508.0269 3.3884 264.2821 372.9878 459.9524 525.1758 538.2205 9.3082 270.2019 378.9076 465.8722 531.0956 544.1403 P(:,:,31) = -243.8616 8.7317 113.9789 198.1767 261.3250 273.9547 -102.8304 149.7629 255.0101 339.2079 402.3562 414.9859 -50.1340 202.4593 307.7065 391.9043 455.0526 467.6823 -10.1521 242.4412 347.6884 431.8862 495.0345 507.6642 18.6687 271.2620 376.5092 460.7069 523.8553 536.4849 24.3193 276.9126 382.1598 466.3576 529.5059 542.1356 P(:,:,32) = -218.1740 26.6307 128.6327 210.2343 271.4355 283.6758 -83.4058 161.3990 263.4010 345.0026 406.2037 418.4440 -33.0496 211.7552 313.7572 395.3588 456.5600 468.8002 5.1567 249.9615 351.9635 433.5651 494.7663 507.0065 32.6976 277.5024 379.5044 461.1060 522.3072 534.5474 38.0974 282.9021 384.9041 466.5057 527.7069 539.9472 P(:,:,33) = -194.3971 43.0851 142.0360 221.1967 280.5673 292.4414 -65.4778 172.0044 270.9553 350.1160 409.4866 421.3607 -17.3070 220.1752 319.1261 398.2868 457.6574 469.5315 19.2411 256.7233 355.6743 434.8350 494.2055 506.0797 45.5868 283.0690 382.0199 461.1806 520.5512 532.4253 50.7521 288.2343 387.1853 466.3460 525.7165 537.5907 P(:,:,34) = -172.3634 58.2216 154.2987 231.1603 288.8066 300.3358 -48.9151 181.6699 277.7470 354.6086 412.2549 423.7841 -2.7885 227.7964 323.8735 400.7352 458.3814 469.9107 32.2086 262.7936 358.8707 435.7323 493.3786 504.9078 57.4362 288.0212 384.0982 460.9599 518.6061 530.1354 62.3824 292.9674 389.0444 465.9061 523.5523 535.0816 P(:,:,35) = -151.9224 72.1546 165.5201 240.2124 296.2317 307.4355 -33.5996 190.4775 283.8429 358.5353 414.5546 425.7584 10.6119 234.6889 328.0544 402.7467 458.7660 469.9698 44.1560 268.2330 361.5985 436.2908 492.3101 503.5140 68.3361 292.4132 385.7786 460.4710 516.4903 527.6941 73.0770 297.1540 390.5195 465.2118 521.2311 532.4349 P(:,:,36) = -132.9394 84.9870 175.7897 248.4318 302.9134 313.8098 -19.4252 198.5012 289.3039 361.9460 416.4276 427.3240 22.9895 240.9159 331.7186 404.3607 458.8423 469.7386 55.1704 273.0968 363.8995 436.5416 491.0232 501.9195 78.3678 296.2942 387.0969 459.7391 514.2207 525.1170 82.9160 300.8424 391.6451 464.2872 518.7688 529.6651 P(:,:,37) = -115.2928 96.8117 185.1885 255.8900 308.9161 319.5213 -6.2963 205.8081 294.1849 364.8864 417.9125 428.5177 34.4303 246.5347 334.9115 405.6130 458.6391 469.2443 65.3304 277.4348 365.8116 436.5131 489.5392 500.1444 87.6046 299.7090 388.0858 458.7873 511.8134 522.4186 91.9717 304.0761 392.4530 463.1545 516.1806 526.7858 P(:,:,38) = -98.8729 107.7125 193.7897 262.6515 314.2979 324.6271 5.8733 212.4587 298.5360 367.3978 419.0441 429.3734 45.0118 251.5972 337.6745 406.5363 458.1826 468.5119 74.7070 281.2924 367.3697 436.2315 487.8778 498.2071 96.1127 302.6981 388.7753 457.6371 509.2835 519.6128 100.3096 306.8949 392.9722 461.8340 513.4803 523.8096 P(:,:,39) = -83.5809 117.7654 201.6597 268.7752 319.1117 329.1790 17.1619 218.5081 302.4024 369.5179 419.8544 429.9218 54.8044 256.1507 340.0450 407.1605 457.4970 467.5643 83.3646 284.7109 368.6052 435.7207 486.0572 496.1246 103.9522 305.2985 389.1928 456.3082 506.6448 516.7121 107.9886 309.3349 393.2292 460.3446 510.6812 520.7485 P(:,:,40) = -69.3270 127.0393 208.8586 274.3141 323.4057 333.2240 27.6399 224.0062 305.8256 371.2810 420.3726 430.1909 63.8716 260.2380 342.0573 407.5128 456.6044 466.4227 91.3614 287.7278 369.5471 435.0026 484.0941 493.9125 111.1773 307.5437 389.3630 454.8185 503.9100 513.7284 115.0625 311.4288 393.2482 458.7036 507.7952 517.6135 P(:,:,41) = -56.0302 135.5966 215.4411 279.3167 327.2234 336.8048 37.3715 228.9983 308.8428 372.7184 420.6252 430.2065 72.2711 263.8980 343.7425 407.6181 455.5248 465.1061 98.7502 290.3770 370.2215 434.0971 482.0038 491.5852 117.8375 309.4644 389.3089 453.1845 501.0912 510.6725 121.5798 313.2067 393.0512 456.9268 504.8335 514.4148 P(:,:,42) = -43.6166 143.4940 221.4568 283.8270 330.6047 339.9602 46.4150 233.5256 311.4884 373.8586 420.6363 429.9918 80.0553 267.1660 345.1288 407.4990 454.2767 463.6322 105.5790 292.6897 370.6524 433.0227 479.8003 489.1559 123.9776 311.0883 389.0511 451.4213 498.1990 507.5545 127.5849 314.6956 392.6583 455.0286 501.8062 511.1618 P(:,:,43) = -32.0193 150.7832 226.9509 287.8851 333.5857 342.7258 54.8232 237.6257 313.7934 374.7275 420.4282 429.5683 87.2719 270.0744 346.2421 407.1763 452.8769 462.0170 111.8914 294.6939 370.8616 431.7958 477.4964 486.6366 129.6383 312.4408 388.6085 449.5427 495.2433 504.3834 133.1178 315.9203 392.0880 453.0222 498.7228 507.8629 P(:,:,44) = -21.1771 157.5112 231.9646 291.5274 336.1994 345.1338 62.6444 241.3326 315.7860 375.3488 420.0208 428.9552 93.9643 272.6525 347.1060 406.6687 451.3408 460.2752 117.7274 296.4156 370.8690 430.4318 475.1038 484.0383 134.8569 313.5451 387.9986 447.5613 492.2334 501.1678 138.2153 316.9036 391.3570 450.9198 495.5918 504.5262 P(:,:,45) = -11.0342 163.7209 236.5355 294.7872 338.4760 347.2137 69.9225 244.6776 317.4923 375.7440 419.4327 428.1705 100.1721 274.9272 347.7418 405.9935 449.6823 458.4200 123.1230 297.8781 370.6928 428.9445 472.6332 481.3710 139.6672 314.4223 387.2369 445.4886 489.1774 497.9151 142.9108 317.6659 390.4806 448.7323 492.4210 501.1588 P(:,:,46) = -1.5400 169.4514 240.6978 297.6950 340.4428 348.9924 76.6979 247.6893 318.9357 375.9328 418.6806 427.2302 105.9315 276.9229 348.1693 405.1664 447.9142 456.4638 128.1116 299.1030 370.3494 427.3466 470.0944 478.6440 144.1001 315.0915 386.3379 443.3350 486.0829 494.6324 147.2348 318.2262 389.4726 446.4698 489.2176 497.7672 P(:,:,47) = 7.3524 174.7387 244.4830 300.2785 342.1251 350.4944 83.0071 250.3935 320.1378 375.9333 417.7799 426.1492 111.2756 278.6619 348.4063 404.2017 446.0483 454.4176 132.7234 300.1098 369.8541 425.6496 467.4962 475.8655 148.1840 315.5704 385.3147 441.1102 482.9568 491.3261 151.2153 318.6017 388.3460 444.1414 485.9880 494.3573 P(:,:,48) = 15.6853 179.6156 247.9198 302.5632 343.5458 351.7423 88.8839 252.8142 321.1184 375.7618 416.7444 424.9409 116.2347 280.1649 348.4691 403.1126 444.0951 452.2916 136.9862 300.9164 369.2207 423.8641 464.8467 473.0432 151.9449 315.8751 384.1794 438.8228 479.8053 488.0019 154.8777 318.8080 387.1122 441.7556 482.7382 490.9347 P(:,:,49) = 23.4981 184.1120 251.0344 304.5724 344.7259 352.7566 94.3592 254.9731 321.8956 375.4336 415.5870 423.6177 120.8365 281.4504 348.3729 401.9109 442.0643 450.0950 140.9254 301.5393 368.4618 421.9998 462.1532 470.1839 155.4064 316.0203 382.9428 436.4808 476.6342 484.6649 158.2456 318.8595 385.7820 439.3199 479.4734 487.5041 P(:,:,50) = 30.8265 188.2556 253.8511 306.3274 345.6847 353.5562 99.4612 256.8903 322.4858 374.9622 414.3195 422.1909 125.1066 282.5358 348.1312 400.6076 439.9649 447.8363 144.5644 301.9935 367.5889 420.0653 459.4226 467.2941 158.5904 316.0195 381.6150 434.0913 473.4486 481.3201 161.3404 318.7695 384.3649 436.8413 476.1986 484.0700 P(:,:,51) = 37.7036 192.0717 256.3918 307.8479 346.4399 354.1583 104.2160 258.5842 322.9043 374.3603 412.9524 420.6708 129.0685 283.4366 347.7567 399.2128 437.8048 445.5232 147.9245 302.2927 366.6127 418.0688 456.6608 464.3793 161.5168 315.8850 380.2051 431.6611 470.2532 477.9716 164.1818 318.5499 382.8700 434.3261 472.9181 480.6365 P(:,:,52) = 44.1597 195.5837 258.6770 309.1517 347.0077 354.5789 108.6476 260.0715 323.1649 373.6395 411.4955 419.0667 132.7435 284.1675 347.2608 397.7355 435.5915 443.1627 151.0256 302.4496 365.5429 416.0175 453.8735 461.4447 164.2041 315.6281 378.7215 429.1961 467.0521 474.6233 166.7880 318.2119 381.3053 431.7799 469.6359 477.2071 P(:,:,53) = 50.2230 198.8130 260.7255 310.2555 347.4030 354.8325 112.7779 261.3679 323.2804 372.8104 409.9579 417.3874 136.1516 284.7416 346.6541 396.1841 433.3316 440.7611 153.8857 302.4757 364.3882 413.9182 451.0657 458.4952 166.6693 315.2593 377.1718 426.7018 463.8493 471.2788 169.1757 317.7657 379.6782 429.2082 466.3557 473.7852 P(:,:,54) = 55.9191 201.7793 262.5543 311.1744 347.6394 354.9324 116.6274 262.4875 323.2626 371.8826 408.3477 415.6407 139.3111 285.1712 345.9463 394.5663 431.0314 438.3244 156.5217 302.3818 363.1569 411.7769 448.2420 455.5350 168.9279 314.7880 375.5631 424.1831 460.6481 467.9412 171.3602 317.2204 377.9954 426.6155 463.0805 470.3735 P(:,:,55) = 61.2720 204.5008 264.1795 311.9224 347.7296 354.8911 120.2148 263.4436 323.1222 370.8652 406.6724 413.8338 142.2388 285.4676 345.1462 392.8892 428.6964 435.8578 158.9488 302.1776 361.8563 409.5992 445.4064 452.5679 170.9942 314.2230 373.9017 421.6446 457.4518 464.6133 173.3559 316.5847 376.2633 424.0063 459.8135 466.9749 P(:,:,56) = 66.3038 206.9945 265.6156 312.5125 347.6852 354.7197 123.5574 264.2481 322.8693 369.7662 404.9388 411.9734 144.9503 285.6410 344.2621 391.1590 426.3317 433.3662 161.1815 301.8722 360.4933 407.3902 442.5629 449.5974 172.8817 313.5724 372.1935 419.0904 454.2631 461.2976 175.1757 315.8664 374.4875 421.3844 456.5571 463.5916 P(:,:,57) = 71.0348 209.2758 266.8762 312.9565 347.5168 354.4288 126.6714 264.9124 322.5128 368.5931 403.1534 410.0654 147.4601 285.7010 343.3014 389.3818 423.9420 430.8541 163.2328 301.4738 359.0742 405.1545 439.7148 446.6268 174.6026 312.8436 370.4440 416.5243 451.0845 457.9966 176.8318 315.0727 372.6731 418.7535 453.3137 460.2258 P(:,:,58) = 75.4840 211.3591 267.9738 313.2655 347.2343 354.0280 129.5715 265.4466 322.0613 367.3530 401.3218 408.1155 149.7814 285.6565 342.2711 387.5628 421.5316 428.3253 165.1150 300.9901 357.6047 402.8964 436.8652 443.6590 176.1682 312.0433 368.6579 413.9496 447.9184 454.7121 178.3353 314.2104 370.8250 416.1167 450.0855 456.8792 P(:,:,59) = 79.6689 213.2578 268.9198 313.4494 346.8467 353.5261 132.2716 265.8604 321.5224 366.0521 399.4493 406.1287 151.9266 285.5154 341.1774 385.7071 419.1043 425.7837 166.8392 300.4281 356.0901 400.6197 434.0169 440.6964 177.5890 311.1778 366.8398 411.3695 444.7667 451.4461 179.6966 313.2854 368.9475 413.4771 446.8743 453.5537 P(:,:,60) = 83.6058 214.9841 269.7251 313.5178 346.3624 352.9313 134.7843 266.1625 320.9035 364.6962 397.5408 404.1097 153.9071 285.2854 340.0263 383.8191 416.6637 423.2326 168.4160 299.7943 354.5352 398.3280 431.1726 437.7415 178.8747 310.2530 364.9939 408.7867 441.6313 448.2002 180.9253 312.3035 367.0445 410.8372 443.6818 450.2507 P(:,:,61) = 87.3098 216.5495 270.3994 313.4793 345.7892 352.2512 137.1215 266.3611 320.2110 363.2909 395.6008 402.0628 155.7336 284.9733 338.8231 381.9030 414.2129 420.6749 169.8550 299.0947 352.9446 396.0245 428.3344 434.7964 180.0344 309.2741 363.1239 406.2038 438.5137 444.9757 182.0302 311.2699 365.1197 408.1996 440.5095 446.9715 P(:,:,62) = 90.7950 217.9646 270.9519 313.3418 345.1342 351.4926 139.2942 266.4638 319.4511 361.8409 393.6333 399.9918 157.4159 284.5855 337.5728 379.9627 411.7551 418.1135 171.1653 298.3349 351.3222 393.7120 425.5044 431.8629 181.0765 308.2460 361.2333 403.6232 435.4156 441.7741 183.0197 310.1892 363.1765 405.5664 437.3588 443.7173 P(:,:,63) = 94.0745 219.2392 271.3912 313.1127 344.4039 350.6622 141.3127 266.4774 318.6294 360.3509 391.6421 397.9004 158.9632 284.1280 336.2800 378.0015 409.2927 415.5510 172.3551 297.5198 349.6718 391.3934 422.6846 428.9428 182.0086 307.1733 359.3253 401.0469 432.3381 438.5963 183.9013 309.0660 361.2180 402.9396 434.2307 440.4890 P(:,:,64) = 97.1604 220.3825 271.7250 312.7991 343.6046 349.7657 143.1864 266.4085 317.7511 358.8251 389.6307 395.7918 160.3841 283.6062 334.9487 376.0228 406.8283 412.9894 173.4323 296.6544 347.9970 389.0710 419.8765 426.0376 182.8380 306.0602 357.4027 398.4768 429.2823 435.4434 184.6822 307.9043 359.2468 400.3209 431.1264 437.2875 P(:,:,65) = 100.0641 221.4030 271.9608 312.4072 342.7419 348.8088 144.9242 266.2632 316.8210 357.2673 387.6021 393.6690 161.6863 283.0252 333.5831 374.0294 404.3641 410.4310 174.4040 295.7429 346.3008 386.7471 417.0818 423.1487 183.5715 304.9104 355.4683 395.9146 426.2493 432.3163 185.3689 306.7078 357.2657 397.7120 428.0467 434.1137 P(:,:,66) = 102.7962 222.3086 272.1054 311.9428 341.8209 347.7966 146.5345 266.0469 315.8437 355.6811 385.5592 391.5349 162.8773 282.3897 332.1865 372.0240 401.9021 407.8777 175.2770 294.7894 344.5862 384.4236 414.3017 420.2774 184.2152 303.7276 353.5244 393.3619 423.2400 429.2156 185.9677 305.4801 355.2769 395.1143 424.9924 430.9681 P(:,:,67) = 105.3666 223.1066 272.1650 311.4116 340.8466 346.7337 148.0248 265.7648 314.8232 354.0698 383.5048 389.3918 163.9640 281.7041 330.7624 370.0091 399.4441 405.3311 176.0575 293.7975 342.8559 382.1025 411.5375 417.4245 184.7750 302.5150 351.5734 390.8201 420.2551 426.1421 186.4842 304.2242 353.2826 392.5292 421.9642 427.8512 P(:,:,68) = 107.7846 223.8040 272.1455 310.8186 339.8235 345.6245 149.4023 265.4218 313.7633 352.4364 381.4413 387.2423 164.9529 280.9723 329.3138 367.9869 396.9918 402.7928 176.7514 292.7708 341.1123 379.7854 408.7903 414.5913 185.2563 301.2757 349.6172 388.2903 417.2952 423.0962 186.9238 302.9432 351.2847 389.9578 418.9627 424.7637 P(:,:,69) = 110.0587 224.4072 272.0524 310.1686 338.7557 344.4731 150.6739 265.0224 312.6676 350.7837 379.3709 385.0883 165.8498 280.1983 327.8435 365.9596 394.5467 400.2642 177.3640 291.7125 339.3577 377.4739 406.0610 411.7784 185.6640 300.0125 347.6577 385.7739 414.3610 420.0784 187.2914 301.6399 349.2851 387.4012 415.9883 421.7058 P(:,:,70) = 112.1971 224.9221 271.8908 309.4658 337.6470 343.2833 151.8457 264.5706 311.5393 349.1143 377.2956 382.9318 166.6604 279.3853 326.3540 363.9290 392.1103 397.7465 177.9006 290.6255 337.5943 375.1692 403.3505 408.9867 186.0031 298.7280 345.6967 383.2717 411.4529 417.0892 187.5917 300.3166 347.2853 384.8603 413.0415 418.6778 P(:,:,71) = 114.2074 225.3543 271.6655 308.7144 336.5011 342.0585 152.9236 264.0704 310.3816 347.4306 375.2173 380.7746 167.3899 278.5367 324.8479 361.8969 389.6836 395.2409 178.3658 289.5126 335.8238 372.8728 400.6595 406.2168 186.2777 297.4246 343.7357 380.7847 408.5714 414.1288 187.8289 298.9758 345.2870 382.3359 410.1226 415.6800 P(:,:,72) = 116.0966 225.7090 271.3808 307.9182 335.3213 340.8019 153.9130 263.5254 309.1972 345.7347 373.1377 378.6184 168.0432 277.6555 323.3273 359.8648 387.2679 392.7485 178.7640 288.3763 334.0482 370.5856 397.9887 403.4693 186.4921 296.1044 341.7762 378.3137 405.7168 411.1974 188.0072 297.6196 343.2914 379.8288 407.2319 412.7126 P(:,:,73) = 117.8714 225.9911 271.0409 307.0808 334.1107 339.5167 154.8193 262.9389 307.9887 344.0286 371.0585 376.4645 168.6248 276.7445 321.7943 357.8342 384.8641 390.2701 179.0994 287.2190 332.2689 368.3088 395.3387 400.7447 186.6500 294.7696 339.8195 375.8593 402.8892 408.2952 188.1304 296.2500 341.2998 377.3397 404.3696 409.7756 P(:,:,74) = 119.5381 226.2051 270.6497 306.2054 332.8721 338.2055 155.6470 262.3140 306.7586 342.3143 368.9811 374.3144 169.1391 275.8062 320.2508 355.8064 382.4732 387.8065 179.3759 286.0429 330.4875 366.0432 392.7100 398.0433 186.7550 293.4221 337.8666 373.4223 400.0891 405.4224 188.2018 294.8688 339.3134 374.8691 401.5358 406.8692 P(:,:,75) = 121.1024 226.3553 270.2107 305.2950 331.6083 336.8709 156.4008 261.6537 305.5091 340.5935 366.9067 372.1693 169.5901 274.8430 318.6984 353.7827 380.0960 385.3586 179.5971 284.8500 328.7054 363.7897 390.1029 395.3656 186.8106 292.0635 335.9189 371.0032 397.3164 402.5791 188.2249 293.4778 337.3332 372.4175 398.7307 403.9934 P(:,:,76) = 122.5699 226.4457 269.7273 304.3526 330.3216 335.5154 157.0849 260.9607 304.2424 338.8676 364.8366 370.0304 169.9815 273.8573 317.1389 351.7642 377.7332 382.9269 179.7663 283.6422 326.9238 361.5491 387.5180 392.7118 186.8197 290.6956 333.9772 368.6025 394.5714 399.7652 188.2026 292.0785 335.3601 369.9854 395.9543 401.1481 P(:,:,77) = 123.9458 226.4801 269.2027 303.3808 329.0144 334.1411 157.7032 260.2376 302.9602 337.1383 362.7719 367.8986 170.3167 272.8511 315.5737 349.7518 375.3854 380.5121 179.8869 282.4212 325.1438 359.3220 384.9555 390.0823 186.7855 289.3198 332.0424 366.2205 391.8541 396.9808 188.1380 290.6724 333.3950 367.5731 393.2067 398.3334 P(:,:,78) = 125.2348 226.4618 268.6398 302.3821 327.6889 332.7502 158.2596 259.4866 301.6645 335.4068 360.7136 365.7749 170.5993 271.8263 314.0042 347.7466 373.0533 378.1147 179.9617 281.1887 323.3666 357.1090 382.4157 387.4771 186.7105 287.9375 330.1155 363.8578 389.1646 394.2259 188.0337 289.2607 331.4387 365.1810 390.4878 395.5491 P(:,:,79) = 126.4417 226.3943 268.0412 301.3588 326.3469 331.3446 158.7574 258.7100 300.3569 333.6745 358.6626 363.6602 170.8321 270.7848 312.4317 345.7492 370.7374 375.7350 179.9935 279.9461 321.5931 354.9106 379.8988 384.8964 186.5975 286.5501 328.1970 361.5146 386.5027 391.5003 187.8922 287.8449 329.4918 362.8093 387.7975 392.7951 P(:,:,80) = 127.5706 226.2806 267.4097 300.3130 324.9905 329.9260 159.1999 257.9099 299.0390 331.9423 356.6198 361.5553 171.0182 269.7282 310.8573 343.7606 368.4381 373.3736 179.9850 278.6950 319.8241 352.7274 377.4049 382.3404 186.4487 285.1586 326.2878 359.1911 383.8685 388.8040 187.7160 286.4259 327.5550 360.4584 385.1358 390.0713 P(:,:,81) = 128.6258 226.1235 266.7475 299.2468 323.6212 328.4961 159.5904 257.0881 297.7122 330.2114 354.5858 359.4607 171.1603 268.6581 309.2821 341.7814 366.1558 371.0307 179.9387 277.4364 318.0605 350.5597 374.9342 379.8090 186.2665 283.7643 324.3883 356.8876 381.2620 386.1369 187.5072 285.0049 325.6290 358.1282 382.5027 387.3775 P(:,:,82) = 129.6108 225.9258 266.0570 298.1620 322.2408 327.0565 159.9316 256.2466 296.3778 328.4828 352.5615 357.3773 171.2610 267.5759 307.7072 339.8122 363.8909 368.7066 179.8568 276.1718 316.3030 348.4080 372.4867 377.3025 186.0531 282.3680 322.4993 354.6043 378.6830 383.4987 187.2679 283.5829 323.7141 355.8191 379.8979 384.7136 P(:,:,83) = 130.5295 225.6900 265.3403 297.0604 320.8506 325.6086 160.2263 255.3869 295.0371 326.7573 350.5474 355.3054 171.3226 266.4831 306.1334 337.8535 361.6437 366.4017 179.7415 274.9021 314.5523 346.2725 370.0626 374.8207 185.8103 280.9709 320.6211 352.3413 376.1314 380.8894 187.0002 282.1607 321.8109 353.5311 377.3213 382.0793 P(:,:,84) = 131.3851 225.4186 264.5992 295.9437 319.4521 324.1537 160.4772 254.5106 293.6913 325.0357 348.5441 353.2458 171.3474 265.3809 304.5615 335.9060 359.4144 364.1161 179.5949 273.6284 312.8090 344.1535 367.6619 372.3636 185.5401 279.5736 318.7542 350.0987 373.6071 378.3087 186.7058 280.7392 319.9198 351.2643 374.7727 379.4744 P(:,:,85) = 132.1809 225.1137 263.8357 294.8133 318.0465 322.6931 160.6866 253.6193 292.3413 323.3189 346.5521 351.1988 171.3377 264.2705 302.9925 333.9701 357.2033 361.8499 179.4189 272.3517 311.0737 342.0513 365.2845 369.9311 185.2443 278.1771 316.8991 347.8767 371.1098 375.7565 186.3864 279.3192 318.0412 349.0188 372.2520 376.8986 P(:,:,86) = 132.9199 224.7775 263.0515 293.6707 316.6351 321.2279 160.8568 252.7143 290.9883 321.6075 344.5719 349.1648 171.2954 263.1529 301.4269 332.0461 355.0105 359.6034 179.2153 271.0729 309.3469 339.9661 362.9305 367.5234 184.9244 276.7820 315.0560 345.6752 368.6396 373.2324 186.0438 277.9013 316.1753 346.7945 369.7589 374.3518 P(:,:,87) = 133.6050 224.4120 262.2482 292.5172 315.2189 319.7593 160.9899 251.7969 289.6331 319.9021 342.6038 347.1442 171.2223 262.0293 299.8655 330.1345 352.8362 357.3766 178.9859 269.7928 307.6290 337.8980 360.5997 365.1401 184.5822 275.3891 313.2253 343.4943 366.1960 370.7364 185.6794 276.4863 314.3225 344.5915 367.2933 371.8336 P(:,:,88) = 134.2389 224.0190 261.4274 291.3541 313.7991 318.2881 161.0882 250.8682 288.2766 318.2033 340.6483 345.1373 171.1204 260.9005 298.3088 328.2355 350.6805 355.1695 178.7321 268.5121 305.9205 335.8472 358.2922 362.7812 184.2189 273.9990 311.4073 341.3340 363.7791 368.2681 185.2947 275.0747 312.4831 342.4098 364.8548 369.3438 P(:,:,89) = 134.8241 223.6003 260.5904 290.1824 312.3765 316.8153 161.1533 249.9295 286.9196 316.5116 338.7057 343.1445 170.9912 259.7674 296.7575 326.3496 348.5436 352.9824 178.4555 267.2316 304.2217 333.8138 356.0078 360.4466 183.8360 272.6122 309.6023 339.1943 361.3884 365.8272 184.8910 273.6671 310.6572 340.2493 362.4433 366.8821 P(:,:,90) = 135.3631 223.1576 259.7386 289.0034 310.9521 315.3418 161.1872 248.9817 285.5628 314.8276 336.7762 341.1659 170.8365 258.6309 295.2120 324.4768 346.4254 350.8152 178.1575 265.9520 302.5330 331.7979 353.7465 358.1362 183.4349 271.2294 307.8104 337.0752 359.0238 363.4136 184.4696 272.2640 308.8451 338.1099 360.0585 364.4483 P(:,:,91) = 135.8581 222.6924 258.8733 287.8181 309.5266 313.8684 161.1917 248.0259 284.2069 313.1516 334.8602 339.2019 170.6576 257.4918 293.6728 322.6175 344.3261 348.6678 177.8395 264.6738 300.8547 329.7995 351.5081 355.8498 183.0166 269.8509 306.0318 334.9766 356.6852 361.0269 184.0317 270.8659 307.0469 335.9916 357.7002 362.0419 P(:,:,92) = 136.3113 222.2062 257.9957 286.6273 308.1010 312.3957 161.1682 247.0630 282.8525 311.4841 332.9578 337.2526 170.4559 256.3508 292.1403 320.7719 342.2456 346.5403 177.5028 263.3976 299.1871 327.8187 349.2924 353.5872 182.5825 268.4773 304.2668 332.8984 354.3721 358.6668 183.5784 269.4732 305.2627 333.8943 355.3680 359.6628 P(:,:,93) = 136.7248 221.7003 257.1068 285.4319 306.6758 310.9246 161.1183 246.0938 281.5003 309.8254 331.0693 335.3181 170.2330 255.2085 290.6149 318.9401 340.1839 344.4327 177.1484 262.1239 297.5304 325.8555 347.0994 351.3482 182.1334 267.1089 302.5154 330.8405 352.0844 356.3332 183.1108 268.0863 303.4927 331.8179 353.0618 357.3105 P(:,:,94) = 137.1005 221.1762 256.2077 284.2329 305.2518 309.4556 161.0436 245.1192 280.1507 308.1759 329.1948 333.3986 169.9899 254.0655 289.0970 317.1222 338.1411 342.3449 176.7776 260.8533 295.8848 323.9100 344.9289 349.1327 181.6706 265.7462 300.7777 328.8029 349.8218 354.0256 182.6299 266.7055 301.7370 329.7622 350.7811 354.9849 P(:,:,95) = 137.4403 220.6349 255.2993 283.0309 303.8295 307.9893 160.9452 244.1399 278.8043 306.5358 327.3345 331.4942 169.7279 252.9225 287.5869 315.3185 336.1171 340.2768 176.3915 259.5861 294.2505 321.9820 342.7807 346.9404 181.1949 264.3895 299.0539 326.7854 347.5841 351.7438 182.1366 265.3312 299.9957 327.7272 348.5259 352.6856 P(:,:,96) = 137.7459 220.0777 254.3827 281.8266 302.4096 306.5262 160.8247 243.1566 277.4615 304.9055 325.4884 329.6050 169.4481 251.7800 286.0849 313.5289 334.1119 338.2285 175.9909 258.3227 292.6277 320.0717 340.6546 344.7712 180.7072 263.0391 297.3440 324.7880 345.3710 349.4875 181.6319 263.9638 298.2687 325.7127 346.2956 350.4122 P(:,:,97) = 138.0189 219.5057 253.4586 280.6209 300.9926 305.0669 160.6831 242.1699 276.1228 303.2851 323.6568 327.7311 169.1516 250.6384 284.5913 311.7536 332.1253 336.1996 175.5768 257.0637 291.0165 318.1788 338.5505 342.6249 180.2084 261.6953 295.6481 322.8104 343.1821 347.2565 181.1165 262.6034 296.5562 323.7185 344.0902 348.1646 P(:,:,98) = 138.2609 218.9199 252.5279 279.4142 299.5789 303.6119 160.5216 241.1806 274.7886 301.6749 321.8396 325.8726 168.8394 249.4984 283.1063 309.9926 330.1574 334.1903 175.1502 255.8092 289.4171 316.3035 336.4682 340.5012 179.6993 260.3583 293.9663 320.8526 341.0173 345.0503 180.5913 261.2503 294.8582 321.7445 341.9093 345.9422 P(:,:,99) = 138.4735 218.3213 251.5912 278.2072 298.1691 302.1615 160.3414 240.1892 273.4591 300.0751 320.0370 324.0294 168.5124 248.3602 281.6301 308.2460 328.2080 332.2004 174.7119 254.5597 287.8296 314.4455 334.4075 338.3999 179.1807 259.0285 292.2984 318.9144 338.8763 342.8687 180.0569 259.9047 293.1746 319.7906 339.7525 343.7449 P(:,:,100) = 138.6580 217.7108 250.6495 277.0004 296.7636 300.7162 160.1435 239.1962 272.1349 298.4858 318.2490 322.2016 168.1715 247.2243 280.1629 306.5139 326.2770 330.2297 174.2626 253.3153 286.2540 312.6049 332.3681 336.3207 178.6533 257.7060 290.6447 316.9956 336.7588 340.7114 179.5141 258.5669 291.5056 317.8565 337.6197 341.5723 P(:,:,101) = 138.8159 217.0893 249.7032 275.7943 295.3626 299.2763 159.9288 238.2022 270.8161 296.9072 316.4756 320.3893 167.8177 246.0911 278.7050 304.7961 324.3645 328.2781 173.8031 252.0765 284.6904 310.7816 330.3499 334.2636 178.1177 256.3911 289.0050 315.0961 334.6645 338.5782 178.9636 257.2370 289.8509 315.9421 335.5104 339.4241 P(:,:,102) = 138.9483 216.4575 248.7530 274.5894 293.9667 297.8422 159.6983 237.2076 269.5031 295.3395 314.7168 318.5923 167.4516 244.9609 277.2564 303.0928 322.4701 326.3456 173.3342 250.8434 283.1390 308.9754 328.3527 332.2281 177.5746 255.0839 287.3794 313.2158 332.5931 336.4686 178.4060 255.9153 288.2108 314.0472 333.4245 337.3000 P(:,:,103) = 139.0565 215.8163 247.7996 273.3862 292.5762 296.4142 159.4530 236.2128 268.1961 293.7827 312.9727 316.8107 167.0741 243.8340 275.8173 301.4039 320.5938 324.4318 172.8565 249.6163 281.5996 307.1862 326.3762 330.2142 177.0247 253.7845 285.7678 311.3544 330.5444 334.3824 177.8419 254.6017 286.5850 312.1716 331.3616 335.1996 P(:,:,104) = 139.1416 215.1664 246.8434 272.1850 291.1913 294.9925 159.1935 235.2183 266.8954 292.2370 311.2432 315.0444 166.6859 242.7108 274.3878 299.7294 318.7356 322.5368 172.3706 248.3954 280.0724 305.4140 324.4202 328.2215 176.4683 252.4932 284.1702 309.5118 328.5180 332.3192 177.2718 253.2966 284.9736 310.3152 329.3214 333.1227 P(:,:,105) = 139.2048 214.5085 245.8851 270.9863 289.8123 293.5775 158.9208 234.2245 265.6011 290.7024 309.5283 313.2935 166.2877 241.5914 272.9680 298.0693 316.8952 320.6604 171.8771 247.1809 278.5574 303.6587 322.4846 326.2498 175.9062 251.2100 282.5865 307.6878 326.5137 330.2789 176.6962 251.9999 283.3765 308.4777 327.3037 331.0689 P(:,:,106) = 139.2471 213.8433 244.9250 269.7904 288.4395 292.1693 158.6356 233.2318 264.3135 289.1789 307.8280 311.5578 165.8801 240.4763 271.5580 296.4235 315.0725 318.8023 171.3767 245.9729 277.0546 301.9200 320.5691 324.2989 175.3388 249.9350 281.0168 305.8822 324.5313 328.2611 176.1157 250.7119 281.7936 306.6590 325.3081 329.0379 P(:,:,107) = 139.2694 213.1713 243.9637 268.5977 287.0731 290.7682 158.3385 232.2404 263.0328 287.6668 306.1422 309.8373 165.4637 239.3656 270.1580 294.7919 313.2674 316.9625 170.8697 244.7716 275.5640 300.1980 318.6734 322.3685 174.7666 248.6685 279.4609 304.0949 322.5703 326.2654 175.5307 249.4325 280.2250 304.8589 323.3344 327.0295 P(:,:,108) = 139.2728 212.4931 243.0016 267.4084 285.7134 289.3745 158.0304 231.2507 261.7591 286.1659 304.4710 308.1320 165.0391 238.2594 268.7679 293.1747 311.4797 315.1407 170.3568 243.5771 274.0856 298.4924 316.7974 320.4584 174.1901 247.4104 277.9188 302.3256 320.6307 324.2917 174.9416 248.1619 278.6704 303.0771 321.3822 325.0432 P(:,:,109) = 139.2582 211.8094 242.0391 266.2228 284.3606 287.9882 157.7117 230.2629 260.4926 284.6764 302.8142 306.4417 164.6069 237.1581 267.3878 291.5715 309.7093 313.3369 169.8385 242.3897 272.6193 296.8031 314.9409 318.5684 173.6096 246.1608 276.3904 300.5742 318.7120 322.3395 174.3489 246.9002 277.1298 301.3136 319.4514 323.0789 P(:,:,110) = 139.2263 211.1206 241.0765 265.0412 283.0148 286.6095 157.3833 229.2775 259.2334 283.1982 301.1717 304.7664 164.1676 236.0618 266.0178 289.9825 307.9561 311.5508 169.3151 241.2093 271.1652 295.1299 313.1035 316.6982 173.0256 244.9198 274.8757 298.8404 316.8140 320.4087 173.7530 245.6473 275.6032 299.5679 317.5415 321.1362 P(:,:,111) = 139.1781 210.4271 240.1142 263.8639 281.6762 285.2386 157.0455 228.2946 257.9817 281.7313 299.5436 303.1061 163.7217 234.9707 264.6578 288.4075 306.2198 309.7822 168.7871 240.0361 269.7232 293.4729 311.2851 314.8476 172.4384 243.6874 273.3745 297.1242 314.9365 318.4989 173.1543 244.4033 274.0904 297.8401 315.6524 319.2148 P(:,:,112) = 139.1143 209.7296 239.1527 262.6911 280.3449 283.8757 156.6991 227.3144 256.7374 280.2759 297.9297 301.4605 163.2696 233.8850 263.3080 286.8465 304.5003 308.0311 168.2549 238.8702 268.2932 291.8317 309.4855 313.0163 171.8484 242.4638 271.8868 295.4253 313.0791 316.6099 172.5530 243.1683 272.5914 296.1298 313.7837 317.3144 P(:,:,113) = 139.0356 209.0284 238.1920 261.5230 279.0212 282.5208 156.3444 226.3372 255.5009 278.8318 296.3300 299.8296 162.8119 232.8047 261.9683 285.2993 302.7975 306.2971 167.7189 237.7117 266.8753 290.2063 307.7044 311.2041 171.2561 241.2488 270.4125 293.7434 311.2416 314.7413 171.9496 241.9424 271.1060 294.4369 311.9351 315.4348 P(:,:,114) = 138.9428 208.3239 237.2327 260.3597 277.7050 281.1741 155.9821 225.3632 254.2720 277.3991 294.7443 298.2134 162.3488 231.7300 260.6388 283.7658 301.1111 304.5802 167.1794 236.5606 265.4694 288.5964 305.9417 309.4108 170.6615 240.0427 268.9515 292.0785 309.4238 312.8929 171.3443 240.7254 269.6342 292.7612 310.1065 313.5756 P(:,:,115) = 138.8365 207.6165 236.2749 259.2016 276.3966 279.8356 155.6125 224.3926 253.0510 275.9777 293.1727 296.6117 161.8809 230.6610 259.3194 282.2461 299.4411 302.8801 166.6369 235.4170 264.0753 287.0020 304.1970 307.6360 170.0652 238.8453 267.5036 290.4303 307.6254 311.0644 170.7374 239.5174 268.1758 291.1025 308.2975 311.7365 P(:,:,116) = 138.7173 206.9067 235.3189 258.0487 275.0960 278.5055 155.2362 223.4256 251.8378 274.5676 291.6149 295.0244 161.4085 229.5979 258.0101 280.7399 297.7872 301.1967 166.0916 234.2809 262.6931 285.4229 302.4703 305.8797 169.4673 237.6567 266.0689 288.7987 305.8460 309.2555 170.1292 238.3185 266.7308 289.4605 306.5079 309.9173 P(:,:,117) = 138.5860 206.1947 234.3649 256.9012 273.8033 277.1838 154.8536 222.4623 250.6325 273.1688 290.0709 293.4514 160.9320 228.5407 256.7109 279.2472 296.1493 299.5298 165.5438 233.1525 261.3228 283.8590 300.7611 304.1416 168.8682 236.4769 264.6472 287.1834 304.0855 307.4660 169.5200 237.1287 265.2989 287.8352 304.7373 308.1178 P(:,:,118) = 138.4430 205.4808 233.4133 255.7592 272.5186 275.8705 154.4650 221.5028 249.4353 271.7812 288.5406 291.8925 160.4517 227.4895 255.4219 277.7678 294.5273 297.8792 164.9939 232.0317 259.9641 282.3100 299.0695 302.4213 168.2681 235.3059 263.2383 285.5842 302.3437 305.6956 168.9100 235.9478 263.8802 286.2262 302.9856 306.3375 P(:,:,119) = 138.2890 204.7655 232.4641 254.6229 271.2420 274.5658 154.0710 220.5475 248.2460 270.4048 287.0239 290.3478 159.9679 226.4444 254.1429 276.3017 292.9209 296.2447 164.4420 230.9185 258.6170 280.7759 297.3950 300.7188 167.6672 234.1437 261.8422 284.0010 300.6202 303.9440 168.2995 234.7760 262.4745 284.6334 301.2525 304.5763 P(:,:,120) = 138.1245 204.0490 231.5175 253.4924 269.9735 273.2697 153.6717 219.5962 247.0648 269.0396 285.5207 288.8169 159.4809 225.4055 252.8740 274.8488 291.3300 294.6262 163.8885 229.8130 257.2816 279.2564 295.7375 299.0338 167.0657 232.9902 260.4588 282.4336 298.9147 302.2110 167.6887 233.6132 261.0817 283.0565 299.5377 302.8339 P(:,:,121) = 137.9499 203.3315 230.5739 252.3677 268.7131 271.9822 153.2677 218.6493 245.8916 267.6855 284.0309 287.3000 158.9911 224.3728 251.6151 273.4090 289.7544 293.0234 163.3337 228.7153 255.9576 277.7515 294.0969 297.3660 166.4640 231.8456 259.0879 280.8818 297.2272 300.4963 167.0777 232.4593 259.7016 281.4955 297.8409 301.1100 P(:,:,122) = 137.7658 202.6134 229.6332 251.2491 267.4610 270.7034 152.8591 217.7067 244.7265 266.3424 282.5543 285.7967 158.4988 223.3463 250.3662 271.9820 288.1939 291.4363 162.7776 227.6252 254.6451 276.2609 292.4728 295.7152 165.8621 230.7097 257.7295 279.3453 295.5572 298.7996 166.4668 231.3144 258.3342 279.9501 296.1620 299.4044 P(:,:,123) = 137.5727 201.8949 228.6958 250.1366 266.2171 269.4332 152.4464 216.7687 243.5696 265.0103 281.0909 284.3070 158.0040 222.3262 249.1272 270.5679 286.6485 289.8646 162.2207 226.5429 253.3438 274.7846 290.8651 294.0812 165.2603 229.5825 256.3834 277.8241 293.9047 297.1208 165.8562 230.1784 256.9793 278.4201 294.5006 297.7167 P(:,:,124) = 137.3709 201.1762 227.7617 249.0302 264.9815 268.1717 152.0299 215.8352 242.4207 263.6892 279.6405 282.8307 157.5073 221.3125 247.8981 269.1665 285.1178 288.3081 161.6630 225.4683 252.0538 273.3223 289.2736 292.4639 164.6587 228.4640 255.0495 276.3179 292.2693 295.4595 165.2460 229.0513 255.6369 276.9053 292.8566 296.0469 P(:,:,125) = 137.1610 200.4576 226.8311 247.9300 263.7542 266.9190 151.6098 214.9064 241.2800 262.3789 278.2030 281.3678 157.0086 220.3052 246.6788 267.7777 283.6018 286.7667 161.1048 224.4014 250.7750 271.8739 287.6980 290.8629 164.0576 227.3542 253.7277 274.8266 290.6507 293.8156 164.6365 227.9331 254.3067 275.4055 291.2297 294.3945 P(:,:,126) = 136.9432 199.7392 225.9042 246.8362 262.5351 265.6749 151.1864 213.9824 240.1474 261.0793 276.7783 279.9181 156.5084 219.3044 245.4693 266.4013 282.1003 285.2401 160.5463 223.3423 249.5072 270.4392 286.1382 289.2780 163.4570 226.2530 252.4179 273.3499 289.0489 292.1887 164.0277 226.8236 252.9886 273.9206 289.6196 292.7594 P(:,:,127) = 136.7181 199.0213 224.9809 245.7487 261.3244 264.4396 150.7600 213.0632 239.0228 259.7906 275.3664 278.4815 156.0068 218.3099 244.2696 265.0373 280.6131 283.7283 159.9876 222.2908 248.2504 269.0182 284.5940 287.7091 162.8572 225.1603 251.1200 271.8877 287.4635 290.5787 163.4198 225.7230 251.6826 272.4503 288.0261 291.1413 P(:,:,128) = 136.4859 198.3040 224.0615 244.6676 260.1221 263.2130 150.3308 212.1489 237.9064 258.5124 273.9669 277.0578 155.5040 217.3220 243.0795 263.6856 279.1401 282.2310 159.4289 221.2470 247.0045 267.6105 283.0651 286.1560 162.2582 224.0763 249.8338 270.4398 285.8944 288.9853 162.8129 224.6310 250.3885 270.9946 286.4491 289.5400 P(:,:,129) = 136.2471 197.5876 223.1461 243.5929 258.9280 261.9951 149.8991 211.2395 236.7981 257.2449 272.5800 275.6470 155.0001 216.3406 241.8991 262.3459 277.6810 280.7481 158.8704 220.2109 245.7694 266.2162 281.5513 284.6183 161.6603 223.0007 248.5593 269.0061 284.3412 287.4082 162.2073 223.5477 249.1063 269.5531 284.8882 287.9552 P(:,:,130) = 136.0020 196.8721 222.2347 242.5248 257.7423 260.7858 149.4650 210.3352 235.6978 255.9878 271.2054 274.2489 154.4955 215.3657 240.7282 261.0183 276.2358 279.2793 158.3122 219.1824 244.5450 264.8350 280.0526 283.0961 161.0635 221.9337 247.2962 267.5863 282.8038 285.8473 161.6029 222.4731 247.8357 268.1257 283.3432 286.3868 P(:,:,131) = 135.7508 196.1579 221.3275 241.4631 256.5649 259.5852 149.0289 209.4359 234.6055 254.7412 269.8429 272.8633 153.9902 214.3972 239.5668 259.7025 274.8043 277.8246 157.7545 218.1615 243.3311 263.4668 278.5686 281.5889 160.4680 220.8750 246.0446 266.1803 281.2820 284.3024 161.0000 221.4070 246.5766 266.7123 281.8140 284.8344 P(:,:,132) = 135.4940 195.4449 220.4244 240.4080 255.3958 258.3933 148.5908 208.5417 233.5212 253.5049 268.4926 271.4901 153.4844 213.4353 238.4149 258.3985 273.3862 276.3838 157.1974 217.1482 242.1278 262.1114 277.0991 280.0967 159.8738 219.8247 244.8042 264.7879 279.7756 282.7731 160.3985 220.3494 245.3290 265.3126 280.3003 283.2979 P(:,:,133) = 135.2317 194.7333 219.5257 239.3595 254.2349 257.2100 148.1511 207.6526 232.4450 252.2788 267.1542 270.1293 152.9784 212.4800 237.2723 257.1061 271.9815 274.9566 156.6409 216.1425 240.9349 260.7687 275.6441 278.6192 159.2811 218.7827 243.5750 263.4089 278.2843 281.2594 159.7987 219.3003 244.0927 263.9265 278.8019 281.7770 P(:,:,134) = 134.9644 194.0234 218.6313 238.3176 253.0823 256.0353 147.7098 206.7688 231.3767 251.0630 265.8277 268.7807 152.4721 211.5311 236.1390 255.8253 270.5900 273.5430 156.0854 215.1443 239.7522 259.4386 274.2033 277.1563 158.6900 217.7490 242.3569 262.0432 276.8079 279.7609 159.2006 218.2596 242.8675 262.5539 277.3186 280.2715 P(:,:,135) = 134.6922 193.3151 217.7413 237.2823 251.9380 254.8692 147.2672 205.8901 230.3163 249.8572 264.5130 267.4441 151.9658 210.5887 235.0149 254.5559 269.2116 272.1427 155.5307 214.1536 238.5798 258.1208 272.7765 275.7077 158.1005 216.7234 241.1496 260.6906 275.3463 278.2775 158.6044 217.2273 241.6535 261.1944 275.8502 278.7813 P(:,:,136) = 134.4155 192.6087 216.8558 236.2536 250.8019 253.7115 146.8234 205.0166 229.2637 248.6615 263.2098 266.1194 151.4596 209.6528 233.9000 253.2977 267.8460 270.7557 154.9772 213.1704 237.4176 256.8153 271.3636 274.2733 157.5128 215.7060 239.9532 259.3509 273.8993 276.8089 158.0100 216.2032 240.4504 259.8481 274.3964 277.3061 P(:,:,137) = 134.1344 191.9041 215.9749 235.2315 249.6739 252.5624 146.3786 204.1483 228.2191 247.4757 261.9181 264.8066 150.9536 208.7234 232.7941 252.0507 266.4932 269.3816 154.4248 212.1946 236.2653 255.5219 269.9643 272.8528 156.9270 214.6967 238.7675 258.0241 272.4665 275.3550 157.4176 215.1873 239.2581 258.5147 272.9571 275.8456 P(:,:,138) = 133.8492 191.2017 215.0985 234.2160 248.5541 251.4217 145.9329 203.2853 227.1822 246.2997 260.6378 263.5054 150.4480 207.8004 231.6973 250.8148 265.1529 268.0205 153.8736 211.2261 235.1230 254.2404 268.5786 271.4462 156.3430 213.6955 237.5923 256.7098 271.0479 273.9156 156.8272 214.1796 238.0765 257.1940 271.5321 274.3997 P(:,:,139) = 133.5602 190.5013 214.2268 233.2071 247.4424 250.2895 145.4865 202.4276 226.1531 245.1335 259.3687 262.2158 149.9428 206.8839 230.6094 249.5897 263.8250 266.6721 153.3238 210.2650 233.9904 252.9708 267.2061 270.0532 155.7611 212.7022 236.4277 255.4080 269.6433 272.4904 156.2389 213.1801 236.9055 255.8859 270.1212 272.9682 P(:,:,140) = 133.2675 189.8032 213.3597 232.2049 246.3388 249.1656 145.0395 201.5752 225.1317 243.9769 258.1108 260.9376 149.4381 205.9738 229.5303 248.3755 262.5095 265.3362 152.7755 209.3111 232.8676 251.7129 265.8468 268.6736 155.1812 211.7168 235.2733 254.1186 268.2525 271.0793 155.6528 212.1885 235.7450 254.5902 268.7241 271.5509 P(:,:,141) = 132.9714 189.1073 212.4973 231.2093 245.2432 248.0500 144.5921 200.7280 224.1180 242.8300 256.8639 259.6707 148.9342 205.0701 228.4601 247.1720 261.2060 264.0128 152.2286 208.3645 231.7545 250.4665 264.5004 267.3072 154.6034 210.7393 234.1293 252.8412 266.8752 269.6820 155.0690 211.2049 234.5949 253.3068 267.3408 270.1476 P(:,:,142) = 132.6721 188.4139 211.6396 230.2202 244.1557 246.9427 144.1444 199.8861 223.1119 241.6925 255.6279 258.4150 148.4310 204.1728 227.3985 245.9791 259.9146 262.7016 151.6833 207.4251 230.6509 249.2315 263.1669 265.9540 154.0278 209.7696 232.9953 251.5759 265.5114 268.2984 154.4874 210.2292 233.4550 252.0356 265.9710 268.7581 P(:,:,143) = 132.3697 187.7229 210.7867 229.2377 243.0760 245.8437 143.6964 199.0496 222.1134 240.5644 254.4027 257.1704 147.9286 203.2818 226.3456 244.7967 258.6349 261.4026 151.1397 206.4929 229.5567 248.0077 261.8460 264.6137 153.4544 208.8076 231.8714 250.3224 264.1607 266.9284 153.9082 209.2614 232.3252 250.7762 264.6145 267.3822 P(:,:,144) = 132.0645 187.0344 209.9385 228.2618 242.0043 244.7528 143.2484 198.2183 221.1224 239.4457 253.1882 255.9367 147.4272 202.3971 225.3013 243.6246 257.3670 260.1155 150.5978 205.5677 228.4718 246.7951 260.5376 263.2861 152.8833 207.8532 230.7573 249.0806 262.8231 265.5716 153.3314 208.3013 231.2054 249.5287 263.2712 266.0197 P(:,:,145) = 131.7566 186.3485 209.0952 227.2925 240.9405 243.6701 142.8003 197.3923 220.1389 238.3362 251.9842 254.7138 146.9268 201.5187 224.2654 242.4627 256.1107 258.8403 150.0577 204.6496 227.3962 245.5936 259.2415 261.9711 152.3145 206.9065 229.6531 247.8504 261.4984 264.2280 152.7570 207.3489 230.0956 248.2929 261.9409 264.6705 P(:,:,146) = 131.4462 185.6653 208.2566 226.3297 239.8845 242.5954 142.3524 196.5715 219.1628 237.2359 250.7906 253.5016 146.4275 200.6466 223.2379 241.3110 254.8657 257.5767 149.5194 203.7385 226.3298 244.4028 257.9576 260.6686 151.7481 205.9673 228.5586 246.6316 260.1864 262.8973 152.1851 206.4042 228.9955 247.0686 260.6233 263.3343 P(:,:,147) = 131.1335 184.9848 207.4229 225.3734 238.8362 241.5288 141.9047 195.7561 218.1941 236.1446 249.6074 252.3000 145.9294 199.7807 222.2188 240.1693 253.6321 256.3247 148.9830 202.8343 225.2724 243.2229 256.6857 259.3783 151.1842 205.0355 227.4736 245.4240 258.8869 261.5794 151.6157 205.4671 227.9052 245.8556 259.3184 262.0110 P(:,:,148) = 130.8186 184.3071 206.5940 224.4235 237.7957 240.4701 141.4573 194.9459 217.2328 235.0623 248.4344 251.1088 145.4325 198.9210 221.2079 239.0375 252.4096 255.0840 148.4485 201.9371 224.2240 242.0535 255.4256 258.1001 150.6226 204.1112 226.3981 244.2276 257.5997 260.2742 151.0489 204.5374 226.8243 244.6539 258.0260 260.7004 P(:,:,149) = 130.5016 183.6322 205.7700 223.4802 236.7628 239.4193 141.0103 194.1409 216.2787 233.9889 247.2715 249.9280 144.9369 198.0675 220.2052 237.9154 251.1981 253.8546 147.9161 201.0467 223.1844 240.8946 254.1773 256.8338 150.0636 203.1942 225.3320 243.0422 256.3248 258.9813 150.4846 203.6153 225.7530 243.4632 256.7459 259.4024 P(:,:,150) = 130.1827 182.9602 204.9507 222.5432 235.7376 238.3764 140.5638 193.3412 215.3318 232.9243 246.1186 248.7575 144.4426 197.2201 219.2107 236.8031 249.9975 252.6364 147.3856 200.1631 222.1536 239.7461 252.9405 255.5794 149.5071 202.2845 224.2751 241.8676 255.0619 257.7008 149.9230 202.7004 224.6910 242.2835 255.4779 258.1167 P(:,:,151) = 129.8621 182.2910 204.1364 221.6127 234.7199 237.3414 140.1178 192.5467 214.3921 231.8684 244.9756 247.5970 143.9498 196.3787 218.2241 235.7004 248.8076 251.4291 146.8573 199.2862 221.1316 238.6079 251.7151 254.3365 148.9531 201.3820 223.2274 240.7037 253.8109 256.4324 149.3640 201.7929 223.6383 241.1146 254.2218 256.8433 P(:,:,152) = 129.5399 181.6248 203.3269 220.6885 233.7098 236.3140 139.6724 191.7574 213.4595 230.8211 243.8423 246.4466 143.4585 195.5434 217.2455 234.6072 247.6284 250.2326 146.3310 198.4160 220.1181 237.4797 250.5009 253.1052 148.4017 200.4867 222.1887 239.5504 252.5716 255.1759 148.8077 200.8926 222.5947 239.9564 252.9776 255.5818 P(:,:,153) = 129.2161 180.9616 202.5222 219.7707 232.7071 235.2943 139.2278 190.9733 212.5339 229.7824 242.7188 245.3060 142.9687 194.7142 216.2748 233.5233 246.4596 249.0469 145.8070 197.5524 219.1131 236.3616 249.2979 251.8852 147.8529 199.5984 221.1590 238.4075 251.3439 253.9312 148.2541 199.9995 221.5602 238.8087 251.7450 254.3323 P(:,:,154) = 128.8910 180.3014 201.7224 218.8592 231.7118 234.2823 138.7839 190.1943 211.6153 228.7522 241.6048 244.1753 142.4804 193.8909 215.3119 232.4487 245.3013 247.8718 145.2851 196.6955 218.1165 235.2533 248.1059 250.6764 147.3068 198.7172 220.1382 237.2750 250.1276 252.6981 147.7031 199.1136 220.5346 237.6714 250.5240 253.0945 P(:,:,155) = 128.5646 179.6442 200.9274 217.9540 230.7239 233.2779 138.3409 189.4205 210.7037 227.7303 240.5002 243.0542 141.9938 193.0735 214.3567 231.3832 244.1531 246.7071 144.7654 195.8450 217.1282 234.1548 246.9247 249.4787 146.7632 197.8429 219.1261 236.1526 248.9226 251.4765 147.1550 198.2346 219.5178 236.5443 249.3143 251.8682 P(:,:,156) = 128.2370 178.9901 200.1373 217.0550 229.7432 232.2809 137.8988 188.6519 209.7990 226.7167 239.4050 241.9427 141.5089 192.2620 213.4091 230.3268 243.0151 245.5528 144.2479 195.0011 216.1482 233.0659 245.7542 248.2918 146.2224 196.9755 218.1227 235.0404 247.7286 250.2663 146.6095 197.3626 218.5098 235.4275 248.1158 250.6534 P(:,:,157) = 127.9084 178.3391 199.3519 216.1622 228.7699 231.2914 137.4576 187.8883 208.9011 225.7114 238.3190 240.8406 141.0256 191.4564 212.4692 229.2794 241.8871 244.4086 143.7328 194.1635 215.1763 231.9866 244.5943 247.1158
VBBV on 5 Apr 2023
you can add another for loop to consider single K for all iterations of T
VBBV on 5 Apr 2023
Edited: VBBV on 5 Apr 2023
%Plot P-V diagram of a given substance using peng robinson equation of state
%Substance properties are defined by the peng robinson constants and
%the critical properties.
clc; clear; close all;
Tc = 765.62; % R
Pc = 550.60; % psi
Vc = 4.0860; % ft3/lbmole
R = 10.732 ; % psi ft3/(lbmole-R)
% peng robinson equation of state constants
% for n-butane
a = 56065.52286;
b = 1.161014;
% alpha
K = [1.250161 1.118431 1.06921 1.031865 1.004945 0.999667];
V = linspace(b*1.2,40*Vc,1000); % vector of volume
% temperature in F
T = [60 180 230 270 300 306];
T = T + 460; % temperature in R
%peng robinson equation
P = zeros(numel(K),numel(T),numel(V));
hold on
for k = 1:length(K)
fprEOSp = @(Tx,Vx)(R*Tx./(Vx-b)-a*K(k)./((Vx.^2+Vx*b)+(Vx*b-b.^2)));
for i= 1:numel(T)
Tx = T(i);
P(k,i,:) = fprEOSp(Tx,V);
Px = reshape(P(k,i,:),1,[]);
plot(V,Px);
end
end
ylim([0 800]);
xlabel('Volume, ft^3');
ylabel('Pressure, psi');

### Categories

Find more on MATLAB in Help Center and File Exchange

R2023a

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!