← Back to branch summary

~ubuntu-branches/ubuntu/hoary/scilab/hoary

« back to all changes in this revision

Viewing changes to demos/npend/Maple/npend.f

  • Committer: Bazaar Package Importer
  • Author(s): Torsten Werner
  • Date: 2005-01-09 22:58:21 UTC
  • mfrom: (1.1.1 upstream)
  • Revision ID: james.westby@ubuntu.com-20050109225821-473xr8vhgugxxx5j
Tags: 3.0-12
http://bugs.debian.org/255869
changed configure.in to build scilab's own malloc.o, closes: #255869

Show diffs side-by-side

added added

removed removed

Lines of Context:
demos/npend/Maple/npend.f
1
 
c      
2
 
c     SUBROUTINE npend
3
 
c   
4
 
c     Copyright INRIA   
5
 
      subroutine npend(neq,t,th,ydot)
6
 
        parameter (n=10)
7
 
        implicit doubleprecision (t)
8
 
        doubleprecision t,th(2*n),ydot(2*n),r(n),j(n),m(n)
9
 
        doubleprecision me3s(n,n),cc3s(n,n),const(n,1)
10
 
        doubleprecision w(3*n),rcond    
11
 
 
12
 
        integer i,k,neq,ierr
13
 
        data g / 9.81/
14
 
        data r / n*1.0/
15
 
        data m / n*1.0/
16
 
        data j / n*0.3/
17
 
      t1 = -th(9)
18
 
      t4 = 2*m(10)
19
 
      t5 = m(9)+t4
20
 
      t9 = 2*r(7)*r(9)*cos(th(7)+t1)*t5
21
 
      t15 = 2*r(3)*r(9)*cos(th(3)+t1)*t5
22
 
      t16 = -th(8)
23
 
      t19 = 2*m(9)
24
 
      t20 = m(8)+t4+t19
25
 
      t24 = 2*r(7)*r(8)*cos(th(7)+t16)*t20
26
 
      t25 = r(2)**2
27
 
      t44 = -th(6)
28
 
      t47 = 2*m(7)
29
 
      t48 = 2*m(8)
30
 
      t49 = m(6)+t19+t4+t47+t48
31
 
      t53 = 2*r(5)*r(6)*cos(th(5)+t44)*t49
32
 
      t59 = 2*r(6)*r(9)*cos(th(6)+t1)*t5
33
 
      t60 = -th(10)
34
 
      t66 = 2*r(5)*m(10)*r(10)*cos(th(5)+t60)
35
 
      t72 = 2*r(1)*r(8)*cos(th(1)+t16)*t20
36
 
      t73 = -th(4)
37
 
      t76 = 2*m(5)
38
 
      t77 = 2*m(6)
39
 
      t78 = t48+m(4)+t19+t76+t47+t4+t77
40
 
      t82 = 2*r(3)*r(4)*cos(th(3)+t73)*t78
41
 
      t88 = 2*r(8)*r(9)*cos(th(8)+t1)*t5
42
 
      t89 = r(5)**2
43
 
      t102 = -th(7)
44
 
      t105 = t4+m(7)+t48+t19
45
 
      t109 = 2*r(6)*r(7)*cos(th(6)+t102)*t105
46
 
      t115 = 2*r(9)*m(10)*r(10)*cos(th(9)+t60)
47
 
      t116 = r(3)**2
48
 
      t138 = 2*r(4)*r(9)*cos(th(4)+t1)*t5
49
 
      t144 = 2*r(3)*m(10)*r(10)*cos(th(3)+t60)
50
 
      t150 = 2*r(4)*m(10)*r(10)*cos(th(4)+t60)
51
 
      t151 = -th(3)
52
 
      t154 = 2*m(4)
53
 
      t155 = t76+t154+t4+t48+t47+m(3)+t19+t77
54
 
      t159 = 2*r(1)*r(3)*cos(th(1)+t151)*t155
55
 
      t160 = r(8)**2
56
 
      t167 = r(6)**2
57
 
      t183 = 2*r(6)*m(10)*r(10)*cos(th(6)+t60)
58
 
      t184 = r(10)**2
59
 
      t192 = 2*r(2)*r(4)*cos(th(2)+t73)*t78
60
 
      t193 = r(9)**2
61
 
      t198 = r(7)**2
