Cody

Problem 1875. GJam: 2013 China Event: Cannon Angle

Solution 319881

Submitted on 16 Sep 2013 by James
  • Size: 19
  • This is the leading solution.
This solution is locked. To view this solution, you need to provide a solution of the same size or smaller.

Test Suite

Test Status Code Input and Output
1   Pass
tic vd=[249 6326 ]; exp=44.5948575; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.5949

2   Pass
vd=[77 41 ]; exp=1.9429163; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.9429

3   Pass
vd=[28 11 ]; exp=3.9516039; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.9516

4   Pass
vd=[105 21 ]; exp=0.5347917; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5348

5   Pass
vd=[147 51 ]; exp=0.6626633; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6627

6   Pass
vd=[300 3210 ]; exp=10.2293872; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 10.2294

7   Pass
vd=[219 61 ]; exp=0.3570850; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3571

8   Pass
vd=[240 41 ]; exp=0.1998405; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1998

9   Pass
vd=[179 3269 ]; exp=44.5041154; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.5041

10   Pass
vd=[300 3480 ]; exp=11.1338133; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 11.1338

11   Pass
vd=[102 81 ]; exp=2.1878908; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.1879

12   Pass
vd=[121 51 ]; exp=0.9781434; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.9781

13   Pass
vd=[138 31 ]; exp=0.4570257; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4570

14   Pass
vd=[210 4500 ]; exp=45.0000000; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 45

15   Pass
vd=[192 3761 ]; exp=44.4745779; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.4746

16   Pass
vd=[37 41 ]; exp=8.5337735; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 8.5338

17   Pass
vd=[32 41 ]; exp=11.5514232; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 11.5514

18   Pass
vd=[199 81 ]; exp=0.5742839; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5743

19   Pass
vd=[300 3360 ]; exp=10.7304643; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 10.7305

20   Pass
vd=[233 71 ]; exp=0.3671784; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3672

21   Pass
vd=[266 81 ]; exp=0.3214027; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3214

22   Pass
vd=[27 21 ]; exp=8.1988972; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 8.1989

23   Pass
vd=[279 21 ]; exp=0.0757408; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0757

24   Pass
vd=[300 5230 ]; exp=17.3572961; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 17.3573

25   Pass
vd=[279 51 ]; exp=0.1839431; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1839

26   Pass
vd=[218 21 ]; exp=0.1240585; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1241

27   Pass
vd=[300 7690 ]; exp=28.4309587; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 28.4310

28   Pass
vd=[89 31 ]; exp=1.0990234; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.0990

29   Pass
vd=[141 2028 ]; exp=44.2618026; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.2618

30   Pass
vd=[178 81 ]; exp=0.7178094; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7178

31   Pass
vd=[132 71 ]; exp=1.1443126; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.1443

32   Pass
vd=[300 690 ]; exp=2.1544417; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.1544

33   Pass
vd=[300 6970 ]; exp=24.6861041; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 24.6861

34   Pass
vd=[41 1 ]; exp=0.1670142; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1670

35   Pass
vd=[123 21 ]; exp=0.3897097; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3897

36   Pass
vd=[48 61 ]; exp=7.5190650; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 7.5191

37   Pass
vd=[77 81 ]; exp=3.8470546; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.8471

38   Pass
vd=[152 71 ]; exp=0.8628902; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.8629

39   Pass
vd=[267 11 ]; exp=0.0433201; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0433

40   Pass
vd=[164 21 ]; exp=0.2192071; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2192

41   Pass
vd=[300 4310 ]; exp=13.9947943; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 13.9948

42   Pass
vd=[256 71 ]; exp=0.3041622; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3042

43   Pass
vd=[192 41 ]; exp=0.3122545; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3123

44   Pass
vd=[300 7110 ]; exp=25.3662790; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 25.3663

45   Pass
vd=[98 31 ]; exp=0.9063600; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.9064

46   Pass
vd=[286 21 ]; exp=0.0720786; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0721

47   Pass
vd=[44 91 ]; exp=13.7142226; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 13.7142

48   Pass
vd=[300 8360 ]; exp=32.7741901; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 32.7742

49   Pass
vd=[300 5060 ]; exp=16.7170592; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 16.7171

50   Pass
vd=[194 81 ]; exp=0.6042720; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6043

51   Pass
vd=[96 61 ]; exp=1.8595639; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.8596

52   Pass
vd=[190 1 ]; exp=0.0077770; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0078

53   Pass
vd=[237 51 ]; exp=0.2549165; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2549

54   Pass
vd=[52 31 ]; exp=3.2254600; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.2255

55   Pass
vd=[250 71 ]; exp=0.3189378; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3189

56   Pass
vd=[29 61 ]; exp=22.6508544; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 22.6509

57   Pass
vd=[300 4030 ]; exp=13.0143792; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 13.0144

58   Pass
vd=[39 31 ]; exp=5.7607904; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 5.7608

59   Pass
vd=[62 61 ]; exp=4.4733355; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.4733

60   Pass
vd=[246 61 ]; exp=0.2829993; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2830

61   Pass
vd=[86 51 ]; exp=1.9374169; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.9374

62   Pass
vd=[252 91 ]; exp=0.4023216; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4023

63   Pass
vd=[275 11 ]; exp=0.0408363; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0408

64   Pass
vd=[49 81 ]; exp=9.6529702; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 9.6530

65   Pass
vd=[127 11 ]; exp=0.1914728; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1915

66   Pass
vd=[223 61 ]; exp=0.3443890; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3444

67   Pass
vd=[160 21 ]; exp=0.2303047; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2303

68   Pass
vd=[196 31 ]; exp=0.2265546; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2266

69   Pass
vd=[55 91 ]; exp=8.5730698; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 8.5731

70   Pass
vd=[110 1234 ]; exp=44.0395169; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.0395

71   Pass
vd=[274 61 ]; exp=0.2281140; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2281

72   Pass
vd=[300 6800 ]; exp=23.8846446; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 23.8846

73   Pass
vd=[300 3310 ]; exp=10.5630642; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 10.5631

74   Pass
vd=[44 31 ]; exp=4.5141266; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.5141

75   Pass
vd=[295 1 ]; exp=0.0032261; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0032

76   Pass
vd=[64 11 ]; exp=0.7540525; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7541

77   Pass
vd=[240 21 ]; exp=0.1023567; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1024

78   Pass
vd=[91 21 ]; exp=0.7120327; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7120

79   Pass
vd=[300 1870 ]; exp=5.8744286; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 5.8744

80   Pass
vd=[300 1730 ]; exp=5.4290641; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 5.4291

81   Pass
vd=[259 51 ]; exp=0.2134486; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2134

82   Pass
vd=[51 91 ]; exp=10.0258563; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 10.0259

83   Pass
vd=[300 4200 ]; exp=13.6075829; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 13.6076

84   Pass
vd=[262 81 ]; exp=0.3312919; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3313

85   Pass
vd=[74 81 ]; exp=4.1674764; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.1675

