Question: Error, (in fprintf) number expected for floating point format

hi .please help me for this error

Error, (in fprintf) number expected for floating point format...

float.mw

T(r, t) := -(1.082744105*(.4426072302*BesselJ(1, 1.511248492*r)+.5201541352*BesselY(1, 1.511248492*r)))*Besselj(0, 3.588021808)*(1.-1.*exp(-3.218475123670897216*t))/Besselj(1, 1.794010904)-(0.7489242453e-1*(-.4693819312*BesselJ(1, 2.197141326*r)+.6841159909*BesselY(1, 2.197141326*r)))*Besselj(0, 9.604121522)*(1.-1.*exp(-23.059787552335899121*t))/Besselj(1, 4.802060761)-(0.1290812036e-1*(-.6766364675*BesselJ(1, 3.111788255*r)-.7242359729*BesselY(1, 3.111788255*r)))*Besselj(0, 15.817923424)*(1.-1.*exp(-62.551675361881970944*t))/Besselj(1, 7.908961712)-(0.5633715672e-1*(.8435152108*BesselJ(1, 4.692507174*r)+.8812474286*BesselY(1, 4.692507174*r)))*Besselj(0, 9.604121522)*(1.-1.*exp(-23.059787552335899121*t))/Besselj(1, 4.802060761)-(0.6254981294e-1*(-.8481299777*BesselJ(1, 5.429681040*r)+1.002082035*BesselY(1, 5.429681040*r)))*Besselj(0, 9.604121522)*(1.-1.*exp(-23.059787552335899121*t))/Besselj(1, 4.802060761)

-1.082744105*(.4426072302*BesselJ(1, 1.511248492*r)+.5201541352*BesselY(1, 1.511248492*r))*Besselj(0, 3.588021808)*(1.-1.*exp(-3.218475123670897216*t))/Besselj(1, 1.794010904)-0.7489242453e-1*(-.4693819312*BesselJ(1, 2.197141326*r)+.6841159909*BesselY(1, 2.197141326*r))*Besselj(0, 9.604121522)*(1.-1.*exp(-23.059787552335899121*t))/Besselj(1, 4.802060761)-0.1290812036e-1*(-.6766364675*BesselJ(1, 3.111788255*r)-.7242359729*BesselY(1, 3.111788255*r))*Besselj(0, 15.817923424)*(1.-1.*exp(-62.551675361881970944*t))/Besselj(1, 7.908961712)-0.5633715672e-1*(.8435152108*BesselJ(1, 4.692507174*r)+.8812474286*BesselY(1, 4.692507174*r))*Besselj(0, 9.604121522)*(1.-1.*exp(-23.059787552335899121*t))/Besselj(1, 4.802060761)-0.6254981294e-1*(-.8481299777*BesselJ(1, 5.429681040*r)+1.002082035*BesselY(1, 5.429681040*r))*Besselj(0, 9.604121522)*(1.-1.*exp(-23.059787552335899121*t))/Besselj(1, 4.802060761)

(1)

``

F53 := proc (t) options operator, arrow; eval(T(r, t), r = 1.5) end proc; G53 := proc (t) options operator, arrow; Re(evalf[20](F53(t))) end proc; Temp3 := cat("T-t (N=", No, ",r=1.5).dat"); fd53 := fopen(Temp3, WRITE, TEXT); i := 0; for i from 0 by 0.5e-1 to 500 do tt := i; y := eval(G53(t), t = tt); fprintf(fd53, "%10.5e %10.5e  \n", tt, y) end do; fclose(fd53)

proc (t) options operator, arrow; eval(T(r, t), r = 1.5) end proc

 

Error, (in fprintf) number expected for floating point format

 

``

``


Download float.mw

Please Wait...