62
 
      t212 = 2*r(3)*r(6)*cos(th(3)+t44)*t49
63
 
      t218 = 2*r(2)*r(3)*cos(th(2)+t151)*t155
64
 
      t224 = 2*r(5)*r(9)*cos(th(5)+t1)*t5
65
 
      t230 = 2*r(5)*r(7)*cos(th(5)+t102)*t105
66
 
      t236 = 2*r(7)*m(10)*r(10)*cos(th(7)+t60)
67
 
      t242 = 2*r(1)*r(4)*cos(th(1)+t73)*t78
68
 
      t248 = 2*r(1)*r(6)*cos(th(1)+t44)*t49
69
 
      t254 = 2*r(2)*r(6)*cos(th(2)+t44)*t49
70
 
      t260 = 2*r(2)*m(10)*r(10)*cos(th(2)+t60)
71
 
      t266 = 2*r(8)*m(10)*r(10)*cos(th(8)+t60)
72
 
      t267 = -th(5)
73
 
      t270 = t77+t48+t19+t47+m(5)+t4
74
 
      t274 = 2*r(1)*r(5)*cos(th(1)+t267)*t270
75
 
      t280 = 2*r(4)*r(5)*cos(th(4)+t267)*t270
76
 
      t286 = 2*r(4)*r(7)*cos(th(4)+t102)*t105
77
 
      t292 = 2*r(6)*r(8)*cos(th(6)+t16)*t20
78
 
      t298 = 2*r(1)*r(9)*cos(th(1)+t1)*t5
79
 
      t304 = 2*r(5)*r(8)*cos(th(5)+t16)*t20
80
 
      t310 = 2*r(3)*r(5)*cos(th(3)+t267)*t270
81
 
      t316 = 2*r(3)*r(7)*cos(th(3)+t102)*t105
82
 
      t322 = 2*r(2)*r(7)*cos(th(2)+t102)*t105
83
 
      t331 = 2*r(1)*r(2)*cos(th(1)-th(2))*(2*m(3)+m(2)+t154+t47+t4+t76+t
84
 
     +77+t19+t48)
85
 
      t332 = r(4)**2
86
 
      t352 = 2*r(4)*r(8)*cos(th(4)+t16)*t20
87
 
      t358 = 2*r(2)*r(5)*cos(th(2)+t267)*t270
88
 
      t364 = 2*r(2)*r(9)*cos(th(2)+t1)*t5
89
 
      t370 = 2*r(2)*r(8)*cos(th(2)+t16)*t20
90
 
      t376 = 2*r(1)*r(7)*cos(th(1)+t102)*t105
91
 
      t382 = 2*r(3)*r(8)*cos(th(3)+t16)*t20
92
 
      t388 = 2*r(4)*r(6)*cos(th(4)+t44)*t49
93
 
      t389 = r(1)**2
94
 
      t417 = 2*r(1)*m(10)*r(10)*cos(th(1)+t60)
95
 
         me3s(9,7) = t9
96
 
         me3s(3,9) = t15
97
 
         me3s(8,7) = t24
98
 
         me3s(2,2) = J(2)+4*t25*m(3)+4*t25*m(4)+4*t25*m(7)+4*t25*m(8)+m(
99
 
     +2)*t25+4*t25*m(9)+4*t25*m(5)+4*t25*m(10)+4*t25*m(6)
100
 
         me3s(6,5) = t53
101
 
         me3s(9,6) = t59
102
 
         me3s(10,5) = t66
103
 
         me3s(8,1) = t72
104
 
         me3s(4,3) = t82
105
 
         me3s(9,8) = t88
106
 
         me3s(5,5) = 4*t89*m(8)+4*t89*m(10)+4*t89*m(7)+4*t89*m(9)+4*t89*
107
 
     +m(6)+J(5)+m(5)*t89
108
 
         me3s(7,6) = t109
109
 
         me3s(10,9) = t115
110
 
         me3s(3,3) = 4*t116*m(4)+4*t116*m(10)+4*t116*m(9)+4*t116*m(5)+4*
111
 
     +t116*m(6)+m(3)*t116+J(3)+4*t116*m(7)+4*t116*m(8)
112
 
         me3s(4,9) = t138
113
 
         me3s(3,10) = t144
114
 
         me3s(10,4) = t150
115
 
         me3s(3,1) = t159