86   Pass
vd=[88 81 ]; exp=2.9417237; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.9417

87   Pass
vd=[255 91 ]; exp=0.3929103; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3929

88   Pass
vd=[72 41 ]; exp=2.2226617; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.2227

89   Pass
vd=[300 7900 ]; exp=29.6707716; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 29.6708

90   Pass
vd=[300 1620 ]; exp=5.0800699; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 5.0801

91   Pass
vd=[253 71 ]; exp=0.3114186; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3114

92   Pass
vd=[265 91 ]; exp=0.3638145; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3638

93   Pass
vd=[248 1 ]; exp=0.0045647; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0046

94   Pass
vd=[300 6090 ]; exp=20.7696393; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 20.7696

95   Pass
vd=[158 71 ]; exp=0.7985813; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7986

96   Pass
vd=[300 2170 ]; exp=6.8338053; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 6.8338

97   Pass
vd=[300 1500 ]; exp=4.7002140; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.7002

98   Pass
vd=[163 31 ]; exp=0.3275780; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3276

99   Pass
vd=[139 31 ]; exp=0.4504729; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4505

100   Pass
vd=[271 21 ]; exp=0.0802786; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0803

101   Pass
vd=[113 41 ]; exp=0.9016072; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.9016

102   Pass
vd=[156 81 ]; exp=0.9346125; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.9346

103   Pass
vd=[100 61 ]; exp=1.7135925; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.7136

104   Pass
vd=[85 91 ]; exp=3.5451226; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.5451

105   Pass
vd=[43 61 ]; exp=9.4316060; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 9.4316

106   Pass
vd=[57 61 ]; exp=5.3012744; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 5.3013

107   Pass
vd=[136 61 ]; exp=0.9260755; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.9261

108   Pass
vd=[300 5200 ]; exp=17.2436039; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 17.2436

109   Pass
vd=[159 71 ]; exp=0.7885653; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7886

110   Pass
vd=[86 61 ]; exp=2.3180653; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.3181

111   Pass
vd=[300 8140 ]; exp=31.2094004; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 31.2094

112   Pass
vd=[300 5190 ]; exp=17.2057754; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 17.2058

113   Pass
vd=[259 6845 ]; exp=45.0000000; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 45

114   Pass
vd=[300 6570 ]; exp=22.8379951; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 22.8380

115   Pass
vd=[285 71 ]; exp=0.2454102; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2454

116   Pass
vd=[300 6460 ]; exp=22.3511092; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 22.3511

117   Pass
vd=[75 573 ]; exp=43.3260450; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 43.3260

118   Pass
vd=[46 41 ]; exp=5.4730826; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 5.4731

119   Pass
vd=[52 21 ]; exp=2.1824862; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.1825

120   Pass
vd=[300 5950 ]; exp=20.1913108; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 20.1913

121   Pass
vd=[68 91 ]; exp=5.5599659; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 5.5600

122   Pass
vd=[64 21 ]; exp=1.4399949; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.4400

123   Pass
vd=[36 11 ]; exp=2.3856596; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.3857

124   Pass
vd=[134 21 ]; exp=0.3283507; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3284

125   Pass
vd=[9 1 ]; exp=3.4745532; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.4746

126   Pass
vd=[143 31 ]; exp=0.4256222; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4256

127   Pass
vd=[116 91 ]; exp=1.9000393; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.9000

128   Pass
vd=[23 41 ]; exp=24.7121091; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 24.7121

129   Pass
vd=[300 8390 ]; exp=33.0022766; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 33.0023

130   Pass
vd=[200 51 ]; exp=0.3579647; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3580

131   Pass
vd=[66 91 ]; exp=5.9068209; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 5.9068

132   Pass
vd=[203 21 ]; exp=0.1430697; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1431

133   Pass
vd=[300 230 ]; exp=0.7175455; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7175

134   Pass
vd=[300 7250 ]; exp=26.0668121; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 26.0668

135   Pass
vd=[81 71 ]; exp=3.0438586; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.0439

136   Pass
vd=[184 91 ]; exp=0.7547006; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7547

137   Pass
vd=[273 41 ]; exp=0.1544470; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1544

138   Pass
vd=[224 91 ]; exp=0.5091983; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5092

139   Pass
vd=[25 1 ]; exp=0.4492173; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4492

140   Pass
vd=[243 51 ]; exp=0.2424831; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2425

141   Pass
vd=[300 4920 ]; exp=16.1968290; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 16.1968

142   Pass
vd=[300 840 ]; exp=2.6239945; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.6240

143   Pass
vd=[300 2740 ]; exp=8.6794304; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 8.6794

144   Pass
vd=[49 91 ]; exp=10.9018740; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 10.9019

145   Pass
vd=[244 81 ]; exp=0.3819768; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3820

146   Pass
vd=[197 31 ]; exp=0.2242603; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2243

147   Pass
vd=[186 31 ]; exp=0.2515707; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2516

148   Pass
vd=[278 71 ]; exp=0.2579249; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2579

149   Pass
vd=[258 1 ]; exp=0.0042177; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0042

150   Pass
vd=[300 2990 ]; exp=9.5002916; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 9.5003

151   Pass
vd=[10 1 ]; exp=2.8120066; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.8120

152   Pass
vd=[300 7390 ]; exp=26.7901102; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 26.7901

153   Pass
vd=[300 430 ]; exp=1.3418485; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.3418

154   Pass
vd=[46 215 ]; exp=42.3568299; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 42.3568

155   Pass
vd=[300 6250 ]; exp=21.4435357; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 21.4435

156   Pass
vd=[109 1 ]; exp=0.0236301; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0236

157   Pass
vd=[152 31 ]; exp=0.3767088; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3767

158   Pass
vd=[290 41 ]; exp=0.1368700; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1369

159   Pass
vd=[277 21 ]; exp=0.0768385; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0768

160   Pass
vd=[136 81 ]; exp=1.2298703; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.2299

161   Pass
vd=[300 910 ]; exp=2.8433535; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.8434

162   Pass
vd=[224 21 ]; exp=0.1175014; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1175

163   Pass
vd=[300 2440 ]; exp=7.7039453; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 7.7039

164   Pass
vd=[119 71 ]; exp=1.4081796; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.4082

165   Pass
vd=[20 31 ]; exp=24.7100694; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 24.7101

166   Pass
vd=[300 5120 ]; exp=16.9419243; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 16.9419

167   Pass
vd=[300 1940 ]; exp=6.0976449; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 6.0976

168   Pass
vd=[119 1445 ]; exp=45.0000000; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 45.0000 - 0.0000i

169   Pass
vd=[255 6635 ]; exp=44.7753101; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.7753

170   Pass
vd=[300 280 ]; exp=0.8735777; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.8736

171   Pass
vd=[219 41 ]; exp=0.2400045; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2400

172   Pass
vd=[262 61 ]; exp=0.2494890; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2495

173   Pass
vd=[300 970 ]; exp=3.0315083; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.0315

174   Pass
vd=[300 7720 ]; exp=28.6029416; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 28.6029

