[tValues unknownsValues] = ...
ode45('Galileo1',[0 3.5]',[0 0]')
{Error using <a href="matlab:matlab.internal.language.introspective.errorDocCallback('odearguments', 'C:\Program Files\MATLAB\R2017b\toolbox\matlab\funfun\private\odearguments.m', 93)" style="font-weight:bold">odearguments</a> (<a href="matlab: opentoline('C:\Program Files\MATLAB\R2017b\toolbox\matlab\funfun\private\odearguments.m',93,0)">line 93</a>)
GALILEO1 must return a column vector.

Error in <a href="matlab:matlab.internal.language.introspective.errorDocCallback('ode45', 'C:\Program Files\MATLAB\R2017b\toolbox\matlab\funfun\ode45.m', 115)" style="font-weight:bold">ode45</a> (<a href="matlab: opentoline('C:\Program Files\MATLAB\R2017b\toolbox\matlab\funfun\ode45.m',115,0)">line 115</a>)
    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options,
    varargin);
} 
ode45('Galileo1',[0 3.5]',[0 0]')
[tValues unknownsValues] = ...
ode45('Galileo1',[0 3.5]',[0 0]')

tValues =

         0
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0001
    0.0001
    0.0001
    0.0003
    0.0004
    0.0005
    0.0006
    0.0013
    0.0019
    0.0026
    0.0032
    0.0064
    0.0096
    0.0128
    0.0160
    0.0320
    0.0480
    0.0640
    0.0800
    0.1600
    0.2400
    0.3201
    0.4001
    0.4876
    0.5751
    0.6626
    0.7501
    0.8376
    0.9251
    1.0126
    1.1001
    1.1876
    1.2751
    1.3626
    1.4501
    1.5376
    1.6251
    1.7126
    1.8001
    1.8876
    1.9751
    2.0626
    2.1501
    2.2376
    2.3251
    2.4126
    2.5001
    2.5876
    2.6751
    2.7626
    2.8501
    2.9376
    3.0251
    3.1126
    3.2001
    3.2751
    3.3500
    3.4250
    3.5000


unknownsValues =

         0         0
    0.0000    0.0001
    0.0000    0.0001
    0.0000    0.0002
    0.0000    0.0002
    0.0000    0.0005
    0.0000    0.0007
    0.0000    0.0010
    0.0000    0.0012
    0.0000    0.0025
    0.0000    0.0037
    0.0000    0.0050
    0.0000    0.0062
    0.0000    0.0125
    0.0000    0.0188
    0.0000    0.0251
    0.0001    0.0313
    0.0002    0.0627
    0.0005    0.0941
    0.0008    0.1255
    0.0013    0.1569
    0.0050    0.3139
    0.0113    0.4709
    0.0201    0.6279
    0.0314    0.7849
    0.1256    1.5699
    0.2826    2.3548
    0.5025    3.1398
    0.7851    3.9248
    1.1661    4.7831
    1.6222    5.6415
    2.1533    6.4999
    2.7596    7.3583
    3.4410    8.2166
    4.1976    9.0750
    5.0292    9.9334
    5.9359   10.7918
    6.9177   11.6501
    7.9747   12.5085
    9.1067   13.3669
   10.3139   14.2253
   11.5961   15.0836
   12.9535   15.9420
   14.3860   16.8004
   15.8936   17.6588
   17.4763   18.5171
   19.1341   19.3755
   20.8670   20.2339
   22.6750   21.0923
   24.5582   21.9506
   26.5164   22.8090
   28.5497   23.6674
   30.6582   24.5258
   32.8417   25.3841
   35.1004   26.2425
   37.4342   27.1009
   39.8431   27.9593
   42.3271   28.8176
   44.8862   29.6760
   47.5204   30.5344
   50.2297   31.3928
   52.6111   32.1283
   55.0477   32.8639
   57.5394   33.5994
   60.0863   34.3350

[tValues unknownsValues] = ...
ode45('Galileo1',[0 3.5]',[0 0]');







plot(tvalues,unknownsValues(:,1))
{Undefined function or variable 'tvalues'.
} 
plot(tValues,unknownsValues(:,1))








[tValues unknownsValues] = ...
R=.1

R =

    0.1000

format compact
more off
[tValues unknownsValues] = ...
ode45(@(tVals,unks) Galileo(tVals,unks,R),[0 3.5]',[0 0]');
plot(tValues,unknownsValues(:,1))
R=.01
R =
    0.0100
[tValues unknownsValues] = ...
ode45(@(tVals,unks) Galileo(tVals,unks,R),[0 3.5]',[0 0]');
plot(tValues,unknownsValues(:,1))