`Produces a table and plot of max VCG for various TCG values

proj maxvt
me plotend

enter pm

title 100m barge
units m

create hull
ends  0 100
top 10
bot 0
out 15
/

write
quit pm

`display
report

lim area from equ to max or ABS 40 > 3.05
lim area from equ to max or fld > 3.05

variable (len) d, v1, v2, v3, t1=0., t2=0.5, t3=1.0

macro mx1
 tcg = {t%1}
 vcg = 0
 solve maxvcg
 set v%1 = {vcg}
/

macro mx
 trim = 0
 heel = 0
 draft = %1
 solve weight, lcg
 .mx1 (3,1) 1
 set d = %1
 \ { d:3} { v1:3} { v2:3} { v3:3}
 me {d}, {v1}, {v2}, {v3}
/

\  Max VCG at various TCG values
\  Draft { t1:1} { t2:1} { t3:1}
me plotstart "Max VCG vs. Draft at various TCG values"
me plotlabel "Draft", "TCG={t1}", "TCG={t2}", "TCG={t3}"

.mx (5,0.25) 3.0

me plotend

print preview
report off