175   Pass
vd=[189 31 ]; exp=0.2436475; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2436

176   Pass
vd=[110 51 ]; exp=1.1836603; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.1837

177   Pass
vd=[298 61 ]; exp=0.1928498; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1928

178   Pass
vd=[231 5445 ]; exp=44.9999996; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 45.0000 - 0.0000i

179   Pass
vd=[69 81 ]; exp=4.7988649; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.7989

180   Pass
vd=[107 11 ]; exp=0.2697430; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2697

181   Pass
vd=[114 21 ]; exp=0.4536767; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4537

182   Pass
vd=[47 225 ]; exp=43.2757298; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 43.2757

183   Pass
vd=[17 1 ]; exp=0.9716372; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.9716

184   Pass
vd=[130 81 ]; exp=1.3460985; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.3461

185   Pass
vd=[300 4630 ]; exp=15.1376561; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 15.1377

186   Pass
vd=[121 91 ]; exp=1.7460565; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.7461

187   Pass
vd=[82 41 ]; exp=1.7129065; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.7129

188   Pass
vd=[199 31 ]; exp=0.2197751; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2198

189   Pass
vd=[93 71 ]; exp=2.3071761; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.3072

190   Pass
vd=[300 5540 ]; exp=18.5513095; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 18.5513

191   Pass
vd=[126 41 ]; exp=0.7251166; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7251

192   Pass
vd=[300 4260 ]; exp=13.8184496; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 13.8184

193   Pass
vd=[300 440 ]; exp=1.3730779; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.3731

194   Pass
vd=[300 3130 ]; exp=9.9634881; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 9.9635

195   Pass
vd=[145 61 ]; exp=0.8146500; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.8147

196   Pass
vd=[300 1170 ]; exp=3.6596870; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.6597

197   Pass
vd=[200 4081 ]; exp=44.4955953; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.4956

198   Pass
vd=[271 41 ]; exp=0.1567351; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1567

199   Pass
vd=[300 1490 ]; exp=4.6685979; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.6686

200   Pass
vd=[291 21 ]; exp=0.0696230; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0696

201   Pass
vd=[139 51 ]; exp=0.7411528; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7412

202   Pass
vd=[53 51 ]; exp=5.1245500; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 5.1246

203   Pass
vd=[42 71 ]; exp=11.6156650; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 11.6157

204   Pass
vd=[243 91 ]; exp=0.4326773; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4327

205   Pass
vd=[300 4430 ]; exp=14.4204070; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 14.4204

206   Pass
vd=[210 11 ]; exp=0.0700282; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0700

207   Pass
vd=[112 61 ]; exp=1.3657683; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.3658

208   Pass
vd=[219 31 ]; exp=0.1814659; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1815

209   Pass
vd=[300 4910 ]; exp=16.1599009; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 16.1599

210   Pass
vd=[44 11 ]; exp=1.5959920; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.5960

211   Pass
vd=[192 11 ]; exp=0.0837741; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0838

212   Pass
vd=[276 71 ]; exp=0.2616766; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2617

213   Pass
vd=[75 91 ]; exp=4.5611459; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.5611

214   Pass
vd=[166 21 ]; exp=0.2139567; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2140

215   Pass
vd=[195 41 ]; exp=0.3027202; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3027

216   Pass
vd=[116 81 ]; exp=1.6909858; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.6910

217   Pass
vd=[282 8114 ]; exp=44.6253583; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.6254

218   Pass
vd=[37 11 ]; exp=2.2581760; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.2582

219   Pass
vd=[300 1910 ]; exp=6.0019356; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 6.0019

220   Pass
vd=[182 1 ]; exp=0.0084757; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0085

221   Pass
vd=[176 21 ]; exp=0.1903338; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1903

222   Pass
vd=[120 1 ]; exp=0.0194965; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0195

223   Pass
vd=[300 1320 ]; exp=4.1319681; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.1320

224   Pass
vd=[298 9061 ]; exp=44.6614758; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.6615

225   Pass
vd=[300 7700 ]; exp=28.4881098; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 28.4881

226   Pass
vd=[79 636 ]; exp=43.5312934; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 43.5313

227   Pass
vd=[300 5280 ]; exp=17.5474818; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 17.5475

228   Pass
vd=[227 31 ]; exp=0.1689006; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1689

229   Pass
vd=[275 41 ]; exp=0.1522086; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1522

230   Pass
vd=[171 1 ]; exp=0.0096012; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0096

231   Pass
vd=[229 91 ]; exp=0.4872031; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4872

232   Pass
vd=[200 1 ]; exp=0.0070187; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0070

233   Pass
vd=[147 41 ]; exp=0.5327126; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5327

234   Pass
vd=[169 31 ]; exp=0.3047300; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3047

235   Pass
vd=[133 41 ]; exp=0.6507837; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6508

236   Pass
vd=[300 1150 ]; exp=3.5967945; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.5968

237   Pass
vd=[147 1 ]; exp=0.0129922; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0130

238   Pass
vd=[144 51 ]; exp=0.6905672; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6906

239   Pass
vd=[30 51 ]; exp=16.8668385; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 16.8668

240   Pass
vd=[217 4805 ]; exp=45.0000000; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 45

241   Pass
vd=[131 1 ]; exp=0.0163597; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0164

242   Pass
vd=[58 41 ]; exp=3.4299252; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.4299

243   Pass
vd=[116 61 ]; exp=1.2731382; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.2731

244   Pass
vd=[43 1 ]; exp=0.1518392; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1518

245   Pass
vd=[167 21 ]; exp=0.2114020; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2114

246   Pass
vd=[295 31 ]; exp=0.1000086; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1000

247   Pass
vd=[217 51 ]; exp=0.3040728; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3041

248   Pass
vd=[300 7220 ]; exp=25.9148699; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 25.9149

249   Pass
vd=[25 51 ]; exp=26.5497777; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 26.5498

250   Pass
vd=[300 3620 ]; exp=11.6073520; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 11.6074

251   Pass
vd=[67 31 ]; exp=1.9402732; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.9403

252   Pass
vd=[300 31 ]; exp=0.0967027; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0967

253   Pass
vd=[300 8610 ]; exp=34.8206771; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 34.8207

254   Pass
vd=[81 91 ]; exp=3.9060385; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.9060

255   Pass
vd=[57 31 ]; exp=2.6826598; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.6827

256   Pass
vd=[300 180 ]; exp=0.5615346; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5615

257   Pass
vd=[73 81 ]; exp=4.2832883; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.2833

258   Pass
vd=[201 21 ]; exp=0.1459311; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1459

259   Pass
vd=[68 11 ]; exp=0.6679330; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6679

260   Pass
vd=[159 51 ]; exp=0.5663996; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5664

261   Pass
vd=[204 1 ]; exp=0.0067462; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0067

262   Pass
vd=[189 81 ]; exp=0.6366722; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6367

263   Pass
vd=[300 590 ]; exp=1.8417362; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.8417

