MaplePrimes - answers and comments on Question, Problem with procedures

use these procedures with ScientificErrorAnalysis

jakubi — Thu, 04 Feb 2010 11:29:05 Z

You have written

dWdt := proc (t::float)...

I.e., you request that the argument of dWdt be of type float. But the ScientificErrorAnalysis:-Quantity object is not of type float:

type(ScientificErrorAnalysis:-Quantity(231.928, 0.5e-2),function);
                                 true

type(ScientificErrorAnalysis:-Quantity(231.928, 0.5e-2),float);
                                false

Hence, the error message.

keeping the float condition

herclau — Fri, 05 Feb 2010 19:59:00 Z

Thanks Jakubi,
I'm clear, you suggest  remove the variable declaration 
procedure flaot?
How else, could declare variables, keeping the float condition 
and that the procedure can be used 
with the package ScientificErrorAnalysis?


What can you say about the first point:
1 - I do not understand because i get different values 
when evaluating different significant figures.


Thank you.

adjustment in the code?

herclau — Sat, 06 Feb 2010 00:54:25 Z

Thanks Jakubi,
I will it. It's that simple? 
Or should I make any further adjustment in the code?

No results found

herclau — Mon, 08 Feb 2010 19:14:17 Z

I've implemented it targeted. No results found.

dWdtSEA := proc (t::{float, specfunc(anything,ScientificErrorAnalysis:-Quantity)})
local DW, DW2, T;
description "Se calcula dWdt según el valor de temperatura";
T := t+273.15; if 273.15 <= T then
DW := add(((1/481)*i-1/481)*C[i, 1]*((1/481)*T+(-1/481)*754.15)^(i-2), i = 2 .. 10)
else
DW2 := add((i-1)*A[i, 1]*((ln(T/273.16)+1.5)/1.5)^(i-2)/(1.5*T), i = 2 .. 13);
DW := Wr(t)*DW2
end if;
evalf[9](DW)
end proc;

t := Quantity(ta, dt, relative);

t := Quantity(ta, dt*ta)

combine(dWdtSEA(t), errors, correlations = false);

Error, (in dWdtSEA) cannot determine if this expression is true or false: 0. <= ScientificErrorAnalysis:-Quantity(ta, dt*ta)

t1 := Quantity(ta, dt);

t1 := Quantity(ta, dt)

combine(dWdtSEA(t1), errors);

Error, (in dWdtSEA) cannot determine if this expression is true or false: 0. <= ScientificErrorAnalysis:-Quantity(ta, dt)

evalf(t); ScientificErrorAnalysis:-GetError(t);
ta
dt*ta

y:=combine(dWdtSEA(t1), errors);
Error, (in dWdtSEA) cannot determine if this expression is true or false: 0. <= ScientificErrorAnalysis:-Quantity(ta, dt)

evalf(y); ScientificErrorAnalysis:-GetError(y);
y
Error, (in ScientificErrorAnalysis:-GetError) expecting a quantity-with-error structure, but got y

Thank you.

boolean

jakubi — Mon, 08 Feb 2010 23:29:31 Z

Boolean expressions do not include ScientificErrorAnalysis:-Quantity objects (actually, ?boolean is not very explicit). For them, as explained above, you have to extract the magnitude by evalf:

if evalf(ScientificErrorAnalysis:-Quantity( 2.3, 0.1))<1 then 1 
else 0 end if;
                                  0

simple example

herclau — Tue, 09 Feb 2010 00:51:49 Z

jakubi,

Could demonstrate a simple example, how to enable this, within a procedure.
If you could be with the procedure dWdt or dWdTSEA better!!!!

HerClau

Gracias!!!!!!