% Modify expression to add input arguments.
% Example:
% a = [1 2 3; 4 5 6];
% foo(a);
pmi
k
k =
0.4000
k^2/100
ans =
0.0016
% Modify expression to add input arguments.
% Example:
% a = [1 2 3; 4 5 6];
% foo(a);
pmi
% Modify expression to add input arguments.
% Example:
% a = [1 2 3; 4 5 6];
% foo(a);
pmi
% Modify expression to add input arguments.
% Example:
% a = [1 2 3; 4 5 6];
% foo(a);
pmi
y(1,:).^2+y(2,:).^2+y(3,:).^2
ans =
Columns 1 through 4
3.0000 3.0000 3.0000 3.0000
Columns 5 through 8
3.0000 3.0000 3.0000 3.0000
Columns 9 through 12
3.0000 3.0000 3.0000 3.0000
Columns 13 through 16
3.0000 3.0000 3.0000 3.0000
Columns 17 through 20
3.0000 3.0000 3.0000 3.0000
Columns 21 through 24
3.0000 3.0000 3.0000 3.0000
Columns 25 through 28
3.0000 3.0000 3.0000 3.0000
Columns 29 through 32
3.0000 3.0000 3.0000 3.0000
Columns 33 through 36
3.0000 3.0000 3.0000 3.0000
Columns 37 through 40
3.0000 3.0000 3.0000 3.0000
Columns 41 through 44
3.0000 3.0000 3.0000 3.0000
Columns 45 through 48
3.0000 3.0000 3.0000 3.0000
Columns 49 through 52
3.0000 3.0000 3.0000 3.0000
Columns 53 through 56
3.0000 3.0000 3.0000 3.0000
Columns 57 through 60
3.0000 3.0000 3.0000 3.0000
Columns 61 through 64
3.0000 3.0000 3.0000 3.0000
Columns 65 through 68
3.0000 3.0000 3.0000 3.0000
Columns 69 through 72
3.0000 3.0000 3.0000 3.0000
Columns 73 through 76
3.0000 3.0000 3.0000 3.0000
Columns 77 through 80
3.0000 3.0000 3.0000 3.0000
Columns 81 through 84
3.0000 3.0000 3.0000 3.0000
Columns 85 through 88
3.0000 3.0000 3.0000 3.0000
Columns 89 through 92
3.0000 3.0000 3.0000 3.0000
Columns 93 through 96
3.0000 3.0000 3.0000 3.0000
Columns 97 through 100
3.0000 3.0000 3.0000 3.0000
Columns 101 through 104
3.0000 3.0000 3.0000 3.0000
Columns 105 through 108
3.0000 3.0000 3.0000 3.0000
Columns 109 through 112
3.0000 3.0000 3.0000 3.0000
Columns 113 through 116
3.0000 3.0000 3.0000 3.0000
Columns 117 through 120
3.0000 3.0000 3.0000 3.0000
Columns 121 through 124
3.0000 3.0000 3.0000 3.0000
Columns 125 through 128
3.0000 3.0000 3.0000 3.0000
Columns 129 through 132
3.0000 3.0000 3.0000 3.0000
Columns 133 through 136
3.0000 3.0000 3.0000 3.0000
Columns 137 through 140
3.0000 3.0000 3.0000 3.0000
Columns 141 through 144
3.0000 3.0000 3.0000 3.0000
Columns 145 through 148
3.0000 3.0000 3.0000 3.0000
Columns 149 through 152
3.0000 3.0000 3.0000 3.0000
Columns 153 through 156
3.0000 3.0000 3.0000 3.0000
Columns 157 through 160
3.0000 3.0000 3.0000 3.0000
Columns 161 through 164
3.0000 3.0000 3.0000 3.0000
Columns 165 through 168
3.0000 3.0000 3.0000 3.0000
Columns 169 through 172
3.0000 3.0000 3.0000 3.0000
Columns 173 through 176
3.0000 3.0000 3.0000 3.0000
Columns 177 through 180
3.0000 3.0000 3.0000 3.0000
Columns 181 through 184
3.0000 3.0000 3.0000 3.0000
Columns 185 through 188
3.0000 3.0000 3.0000 3.0000
Columns 189 through 192
3.0000 3.0000 3.0000 3.0000
Columns 193 through 196
3.0000 3.0000 3.0000 3.0000
Columns 197 through 200
3.0000 3.0000 3.0000 3.0000
Columns 201 through 204
3.0000 3.0000 3.0000 3.0000
Columns 205 through 208
3.0000 3.0000 3.0000 3.0000
Columns 209 through 212
3.0000 3.0000 3.0000 3.0000
Columns 213 through 216
3.0000 3.0000 3.0000 3.0000
Columns 217 through 220
3.0000 3.0000 3.0000 3.0000
Columns 221 through 224
3.0000 3.0000 3.0000 3.0000
Columns 225 through 228
3.0000 3.0000 3.0000 3.0000
Columns 229 through 232
3.0000 3.0000 3.0000 3.0000
Columns 233 through 236
3.0000 3.0000 3.0000 3.0000
Columns 237 through 240
3.0000 3.0000 3.0000 3.0000
Columns 241 through 244
3.0000 3.0000 3.0000 3.0000
Columns 245 through 248
3.0000 3.0000 3.0000 3.0000
Columns 249 through 252
3.0000 3.0000 3.0000 3.0000
Columns 253 through 256
3.0000 3.0000 3.0000 3.0000
Columns 257 through 260
3.0000 3.0000 3.0000 3.0000
Columns 261 through 264
3.0000 3.0000 3.0000 3.0000
Columns 265 through 268
3.0000 3.0000 3.0000 3.0000
Columns 269 through 272
3.0000 3.0000 3.0000 3.0000
Columns 273 through 276
3.0000 3.0000 3.0000 3.0000
Columns 277 through 280
3.0000 3.0000 3.0000 3.0000
Columns 281 through 284
3.0000 3.0000 3.0000 3.0000
Columns 285 through 288
3.0000 3.0000 3.0000 3.0000
Columns 289 through 292
3.0000 3.0000 3.0000 3.0000
Columns 293 through 296
3.0000 3.0000 3.0000 3.0000
Columns 297 through 300
3.0000 3.0000 3.0000 3.0000
Column 301
3.0000
format long e
y(1,:).^2+y(2,:).^2+y(3,:).^2
ans =
Column 1
3.000000000000000e+00
Column 2
3.000000000000002e+00
Column 3
2.999999999999998e+00
Column 4
2.999999999999983e+00
Column 5
2.999999999999982e+00
Column 6
2.999999999999983e+00
Column 7
2.999999999999965e+00
Column 8
2.999999999999964e+00
Column 9
2.999999999999964e+00
Column 10
2.999999999999964e+00
Column 11
2.999999999999964e+00
Column 12
2.999999999999964e+00
Column 13
2.999999999999965e+00
Column 14
2.999999999999964e+00
Column 15
2.999999999999964e+00
Column 16
2.999999999999964e+00
Column 17
2.999999999999964e+00
Column 18
2.999999999999964e+00
Column 19
2.999999999999963e+00
Column 20
2.999999999999963e+00
Column 21
2.999999999999963e+00
Column 22
2.999999999999963e+00
Column 23
2.999999999999963e+00
Column 24
2.999999999999963e+00
Column 25
2.999999999999963e+00
Column 26
2.999999999999963e+00
Column 27
2.999999999999963e+00
Column 28
2.999999999999964e+00
Column 29
2.999999999999963e+00
Column 30
2.999999999999963e+00
Column 31
2.999999999999963e+00
Column 32
2.999999999999962e+00
Column 33
2.999999999999963e+00
Column 34
2.999999999999963e+00
Column 35
2.999999999999963e+00
Column 36
2.999999999999963e+00
Column 37
2.999999999999963e+00
Column 38
2.999999999999962e+00
Column 39
2.999999999999963e+00
Column 40
2.999999999999962e+00
Column 41
2.999999999999962e+00
Column 42
2.999999999999963e+00
Column 43
2.999999999999963e+00
Column 44
2.999999999999963e+00
Column 45
2.999999999999963e+00
Column 46
2.999999999999931e+00
Column 47
2.999999999999916e+00
Column 48
2.999999999999908e+00
Column 49
2.999999999999904e+00
Column 50
2.999999999999901e+00
Column 51
2.999999999999901e+00
Column 52
2.999999999999901e+00
Column 53
2.999999999999900e+00
Column 54
2.999999999999901e+00
Column 55
2.999999999999901e+00
Column 56
2.999999999999900e+00
Column 57
2.999999999999900e+00
Column 58
2.999999999999900e+00
Column 59
2.999999999999900e+00
Column 60
2.999999999999900e+00
Column 61
2.999999999999901e+00
Column 62
2.999999999999901e+00
Column 63
2.999999999999901e+00
Column 64
2.999999999999901e+00
Column 65
2.999999999999902e+00
Column 66
2.999999999999903e+00
Column 67
2.999999999999905e+00
Column 68
2.999999999999910e+00
Column 69
2.999999999999919e+00
Column 70
2.999999999999939e+00
Column 71
2.999999999999938e+00
Column 72
2.999999999999939e+00
Column 73
2.999999999999938e+00
Column 74
2.999999999999938e+00
Column 75
2.999999999999938e+00
Column 76
2.999999999999938e+00
Column 77
2.999999999999938e+00
Column 78
2.999999999999938e+00
Column 79
2.999999999999938e+00
Column 80
2.999999999999939e+00
Column 81
2.999999999999938e+00
Column 82
2.999999999999939e+00
Column 83
2.999999999999938e+00
Column 84
2.999999999999938e+00
Column 85
2.999999999999937e+00
Column 86
2.999999999999938e+00
Column 87
2.999999999999937e+00
Column 88
2.999999999999937e+00
Column 89
2.999999999999936e+00
Column 90
2.999999999999938e+00
Column 91
2.999999999999937e+00
Column 92
2.999999999999937e+00
Column 93
2.999999999999938e+00
Column 94
2.999999999999938e+00
Column 95
2.999999999999938e+00
Column 96
2.999999999999938e+00
Column 97
2.999999999999938e+00
Column 98
2.999999999999939e+00
Column 99
2.999999999999939e+00
Column 100
2.999999999999938e+00
Column 101
2.999999999999938e+00
Column 102
2.999999999999938e+00
Column 103
2.999999999999937e+00
Column 104
2.999999999999938e+00
Column 105
2.999999999999938e+00
Column 106
2.999999999999938e+00
Column 107
2.999999999999939e+00
Column 108
2.999999999999940e+00
Column 109
2.999999999999955e+00
Column 110
2.999999999999956e+00
Column 111
2.999999999999956e+00
Column 112
2.999999999999926e+00
Column 113
2.999999999999925e+00
Column 114
2.999999999999925e+00
Column 115
2.999999999999925e+00
Column 116
2.999999999999925e+00
Column 117
2.999999999999914e+00
Column 118
2.999999999999914e+00
Column 119
2.999999999999914e+00
Column 120
2.999999999999914e+00
Column 121
2.999999999999914e+00
Column 122
2.999999999999914e+00
Column 123
2.999999999999914e+00
Column 124
2.999999999999914e+00
Column 125
2.999999999999914e+00
Column 126
2.999999999999914e+00
Column 127
2.999999999999914e+00
Column 128
2.999999999999914e+00
Column 129
2.999999999999914e+00
Column 130
2.999999999999914e+00
Column 131
2.999999999999915e+00
Column 132
2.999999999999915e+00
Column 133
2.999999999999915e+00
Column 134
2.999999999999915e+00
Column 135
2.999999999999915e+00
Column 136
2.999999999999915e+00
Column 137
2.999999999999915e+00
Column 138
2.999999999999915e+00
Column 139
2.999999999999914e+00
Column 140
2.999999999999914e+00
Column 141
2.999999999999914e+00
Column 142
2.999999999999914e+00
Column 143
2.999999999999914e+00
Column 144
2.999999999999914e+00
Column 145
2.999999999999914e+00
Column 146
2.999999999999914e+00
Column 147
2.999999999999914e+00
Column 148
2.999999999999914e+00
Column 149
2.999999999999914e+00
Column 150
2.999999999999914e+00
Column 151
2.999999999999913e+00
Column 152
2.999999999999914e+00
Column 153
2.999999999999913e+00
Column 154
2.999999999999913e+00
Column 155
2.999999999999913e+00
Column 156
2.999999999999885e+00
Column 157
2.999999999999872e+00
Column 158
2.999999999999865e+00
Column 159
2.999999999999862e+00
Column 160
2.999999999999861e+00
Column 161
2.999999999999860e+00
Column 162
2.999999999999859e+00
Column 163
2.999999999999859e+00
Column 164
2.999999999999859e+00
Column 165
2.999999999999858e+00
Column 166
2.999999999999858e+00
Column 167
2.999999999999859e+00
Column 168
2.999999999999859e+00
Column 169
2.999999999999859e+00
Column 170
2.999999999999859e+00
Column 171
2.999999999999859e+00
Column 172
2.999999999999859e+00
Column 173
2.999999999999859e+00
Column 174
2.999999999999859e+00
Column 175
2.999999999999858e+00
Column 176
2.999999999999859e+00
Column 177
2.999999999999860e+00
Column 178
2.999999999999861e+00
Column 179
2.999999999999863e+00
Column 180
2.999999999999868e+00
Column 181
2.999999999999877e+00
Column 182
2.999999999999894e+00
Column 183
2.999999999999895e+00
Column 184
2.999999999999894e+00
Column 185
2.999999999999894e+00
Column 186
2.999999999999894e+00
Column 187
2.999999999999895e+00
Column 188
2.999999999999895e+00
Column 189
2.999999999999895e+00
Column 190
2.999999999999894e+00
Column 191
2.999999999999895e+00
Column 192
2.999999999999895e+00
Column 193
2.999999999999895e+00
Column 194
2.999999999999895e+00
Column 195
2.999999999999894e+00
Column 196
2.999999999999894e+00
Column 197
2.999999999999894e+00
Column 198
2.999999999999894e+00
Column 199
2.999999999999893e+00
Column 200
2.999999999999893e+00
Column 201
2.999999999999893e+00
Column 202
2.999999999999893e+00
Column 203
2.999999999999893e+00
Column 204
2.999999999999893e+00
Column 205
2.999999999999893e+00
Column 206
2.999999999999893e+00
Column 207
2.999999999999893e+00
Column 208
2.999999999999894e+00
Column 209
2.999999999999894e+00
Column 210
2.999999999999894e+00
Column 211
2.999999999999894e+00
Column 212
2.999999999999895e+00
Column 213
2.999999999999895e+00
Column 214
2.999999999999895e+00
Column 215
2.999999999999895e+00
Column 216
2.999999999999895e+00
Column 217
2.999999999999896e+00
Column 218
2.999999999999896e+00
Column 219
2.999999999999896e+00
Column 220
2.999999999999897e+00
Column 221
2.999999999999909e+00
Column 222
2.999999999999909e+00
Column 223
2.999999999999909e+00
Column 224
2.999999999999896e+00
Column 225
2.999999999999893e+00
Column 226
2.999999999999872e+00
Column 227
2.999999999999872e+00
Column 228
2.999999999999872e+00
Column 229
2.999999999999857e+00
Column 230
2.999999999999857e+00
Column 231
2.999999999999856e+00
Column 232
2.999999999999857e+00
Column 233
2.999999999999857e+00
Column 234
2.999999999999857e+00
Column 235
2.999999999999856e+00
Column 236
2.999999999999857e+00
Column 237
2.999999999999856e+00
Column 238
2.999999999999856e+00
Column 239
2.999999999999856e+00
Column 240
2.999999999999855e+00
Column 241
2.999999999999855e+00
Column 242
2.999999999999856e+00
Column 243
2.999999999999856e+00
Column 244
2.999999999999856e+00
Column 245
2.999999999999855e+00
Column 246
2.999999999999855e+00
Column 247
2.999999999999855e+00
Column 248
2.999999999999855e+00
Column 249
2.999999999999855e+00
Column 250
2.999999999999855e+00
Column 251
2.999999999999855e+00
Column 252
2.999999999999854e+00
Column 253
2.999999999999854e+00
Column 254
2.999999999999854e+00
Column 255
2.999999999999855e+00
Column 256
2.999999999999855e+00
Column 257
2.999999999999856e+00
Column 258
2.999999999999856e+00
Column 259
2.999999999999855e+00
Column 260
2.999999999999856e+00
Column 261
2.999999999999856e+00
Column 262
2.999999999999856e+00
Column 263
2.999999999999856e+00
Column 264
2.999999999999856e+00
Column 265
2.999999999999856e+00
Column 266
2.999999999999856e+00
Column 267
2.999999999999856e+00
Column 268
2.999999999999826e+00
Column 269
2.999999999999810e+00
Column 270
2.999999999999803e+00
Column 271
2.999999999999799e+00
Column 272
2.999999999999798e+00
Column 273
2.999999999999797e+00
Column 274
2.999999999999797e+00
Column 275
2.999999999999797e+00
Column 276
2.999999999999797e+00
Column 277
2.999999999999797e+00
Column 278
2.999999999999796e+00
Column 279
2.999999999999795e+00
Column 280
2.999999999999795e+00
Column 281
2.999999999999795e+00
Column 282
2.999999999999796e+00
Column 283
2.999999999999796e+00
Column 284
2.999999999999795e+00
Column 285
2.999999999999796e+00
Column 286
2.999999999999796e+00
Column 287
2.999999999999796e+00
Column 288
2.999999999999797e+00
Column 289
2.999999999999799e+00
Column 290
2.999999999999803e+00
Column 291
2.999999999999811e+00
Column 292
2.999999999999828e+00
Column 293
2.999999999999860e+00
Column 294
2.999999999999860e+00
Column 295
2.999999999999860e+00
Column 296
2.999999999999860e+00
Column 297
2.999999999999860e+00
Column 298
2.999999999999861e+00
Column 299
2.999999999999861e+00
Column 300
2.999999999999860e+00
Column 301
2.999999999999860e+00
f = @(y) [-2*y(2)*y(1);y(1)^2+y(3)^2-y(2)^2-1;-2*(y(2)+y(1))*y(3)];
y0=[0.5,2,10];
tspan=[0,15];
[tout,yout]=rk23(f,tspan,y0);
??? Error using ==>
@(y)[-2*y(2)*y(1);y(1)^2+y(3)^2-y(2)^2-1;-2*(y(2)+y(1))*y(3)]
Too many input arguments.
Error in ==> rk23 at 35
f(:,1) = odefun(tout(n),yout(:,n));
f = @(t,y) [-2*y(2)*y(1);y(1)^2+y(3)^2-y(2)^2-1;-2*(y(2)+y(1))*y(3)];
[tout,yout]=rk23(f,tspan,y0);
size(yout)
ans =
416 3
E=(yout(:,1).^2+yout(:,2).^2+yout(:,3).^2)./(2*yout(:,1));
E=(yout(:,1).^2+yout(:,2).^2+yout(:,3).^2+1)./(2*yout(:,1));
plot(tout,E)
plot(tout,E,'.')
options.RelTol=1e-2;
options.AbsTol=1e-2;
[tout,yout]=rk23(f,tspan,y0,options);
Warning: time step rejected
> In rk23 at 46
Warning: time step rejected
> In rk23 at 46
Warning: time step rejected
> In rk23 at 46
Warning: time step rejected
> In rk23 at 46
Warning: time step rejected
> In rk23 at 46
Warning: time step rejected
> In rk23 at 46
Warning: time step rejected
> In rk23 at 46
Warning: time step rejected
> In rk23 at 46
plot(tout,E,'.')
??? Error using ==> plot
Vectors must be the same lengths.
E=(yout(:,1).^2+yout(:,2).^2+yout(:,3).^2+1)./(2*yout(:,1));
plot(tout,E,'.')
plot(tout,yout(:,1),'.')
plot(tout,yout(:,2),'.')
plot(tout,yout(:,3),'.')
plot(tout,E,'.')
diary off