116
 
         me3s(8,8) = 4*t160*m(9)+J(8)+m(8)*t160+4*t160*m(10)
117
 
         me3s(6,6) = m(6)*t167+4*t167*m(7)+4*t167*m(8)+J(6)+4*t167*m(9)+
118
 
     +4*t167*m(10)
119
 
         me3s(10,6) = t183
120
 
         me3s(10,10) = m(10)*t184+J(10)
121
 
         me3s(4,2) = t192
122
 
         me3s(9,9) = 4*t193*m(10)+J(9)+m(9)*t193
123
 
         me3s(7,7) = 4*t198*m(9)+4*t198*m(10)+4*t198*m(8)+J(7)+m(7)*t198
124
 
         me3s(6,3) = t212
125
 
         me3s(3,2) = t218
126
 
         me3s(10,3) = t144
127
 
         me3s(9,5) = t224
128
 
         me3s(7,5) = t230
129
 
         me3s(10,7) = t236
130
 
         me3s(1,4) = t242
131
 
         me3s(4,10) = t150
132
 
         me3s(5,6) = t53
133
 
         me3s(6,1) = t248
134
 
         me3s(6,2) = t254
135
 
         me3s(10,2) = t260
136
 
         me3s(8,9) = t88
137
 
         me3s(6,7) = t109
138
 
         me3s(10,8) = t266
139
 
         me3s(1,5) = t274
140
 
         me3s(9,10) = t115
141
 
         me3s(5,7) = t230
142
 
         me3s(5,4) = t280
143
 
         me3s(7,8) = t24
144
 
         me3s(7,4) = t286
145
 
         me3s(9,4) = t138
146
 
         me3s(1,3) = t159
147
 
         me3s(8,10) = t266
148
 
         me3s(6,8) = t292
149
 
         me3s(9,1) = t298
150
 
         me3s(1,6) = t248
151
 
         me3s(5,8) = t304
152
 
         me3s(5,3) = t310
153
 
         me3s(7,9) = t9
154
 
         me3s(7,3) = t316
155
 
         me3s(9,3) = t15
156
 
         me3s(2,7) = t322
157
 
         me3s(1,2) = t331
158
 
         me3s(4,4) = 4*t332*m(10)+4*t332*m(9)+4*t332*m(6)+4*t332*m(5)+4*
159
 
     +t332*m(7)+J(4)+m(4)*t332+4*t332*m(8)
160
 
         me3s(4,1) = t242
161
 
         me3s(6,9) = t59
162
 
         me3s(8,4) = t352
163
 
         me3s(4,5) = t280
164
 
         me3s(3,4) = t82
165
 
         me3s(5,9) = t224
166
 
         me3s(5,2) = t358
167
 
         me3s(7,10) = t236
168
 
         me3s(7,2) = t322
169
 
         me3s(9,2) = t364
170
 
         me3s(2,8) = t370
171
 
         me3s(1,7) = t376
172
 
         me3s(6,10) = t183
173
 
         me3s(8,3) = t382
174
 
         me3s(4,6) = t388
175
 
         me3s(3,5) = t310
176
 
         me3s(2,9) = t364
177
 
         me3s(7,1) = t376
178
 
         me3s(1,8) = t72
179
 
         me3s(5,10) = t66
180
 
         me3s(2,4) = t192
181
 
         me3s(8,2) = t370
182
 
         me3s(4,7) = t286
183
 
         me3s(3,6) = t212
184
 
         me3s(2,1) = t331
185
 
         me3s(1,1) = 4*t389*m(8)+4*t389*m(10)+m(1)*t389+4*t389*m(6)+4*t3
186
 
     +89*m(9)+4*t389*m(7)+4*t389*m(2)+J(1)+4*t389*m(5)+4*t389*m(4)+4*t38
187
 
     +9*m(3)
188
 
         me3s(2,10) = t260
189
 
         me3s(1,9) = t298
190
 
         me3s(10,1) = t417
191
 
         me3s(2,5) = t358
192
 
         me3s(4,8) = t352
193
 
         me3s(3,7) = t316
194
 
         me3s(8,5) = t304
195
 
         me3s(2,3) = t218
196
 
         me3s(1,10) = t417
197
 
         me3s(5,1) = t274
198
 
         me3s(2,6) = t254
199
 
         me3s(3,8) = t382
200
 
         me3s(8,6) = t292