264   Pass
vd=[149 61 ]; exp=0.7714868; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7715

265   Pass
vd=[300 8980 ]; exp=38.9553303; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 38.9553

266   Pass
vd=[216 11 ]; exp=0.0661918; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0662

267   Pass
vd=[119 91 ]; exp=1.8053178; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.8053

268   Pass
vd=[125 81 ]; exp=1.4560313; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.4560

269   Pass
vd=[300 6750 ]; exp=23.6536081; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 23.6536

270   Pass
vd=[300 790 ]; exp=2.4674046; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.4674

271   Pass
vd=[300 4520 ]; exp=14.7419107; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 14.7419

272   Pass
vd=[300 6830 ]; exp=24.0242561; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 24.0243

273   Pass
vd=[94 51 ]; exp=1.6213062; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.6213

274   Pass
vd=[75 21 ]; exp=1.0483648; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.0484

275   Pass
vd=[300 5630 ]; exp=18.9049816; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 18.9050

276   Pass
vd=[86 71 ]; exp=2.6991243; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.6991

277   Pass
vd=[43 31 ]; exp=4.7284323; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.7284

278   Pass
vd=[193 31 ]; exp=0.2336526; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2337

279   Pass
vd=[259 31 ]; exp=0.1297425; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1297

280   Pass
vd=[287 21 ]; exp=0.0715772; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0716

281   Pass
vd=[300 6870 ]; exp=24.2115915; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 24.2116

282   Pass
vd=[57 41 ]; exp=3.5519442; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.5519

283   Pass
vd=[239 11 ]; exp=0.0540650; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0541

284   Pass
vd=[192 71 ]; exp=0.5407548; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5408

285   Pass
vd=[300 8470 ]; exp=33.6315148; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 33.6315

286   Pass
vd=[77 61 ]; exp=2.8933815; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.8934

287   Pass
vd=[105 1125 ]; exp=45.0000000; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 45

288   Pass
vd=[300 9183 ]; exp=44.6530547; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.6531

289   Pass
vd=[300 7170 ]; exp=25.6638839; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 25.6639

290   Pass
vd=[300 5560 ]; exp=18.6296154; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 18.6296

291   Pass
vd=[300 980 ]; exp=3.0628800; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.0629

292   Pass
vd=[145 91 ]; exp=1.2154985; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.2155

293   Pass
vd=[300 5820 ]; exp=19.6630507; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 19.6631

294   Pass
vd=[32 61 ]; exp=17.8587444; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 17.8587

295   Pass
vd=[300 5270 ]; exp=17.5093740; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 17.5094

296   Pass
vd=[237 31 ]; exp=0.1549479; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1549

297   Pass
vd=[168 1 ]; exp=0.0099472; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0099

298   Pass
vd=[41 71 ]; exp=12.2256701; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 12.2257

299   Pass
vd=[207 11 ]; exp=0.0720728; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0721

300   Pass
vd=[250 61 ]; exp=0.2740155; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2740

301   Pass
vd=[300 51 ]; exp=0.1590921; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1591

302   Pass
vd=[161 51 ]; exp=0.5524132; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5524

303   Pass
vd=[265 1 ]; exp=0.0039979; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0040

304   Pass
vd=[246 21 ]; exp=0.0974246; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0974

305   Pass
vd=[195 61 ]; exp=0.4503988; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4504

306   Pass
vd=[160 2612 ]; exp=44.6077195; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.6077

307   Pass
vd=[300 6150 ]; exp=21.0206716; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 21.0207

308   Pass
vd=[83 21 ]; exp=0.8559461; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.8559

309   Pass
vd=[160 61 ]; exp=0.6690338; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6690

310   Pass
vd=[156 61 ]; exp=0.7037899; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7038

311   Pass
vd=[162 1 ]; exp=0.0106977; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0107

312   Pass
vd=[166 61 ]; exp=0.6215362; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6215

313   Pass
vd=[107 51 ]; exp=1.2510059; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.2510

314   Pass
vd=[215 41 ]; exp=0.2490182; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2490

315   Pass
vd=[212 91 ]; exp=0.5684822; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5685

316   Pass
vd=[84 61 ]; exp=2.4300259; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.4300

317   Pass
vd=[300 6630 ]; exp=23.1071658; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 23.1072

318   Pass
vd=[119 61 ]; exp=1.2097168; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.2097

319   Pass
vd=[7 1 ]; exp=5.7684795; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 5.7685

320   Pass
vd=[300 6730 ]; exp=23.5617578; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 23.5618

321   Pass
vd=[300 2860 ]; exp=9.0724806; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 9.0725

322   Pass
vd=[274 21 ]; exp=0.0785303; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0785

323   Pass
vd=[300 1510 ]; exp=4.7318359; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.7318

324   Pass
vd=[42 21 ]; exp=3.3498826; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.3499

325   Pass
vd=[300 5220 ]; exp=17.3193640; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 17.3194

326   Pass
vd=[300 8670 ]; exp=35.3730578; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 35.3731

327   Pass
vd=[18 31 ]; exp=34.8306447; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 34.8306

328   Pass
vd=[300 8540 ]; exp=34.2104815; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 34.2105

329   Pass
vd=[256 21 ]; exp=0.0899619; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0900

330   Pass
vd=[91 71 ]; exp=2.4099422; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.4099

331   Pass
vd=[247 51 ]; exp=0.2346928; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2347

332   Pass
vd=[56 21 ]; exp=1.8813698; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.8814

333   Pass
vd=[135 91 ]; exp=1.4023810; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.4024

334   Pass
vd=[300 3930 ]; exp=12.6682378; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 12.6682

335   Pass
vd=[97 41 ]; exp=1.2237457; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.2237

336   Pass
vd=[280 8000 ]; exp=45.0000000; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 45

337   Pass
vd=[280 41 ]; exp=0.1468211; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1468

338   Pass
vd=[147 21 ]; exp=0.2728412; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2728

339   Pass
vd=[23 11 ]; exp=5.8790664; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 5.8791

340   Pass
vd=[139 1 ]; exp=0.0145308; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0145

341   Pass
vd=[78 620 ]; exp=43.5307164; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 43.5307

342   Pass
vd=[250 41 ]; exp=0.1841728; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1842

343   Pass
vd=[69 41 ]; exp=2.4205902; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.4206

344   Pass
vd=[300 2730 ]; exp=8.6467533; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 8.6468

345   Pass
vd=[141 91 ]; exp=1.2854868; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.2855

346   Pass
vd=[300 3800 ]; exp=12.2211854; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 12.2212

347   Pass
vd=[290 11 ]; exp=0.0367211; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0367

348   Pass
vd=[276 21 ]; exp=0.0773963; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0774

349   Pass
vd=[105 11 ]; exp=0.2801172; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2801

350   Pass
vd=[238 81 ]; exp=0.4014801; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4015

351   Pass
vd=[200 91 ]; exp=0.6387576; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6388

352   Pass
vd=[293 41 ]; exp=0.1340815; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1341

