exactInt=4.1302e4
exactInt =
41302
intSpline=integral(@(t) spline(tMeasured,Tmeasured,t),0,2)*60*sigma
{Undefined function or variable 'Tmeasured'.
Error in @(t)spline(tMeasured,Tmeasured,t)
Error in integralCalc/iterateScalarValued (line 314)
fx = FUN(t);
Error in integralCalc/vadapt (line 132)
[q,errbnd] = iterateScalarValued(u,tinterval,pathlen);
Error in integralCalc (line 75)
[q,errbnd] = vadapt(@AtoBInvTransform,interval);
Error in integral (line 88)
Q = integralCalc(fun,a,b,opstruct);}
sigma=5.670373e-8
sigma =
5.6704e-08
T0=273.15
T0 =
273.1500
intExact=4.1302e4
intExact =
41302
intSpline=integral(@(t) spline(tMeasured,Tmeasured,t),0,2)
{Undefined function or variable 'Tmeasured'.
Error in @(t)spline(tMeasured,Tmeasured,t)
Error in integralCalc/iterateScalarValued (line 314)
fx = FUN(t);
Error in integralCalc/vadapt (line 132)
[q,errbnd] = iterateScalarValued(u,tinterval,pathlen);
Error in integralCalc (line 75)
[q,errbnd] = vadapt(@AtoBInvTransform,interval);
Error in integral (line 88)
Q = integralCalc(fun,a,b,opstruct);}
intSpline=integral(@(t) spline(tMeasured,TMeasured,t),0,2)
intSpline =
11.8333
intSpline=integral(@(t) spline(tMeasured,TMeasured,t),0,2)*60*sigma
intSpline =
4.0260e-05
intSpline=integral(@(t) (spline(tMeasured,TMeasured,t)+T0)^4,0,2)*60*sigma
{Error using ^
Inputs must be a scalar and a square matrix.
To compute elementwise POWER, use POWER (.^) instead.
Error in @(t)(spline(tMeasured,TMeasured,t)+T0)^4
Error in integralCalc/iterateScalarValued (line 314)
fx = FUN(t);
Error in integralCalc/vadapt (line 132)
[q,errbnd] = iterateScalarValued(u,tinterval,pathlen);
Error in integralCalc (line 75)
[q,errbnd] = vadapt(@AtoBInvTransform,interval);
Error in integral (line 88)
Q = integralCalc(fun,a,b,opstruct);}
intSpline=integral(@(t) (spline(tMeasured,TMeasured,t)+T0).^4,0,2)*60*sigma
intSpline =
4.1307e+04
intExact=4.1302e4
intExact =
41302
intSpline=integral(@(t) (interp1(tMeasured,TMeasured,t)+T0).^4,0,2)*60*sigma
intSpline =
4.1434e+04
intParFit=integral(@(t) (polyval(CoefPar,t)+T0).^4,0,2)*60*sigma
intParFit =
4.1352e+04
intQuartFit=integral(@(t) (polyval(CoefQuart,t)+T0).^4,0,2)*60*sigma
intQuartFit =
4.1297e+04
intQuartFit=quad(@(t) (polyval(CoefQuart,t)+T0).^4,0,2)*60*sigma
intQuartFit =
4.1297e+04
intQuartFit=quad(@(t) (polyval(CoefQuart,t)+T0)^4,0,2)*60*sigma
{Error using ^
Inputs must be a scalar and a square matrix.
To compute elementwise POWER, use POWER (.^) instead.
Error in @(t)(polyval(CoefQuart,t)+T0)^4
Error in quad (line 67)
y = f(x, varargin{:});}
tNoAir=3.5
tNoAir =
3.5000
sNoAir=60;
unknowns0=[0 0]';
r=0.2
r =
0.2000
[tValues, unknownsValues]=ode45(@(t,y) Galileo(t,y,r),[0 3.5],unknowns0);
tValues
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
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.1396
0.7851 3.9244
1.1660 4.7826
1.6220 5.6406
2.1531 6.4984
2.7592 7.3561
3.4404 8.2137
4.1966 9.0710
5.0278 9.9281
5.9340 10.7850
6.9152 11.6417
7.9713 12.4980
9.1024 13.3541
10.3083 14.2098
11.5891 15.0653
12.9447 15.9203
14.3751 16.7750
15.8803 17.6293
17.4602 18.4832
19.1148 19.3366
20.8441 20.1896
22.6480 21.0421
24.5265 21.8941
26.4795 22.7456
28.5069 23.5965
30.6088 24.4469
32.7851 25.2967
35.0357 26.1460
37.3606 26.9946
39.7598 27.8426
42.2331 28.6899
44.7805 29.5366
47.4020 30.3825
50.0974 31.2278
52.4660 31.9515
54.8889 32.6747
57.3659 33.3974
59.8971 34.1194
[tValues, unknownsValues]=ode45(@(t,y) Galileo(t,y,r),[0 3.5],unknowns0);
[tValues, unknownsValues]=ode45(@(t,y) Galileo(t,y,r),[0 3.5],unknowns0);
r=0.1
r =
0.1000
[tValues, unknownsValues]=ode45(@(t,y) Galileo(t,y,r),[0 3.5],unknowns0);
unknownsValues(end,1)
ans =
59.7099
r=0.1
r =
0.1000
r=0.05
r =
0.0500
[tValues, unknownsValues]=ode45(@(t,y) Galileo(t,y,r),[0 3.5],unknowns0);
unknownsValues(end,1)
ans =
59.3411