diff | fsolve | plot |
[...] | := | ? |
Jump to section:
The Otsego Electric Company must run a power line from Power Station 23 in Section 34 to the saw mill under construction in Section 26. The saw mill is 2.3 km north and 5.2 km east of the power station; the intervening private land is 1,300 m wide (east-west). The power line must cross private land. The cost of an easement is $0.75 per meter through the state forest and $2.25 per meter on the private land.Determine the minimum cost route from the power station to the saw mill. Find the optimal size of the wire. Finally, determine the total cost of the wire for the route.
AWG | Cost ($/m) | Area (mcm) |
000 | 0.25 | 167 |
According to the theory of electricity, the resistance of the wire causes a power loss due to energy dissipation in the form of heat. The cost per meter of this loss, C_heat, is inversely proportional to the cross section area of the wire, A. From tests run by Otsego's Engineering Department, the following data are known:
AWG | Cost ($/m) | Area (mcm) |
00 | 0.2385 | 105 |
Let x be the distance on private land and let w be the distance through the
state forest, both measured in meters. Then, the cost of easement equation is
C_ease := $0.75 w + $2.25 x
Before invoking fsolve to find the critical points of C_ease as a function of x, examine a plot of C_ease over the x domain 1300..3600 (why this domain?) to determine the approximate coordinates of the minimum. Giving fsolve a sensible search interval may greatly speed your computations and, in some instances, is necessary.
Once we have found the distance d meters in the diagram above, we can plot the
path of the power line as follows.
We use d := 200 meters. It's not the correct answer! It's just for illustration.
d := 200:
Line := [0, 0], [3.9, 2.3 - d/1000], [5.2, 2.3];
Connect back to the origin to show the straight path, too.
Line := Line, [0, 0]:
Enter the following plot command to see the path.
plot([Line], x=0..6, y=-1..3, title=`Minimum Cost Path`);
Go to: Project Description.
Last update: Nov 12, 1995 (WmCB)