353   Pass
vd=[299 9122 ]; exp=44.6851269; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.6851

354   Pass
vd=[21 1 ]; exp=0.6366722; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6367

355   Pass
vd=[36 91 ]; exp=21.7406299; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 21.7406

356   Pass
vd=[270 7438 ]; exp=44.5863296; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.5863

357   Pass
vd=[300 5870 ]; exp=19.8652679; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 19.8653

358   Pass
vd=[125 21 ]; exp=0.3773380; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3773

359   Pass
vd=[132 81 ]; exp=1.3055882; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.3056

360   Pass
vd=[300 4510 ]; exp=14.7060881; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 14.7061

361   Pass
vd=[111 91 ]; exp=2.0753631; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.0754

362   Pass
vd=[176 11 ]; exp=0.0996981; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0997

363   Pass
vd=[76 81 ]; exp=3.9496009; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.9496

364   Pass
vd=[183 51 ]; exp=0.4275657; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4276

365   Pass
vd=[144 41 ]; exp=0.5551429; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5551

366   Pass
vd=[111 11 ]; exp=0.2506519; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2507

367   Pass
vd=[232 21 ]; exp=0.1095376; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1095

368   Pass
vd=[153 2388 ]; exp=44.3197033; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.3197

369   Pass
vd=[264 41 ]; exp=0.1651571; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1652

370   Pass
vd=[203 1 ]; exp=0.0068128; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0068

371   Pass
vd=[206 1 ]; exp=0.0066158; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0066

372   Pass
vd=[190 41 ]; exp=0.3188632; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3189

373   Pass
vd=[67 51 ]; exp=3.1962494; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.1962

374   Pass
vd=[59 21 ]; exp=1.6946784; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.6947

375   Pass
vd=[300 8380 ]; exp=32.9257979; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 32.9258

376   Pass
vd=[108 51 ]; exp=1.2279321; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.2279

377   Pass
vd=[267 7274 ]; exp=44.7042055; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.7042

378   Pass
vd=[298 21 ]; exp=0.0663905; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0664

379   Pass
vd=[300 9070 ]; exp=40.4879078; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 40.4879

380   Pass
vd=[247 31 ]; exp=0.1426554; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1427

381   Pass
vd=[300 1580 ]; exp=4.9533543; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.9534

382   Pass
vd=[300 310 ]; exp=0.9672092; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.9672

383   Pass
vd=[246 91 ]; exp=0.4221877; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4222

384   Pass
vd=[110 81 ]; exp=1.8807470; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.8807

385   Pass
vd=[189 1 ]; exp=0.0078595; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0079

386   Pass
vd=[299 91 ]; exp=0.2857755; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2858

387   Pass
vd=[300 4820 ]; exp=15.8288860; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 15.8289

388   Pass
vd=[184 41 ]; exp=0.3399986; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3400

389   Pass
vd=[173 81 ]; exp=0.7599105; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7599

390   Pass
vd=[201 61 ]; exp=0.4239087; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4239

391   Pass
vd=[300 1820 ]; exp=5.7152103; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 5.7152

392   Pass
vd=[130 41 ]; exp=0.6811720; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6812

393   Pass
vd=[65 71 ]; exp=4.7395084; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.7395

394   Pass
vd=[192 91 ]; exp=0.6931066; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6931

395   Pass
vd=[136 1887 ]; exp=44.4506934; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.4507

396   Pass
vd=[171 51 ]; exp=0.4896861; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4897

397   Pass
vd=[236 91 ]; exp=0.4587274; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4587

398   Pass
vd=[173 31 ]; exp=0.2908008; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2908

399   Pass
vd=[35 41 ]; exp=9.5737143; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 9.5737

400   Pass
vd=[300 9020 ]; exp=39.5832887; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 39.5833

401   Pass
vd=[296 91 ]; exp=0.2915978; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2916

402   Pass
vd=[33 41 ]; exp=10.8258185; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 10.8258

403   Pass
vd=[300 4110 ]; exp=13.2927675; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 13.2928

404   Pass
vd=[277 1 ]; exp=0.0036590; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0037

405   Pass
vd=[254 31 ]; exp=0.1349008; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1349

406   Pass
vd=[82 81 ]; exp=3.3899244; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.3899

407   Pass
vd=[54 21 ]; exp=2.0235395; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.0235

408   Pass
vd=[54 91 ]; exp=8.9040499; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 8.9040

409   Pass
vd=[85 71 ]; exp=2.7632032; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.7632

410   Pass
vd=[14 1 ]; exp=1.4329920; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.4330

411   Pass
vd=[170 81 ]; exp=0.7869742; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7870

412   Pass
vd=[42 51 ]; exp=8.2296248; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 8.2296

413   Pass
vd=[35 125 ]; exp=45.0000000; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 45

414   Pass
vd=[300 3390 ]; exp=10.8310888; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 10.8311

415   Pass
vd=[122 81 ]; exp=1.5285870; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.5286

416   Pass
vd=[74 31 ]; exp=1.5901570; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.5902

417   Pass
vd=[300 8170 ]; exp=31.4129069; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 31.4129

418   Pass
vd=[300 1660 ]; exp=5.2068859; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 5.2069

419   Pass
vd=[134 41 ]; exp=0.6411051; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6411

420   Pass
vd=[184 3454 ]; exp=44.4258154; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.4258

421   Pass
vd=[143 81 ]; exp=1.1123483; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.1123

422   Pass
vd=[67 91 ]; exp=5.7294046; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 5.7294

423   Pass
vd=[208 41 ]; exp=0.2660616; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2661

424   Pass
vd=[169 81 ]; exp=0.7963175; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7963

425   Pass
vd=[213 11 ]; exp=0.0680695; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0681

426   Pass
vd=[152 61 ]; exp=0.7413270; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7413

427   Pass
vd=[259 21 ]; exp=0.0878899; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0879

428   Pass
vd=[174 71 ]; exp=0.6584409; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6584

429   Pass
vd=[242 41 ]; exp=0.1965510; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1966

430   Pass
vd=[265 31 ]; exp=0.1239339; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1239

431   Pass
vd=[300 6930 ]; exp=24.4952155; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 24.4952

432   Pass
vd=[177 1 ]; exp=0.0089613; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0090

433   Pass
vd=[300 8630 ]; exp=35.0015489; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 35.0015

434   Pass
vd=[281 8057 ]; exp=44.7766386; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.7766

435   Pass
vd=[300 4450 ]; exp=14.4916789; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 14.4917

436   Pass
vd=[257 51 ]; exp=0.2167838; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2168

437   Pass
vd=[300 7210 ]; exp=25.8644501; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 25.8645

438   Pass
vd=[300 8110 ]; exp=31.0086245; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 31.0086

439   Pass
vd=[120 31 ]; exp=0.6044357; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6044

440   Pass
vd=[300 4010 ]; exp=12.9449893; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 12.9450

441   Pass
vd=[85 737 ]; exp=44.2615751; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.2616

442   Pass
vd=[287 1 ]; exp=0.0034084; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0034

