# File diffeq.rb, line 163
def atomall
 k = 1
 while k <= 1  do 
 $y[k] = $y_higher[1][k]
 k += 1
end
 while k <= $max_taylor  do 
 $y[k] = $zero
 k = k + 1
end
# END OUTFILE1
# BEGIN ATOMHDR1
# emit pre sin $eq_no = 1 i = 1
# emit pre sin 1 $eq_no = 1
$tmp1[1] = $x[1].sin
$tmp1_g[1] = $x[1].cos
# emit pre add $eq_no = 1 i = 1
$tmp2[1] = (const(0.0)[1] + ($tmp1[1]))
# emit pre assign $eq_no = 1 i = 1 $min_hdrs = 4
if 2 <= $max_taylor then
$y[2] = $tmp2[1] * ($h ** 1) / factorial_2(0,1)
end
# END ATOMHDR1
# BEGIN ATOMHDR2
# emit pre sin $eq_no = 1 i = 2
# emit pre sin  2 $eq_no = 1
$tmp1[2] = att(1,$tmp1_g,$x,1)
$tmp1_g[2] = $minus_one*att(1,$tmp1,$x,1)
# emit pre add $eq_no = 1 i = 2
$tmp2[2] = (const(0.0)[2] + ($tmp1[2]))
# emit pre assign $eq_no = 1 i = 2 $min_hdrs = 4
if 3 <= $max_taylor then
$y[3] = $tmp2[2] * ($h ** 1) / factorial_2(1,2)
end
# END ATOMHDR2
# BEGIN ATOMHDR3
# emit pre sin $eq_no = 1 i = 3
# emit pre sin  3 $eq_no = 1
$tmp1[3] = att(2,$tmp1_g,$x,1)
$tmp1_g[3] = $minus_one*att(2,$tmp1,$x,1)
# emit pre add $eq_no = 1 i = 3
$tmp2[3] = (const(0.0)[3] + ($tmp1[3]))
# emit pre assign $eq_no = 1 i = 3 $min_hdrs = 4
if 4 <= $max_taylor then
$y[4] = $tmp2[3] * ($h ** 1) / factorial_2(2,3)
end
# END ATOMHDR3
# BEGIN ATOMHDR4
# emit pre sin $eq_no = 1 i = 4
# emit pre sin  4 $eq_no = 1
$tmp1[4] = att(3,$tmp1_g,$x,1)
$tmp1_g[4] = $minus_one*att(3,$tmp1,$x,1)
# emit pre add $eq_no = 1 i = 4
$tmp2[4] = (const(0.0)[4] + ($tmp1[4]))
# emit pre assign $eq_no = 1 i = 4 $min_hdrs = 4
if 5 <= $max_taylor then
$y[5] = $tmp2[4] * ($h ** 1) / factorial_2(3,4)
end
# END ATOMHDR4
# BEGIN OUTFILE3
k = 4 + 0 + 1
while k <= $max_taylor do
# END OUTFILE3
# BEGIN OUTFILE4
# emit sin $eq_no = 1
pr = 0
k2 =  k - pr 
$tmp1[k2] = att(k2-1,$tmp1_g,$x,1)
$tmp1_g[k2] = $minus_one*att(k2-1,$tmp1,$x,1)
# emit add $eq_no = 1
pr = 0
k2 =  k - pr 
$tmp2[k2] = (const(0.0)[k2]+($tmp1[k2]))
# emit assign $eq_no = 1
pr = 0
k2 =  k - pr 
if k2 + 1 <= $max_taylor then
$y[k2 + 1 ] = $tmp2[k2] * ($h ** 1)  / factorial_2( k2-1 ,1+k2-1)
end
k = k + 1
end
k = 1
while k <= $max_taylor do
$y_higher[1 + 1][k] = $tmp2[k]
k = k + 1
end
# END OUTFILE4
# BEGIN OUTFILE5
end