201
 
         me3s(6,4) = t388
202
 
      t1 = -th(9)
203
 
      t4 = 2*m(10)
204
 
      t5 = m(9)+t4
205
 
      t8 = r(7)*r(9)*sin(th(7)+t1)*t5
206
 
      t14 = r(3)*r(9)*sin(th(3)+t1)*t5
207
 
      t16 = -th(8)
208
 
      t19 = 2*m(9)
209
 
      t20 = m(8)+t4+t19
210
 
      t23 = r(7)*r(8)*sin(th(7)+t16)*t20
211
 
      t25 = -th(6)
212
 
      t28 = 2*m(7)
213
 
      t29 = 2*m(8)
214
 
      t30 = m(6)+t19+t4+t28+t29
215
 
      t33 = r(5)*r(6)*sin(th(5)+t25)*t30
216
 
      t39 = r(6)*r(9)*sin(th(6)+t1)*t5
217
 
      t41 = -th(10)
218
 
      t46 = r(5)*m(10)*r(10)*sin(th(5)+t41)
219
 
      t52 = r(1)*r(8)*sin(th(1)+t16)*t20
220
 
      t54 = -th(4)
221
 
      t57 = 2*m(5)
222
 
      t58 = 2*m(6)
223
 
      t59 = t29+m(4)+t19+t57+t28+t4+t58
224
 
      t62 = r(3)*r(4)*sin(th(3)+t54)*t59
225
 
      t68 = r(8)*r(9)*sin(th(8)+t1)*t5
226
 
      t70 = -th(7)
227
 
      t73 = t4+m(7)+t29+t19
228
 
      t76 = r(6)*r(7)*sin(th(6)+t70)*t73
229
 
      t82 = r(9)*m(10)*r(10)*sin(th(9)+t41)
230
 
      t88 = r(4)*r(9)*sin(th(4)+t1)*t5
231
 
      t94 = r(3)*m(10)*r(10)*sin(th(3)+t41)
232
 
      t100 = r(4)*m(10)*r(10)*sin(th(4)+t41)
233
 
      t102 = -th(3)
234
 
      t105 = 2*m(4)
235
 
      t106 = t57+t105+t4+t29+t28+m(3)+t19+t58
236
 
      t109 = r(1)*r(3)*sin(th(1)+t102)*t106
237
 
      t115 = r(6)*m(10)*r(10)*sin(th(6)+t41)
238
 
      t121 = r(2)*r(4)*sin(th(2)+t54)*t59
239
 
      t127 = r(3)*r(6)*sin(th(3)+t25)*t30
240
 
      t133 = r(2)*r(3)*sin(th(2)+t102)*t106
241
 
      t140 = r(5)*r(9)*sin(th(5)+t1)*t5
242
 
      t146 = r(5)*r(7)*sin(th(5)+t70)*t73
243
 
      t152 = r(7)*m(10)*r(10)*sin(th(7)+t41)
244
 
      t158 = r(1)*r(4)*sin(th(1)+t54)*t59
245
 
      t166 = r(1)*r(6)*sin(th(1)+t25)*t30
246
 
      t172 = r(2)*r(6)*sin(th(2)+t25)*t30
247
 
      t178 = r(2)*m(10)*r(10)*sin(th(2)+t41)
248
 
      t186 = r(8)*m(10)*r(10)*sin(th(8)+t41)
249
 
      t188 = -th(5)
250
 
      t191 = t58+t29+t19+t28+m(5)+t4
251
 
      t194 = r(1)*r(5)*sin(th(1)+t188)*t191
252
 
      t202 = r(4)*r(5)*sin(th(4)+t188)*t191
253
 
      t209 = r(4)*r(7)*sin(th(4)+t70)*t73
254
 
      t218 = r(6)*r(8)*sin(th(6)+t16)*t20
255
 
      t224 = r(1)*r(9)*sin(th(1)+t1)*t5
256
 
      t231 = r(5)*r(8)*sin(th(5)+t16)*t20
257
 
      t237 = r(3)*r(5)*sin(th(3)+t188)*t191
258
 
      t244 = r(3)*r(7)*sin(th(3)+t70)*t73
259
 
      t251 = r(2)*r(7)*sin(th(2)+t70)*t73
260
 
      t260 = r(1)*r(2)*sin(th(1)-th(2))*(2*m(3)+m(2)+t105+t28+t4+t57+t58
261
 
     ++t19+t29)