443   Pass
vd=[300 6900 ]; exp=24.3530039; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 24.3530

444   Pass
vd=[107 81 ]; exp=1.9878552; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.9879

445   Pass
vd=[73 91 ]; exp=4.8168446; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.8168

446   Pass
vd=[300 5830 ]; exp=19.7034001; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 19.7034

447   Pass
vd=[125 51 ]; exp=0.9165221; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.9165

448   Pass
vd=[174 31 ]; exp=0.2874678; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2875

449   Pass
vd=[258 91 ]; exp=0.3838254; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3838

450   Pass
vd=[88 51 ]; exp=1.8502291; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.8502

451   Pass
vd=[300 890 ]; exp=2.7806630; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.7807

452   Pass
vd=[131 71 ]; exp=1.1618593; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.1619

453   Pass
vd=[300 5890 ]; exp=19.9464871; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 19.9465

454   Pass
vd=[275 81 ]; exp=0.3007089; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3007

455   Pass
vd=[201 31 ]; exp=0.2154232; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2154

456   Pass
vd=[103 61 ]; exp=1.6151178; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.6151

457   Pass
vd=[169 11 ]; exp=0.1081282; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1081

458   Pass
vd=[234 81 ]; exp=0.4153242; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4153

459   Pass
vd=[127 41 ]; exp=0.7137400; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7137

460   Pass
vd=[45 61 ]; exp=8.5850615; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 8.5851

461   Pass
vd=[268 71 ]; exp=0.2775328; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2775

462   Pass
vd=[120 81 ]; exp=1.5800158; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.5800

463   Pass
vd=[300 2520 ]; exp=7.9631306; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 7.9631

464   Pass
vd=[39 91 ]; exp=17.9483177; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 17.9483

465   Pass
vd=[113 31 ]; exp=0.6816548; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6817

466   Pass
vd=[152 91 ]; exp=1.1060655; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.1061

467   Pass
vd=[300 1260 ]; exp=3.9429273; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.9429

468   Pass
vd=[300 6980 ]; exp=24.7340581; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 24.7341

469   Pass
vd=[212 61 ]; exp=0.3810567; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3811

470   Pass
vd=[8 1 ]; exp=4.4040343; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.4040

471   Pass
vd=[153 61 ]; exp=0.7316660; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7317

472   Pass
vd=[55 1 ]; exp=0.0928099; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0928

473   Pass
vd=[300 3080 ]; exp=9.7977560; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 9.7978

474   Pass
vd=[108 81 ]; exp=1.9511562; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.9512

475   Pass
vd=[81 41 ]; exp=1.7555143; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.7555

476   Pass
vd=[242 71 ]; exp=0.3403741; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3404

477   Pass
vd=[46 81 ]; exp=11.0165379; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 11.0165

478   Pass
vd=[221 1 ]; exp=0.0057482; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0057

479   Pass
vd=[249 71 ]; exp=0.3215048; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3215

480   Pass
vd=[300 7440 ]; exp=27.0544802; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 27.0545

481   Pass
vd=[45 21 ]; exp=2.9165098; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.9165

482   Pass
vd=[289 21 ]; exp=0.0705899; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0706

483   Pass
vd=[300 7520 ]; exp=27.4846425; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 27.4846

484   Pass
vd=[235 1 ]; exp=0.0050837; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0051

485   Pass
vd=[299 1234 ]; exp=3.8870928; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.8871

486   Pass
vd=[250 1 ]; exp=0.0044920; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0045

487   Pass
vd=[176 51 ]; exp=0.4622558; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4623

488   Pass
vd=[300 8510 ]; exp=33.9588221; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 33.9588

489   Pass
vd=[151 31 ]; exp=0.3817151; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3817

490   Pass
vd=[61 1 ]; exp=0.0754501; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0755

491   Pass
vd=[132 31 ]; exp=0.4995219; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4995

492   Pass
vd=[238 21 ]; exp=0.1040843; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1041

493   Pass
vd=[238 91 ]; exp=0.4510494; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4510

494   Pass
vd=[101 31 ]; exp=0.8533003; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.8533

495   Pass
vd=[252 1 ]; exp=0.0044210; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0044

496   Pass
vd=[101 41 ]; exp=1.1286835; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.1287

497   Pass
vd=[152 51 ]; exp=0.6197772; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6198

498   Pass
vd=[300 8830 ]; exp=37.0236469; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 37.0236

499   Pass
vd=[240 31 ]; exp=0.1510984; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1511

500   Pass
vd=[134 51 ]; exp=0.7975086; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7975

501   Pass
vd=[300 1020 ]; exp=3.1884039; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.1884

502   Pass
vd=[257 91 ]; exp=0.3868183; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3868

503   Pass
vd=[128 11 ]; exp=0.1884927; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1885

504   Pass
vd=[300 6160 ]; exp=21.0627029; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 21.0627

505   Pass
vd=[213 31 ]; exp=0.1918335; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1918

506   Pass
vd=[255 51 ]; exp=0.2201977; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2202

507   Pass
vd=[126 11 ]; exp=0.1945242; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1945

508   Pass
vd=[280 51 ]; exp=0.1826315; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1826

509   Pass
vd=[150 41 ]; exp=0.5116148; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5116

510   Pass
vd=[300 7480 ]; exp=27.2684278; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 27.2684

511   Pass
vd=[300 4140 ]; exp=13.3975114; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 13.3975

512   Pass
vd=[133 51 ]; exp=0.8095495; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.8095

513   Pass
vd=[168 41 ]; exp=0.4078483; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4078

514   Pass
vd=[281 71 ]; exp=0.2524469; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2524

515   Pass
vd=[147 2205 ]; exp=45.0000000; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 45

516   Pass
vd=[113 1 ]; exp=0.0219868; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0220

517   Pass
vd=[242 61 ]; exp=0.2924323; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2924

518   Pass
vd=[56 91 ]; exp=8.2607465; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 8.2607

519   Pass
vd=[260 21 ]; exp=0.0872152; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0872

520   Pass
vd=[22 31 ]; exp=19.4398009; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 19.4398

521   Pass
vd=[165 71 ]; exp=0.7322451; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7322

522   Pass
vd=[300 2190 ]; exp=6.8980297; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 6.8980

523   Pass
vd=[300 7570 ]; exp=27.7582308; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 27.7582

524   Pass
vd=[141 41 ]; exp=0.5790205; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5790

525   Pass
vd=[300 4840 ]; exp=15.9022391; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 15.9022

526   Pass
vd=[142 11 ]; exp=0.1531570; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1532

527   Pass
vd=[300 7750 ]; exp=28.7765421; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 28.7765

528   Pass
vd=[232 81 ]; exp=0.4225163; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4225

529   Pass
vd=[53 21 ]; exp=2.1007558; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.1008

530   Pass
vd=[123 81 ]; exp=1.5038100; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.5038

531   Pass
vd=[210 51 ]; exp=0.3246830; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3247

