## 45 Reputation

8 years, 274 days

## semi circular plate...

Thanks for all, and especially @Rouben Rostamian  . Your code is Ok:).

I have to apply it to a semicircular plate, whose boundary conditions are that the temperature is 0 along the diameter and  100 along the contour of the semicircle. In the The boundary conditions match the exact solution: from the code I wrote this

g:= (x,y)-> `if`(y=sqrt(1-x^2), 100, `if`(x=sqrt(1-y^2), 100,`if`(x=0, 0, 'procname'(args))));

with

Here is our solution: from the code:
x, y, w := algo_12_1([-1,1],[0,1],[6,3],1e-10,100):
Error, (in algo_12_1) unable to store '.2500000000*g(-1., HFloat(0.6666666666))+.2500000000*g(HFloat(-0.6666666667), 1.)' when datatype=float[8]
but this is not well :(.

Regards

## try more conditons...

This is not as well for more conditons.:

g:= (x,y)-> `if`(x=0.5, 200*y, 'procname'(args)) or  (y=0.5, 200*x, 'procname'(args));

regards.

## more conditions...

Thank you for all,

To insert more than one condition?, for example four conditions for the same function g:

g(0,y):=0, g(x,0):=0

g(x,0.5)=200*x ; g(0.5,y)=200*y

Thank you so much.

## @acer  I'm sorry, I do not understa...

I'm sorry, I do not understand you. I sent the file in which I am working. How do I make the succession of matrix?.

example.mw

## @acer @Preben Alsholm Thank you for your...

@acer @Preben Alsholm

Thank you for your comments. I have a sequence of matrices 5x5, and have to work with the elements of each matrix. Given a matrix A [k] 5x5 , with k: = 0, I have to work with the same elements as follows:

j := 1;
1
l := 2;
2
A[k+1](j,l):=A[k](j,l)-v[l] z[j]-v[j] z[l];

v,z are vectors.

But it fails. Can I do it ?, otherwise fails Why ??
Error, Got internal error in Typesetting:-Parse : "'_Inert_DELAYSUM' is not a valid inert form"
Regards;

## @Kitonum Thanks!!!...

@Kitonum Thanks!!!

## @Mac Dude  Hi, I resolved this with...

Hi,

I resolved this with w[n].

Regards.

Thanks.

Best Regards.

## Thanks...

@Adri van der Meer , This is correct.

## @Carl Love  when I run the command ...

when I run the command num_n (1) for example, the system returns x =1; u (1) = ..., u '(1) = ....

The question is:
Through num_n (), How I can extract only u (x)?

Regards.