262
 
      t268 = r(4)*r(8)*sin(th(4)+t16)*t20
263
 
      t277 = r(2)*r(5)*sin(th(2)+t188)*t191
264
 
      t285 = r(2)*r(9)*sin(th(2)+t1)*t5
265
 
      t291 = r(2)*r(8)*sin(th(2)+t16)*t20
266
 
      t297 = r(1)*r(7)*sin(th(1)+t70)*t73
267
 
      t304 = r(3)*r(8)*sin(th(3)+t16)*t20
268
 
      t310 = r(4)*r(6)*sin(th(4)+t25)*t30
269
 
      t328 = r(1)*m(10)*r(10)*sin(th(1)+t41)
270
 
         cc3s(9,7) = -2*t8
271
 
         cc3s(3,9) = 2*t14
272
 
         cc3s(8,7) = -2*t23
273
 
         cc3s(2,2) = 0
274
 
         cc3s(6,5) = -2*t33
275
 
         cc3s(9,6) = -2*t39
276
 
         cc3s(10,5) = -2*t46
277
 
         cc3s(8,1) = -2*t52
278
 
         cc3s(4,3) = -2*t62
279
 
         cc3s(9,8) = -2*t68
280
 
         cc3s(5,5) = 0
281
 
         cc3s(7,6) = -2*t76
282
 
         cc3s(10,9) = -2*t82
283
 
         cc3s(3,3) = 0
284
 
         cc3s(4,9) = 2*t88
285
 
         cc3s(3,10) = 2*t94
286
 
         cc3s(10,4) = -2*t100
287
 
         cc3s(3,1) = -2*t109
288
 
         cc3s(8,8) = 0
289
 
         cc3s(6,6) = 0
290
 
         cc3s(10,6) = -2*t115
291
 
         cc3s(10,10) = 0
292
 
         cc3s(4,2) = -2*t121
293
 
         cc3s(9,9) = 0
294
 
         cc3s(7,7) = 0
295
 
         cc3s(6,3) = -2*t127
296
 
         cc3s(3,2) = -2*t133
297
 
         cc3s(10,3) = -2*t94
298
 
         cc3s(9,5) = -2*t140
299
 
         cc3s(7,5) = -2*t146
300
 
         cc3s(10,7) = -2*t152
301
 
         cc3s(1,4) = 2*t158
302
 
         cc3s(4,10) = 2*t100
303
 
         cc3s(5,6) = 2*t33
304
 
         cc3s(6,1) = -2*t166
305
 
         cc3s(6,2) = -2*t172
306
 
         cc3s(10,2) = -2*t178
307
 
         cc3s(8,9) = 2*t68
308
 
         cc3s(6,7) = 2*t76
309
 
         cc3s(10,8) = -2*t186
310
 
         cc3s(1,5) = 2*t194
311
 
         cc3s(9,10) = 2*t82
312
 
         cc3s(5,7) = 2*t146
313
 
         cc3s(5,4) = -2*t202
314
 
         cc3s(7,8) = 2*t23
315
 
         cc3s(7,4) = -2*t209
316
 
         cc3s(9,4) = -2*t88
317
 
         cc3s(1,3) = 2*t109
318
 
         cc3s(8,10) = 2*t186
319
 
         cc3s(6,8) = 2*t218
320
 
         cc3s(9,1) = -2*t224
321
 
         cc3s(1,6) = 2*t166
322
 
         cc3s(5,8) = 2*t231
323
 
         cc3s(5,3) = -2*t237
324
 
         cc3s(7,9) = 2*t8
325
 
         cc3s(7,3) = -2*t244
326
 
         cc3s(9,3) = -2*t14
327
 
         cc3s(2,7) = 2*t251
328
 
         cc3s(1,2) = 2*t260
329
 
         cc3s(4,4) = 0
330
 
         cc3s(4,1) = -2*t158
331
 
         cc3s(6,9) = 2*t39
332
 
         cc3s(8,4) = -2*t268
333
 
         cc3s(4,5) = 2*t202
334
 
         cc3s(3,4) = 2*t62
335
 
         cc3s(5,9) = 2*t140
336
 
         cc3s(5,2) = -2*t277
337
 
         cc3s(7,10) = 2*t152
338
 
         cc3s(7,2) = -2*t251