532   Pass
vd=[273 1 ]; exp=0.0037670; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0038

533   Pass
vd=[50 81 ]; exp=9.2565056; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 9.2565

534   Pass
vd=[298 1 ]; exp=0.0031614; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0032

535   Pass
vd=[166 11 ]; exp=0.1120718; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1121

536   Pass
vd=[300 6260 ]; exp=21.4861399; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 21.4861

537   Pass
vd=[253 81 ]; exp=0.3552825; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3553

538   Pass
vd=[144 91 ]; exp=1.2324494; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.2324

539   Pass
vd=[288 8463 ]; exp=44.6385985; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.6386

540   Pass
vd=[268 7328 ]; exp=44.5316039; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.5316

541   Pass
vd=[297 61 ]; exp=0.1941507; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1942

542   Pass
vd=[62 41 ]; exp=2.9999445; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.9999

543   Pass
vd=[223 5074 ]; exp=44.6458418; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.6458

544   Pass
vd=[248 91 ]; exp=0.4154052; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4154

545   Pass
vd=[106 61 ]; exp=1.5249020; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.5249

546   Pass
vd=[211 41 ]; exp=0.2585494; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2585

547   Pass
vd=[120 21 ]; exp=0.4094400; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4094

548   Pass
vd=[210 1 ]; exp=0.0063662; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0064

549   Pass
vd=[300 4660 ]; exp=15.2461386; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 15.2461

550   Pass
vd=[172 71 ]; exp=0.6738453; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6738

551   Pass
vd=[214 91 ]; exp=0.5579047; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5579

552   Pass
vd=[131 11 ]; exp=0.1799582; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1800

553   Pass
vd=[50 255 ]; exp=44.1896883; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.1897

554   Pass
vd=[227 51 ]; exp=0.2778715; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2779

555   Pass
vd=[300 4330 ]; exp=14.0654935; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 14.0655

556   Pass
vd=[115 1 ]; exp=0.0212287; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0212

557   Pass
vd=[80 91 ]; exp=4.0049368; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.0049

558   Pass
vd=[228 21 ]; exp=0.1134147; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1134

559   Pass
vd=[104 11 ]; exp=0.2855301; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2855

560   Pass
vd=[260 61 ]; exp=0.2533422; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2533

561   Pass
vd=[290 1 ]; exp=0.0033383; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0033

562   Pass
vd=[105 81 ]; exp=2.0644344; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.0644

563   Pass
vd=[57 71 ]; exp=6.1830728; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 6.1831

564   Pass
vd=[43 188 ]; exp=42.5787456; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 42.5787

565   Pass
vd=[214 71 ]; exp=0.4352775; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4353

566   Pass
vd=[284 11 ]; exp=0.0382891; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0383

567   Pass
vd=[267 41 ]; exp=0.1614665; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1615

568   Pass
vd=[270 31 ]; exp=0.1193862; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1194

569   Pass
vd=[131 51 ]; exp=0.8344642; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.8345

570   Pass
vd=[70 71 ]; exp=4.0817971; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.0818

571   Pass
vd=[21 41 ]; exp=32.8296766; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 32.8297

572   Pass
vd=[282 91 ]; exp=0.3212706; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3213

573   Pass
vd=[24 11 ]; exp=5.3933350; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 5.3933

574   Pass
vd=[226 21 ]; exp=0.1154310; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1154

575   Pass
vd=[227 11 ]; exp=0.0599322; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0599

576   Pass
vd=[300 4800 ]; exp=15.7556484; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 15.7556

577   Pass
vd=[134 81 ]; exp=1.2668806; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.2669

578   Pass
vd=[256 61 ]; exp=0.2613212; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2613

579   Pass
vd=[300 8860 ]; exp=37.3715407; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 37.3715

580   Pass
vd=[230 81 ]; exp=0.4298969; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4299

581   Pass
vd=[182 91 ]; exp=0.7713825; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7714

582   Pass
vd=[142 2057 ]; exp=44.3369812; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.3370

583   Pass
vd=[32 71 ]; exp=21.4019868; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 21.4020

584   Pass
vd=[142 61 ]; exp=0.8494455; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.8494

585   Pass
vd=[300 3170 ]; exp=10.0963242; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 10.0963

586   Pass
vd=[191 61 ]; exp=0.4694629; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4695

587   Pass
vd=[190 71 ]; exp=0.5522004; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5522

588   Pass
vd=[205 71 ]; exp=0.4743394; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4743

589   Pass
vd=[248 71 ]; exp=0.3241029; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3241

590   Pass
vd=[166 31 ]; exp=0.3158443; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3158

591   Pass
vd=[145 51 ]; exp=0.6810732; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6811

592   Pass
vd=[300 6050 ]; exp=20.6033628; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 20.6034

593   Pass
vd=[254 71 ]; exp=0.3089712; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3090

594   Pass
vd=[185 71 ]; exp=0.5824565; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5825

595   Pass
vd=[300 3380 ]; exp=10.7975318; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 10.7975

596   Pass
vd=[296 71 ]; exp=0.2275088; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2275

597   Pass
vd=[38 81 ]; exp=16.6740056; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 16.6740

598   Pass
vd=[24 1 ]; exp=0.4874355; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4874

599   Pass
vd=[280 31 ]; exp=0.1110109; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1110

600   Pass
vd=[300 8690 ]; exp=35.5640799; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 35.5641

601   Pass
vd=[300 4160 ]; exp=13.4674481; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 13.4674

602   Pass
vd=[252 61 ]; exp=0.2696832; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2697

603   Pass
vd=[98 71 ]; exp=2.0773304; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.0773

604   Pass
vd=[209 81 ]; exp=0.5206370; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5206

605   Pass
vd=[300 2710 ]; exp=8.5814338; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 8.5814

606   Pass
vd=[57 21 ]; exp=1.8158465; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.8158

607   Pass
vd=[102 51 ]; exp=1.3767520; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.3768

608   Pass
vd=[300 1770 ]; exp=5.5561708; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 5.5562

609   Pass
vd=[300 8160 ]; exp=31.3447605; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 31.3448

610   Pass
vd=[74 71 ]; exp=3.6499694; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.6500

611   Pass
vd=[51 41 ]; exp=4.4432920; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.4433

612   Pass
vd=[297 91 ]; exp=0.2896374; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2896

613   Pass
vd=[48 71 ]; exp=8.7887794; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 8.7888

614   Pass
vd=[140 41 ]; exp=0.5873229; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5873

615   Pass
vd=[284 61 ]; exp=0.2123322; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2123

616   Pass
vd=[229 41 ]; exp=0.2195007; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2195

617   Pass
vd=[294 81 ]; exp=0.2630966; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2631

618   Pass
vd=[138 71 ]; exp=1.0469250; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.0469

619   Pass
vd=[84 11 ]; exp=0.4376931; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4377

620   Pass
vd=[220 11 ]; exp=0.0638067; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0638

621   Pass
vd=[300 650 ]; exp=2.0293307; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.0293

622   Pass
vd=[261 31 ]; exp=0.1277617; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1278

623   Pass
vd=[259 11 ]; exp=0.0460375; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0460

624   Pass
vd=[150 21 ]; exp=0.2620364; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2620

625   Pass
vd=[300 3580 ]; exp=11.4717190; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 11.4717

626   Pass
vd=[295 71 ]; exp=0.2290539; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2291

627   Pass
vd=[300 8330 ]; exp=32.5500293; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 32.5500

628   Pass
vd=[127 21 ]; exp=0.3655463; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3655

629   Pass
vd=[196 91 ]; exp=0.6651000; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.6651

630   Pass
vd=[300 7230 ]; exp=25.9654030; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 25.9654

631   Pass
vd=[55 51 ]; exp=4.7550987; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.7551

632   Pass
vd=[118 61 ]; exp=1.2303199; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.2303

633   Pass
vd=[172 61 ]; exp=0.5789235; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5789

634   Pass
vd=[300 3990 ]; exp=12.8756810; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 12.8757

635   Pass
vd=[49 41 ]; exp=4.8168001; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.8168

636   Pass
vd=[300 8400 ]; exp=33.0792168; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 33.0792

637   Pass
vd=[272 51 ]; exp=0.1935327; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1935

638   Pass
vd=[246 51 ]; exp=0.2366048; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2366

639   Pass
vd=[209 71 ]; exp=0.4563550; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4564

640   Pass
vd=[113 51 ]; exp=1.1216128; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.1216

641   Pass
vd=[58 11 ]; exp=0.9181841; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.9182

642   Pass
vd=[90 21 ]; exp=0.7279469; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7279

643   Pass
vd=[300 5580 ]; exp=18.7080845; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 18.7081

644   Pass
vd=[26 1 ]; exp=0.4153242; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4153

645   Pass
vd=[300 6020 ]; exp=20.4792082; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 20.4792

646   Pass
vd=[153 71 ]; exp=0.8516442; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.8516

647   Pass
vd=[231 71 ]; exp=0.3735643; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3736

648   Pass
vd=[296 41 ]; exp=0.1313774; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1314

649   Pass
vd=[218 51 ]; exp=0.3012894; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3013

650   Pass
vd=[156 71 ]; exp=0.8191945; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.8192

651   Pass
vd=[244 91 ]; exp=0.4291377; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.4291

652   Pass
vd=[97 61 ]; exp=1.8213680; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.8214

653   Pass
vd=[33 21 ]; exp=5.4466530; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 5.4467

654   Pass
vd=[48 21 ]; exp=2.5623281; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.5623

655   Pass
vd=[300 5180 ]; exp=17.1679811; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 17.1680

656   Pass
vd=[300 3630 ]; exp=11.6413032; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 11.6413

657   Pass
vd=[164 81 ]; exp=0.8456275; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.8456

658   Pass
vd=[180 3306 ]; exp=44.7534375; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.7534

659   Pass
vd=[49 61 ]; exp=7.2085690; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 7.2086

660   Pass
vd=[89 21 ]; exp=0.7444009; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7444

661   Pass
vd=[300 2110 ]; exp=6.6413404; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 6.6413

662   Pass
vd=[199 4040 ]; exp=44.3892211; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.3892

663   Pass
vd=[297 31 ]; exp=0.0986662; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0987

664   Pass
vd=[300 5920 ]; exp=20.0686784; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 20.0687

665   Pass
vd=[300 7180 ]; exp=25.7138602; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 25.7139

666   Pass
vd=[300 5210 ]; exp=17.2814667; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 17.2815

667   Pass
vd=[212 11 ]; exp=0.0687132; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0687

668   Pass
vd=[246 1 ]; exp=0.0046393; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0046

669   Pass
vd=[250 11 ]; exp=0.0494119; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0494

670   Pass
vd=[47 1 ]; exp=0.1270938; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1271

671   Pass
vd=[112 41 ]; exp=0.9177847; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.9178

672   Pass
vd=[300 5110 ]; exp=16.9043649; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 16.9044

673   Pass
vd=[300 6410 ]; exp=22.1324926; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 22.1325

674   Pass
vd=[300 7030 ]; exp=24.9752527; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 24.9753

675   Pass
vd=[300 7970 ]; exp=30.1045540; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 30.1046

676   Pass
vd=[255 1 ]; exp=0.0043176; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0043

677   Pass
vd=[222 5028 ]; exp=44.4345459; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 44.4345

678   Pass
vd=[247 61 ]; exp=0.2807124; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2807

679   Pass
vd=[40 163 ]; exp=43.3665978; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 43.3666

680   Pass
vd=[88 21 ]; exp=0.7614192; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.7614

681   Pass
vd=[91 51 ]; exp=1.7300957; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.7301

682   Pass
vd=[145 21 ]; exp=0.2804200; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2804

683   Pass
vd=[170 61 ]; exp=0.5926273; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.5926

684   Pass
vd=[95 11 ]; exp=0.3421957; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.3422

685   Pass
vd=[137 1 ]; exp=0.0149581; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0150

686   Pass
vd=[220 51 ]; exp=0.2958362; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2958

687   Pass
vd=[132 51 ]; exp=0.8218652; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.8219

688   Pass
vd=[300 6850 ]; exp=24.1177523; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 24.1178

689   Pass
vd=[271 61 ]; exp=0.2331926; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2332

690   Pass
vd=[184 31 ]; exp=0.2570695; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2571

691   Pass
vd=[300 3250 ]; exp=10.3626812; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 10.3627

692   Pass
vd=[11 1 ]; exp=2.3227865; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 2.3228

693   Pass
vd=[148 81 ]; exp=1.0384257; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.0384

694   Pass
vd=[300 5400 ]; exp=18.0076051; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 18.0076

695   Pass
vd=[287 31 ]; exp=0.1056617; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1057

696   Pass
vd=[183 11 ]; exp=0.0922168; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0922

697   Pass
vd=[278 81 ]; exp=0.2942536; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2943

698   Pass
vd=[291 1 ]; exp=0.0033154; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0033

699   Pass
vd=[64 71 ]; exp=4.8902189; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 4.8902

700   Pass
vd=[290 61 ]; exp=0.2036368; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2036

701   Pass
vd=[124 51 ]; exp=0.9313695; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.9314

702   Pass
vd=[230 21 ]; exp=0.1114509; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1115

703   Pass
vd=[109 51 ]; exp=1.2054913; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.2055

704   Pass
vd=[84 91 ]; exp=3.6304847; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 3.6305

705   Pass
vd=[289 81 ]; exp=0.2722793; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.2723

706   Pass
vd=[138 1 ]; exp=0.0147421; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.0147

707   Pass
vd=[130 91 ]; exp=1.5124296; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 1.5124

708   Pass
vd=[182 21 ]; exp=0.1779910; theta=Cannon(vd(1),vd(2)); assert(abs(theta-exp)<2e-6)

ans = 0.1780

709   Pass
vd=[300 8780 ]; exp=36.4745305

exp = 36.4745