339
 
         cc3s(9,2) = -2*t285
340
 
         cc3s(2,8) = 2*t291
341
 
         cc3s(1,7) = 2*t297
342
 
         cc3s(6,10) = 2*t115
343
 
         cc3s(8,3) = -2*t304
344
 
         cc3s(4,6) = 2*t310
345
 
         cc3s(3,5) = 2*t237
346
 
         cc3s(2,9) = 2*t285
347
 
         cc3s(7,1) = -2*t297
348
 
         cc3s(1,8) = 2*t52
349
 
         cc3s(5,10) = 2*t46
350
 
         cc3s(2,4) = 2*t121
351
 
         cc3s(8,2) = -2*t291
352
 
         cc3s(4,7) = 2*t209
353
 
         cc3s(3,6) = 2*t127
354
 
         cc3s(2,1) = -2*t260
355
 
         cc3s(1,1) = 0
356
 
         cc3s(2,10) = 2*t178
357
 
         cc3s(1,9) = 2*t224
358
 
         cc3s(10,1) = -2*t328
359
 
         cc3s(2,5) = 2*t277
360
 
         cc3s(4,8) = 2*t268
361
 
         cc3s(3,7) = 2*t244
362
 
         cc3s(8,5) = -2*t231
363
 
         cc3s(2,3) = 2*t133
364
 
         cc3s(1,10) = 2*t328
365
 
         cc3s(5,1) = -2*t194
366
 
         cc3s(2,6) = 2*t172
367
 
         cc3s(3,8) = 2*t304
368
 
         cc3s(8,6) = -2*t218
369
 
         cc3s(6,4) = -2*t310
370
 
      t2 = 2*m(10)
371
 
      t3 = 2*m(9)
372
 
      t9 = 2*m(5)
373
 
      t10 = 2*m(4)
374
 
      t11 = 2*m(8)
375
 
      t12 = 2*m(7)
376
 
      t13 = 2*m(6)
377
 
      t39 = 2*m(3)
378
 
         const(8,1) = g*cos(th(8))*r(8)*(m(8)+t2+t3)
379
 
         const(3,1) = g*r(3)*cos(th(3))*(t9+t10+t2+t11+t12+m(3)+t3+t13)
380
 
         const(6,1) = g*r(6)*cos(th(6))*(m(6)+t3+t2+t12+t11)
381
 
         const(9,1) = g*r(9)*cos(th(9))*(m(9)+t2)
382
 
         const(4,1) = g*r(4)*cos(th(4))*(t11+m(4)+t3+t9+t12+t2+t13)
383
 
         const(7,1) = g*r(7)*cos(th(7))*(t2+m(7)+t11+t3)
384
 
         const(2,1) = g*r(2)*cos(th(2))*(t39+m(2)+t10+t12+t2+t9+t13+t3+t
385
 
     +11)
386
 
         const(1,1) = g*r(1)*cos(th(1))*(m(1)+t3+2*m(2)+t10+t39+t11+t12+
387
 
     +t13+t9+t2)
388
 
         const(10,1) = m(10)*g*r(10)*cos(th(10))
389
 
         const(5,1) = g*r(5)*cos(th(5))*(t13+t11+t3+t12+m(5)+t2)
390
 
c         
391
 
        do 1000, i =1,n ,1
392
 
          ydot(i) = th(i+n)
393
 
 1000   continue
394
 
c       
395
 
c         
396
 
        do 1001, i =1,n ,1
397
 
c           
398
 
          do 1002, k =1,n ,1
399
 
            Const(i,1) =  Const(i,1)+CC3S(i,k)*(th(k+n)**2)
400
 
 1002     continue
401
 
c         
402
 
           Const(i,1) = -Const(i,1)
403
 
 1001   continue
404
 
c       
405
 
c        we must solve  M z =const 
406
 
c        which gives ydot((n+1)..2*n) 
407
 
        call dlslv(me3s,n,n,Const,n,1,w, rcond,ierr,1)
408
 
        if (ierr.ne.0) then
409
 
          write(6,2000) 
410
 
 2000     format('Matrice mal conditionnee')
411
 
        endif
412
 
c         
413
 
        do 1003, i =1,n ,1
414
 
          ydot(n+i) = const(i,1)
415
 
 1003   continue
416
 
c       
417
 
        return
418
 
      end